*** Hypoellipse: Unix/Non-Xpick/Y2K version 3.9 11/1/2001 *** Configured for up to 1024 stations in station list and up to 1501 records per earthquake. Run on 04/12/03 at 11:42 jump hypoellipse.controls 12/3/04 cent va 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 12/3/04 cent va begin station list -1 20000101 station list code = -1 set up for events starting on 20000101 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 cnyrmody hr 12/3/04 cent va 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). 51 1000.0000 1000.0000 for epicentral distance beyond this, use first travel-time table. 52 2800.0000 2800.0000 Wood Anderson static magnification assumed for local magnitude determination. 53 1.0000 1.0000 if .eq. 1 stations with 4-letter codes ending e or n treated as horizontals. 54 200.0000 200.0000 if 1st computed trial location > this (km) from closest station, start at closest station. 55 19.0000 20.0000****assumed century for events without summary record. 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 4 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 -------------------------------------------------------------- ---------------------------------------------------------------- ---------------------------------------------------------------- 04/12/03 01:27 12/3/04 cent va -az/dp--step---se =az/dp==step===se -az/dp--step---se 259/ 3 -.012 2.49 351/38 -22.0 32.0 165/52 -1.52 7.46 horizontal and vertical single variable standard deviations (68% - one degree of freedom; max 99 km) seh = 1.35 seh = 13.71 sez = 10.96 quality = d az = -99. az = -9. se of orig = 0.99; # of iterations = 7; dmax = 400.00; sequence number = event type = " "; processing status = " " s minus p interval for closest station = 6.82 date origin lat long depth mag no d1 gap d rms avwt se 20041203 127 13.99 37n52.70 77w57.78 10.00* 9 52 197 1 0.3762 1.00 0.50 37.8784 77.9631 seh sez q sqd adj in nr avr aar nm avxm mdxm sdxm nf avfm mdfm sdfm vpvs 13.7 11.0 d d d 0.11 11 16 .0000 .302 0 0.0 0 0.0 0.000 -- 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 cvva z ep 22.50 -0.27 0.51 52.5 288 101 1 8.51 8.77 .005 -0.27 cvva cvva z s es 29.32 0.14 0.59 101 2 15.33 15.18 .008 0.14 cvva s cbn z ep 24.23 -0.27 0.51 63.2 55 99 1 10.24 10.50 .002 -0.27 cbn cbn z s es 32.49 0.32 0.58 99 2 18.50 18.17 .003 0.32 cbn s bla z ep 3 47.80 -0.42 1.51 229.4 252 48 1 33.81 34.16 .070 -0.42 bla bla z s es 3 74.59 1.49 1.75 48 2 60.60 58.99 .120 1.49 bla s fwv z ep 3 51.88 0.72 1.77 253.3 263 48 1 37.89 37.08 .085 0.72 fwv fwv z s es 3 81.64 3.46 m ----- 48 2 67.65 64.04 .147 3.46 fwv s eln z ep 3 51.39 0.03 1.76 255.1 256 48 1 37.40 37.30 .070 0.03 eln eln z s es 3 81.32 2.79 m ----- 48 2 67.33 64.42 .120 2.79 eln s mcwv z ep 4 4.00 12.50 ----- 256.5 321 48 1 50.01 37.48 .031 12.50 mcwv mcwv z s es 4 25.28 6.51 ----- 48 2 71.29 64.73 .053 6.51 mcwv s pwv z ep 3 55.38 1.08 2.11 278.9 258 48 1 41.39 40.21 .090 1.08 pwv pwv z s es 3 88.69 5.08 m ----- 48 2 74.70 69.46 .156 5.08 pwv s sspa z ep 4 1.59 4.03 ----- 306.2 1 48 1 47.60 43.55 .017 4.03 sspa sspa z s es 4 38.44 9.18 ----- 48 2 84.45 75.24 .030 9.18 sspa 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 z n nw se ne sw e ave. of end points 0.20 0.29 0.59 0.67 0.93 1.06 1.54 number rms min drms ave drms quality 9 0.38 0.27 0.81 b jump back to main input stream. sound bell  completed reading input phase file average rms of all events = 0.86874 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 0.0 1.0 1.0 percentage/ 0.0 0.0 0.0 100.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 blaz 1 0.1-0.424 0.000 1 1.0-0.424 0.000 1 0.1 1.485 0.000 1 0.7 1.485 0.000 bla cbnz 1 1.0-0.270 0.000 1 8.0-0.270 0.000 1 0.7 0.321 0.000 1 6.0 0.321 0.000 cbn cvvaz 1 1.0-0.274 0.000 1 8.0-0.274 0.000 1 0.7 0.137 0.000 1 6.0 0.137 0.000 cvva elnz 1 0.1 0.027 0.000 1 0.7 0.027 0.000 1 0.1 2.786 0.000 1 0.5 2.786 0.000 eln fwvz 1 0.1 0.721 0.000 1 0.7 0.721 0.000 1 0.1 3.459 0.000 1 0.6 3.459 0.000 fwv pwvz 1 0.1 1.081 0.000 1 0.5 1.081 0.000 1 0.1 5.076 0.000 1 0.4 5.076 0.001 pwv s-p residuals x-mag res f-mag res station n wt ave sd n ave sd n ave sd irelo = 0 nreloc = 0