*** Hypoellipse - Unix version 2.5 October 1993 *** Configured for up to 556 stations in station list and up to 1501 records per earthquake. run on 98/04/22 at 11:36 jump hypoellipse.controls cva eq 4/21/98 220 do not scan for missing stations. code = 1 global minimum search turned off subroutine input1 found end of file on unit 12 list of stations available for these solutions cva eq 4/21/98 220 begin station list +1 010101 non indexed station list. code = 1 set up for events starting on 10101 0 name latitude longitude elev p thickness p p pdy1 sdy1 pdy2 sdy2 pdy3 sdy3 pdy4 sdy4 pdy5 sdy5 calr xmgc mgwt fmgc wt * continuation record * thk 1 2 mod dly sys ps polarity stawt teldy code altdy yrmody hrmn jump to station list: hypoellipse.sta bla 37n 12.68 80w 25.21 634 1 0.000.00 1 1 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 0.00 0.00 1 1 1.00 1010 * 1.0 0.00 0.00 999999 9999 cvl 37n 58.88 78w 27.65 167 1 0.000.00 1 1 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 0.00 0.00 1 1 1.00 1010 * 1.0 0.00 0.00 999999 9999 ghv 37n 47.65 78w 6.44 107 1 0.000.00 1 1 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 0.00 0.00 1 1 1.00 1010 * 1.0 0.00 0.00 999999 9999 na12 37n 59.29 77w 52.62 134 1 0.000.00 1 1 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 0.00 0.00 1 1 1.00 1010 * 1.0 0.00 0.00 999999 9999 cva eq 4/21/98 220 test variables description standard reset to 1 1.7800 1.7000****ratio of p-wave velocity to s-wave velocity. 2 5.0000 0.0000****lt 0 no elev cor/ =0 use 1st vel/ gt 0 use this. 3 0.0000 0.0000****first trial latitude in degrees. 4 0.0000 0.0000****first trial longitude in degrees. 5 -99.0000 5.0000****first trial depth in kilometers, unless = -99. 6 0.0000 -10.0000****sphere rad for aux rms values. if neg cont iteration at most neg point. 7 10.0000 99.0000****minimum number of first motions required to plot. 8 0.0000 0.0000 elevation of top of layered models (km). 9 0.0000 1.0000****if 0 allow neg depths in summary and archive files. 10 50.0000 1.0000****apply distance weighting on this iteration. 11 50.0000 50.0000****xnear = greatest distance with weight of 1.0 12 100.0000 400.0000****xfar = least distnace with weight of 0.0 13 50.0000 50.0000****apply azimuthal weighting on this iteration. 14 50.0000 4.0000****weight out large residuals on this iteration. 15 10.0000 2.0000****give zero weight to residuals gt this. 16 50.0000 50.0000****apply boxcar weighting on this iteration. 17 2.0000 2.0000****give zero weight to residuals gt this*stand. dev. 18 50.0000 50.0000****begin jeffreys weighting on this iteration. 19 0.0500 0.0500****use jeffreys weighting only if rms gt this. 20 0.0500 0.0500****mu of jeffreys weighting funct. 21 9.0000 9.0000 maximum number of iterations. 22 35.0000 35.0000 limit change in focal depth to this amount (km). 23 0.7000 0.5000****if delz would make z neg, set delz = -this*z (km). 24 35.0000 35.0000 limit change in epicenter to this. (km). 25 40.0000 40.0000 fix depth if epicentral change gt this. (km). 26 0.0025 0.0025 stop iterating if square of adjustment lt this. 27 20.0000 0.0000****global opt: if deep solution z > this, continue with z 1/2 way to surface. 28 0.0000 400.0000****for fixed hypo on plane, set = plunge azimuth. if neg. continue as free sol. 29 -0.1000 0.3000****set std err of res=+this if degrees of freedom =0 or =-this if this lt 0. 30 0.0000 0.0000 dip of plunge vector for epi. fixed on plane. see test(28) & (47) also. 31 -1.1500 -3.4500****duration magnitude c1, constant. 32 2.0000 2.8500****duration magnitude c2, *log((f - p)*fmgc). 33 0.0000 0.0000****duration magnitude c3, *delta. 34 0.0000 0.0000 if not 0, scale the normal equations. 35 0.0010 0.0010****minimum damping of normal equations. 36 100.0000 100.0000****maximum first trial depth if computed from p-arrival times. 37 3.0000 3.0000 if termination occurs before this iteration, set iteration number to this and continue. 38 0.0000 0.0000****if this =1, run all with and then without s/ =2,run with s/ =3, run without s/ =4, fix hypo / neg, use s to fix origin. 39 1.0000 0.7500****multiply the s and s-p weights by this factor. 40 0.0070 0.0000****duration magnitude c4, *depth. 41 0.0000 0.0000 if this =1, print opt. ge 1, & summary opt. =+ or -1, then write sum. record each itteration. 42 75.0000 75.0000****global opt: deep starting z wrt top of model. 43 0.0000 0.0000 duration magnitude c5, *(log((f - p)*fmgc)**2). 44 0.0000 0.0000 if =1 rerun debug eqs with critical sta/ =2 continue iter with crit sta. 45 0.1379 0.1379 x scale factor for focal mechanism plot. 46 0.0000 0.0000 xfar set ge dist of test(46)th station + 10. if lt 0 then fill gap. 47 0.0000 0.0000 weight for fix on plane. see test(28) and (30). 48 6.5000 6.5000 half-space velocity for first trial location. 49 0.0000 0.0000****if .ne. 0 calculate vp/vs ratio; if abs val >1 make wadati plot; if neg, use wadati origin in solution. 50 0.0000 0.0000 for exploring rms space, compute this number of fixed depth solutions (up to 22). weight option - relative standard errors for code: 0 1 2 3 1.000 1.153 1.414 2.000 printer option 0 summary option 2 magnitude option 1 tabulation option 4 no event output -2 one line/eq -1 final solution 0 no sum records 0 use xmag 0 no summary 0 one line per iter 1 summary records 1 use fmag 1 a 1 sta res each iter 2 sum + archive file 2 use (xmag+fmag)/2 2 a + b 2 regres each iter 3 archive file 3 prefer fmag /xmag 3 a,b + c 3 "corrected" input 5 prefer xmag /fmag 4 a,b,c + d 4 if neg use fms not fmp positive/q from std errors negative/q from sol+sta u of a cal data file: pub1:[alaska.data]uofacal.dat make compensating change in layer below variable layer. velocity model 1 layer velocity depth thickness vpvs km/sec km km 1 6.090 0.000 15.000 1.730 2 6.500 15.000 21.000 1.720 3 8.180 36.000 1000.000 1.730 the next model is for s only: velocity model 2 layer velocity depth thickness vpvs km/sec km km 4 3.520 0.000 15.000 0.000 5 3.779 15.000 21.000 0.000 6 4.728 36.000 1000.000 0.000 jump to hypoellipse.arrivals ---------------------------------------------------------------- ---------------------------------------------------------------- ---------------------------------------------------------------- 98/04/21 23:28 cva eq 4/21/98 220 -az/dp--step---se =az/dp==step===se -az/dp--step---se 11/ 9 -.002 3.43 277/24 -.001 1.70 119/64 .0366 8.65 se of orig = 0.44; # of iterations = 4; dmax = 400.00; sequence number = event type = " "; processing status = " " s minus p interval for closest station = 3.15 date origin lat long depth mag no d1 gap d rms avwt se 980421 2328 26.60 38n10.24 78w34.16 1.97 8 23 231 1 0.1627 1.00 0.30 seh sez q sqd adj in nr avr aar nm avxm mdxm sdxm nf avfm mdfm sdfm 2.3 4.2 d c d 0.10 10 8 .0000 .151 0 0.0 0 0.0 -- travel times and delays -- stn c pha remk p p-sec s-sec resid std-er dist azm ain tc c vthk ttob-ttcal-dlay-edly=resid rmk stn pha sources cvl ipd d 30.30 -0.11 0.30 23.1 156 95 1 3.70 3.80 .002 -0.11 cvl cvl s is 33.45 0.27 0.35 95 2 6.85 6.58 .004 0.27 cvl s ghv ipc c 36.00 -0.17 0.30 58.3 136 92 1 9.40 9.57 .001 -0.17 ghv ghv s is 42.98 -0.18 0.35 92 2 16.38 16.56 .001 -0.18 ghv s na12 epd d 37.20 0.08 0.30 64.0 108 92 1 10.60 10.52 .001 0.08 na12 na12 s es 44.90 0.10 0.35 92 2 18.30 18.20 .001 0.10 na12 s bla epn n 57.70 0.21 0.39 194.9 237 48 1 31.10 30.82 .070 0.21 bla bla s esn 79.80 -0.14 0.45 48 2 53.20 53.21 .120 -0.14 bla s -- magnitude data -- stn c source sys c10 amx gr ink amf per unit/mm gnd mot u xmgc xmag fmp fmag quality evaluation diagonals in order of strength se ne z n sw nw e ave. of end points 0.41 0.44 0.46 0.50 0.85 0.98 1.06 number rms min drms ave drms quality 8 0.16 0.37 0.70 a jump back to main input stream. sound bell  completed reading input phase file average rms of all events = 0.16295 average vp/vs ratio = 0.00 for 0 events. standard deviation of ratio = 0.00 ***** class/ a b c d total ***** number/ 0.0 0.0 1.0 0.0 1.0 percentage/ 0.0 0.0 100.0 0.0 include only class d and better in the following statistics. ----------- p residuals -------------- ----------- s residuals ------------- no event wting event wting no event wting event wting station n wt ave sd n wt ave sd n wt ave sd n wt ave sd station bla 1 0.6 0.208 0.000 1 2.5 0.208 0.000 1 0.5-0.135 0.000 1 1.9-0.135 0.000 bla cvl 1 1.0-0.106 0.000 1 4.0-0.106 0.000 1 0.8 0.266 0.000 1 3.0 0.266 0.000 cvl ghv 1 1.0-0.174 0.000 1 4.0-0.174 0.000 1 0.7-0.181 0.000 1 3.0-0.181 0.000 ghv na12 1 1.0 0.078 0.000 1 4.0 0.078 0.000 1 0.7 0.099 0.000 1 3.0 0.099 0.000 na12 s-p residuals x-mag res f-mag res station n wt ave sd n ave sd n ave sd irelo = 0 nreloc = 0