*** 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 97/12/29 at 17:14 c* jump vtso.controls 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 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: 218.stations 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 eln 37n 13.70 80w 45.10 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 wmv 37n 6.51 80w 58.23 1157 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 mcwv 39n 39.49 79w 50.74 280 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 ceh 35n 53.46 79w 5.58 152 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 sltn 36n 26.59 82w 7.23 1280 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 subroutine input1 found end of file on unit 12 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 -3.0000****minimum number of first motions required to plot. 8 0.0000 1000.0000****elevation of top of layered models (km). 9 0.0000 360.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.7000 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 20.0000 global opt: if deep solution z > this, continue with z 1/2 way to surface. 28 0.0000 0.0000 for fixed hypo on plane, set = plunge azimuth. if neg. continue as free sol. 29 -0.1000 -0.1000 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 1.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 5.630 0.000 5.700 1.640 2 6.050 5.700 9.000 1.720 3 6.530 14.700 36.000 1.700 4 8.180 50.700 1000.000 1.710 the next model is for s only: velocity model 2 layer velocity depth thickness vpvs km/sec km km 5 3.433 0.000 5.700 0.000 6 3.517 5.700 9.000 0.000 7 3.841 14.700 36.000 0.000 8 4.784 50.700 1000.000 0.000 jump to 218.arrivals ---------------------------------------------------------------- ---------------------------------------------------------------- ---------------------------------------------------------------- 97/12/21 07:44 -az/dp--step---se =az/dp==step===se -az/dp--step---se 116/ 3 -.020 4.42 207/ 6 -.036 2.15 359/84 1153. 99.0 se of orig = 0.21; # of iterations = 10; dmax = 400.00; sequence number = event type = " "; processing status = " " s minus p interval for closest station = 9.22 date origin lat long depth mag no d1 gap d rms avwt se 971221 744 46.71 37n59.13 79w57.20 -9.50 5 95 343 1 0.7319 1.00 0.10 seh sez q sqd adj in nr avr aar nm avxm mdxm sdxm nf avfm mdfm sdfm 8.5 99.0 d d d 2.50 10 12 0.000 .677 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 bla EP 58.00 -0.41 0.10 95.3 206 85 1 11.29 11.70 -0.41 bla bla s ES 67.22 0.50 0.12 85 2 20.51 20.01 0.50 bla s eln EP 59.42 -0.75 0.10 109.7 220 85 1 12.71 13.45 -0.75 eln eln s ES 69.23 -0.49 0.12 85 2 22.52 23.00 -0.49 eln s wmv EP 04.09 1.15 0.10 132.5 223 86 1 17.38 16.23 1.15 wmv wmv s ES 17.16 2.70 m ----- 86 2 30.45 27.75 2.70 wmv s mcwv EP 4 00.00 -9.47 ----- 185.9 3 87 1 13.29 22.76 -9.47 mcwv mcwv s ES 40.36 14.73 m ----- 87 2 53.65 38.91 14.73 mcwv s ceh EP 4 00.00 -16.65 ----- 244.7 161 88 1 13.29 29.94 -16.65 ceh ceh s ES 46.60 8.69 m ----- 88 2 59.89 51.20 8.69 ceh s sltn EP 4 25.52 7.32 ----- 257.5 229 88 1 38.81 31.49 7.32 sltn sltn s ES 50.86 10.30 m ----- 88 2 64.15 53.85 10.30 sltn 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 nw se e ne sw n ave. of end points 0.00 0.02 0.03 0.04 0.06 0.08 0.09 number rms min drms ave drms quality 5 0.73 0.00 0.05 d jump back to main input stream. sound bell  completed reading input phase file average rms of all events = 4.72975 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 bla 1 1.0-0.413 0.000 1 4.8-0.413 0.000 1 0.7 0.499 0.000 1 3.6 0.499 0.000 bla eln 1 0.9-0.745 0.000 1 4.6-0.745 0.000 1 0.7-0.487 0.000 1 3.5-0.487 0.000 eln wmv 1 0.9 1.149 0.000 1 4.3 1.149 0.000 1 0.7 2.697 0.000 1 3.3 2.697 0.000 wmv mcwv 0 0.0 0.000 0.000 0 0.0 0.000 0.000 1 0.514.734 0.000 1 2.514.734 0.000 mcwv ceh 0 0.0 0.000 0.000 0 0.0 0.000 0.000 1 0.3 8.690 0.000 1 1.5 8.690 0.000 ceh sltn 0 0.0 0.000 0.000 0 0.0 0.000 0.000 1 0.310.298 0.000 1 1.310.298 0.000 sltn s-p residuals x-mag res f-mag res station n wt ave sd n ave sd n ave sd irelo = 0 nreloc = 0