1 | C--COMMON BLOCK INCLUDE FILE FOR HYPOINVERSE |
---|
2 | LOGICAL LSUM,LARC,LPRT,LERR,LEJCT,LREP,DONE,LCUSP,LPRALL |
---|
3 | LOGICAL LAPP,LJUNK,LMULT,LTBIG,LMAG,LBSTA,LBCRU,SUBMOD,LP153 |
---|
4 | LOGICAL L2000,HYPOFIX,ALLFIX, LASK |
---|
5 | |
---|
6 | C--SHORT INTEGER AND LOGICAL DECLARATIONS |
---|
7 | INCLUDE 'integer.for' !SUN & VAX VERSION. USES *2 LENGTHS |
---|
8 | |
---|
9 | C--GENERAL PARAMETERS |
---|
10 | CHARACTER CM*3,INST*132,CTEMP*1,IE*1,IS*1,TERM*80 |
---|
11 | C CHARACTER CPADD*2 |
---|
12 | CHARACTER RUNLAB*1,CP1*1,CP2*1,CP3*1,FIXCHR*1, CDOMAN*2,CPVERS*2 |
---|
13 | COMMON /D/ SUBMOD !F=NORMAL PROGRAM, T=SUBROUTINE VERSION |
---|
14 | C COMMON /C/ CPADD !DUMMY PADDING IF NEEDED |
---|
15 | COMMON /C/ CM !CURRENT COMMAND NAME |
---|
16 | COMMON /C/ INST !CURRENT COMMAND PARAMETERS OR PHASE CARD |
---|
17 | COMMON /C/ CTEMP !CHARACTER TEMPORARY |
---|
18 | COMMON /C/ TERM !TERMINATOR (INSTRUCTION) CARD |
---|
19 | COMMON /C/ IS !CHARACTER TEMPORARY, 'S' FOR SOUTH LAT |
---|
20 | COMMON /C/ IE !CHARACTER TEMPORARY, 'E' FOR EAST LON |
---|
21 | COMMON /D/ PI !PI |
---|
22 | COMMON /D/ RDEG !DEGREES PER RADIAN |
---|
23 | COMMON /D/ INP !COMMAND INPUT UNIT NUMBER |
---|
24 | COMMON /D/ ISTAT !A STATUS FLAG SET BY HYCMD FOR HYP |
---|
25 | COMMON /D/ ISTAT2 !A STATUS FLAG SET BY HYOPEN FOR HYP |
---|
26 | COMMON /D/ IRES !STATUS FLAG FROM HYPOINV & CUSP SUBROUTINES |
---|
27 | COMMON /D/ IPADD !DUMMY PADDING |
---|
28 | COMMON /D/ LTBIG !BIG EVENT FLAG (HYPHS TO HYLST) |
---|
29 | COMMON /C/ RUNLAB !1-LETTER LABEL FOR ENTIRE RUN (LAB) |
---|
30 | COMMON /D/ LP153 !T=PASS COL 153 THRU, F=RUNLAB (LAB) |
---|
31 | COMMON /C/ CP1 !PASS THRU LABEL, IE AUTHORITY |
---|
32 | COMMON /C/ CP2 !PASS THRU LABEL, IE VERSION |
---|
33 | COMMON /C/ CP3 !PASS THRU LABEL, IE REVISION |
---|
34 | COMMON /D/ HYPOFIX !INDIVIDUAL FLAG TO FIX ONE HYPOCENTER |
---|
35 | COMMON /D/ ALLFIX !INDIVIDUAL FLAG TO FIX ONE HYPO & ORIGIN |
---|
36 | COMMON /C/ FIXCHR !CHARACTER USED TO SIGNAL HYPOCENTER FIX |
---|
37 | COMMON /C/ CDOMAN !CODE FOR PROCESSING DOMAIN (IE NC) (VER) |
---|
38 | COMMON /C/ CPVERS !CODE FOR PROCESSING VERSION (IE 01) (VER) |
---|
39 | |
---|
40 | C--FILES AND OUTPUT CONTROLS |
---|
41 | C--DEFINING COMMAND NAMES ARE IN PARENTHESES |
---|
42 | CHARACTER*60 STAFIL,INFILE,PHSFIL,SUMFIL,ARCFIL,PRTFIL,DELFIL |
---|
43 | CHARACTER*60 ATNFIL,FMCFIL,XMCFIL,MAGFIL,BSTAFL,BCRUFL,TERMIN |
---|
44 | CHARACTER*60 CALFIL |
---|
45 | CHARACTER FORID*30, GREETING*70 |
---|
46 | COMMON /C/ SUMFIL !HYPOINVERSE SUMMARY OUTPUT FILENAME (SUM) |
---|
47 | COMMON /C/ ARCFIL !ARCHIVE OUTPUT FILENAME (ARC) |
---|
48 | COMMON /C/ PRTFIL !PRINT OUTPUT FILENAME (PRT) |
---|
49 | COMMON /D/ LSUM !TRUE TO OUTPUT SUMMARY DATA (SUM) |
---|
50 | COMMON /D/ IH71S !SUMMARY FORMAT 1=HI 2=H71 (H71) |
---|
51 | COMMON /D/ IH71T !TERMINATOR FORMAT 1=HI 2=H71 3=HDR TRIAL (H71) |
---|
52 | COMMON /D/ LARC !TRUE TO OUTPUT ARCHIVE DATA (ARC) |
---|
53 | COMMON /D/ LPRT !TRUE TO OUTPUT A PRINTOUT FILE (PRT) |
---|
54 | COMMON /D/ LMAG !TRUE TO OUTPUT MAGNITUDE DATA FILE (MFL) |
---|
55 | COMMON /D/ LERR !TRUE TO OUTPUT ERROR MESSAGES TO TERMINAL (ERF) |
---|
56 | COMMON /D/ LCUSP !TRUE IF READING FROM CUSP MEM FILES (COP) |
---|
57 | COMMON /C/ FORID !CUSP ID READING FORMAT (FID) |
---|
58 | COMMON /C/ STAFIL !STATION INPUT FILENAME (STA) |
---|
59 | COMMON /C/ BSTAFL !BINARY STATION SNAPSHOT (WST,RST) |
---|
60 | COMMON /D/ LBSTA !T IF BINARY STATIONS, F IF ASCII (RST,STA,DEL) |
---|
61 | COMMON /C/ BCRUFL !BINARY CRUST MODELS SNAPSHOT (WCR,RCR) |
---|
62 | COMMON /D/ LBCRU !T IF BINARY CRUST MODELS, F IF ASCII (RCR,CRT) |
---|
63 | COMMON /C/ DELFIL !STATION DELAY FILENAME (DEL) |
---|
64 | COMMON /C/ ATNFIL !STATION ATTENUATION FILENAME (ATE) |
---|
65 | COMMON /C/ CALFIL !STATION CAL FACTOR FILENAME (CAL) |
---|
66 | COMMON /C/ FMCFIL !STATION FMAG CORRECTION FILENAME (FMC) |
---|
67 | COMMON /C/ XMCFIL !STATION XMAG CORRECTION FILENAME (XMC) |
---|
68 | COMMON /C/ PHSFIL !PHASE DATA INPUT FILENAME (PHS) |
---|
69 | COMMON /C/ MAGFIL !MAGNITUDE DATA OUTPUT FILE (MFL) |
---|
70 | COMMON /C/ INFILE(0:4) !COMMAND INPUT FILENAME (@) |
---|
71 | COMMON /D/ ISTFMT !STATION FORMAT NO. 1=HI 2=H71 3=HI2 (STF) |
---|
72 | COMMON /D/ JCP !PHASE DATA FORMAT CODE (COP) |
---|
73 | COMMON /D/ JCPO !1=MODIFY CUSP MEM FILES, 0=NO (COP) |
---|
74 | COMMON /D/ JCA !ARCHIVE DATA FORMAT CODE (CAR) |
---|
75 | C THE CODES ARE: 1=FULL 2=CONDENSED 3=CONDENSED PICKER STYLE |
---|
76 | COMMON /D/ JST !0=PRINT EVENTS ONLY, 1=ADD PARAMS & FILES TO START |
---|
77 | C OF PRINTOUT, 2=ADD STA LIST & CRUST MODELS (LST) |
---|
78 | COMMON /D/ JST2 !STATION LIST QUANTITY CONTROL 0,1,2 (LST) |
---|
79 | COMMON /D/ JST3 !CRUST LIST QUANTITY CONTROL 0,1 (LST) |
---|
80 | COMMON /D/ KPRINT !CONTROLS AMOUNT OF PRINTOUT (KPR) |
---|
81 | COMMON /D/ LEJCT !TRUE TO PAGE EJECT PRINTOUT BEFORE EACH EVENT (TOP) |
---|
82 | COMMON /D/ LREP !TRUE TO REPORT EACH EVENT ON TERMINAL AS LOCATED (REP) |
---|
83 | COMMON /D/ LPRALL !TRUE TO PRINT UNWEIGHTED STATIONS (REP) |
---|
84 | COMMON /D/ LAPP(3) !INDICATES WHETHER OUTPUT FILES ARE APPENDED TO (APP) |
---|
85 | C 1=PRINT FILE 2=SUMMARY 3=ARCHIVE |
---|
86 | COMMON /C/ TERMIN !THE CURRENT TERMINAL DEVICE CODE |
---|
87 | COMMON /D/ L2000 !T=SUPPORT YR2000 FORMATS, F=OLD FORMATS |
---|
88 | COMMON /D/ ICENT !DEFAULT CENTURY OLD PHASE FORMAT |
---|
89 | COMMON /D/ IAMPU !DEFAULT AMP UNITS CODE |
---|
90 | COMMON /C/ GREETING !GREETING MESSAGE FOR TERMINAL & PRT FILE |
---|
91 | |
---|
92 | C--HYPOINVERSE CRUSTAL MODELS ******************************************** |
---|
93 | C--DATA PRESENT FOR EACH OF THE LH=LM MODELS |
---|
94 | C--LH IS THE MAX NO. OF LAYER MODELS WHICH MUST = LM, THE MAX NO. OF MODELS |
---|
95 | PARAMETER (NLYR=20,LH=36) |
---|
96 | CHARACTER MODNAM*30 |
---|
97 | COMMON /MC/ MODNAM(LH) !MODEL NAME OR LABEL |
---|
98 | COMMON /M/ LAY(LH) !NUMBER OF LAYERS OR V-D POINTS |
---|
99 | COMMON /M/ D(NLYR,LH) !DEPTH TO LAYER TOP OR VELOCITY POINT |
---|
100 | COMMON /M/ VEL(NLYR,LH) !LAYER OR POINT VELOCITY |
---|
101 | COMMON /M/ THK(NLYR,LH) !THICKNESS OF HOMOGENEOUS LAYER |
---|
102 | COMMON /M/ VSQ(NLYR,LH) !SQUARED VELOCITY OF HOMOGENEOUS LAYER |
---|
103 | COMMON /M/ MODTYP(LH) !MDL TYPE (-1=UNDEF, 0=GRAD, 1=HOMO LAYER) |
---|
104 | |
---|
105 | C--DATA USED ONLY FOR LINEAR GRADIENT TRAVEL TIME TABLES |
---|
106 | PARAMETER (LN=36) !THE NUMBER OF LINEAR-GRAD MODELS ALLOWED |
---|
107 | LOGICAL GD1,GD2,GZ1,GZ2 |
---|
108 | COMMON /M/ REDV(LN) !ONE OVER THE REDUCING VELOCITY |
---|
109 | COMMON /M/ NZ(LN) !NUMBER OF DEPTH GRID POINTS |
---|
110 | COMMON /M/ NZ1(LN),DZ1(LN),NZ2(LN),DZ2(LN) !DEPTH GRID PARAMS |
---|
111 | COMMON /M/ ND(LN) !NUMBER OF DISTANCE GRID POINTS |
---|
112 | COMMON /M/ ND1(LN),DD1(LN),ND2(LN),DD2(LN) !DIST GRID PARAMS |
---|
113 | COMMON /M/ GD1(LN),GD2(LN),GZ1(LN),GZ2(LN) !GRID FLAGS |
---|
114 | |
---|
115 | C--DECLARE ARRAY STORAGE FOR TRAVEL TIME TABLE. THE ARRAY KT HOLDS |
---|
116 | C CONDENSED TRAVEL TIMES TO 28 DEPTH & 42 DISTANCE POINTS. |
---|
117 | C KDHR HOLDS THE DISTANCES AT WHICH HORIZONTAL RAYS EMERGE. |
---|
118 | COMMON /M/ KDHR(LN,28), KT(LN,28,42) |
---|
119 | |
---|
120 | C--PARAMETERS FOR MULTIPLE CRUSTAL MODELS |
---|
121 | C--LM IS THE MAXIMUM NUMBER OF MODELS OF EITHER GRADIENT OR LAYER TYPE. |
---|
122 | C LM MUST = LH (MAX NO. OF LAYER MODELS). LN (MAX NO. OF GRADIENT MODELS) |
---|
123 | C MAY BE LESS THAN OR = LM. LH & LN ARE PARAMETERS IN HYPCOMC. |
---|
124 | PARAMETER (LM=36) !MAX NUMBER OF CRUST MODELS |
---|
125 | CHARACTER CRODE*3, CRUFIL*60 |
---|
126 | COMMON /D/ MOD !CURRENT CRUST MODEL NUMBER |
---|
127 | COMMON /C/ CRUFIL(LM) !CRUST MODEL INPUT FILENAMES (CRT,CRH) |
---|
128 | COMMON /C/ CRODE(LM) !CRUST MODEL CODE FROM NAME |
---|
129 | COMMON /D/ MODALT(LM) !=0 NO ALTERNATE MOD, =N USE MODEL N W/SOME STAS |
---|
130 | COMMON /D/ LMULT !T TO PROCESS W/ REGION-DEPENDENT MODELS (MUL) |
---|
131 | COMMON /D/ MODDEF !DEFAULT MODEL FOR UNLISTED REGIONS (MUL) |
---|
132 | COMMON /D/ MAXMOD !HIGHEST MODEL NUMBER DEFINED (CRT,CRH) |
---|
133 | COMMON /D/ NMOD !NO. OF POS WEIGHTED MODELS FOR AN EVENT (1-3) |
---|
134 | COMMON /D/ MODS (3) !MODEL NUMBERS FOR THIS EVENT |
---|
135 | COMMON /D/ WMOD (3) !WEIGHTS OF EACH MODEL |
---|
136 | |
---|
137 | PARAMETER (NODMAX=124) !MAX NO. OF GEOGRAPHIC NODES FOR DEFINING MODELS |
---|
138 | COMMON /D/ NNODE !NUMBER OF NODES CURRENTLY DEFINED (NOD) |
---|
139 | COMMON /D/ HLAT(NODMAX) !NODE CIRCLE CENTER LAT (NOD) |
---|
140 | COMMON /D/ HLON(NODMAX) !NODE CIRCLE CENTER LON (NOD) |
---|
141 | COMMON /D/ RAD1(NODMAX) !NODE CIRCLE INNER RADIUS (NOD) |
---|
142 | COMMON /D/ RAD2(NODMAX) !NODE CIRCLE OUTER RADIUS |
---|
143 | COMMON /D/ DRAD(NODMAX) !NODE CIRCLE TRANSITION WIDTH (NOD) |
---|
144 | COMMON /D/ MODH(NODMAX) !MODEL NO. FOR THIS NODE (NOD) |
---|
145 | COMMON /D/ HDIS(NODMAX) !EPICENTRAL DIST FROM NODE |
---|
146 | |
---|
147 | C--DATA FOR ALL AVAILABLE STATIONS STORED AS LOOKUP TABLE (SUBSCRIPT J) |
---|
148 | C USE 2-BYTE INTEGERS WHERE POSSIBLE |
---|
149 | CHARACTER STANAM*5,STRMK*1,JCOMP1*1,JCOMP3*3,JNET*2,JSLOC*2 |
---|
150 | CHARACTER JCOMPA*3,JSLOC2*2 |
---|
151 | LOGICAL LLOC2 |
---|
152 | C INTEGER*4 JCEXP,JFEXP,JXEXP !LONGER LENGTH REQUIRED |
---|
153 | PARAMETER (MAXSTA=9000) !MAX NO. OF STATIONS STORED IN LOOKUP TABLE |
---|
154 | COMMON /D/ JSTA !NUMBER OF STATIONS ACTUALLY READ INTO TABLE |
---|
155 | C JSTA <= MAXSTA |
---|
156 | COMMON /C/ STANAM (MAXSTA) !5-LETTER STATION NAME |
---|
157 | COMMON /C/ JNET (MAXSTA) !2-LETTER SEISMIC NET CODE |
---|
158 | COMMON /C/ JCOMP1 (MAXSTA) !1-LETTER COMPONENT CODE |
---|
159 | COMMON /C/ JCOMP3 (MAXSTA) !3-LETTER COMPONENT CODE |
---|
160 | COMMON /C/ JCOMPA (MAXSTA) !ALTERNATIVE 3-LETTER COMPONENT CODE |
---|
161 | COMMON /C/ JSLOC (MAXSTA) !2-LETTER LOCATION CODE |
---|
162 | COMMON /C/ JSLOC2 (MAXSTA) !EQUIVALENT LOCATION CODE (' '='--') |
---|
163 | COMMON /D/ LLOC2 !T IF 2ND LOCATION CODE MATCHED |
---|
164 | COMMON /D/ JLATD (MAXSTA) !LATITUDE (DEG) |
---|
165 | COMMON /D/ JLATM (MAXSTA) !LATITUDE (IN .01 MIN) |
---|
166 | COMMON /D/ JLOND (MAXSTA) !LONGITUDE (DEG) |
---|
167 | COMMON /D/ JLONM (MAXSTA) !LONGITUDE (IN .01 MIN) |
---|
168 | COMMON /D/ JPER (MAXSTA) !STANDARD PERIOD (PERIOD IN SEC * 10) |
---|
169 | COMMON /D/ JTYPE (MAXSTA) !RESPONSE TYPE 0=WA 1,3=L4C 4=HS1 ETC |
---|
170 | COMMON /D/ JCAL (MAXSTA) !STATION CAL FACTOR *1000 |
---|
171 | COMMON /D/ JCEXP (MAXSTA) !CAL EXPIRATION DATE (IN YMDH) |
---|
172 | COMMON /D/ JLMOD (MAXSTA) !T TO USE ALTERNATE MODEL FOR THIS STA |
---|
173 | COMMON /D/ JFCOR (MAXSTA) !DUR MAG CORRECTION *100 |
---|
174 | COMMON /D/ JFEXP (MAXSTA) !FCOR EXPIRATION DATE (IN YMDH) |
---|
175 | COMMON /D/ JXCOR (MAXSTA) !AMP MAG CORRECTION *100 |
---|
176 | COMMON /D/ JXEXP (MAXSTA) !XCOR EXPIRATION DATE (IN YMDH) |
---|
177 | COMMON /D/ JPSWT (MAXSTA) !P & S STATION WEIGHTS |
---|
178 | COMMON /D/ JXWT (MAXSTA) !AMPLITUDE MAGNITUDE WEIGHT *10 |
---|
179 | COMMON /D/ JFWT (MAXSTA) !DURATION MAGNITUDE WEIGHT *10 |
---|
180 | COMMON /D/ JFGWT (MAXSTA) !DURATION GAIN CORRECTION WEIGHT *1 |
---|
181 | COMMON /D/ JPD (LM,MAXSTA) !P DELAYS IN .01 SEC |
---|
182 | COMMON /C/ STRMK (MAXSTA) !1-LETTER STATION REMARK |
---|
183 | C--THE FOLLOWING 4 VARIABLES DEPEND ON COMP TYPE |
---|
184 | COMMON /D/ JFM1 (MAXSTA) !USE THIS STATION FOR 1ST DUR MAG |
---|
185 | COMMON /D/ JFM2 (MAXSTA) !USE THIS STATION FOR 2ND DUR MAG |
---|
186 | COMMON /D/ JXM1 (MAXSTA) !USE THIS STATION FOR 1ST AMP MAG |
---|
187 | COMMON /D/ JXM2 (MAXSTA) !USE THIS STATION FOR 2ND AMP MAG |
---|
188 | COMMON /D/ JPM1 (MAXSTA) !USE THIS STATION FOR 1ST P AMP MAG |
---|
189 | COMMON /D/ JPM2 (MAXSTA) !USE THIS STATION FOR 2ND P AMP MAG |
---|
190 | |
---|
191 | C************** DATA FOR A PARTICULAR EVENT ********************* |
---|
192 | C--DATA PRESENT FOR EACH STATION REPORTING AN EVENT (SUBSCRIPT K) |
---|
193 | C USE 2-BYTE INTEGERS WHERE POSSIBLE |
---|
194 | PARAMETER (MAXPHS=1500) !MAX NO. OF STATIONS FOR EACH EVENT |
---|
195 | CHARACTER KPRK*3,KSRK*2,KRMK*1,KSOU*1,KDEV*3,KRMK6*6,PARMK*1 |
---|
196 | C CHARACTER KPSWT*1 |
---|
197 | CHARACTER KSHAD*104 !A MAX LENGTH OF 104 IS TESTED IN HYPHS |
---|
198 | COMMON /D/ KSTA !NO. OF STATIONS (PHASE CARDS) FOR THIS EVENT |
---|
199 | C KSTA <= MAXPHS |
---|
200 | COMMON /D/ KINDX (MAXPHS) !POINTER TO STATION IN LOOKUP TABLE |
---|
201 | COMMON /D/ KP (MAXPHS) !P TIME IN .01" |
---|
202 | COMMON /D/ KS (MAXPHS) !S TIME IN .01" |
---|
203 | COMMON /D/ DIS (MAXPHS) !DISTANCE TO EPICENTER IN KM |
---|
204 | COMMON /D/ KAZEM (MAXPHS) !CONDENSED AZIMUTH & EMERGENCE ANGLE |
---|
205 | COMMON /D/ KWT (MAXPHS) !CONDENSED ASSIGNED P & S WEIGHT CODES |
---|
206 | COMMON /D/ KFWT (MAXPHS) !ASSIGNED DUR MAGNITUDE WEIGHT CODE |
---|
207 | COMMON /D/ KXWT (MAXPHS) !ASSIGNED AMP MAGNITUDE WEIGHT CODE |
---|
208 | COMMON /D/ KPAWT (MAXPHS) !ASSIGNED P AMP MAG WEIGHT CODE |
---|
209 | COMMON /D/ PAWT (MAXPHS) !P AMP MAGNITUDE WEIGHT |
---|
210 | COMMON /D/ KFMP (MAXPHS) !CODA DURATION IN SECONDS |
---|
211 | COMMON /D/ KXMAG (MAXPHS) !CALCULATED AMP MAG FOR THIS STA *100 |
---|
212 | COMMON /D/ KFMAG (MAXPHS) !CALCULATED DUR MAG FOR THIS STA *100 |
---|
213 | COMMON /D/ KPMAG (MAXPHS) !CALCULATED P MAG FOR THIS STA *100 |
---|
214 | COMMON /D/ AMPK (MAXPHS) !MAX PEAK-TO-PEAK AMP IN MM OR COUNTS |
---|
215 | COMMON /D/ KAMPU (MAXPHS) !AMP UNITS CODE 0=PP-MM 1=0-P 2=COUNTS |
---|
216 | COMMON /D/ KAMPTYP (MAXPHS) !AMP TYP 0=? 1=WA 2=VEL 3=ACCEL 4=NOUSE |
---|
217 | COMMON /D/ KIMTYP (MAXPHS) !INST TYP 0=NOMAG 1=ML 2=MX |
---|
218 | COMMON /D/ KPAMP (MAXPHS) !PP P AMP IN COUNTS, 1ST 3 PEAKS |
---|
219 | COMMON /D/ KPER (MAXPHS) !PERIOD AT AMP READING IN .01 SEC. |
---|
220 | COMMON /D/ KCAL (MAXPHS) !CAL FACTOR FOR THIS AMP ONLY IN .01 |
---|
221 | COMMON /C/ KPRK (MAXPHS) !P REMARK INCL FIRST MOTION |
---|
222 | COMMON /C/ KSRK (MAXPHS) !S REMARK |
---|
223 | COMMON /C/ PARMK (MAXPHS) !P AMPLITUDE DESCRIPTOR |
---|
224 | COMMON /D/ PNORM (MAXPHS) !PMAG L1 NORM FOR STA TB MINIMIZED |
---|
225 | COMMON /D/ PNORM2 (MAXPHS) !PMAG2 L1 NORM FOR STA TB MINIMIZED |
---|
226 | COMMON /C/ KRMK (MAXPHS) !ANY 1 LETTER STATION REMARK |
---|
227 | COMMON /C/ KRMK6 (MAXPHS) !ORIGINAL 6-LETTER REMARK |
---|
228 | COMMON /C/ KSHAD (MAXPHS) !OPTIONAL SHADOW PHASE RECORD |
---|
229 | COMMON /D/ KLSHA (MAXPHS) !LENGTH OF KSHAD |
---|
230 | COMMON /C/ KSOU (MAXPHS) !1-LETTER DATA SOURCE CODE |
---|
231 | COMMON /C/ KDEV (MAXPHS) !3-LETTER DIGITIZER DEVICE CODE |
---|
232 | |
---|
233 | C--DATA PRESENT FOR EACH PHASE (P OR S) REPORTED FOR THIS EVENT (SUBSCRIPT IM) |
---|
234 | PARAMETER (MMAX=1200) !MAX NO. OF PHASES (P OR S) PER EVENT |
---|
235 | COMMON /D/ M !NO. OF ARRIVAL TIMES FOR THIS EVENT |
---|
236 | COMMON /D/ IND (MMAX) !POINTER TO STATION IN TABLE & S FLAG |
---|
237 | COMMON /D/ A(MMAX,5) !TRAVEL TIMES & DERIVATIVES |
---|
238 | COMMON /D/ R (MMAX) !CURRENT TRAVEL TIME RESIDUAL |
---|
239 | COMMON /D/ W (MMAX) !COMPUTED WEIGHT OF READING |
---|
240 | COMMON /D/ IMPORT (MMAX) !COMPUTED IMPORTANCE OF READING IN .001 |
---|
241 | COMMON /D/ MTCAL (MMAX) !CALC TT IN .01. MULT S TIMES BY POS |
---|
242 | C--THESE VALUES ARE PASSED BY THE ARRAY A. |
---|
243 | C 1 AIN !ANGLE OF EMERGENCE AT THE SOURCE, MEASURED IN DEGREES FROM NADIR. |
---|
244 | C 2 T !TRAVEL TIME |
---|
245 | C 3 DTDR !TT DERIVATIVE WRT DISTANCE |
---|
246 | C 4 DTDZ !TT DERIVATIVE WRT DEPTH |
---|
247 | |
---|
248 | C--HYPOCENTER DATA FOR THIS EVENT |
---|
249 | COMMON /D/ KYEAR !TWO DIGIT YEAR |
---|
250 | COMMON /D/ KYEAR2 !4 DIGIT YEAR |
---|
251 | COMMON /D/ KMONTH !MONTH |
---|
252 | COMMON /D/ KDAY !DAY |
---|
253 | COMMON /D/ KHOUR !HOUR |
---|
254 | COMMON /D/ KMIN !MINUTE |
---|
255 | COMMON /D/ T1 !CURRENT (OR TRIAL) ORIGIN TIME, SEC |
---|
256 | COMMON /D/ Z1 !CURRENT (OR TRIAL) DEPTH |
---|
257 | COMMON /D/ CLAT !CURRENT (OR TRIAL) LAT, DECIMAL DEGREES |
---|
258 | COMMON /D/ CLON !CURRENT (OR TRIAL) LON, DECIMAL DEG, POS WEST |
---|
259 | COMMON /D/ LAT !EVENT LAT, DEG |
---|
260 | COMMON /D/ XLTM !EVENT LAT, MIN |
---|
261 | COMMON /D/ LON !EVENT LON, DEG |
---|
262 | COMMON /D/ XLNM !EVENT LON, MIN |
---|
263 | |
---|
264 | C--ERROR DATA FOR THIS EVENT |
---|
265 | COMMON /D/ ERH !HORIZONTAL LOCATION ERROR |
---|
266 | COMMON /D/ ERZ !DEPTH ERROR |
---|
267 | COMMON /D/ SERR(3) !MAGNITUDES OF 3 STD. ERRORS |
---|
268 | COMMON /D/ IAZ(3) !AZIMUTHS OF 3 STD. ERRORS |
---|
269 | COMMON /D/ IDIP(3) !DIPS OF 3 STD. ERRORS |
---|
270 | |
---|
271 | C--FIRST MAGNITUDE DATA FOR EVENT |
---|
272 | CHARACTER BMTYP*1,BMTYPX*1 |
---|
273 | COMMON /D/ FMAG !MEDIAN DURATION MAG |
---|
274 | COMMON /D/ XMAG !MEDIAN AMPLITUDE MAG |
---|
275 | C COMMON /D/ MFMAG !TOTAL OF DURATION MAG WEIGHTS *100 |
---|
276 | COMMON /D/ NFMAG !NUMBER OF DURATION MAGS W/ WEIGHTS>0 |
---|
277 | C COMMON /D/ MXMAG !TOTAL OF AMPLITUDE MAG WEIGHTS *100 |
---|
278 | COMMON /D/ NXMAG !NUMBER OF AMPLITUDE MAGS W/ WEIGHTS>0 |
---|
279 | COMMON /D/ FMMAD !MEDIAN ABSOL DIFFERENCE OF DURATION MAGS |
---|
280 | COMMON /D/ XMMAD !MEDIAN ABSOL DIFFERENCE OF AMPLITUDE MAGS |
---|
281 | COMMON /D/ BMAG !PASS THRU (BERKELEY) MAGNITUDE |
---|
282 | C COMMON /D/ MBMAG !TOTAL OF PASS THRU MAG WEIGHTS *100 |
---|
283 | COMMON /D/ NBMAG !NUMBER OF PASS THRU MAGS W/ WEIGHTS>0 |
---|
284 | COMMON /C/ BMTYP !PASS THRU MAGNITUDE LABEL CODE |
---|
285 | COMMON /D/ BMAGX !PASS THRU X-MAGNITUDE |
---|
286 | C COMMON /D/ MBMAGX !NUMBER OF PASS THRU XMAG READINGS *100 |
---|
287 | COMMON /D/ NBMAGX !NUMBER OF PASS THRU XMAG READINGS |
---|
288 | COMMON /C/ BMTYPX !PASS THRU X-MAGNITUDE LABEL CODE |
---|
289 | |
---|
290 | C--SECOND MAGNITUDE DATA FOR EVENT |
---|
291 | COMMON /D/ FMAG2 !MEDIAN DURATION MAG |
---|
292 | COMMON /D/ XMAG2 !MEDIAN AMPLITUDE MAG |
---|
293 | C COMMON /D/ MFMAG2 !TOTAL OF DURATION MAG WEIGHTS *100 |
---|
294 | COMMON /D/ NFMAG2 !NUMBER OF DURATION MAGS W/ WEIGHTS>0 |
---|
295 | C COMMON /D/ MXMAG2 !TOTAL OF AMPLITUDE MAG WEIGHTS *100 |
---|
296 | COMMON /D/ NXMAG2 !NUMBER OF AMPLITUDE MAGS W/ WEIGHTS>0 |
---|
297 | COMMON /D/ FMMAD2 !MEDIAN ABSOL DIFFERENCE OF DURATION MAGS |
---|
298 | COMMON /D/ XMMAD2 !MEDIAN ABSOL DIFFERENCE OF AMPLITUDE MAGS |
---|
299 | |
---|
300 | C--MISCELLANEOUS DATA FOR THIS EVENT |
---|
301 | CHARACTER REMK*3,RMK1*1,RMK2*1,FULNAM*25,SHADO*95,SHAD1*95 |
---|
302 | CHARACTER*1 SOUCOD,FMSOU,XMSOU,FMSOU2,XMSOU2 |
---|
303 | COMMON /D/ RMS !RMS RESIDUAL AFTER RESIDUAL WEIGHTING |
---|
304 | COMMON /D/ RMSWT !RMS RESIDUAL BEFORE RESIDUAL WEIGHTING |
---|
305 | COMMON /C/ REMK !3-LETTER REMARK BASED ON LOCATION |
---|
306 | COMMON /C/ RMK1 !FIRST 1-LETTER AUXILIARY REMARK |
---|
307 | COMMON /C/ RMK2 !SECOND 1-LETTER AUXILIARY REMARK |
---|
308 | COMMON /C/ FULNAM !FULL 25-LETTER REGION NAME |
---|
309 | COMMON /D/ DONE !FLAG SET WHEN EVENT IS LOCATED |
---|
310 | COMMON /D/ V(4,4) !EIGENVECTOR MATRIX |
---|
311 | COMMON /D/ EIGVAL(4) !EIGENVALUES |
---|
312 | COMMON /D/ COVAR(4,4) !COVARIANCE MATRIX |
---|
313 | COMMON /D/ IDNO !EVENT ID NO. INPUT FOR EVENT |
---|
314 | COMMON /D/ MEMDSK !0=GET CUSP EVENT FROM MEMORY 1=DISK (LOC) |
---|
315 | COMMON /D/ INUM !EVENT NO. WITHIN THIS LOCATION RUN |
---|
316 | COMMON /D/ LENSHA !LENGTH OF OPTIONAL REC AFTER TERMINATOR CARD |
---|
317 | COMMON /C/ SHADO !OPTIONAL RECORD AFTER TERMINATOR CARD |
---|
318 | |
---|
319 | PARAMETER (MSHA=4) !MAX NUMBER OF SHADOW RECS AFTER HEADER |
---|
320 | COMMON /D/ NSHA1 !ACTUAL NUMBER OF SHADOW RECS AFTER HEADER |
---|
321 | COMMON /D/ LSHA1 (MSHA) !LENGTH OF OPTIONAL REC AFTER ARCHIVE HEADER |
---|
322 | COMMON /C/ SHAD1 (MSHA) !OPTIONAL RECORD AFTER ARCHIVE HEADER |
---|
323 | COMMON /D/ ITR !CURRENT ITERATION NUMBER |
---|
324 | COMMON /D/ DMIN !DISTANCE TO CLOSEST STATION |
---|
325 | COMMON /D/ KEND !PROGRAM TERMINATION FLAG |
---|
326 | COMMON /D/ MAXGAP !MAXIMUM AZIMUTHAL GAP |
---|
327 | COMMON /D/ NFRM !NUMBER OF FIRST MOTIONS |
---|
328 | COMMON /D/ NWR !NUMBER OF WEIGHTED READINGS (P & S) |
---|
329 | COMMON /D/ NWS !NUMBER OF WEIGHTED S READINGS |
---|
330 | COMMON /D/ NVR !NUMBER OF VALID READINGS (>0 ASSIGNED WEIGHTS) |
---|
331 | COMMON /D/ RR !LENGTH OF ADJUSTMENT VECTOR IN KM |
---|
332 | COMMON /C/ SOUCOD !MOST COMMON PHASE SOURCE CODE |
---|
333 | COMMON /C/ FMSOU !MOST COMMON 1ST FMAG SOURCE CODE |
---|
334 | COMMON /C/ XMSOU !MOST COMMON 1ST XMAG SOURCE CODE |
---|
335 | COMMON /C/ FMSOU2 !MOST COMMON 2ND FMAG SOURCE CODE |
---|
336 | COMMON /C/ XMSOU2 !MOST COMMON 2ND XMAG SOURCE CODE |
---|
337 | |
---|
338 | C***************** PARAMETERS WHICH APPLY TO ALL EVENTS ****************** |
---|
339 | C--THE DEFINING COMMANDS ARE GIVEN IN PARENTHESES |
---|
340 | |
---|
341 | CC--STATION DELAY PARAMETERS (DLY) [MADE OBSOLETE BY MULTIPLE MODEL OPTION] |
---|
342 | C COMMON /D/ KDLY !DELAY MODEL CONTROL: 1=MODEL 1, 2=MODEL 2, 3=LINE |
---|
343 | CC SEPARATES 2 MODELS, 4=CIRCLE SEPS. 2 MODELS. |
---|
344 | C COMMON /D/ DLYBAL !FACTOR WHICH BALENCES BETWEEN 2 DELAY MODELS |
---|
345 | CC (0 FOR MOD1, 1 FOR MOD2) |
---|
346 | C COMMON /D/ DLYAZ !AZIMUTH OF LINE OR RADIUS OF CIRCLE WHICH |
---|
347 | CC SEPARATES THE 2 DELAY MODELS |
---|
348 | C COMMON /D/ DLYWD !WIDTH OF TRANSITION ZONE BETWEEN MODELS |
---|
349 | C COMMON /D/ DLYLON !LON OF CIRCLE CENTER OR POINT ON LINE |
---|
350 | C COMMON /D/ DLYLAT !LAT OF CIRCLE CENTER OR POINT ON LINE |
---|
351 | CC BOTH IN DECIMAL DEGREES, POS N & W |
---|
352 | |
---|
353 | C--TERMINATING LOCATION UPON CONVERGENCE (CON) |
---|
354 | COMMON /D/ ITRLIM !MAXIMUM NUMBER OF ITERATIONS |
---|
355 | COMMON /D/ DQUIT !STOP WHEN HYPO ADJUSTMENT IS < THIS |
---|
356 | COMMON /D/ DRQT !STOP WHEN CHANGE IN RMS < THIS |
---|
357 | COMMON /D/ MINSTA !MIN NO. OF PHASE CARDS TO ATTEMPT LOC (MIN) |
---|
358 | |
---|
359 | C--ITERATION & DAMPING CONTROLS (DAM) |
---|
360 | COMMON /D/ DXFIX !KEEP DEPTH FIXED UNTIL ADJUST IS < THIS |
---|
361 | COMMON /D/ DZMAX !DAMP DEPTH ADJUST IF > THIS |
---|
362 | COMMON /D/ DXMAX !DAMP DISTANCE ADJUST IF > THIS |
---|
363 | COMMON /D/ D2FAR !STOP WHEN STA-DIST(2) > THIS |
---|
364 | COMMON /D/ DZAIR !FRACTION OF DEPTH ADJ. IF QUAKE GOES INTO AIR |
---|
365 | COMMON /D/ DAMP !DAMPING FACTOR FOR ALL HYPO ADJUSTMENTS |
---|
366 | COMMON /D/ EIGTOL !MIN EIGENVALUE REQUIRED FOR ADJUSTMENT |
---|
367 | COMMON /D/ RBACK !LIMIT ON RMS INCREASE FOR ONE ITERATION |
---|
368 | COMMON /D/ BACFAC !FRACTION OF STEP TO MOVE BACK IF |
---|
369 | C INCREASE IN RMS > RBACK |
---|
370 | |
---|
371 | C--DURATION MAG CONSTANTS |
---|
372 | LOGICAL LATEN,LCOWT, LNOFMC,LNOXMC |
---|
373 | CHARACTER CLA0*3, CDUG*3 |
---|
374 | C INTEGER*4 ICDATE,IFDATE,IXDATE !LONGER LENGTH REQUIRED |
---|
375 | COMMON /D/ LATEN !TRUE TO CONVERT STA ATTEN TO CAL FACTOR (ATN) |
---|
376 | COMMON /D/ FMA1,FMB1,FMZ1,FMD1,FMF1 !LOW DUR MAG CONSTANTS (DUR) |
---|
377 | COMMON /D/ FMA2,FMB2,FMZ2,FMD2,FMF2 !HIGH DUR MAG CONSTANTS (DUR) |
---|
378 | COMMON /D/ FMBRK !DURATION SEPARATING LO & HI CONSTANTS (DUR) |
---|
379 | COMMON /D/ FMGN !COEFF OF GAIN CORR TO OLD FMAGS (DUR) |
---|
380 | COMMON /D/ IDUG !NUMBER OF GAIN CORR COMPS, -1=ALL (DUG) |
---|
381 | COMMON /C/ CDUG(10) !COMPONENTS TO CORRECT (DUG) |
---|
382 | COMMON /D/ DCOFM1 !COEFF OF ADDL LINEAR DIST TERM (DU2) |
---|
383 | COMMON /D/ DBRKM1 !START DIST OF ADDL LINEAR DIST TERM (DU2) |
---|
384 | COMMON /D/ DCOFM2 !COEFF OF ADDL LINEAR DIST TERM (DU2) |
---|
385 | COMMON /D/ DBRKM2 !START DIST OF ADDL LINEAR DIST TERM (DU2) |
---|
386 | COMMON /D/ ZCOFM !COEFF OF ADDL LINEAR DEPTH TERM (DU2) |
---|
387 | COMMON /D/ ZBRKM !START DEPTH OF ADDL LINEAR DEPTH TERM (DU2) |
---|
388 | COMMON /D/ FMA1B,FMB1B,FMZ1B,FMD1B,FMF1B !LOW 2nd DUR MAG (DUB) |
---|
389 | COMMON /D/ FMA2B,FMB2B,FMZ2B,FMD2B,FMF2B !HIGH 2nd DUR MAG (DUB) |
---|
390 | COMMON /D/ FMBRKB !2ND DURATION SEPARATING LO & HI CONSTS (DUB) |
---|
391 | COMMON /D/ FMGNB !COEFF OF GAIN CORR FOR 2nd FMAGS (DUB) |
---|
392 | COMMON /D/ MAGSEL !FMAG1: 1=CODA MAGS 2=TAU MAGS 3=2nd CODA (MAG) |
---|
393 | COMMON /D/ MAGSL2 !FMAG2: 1=CODA MAGS 2=TAU MAGS 3=2nd CODA (MAG) |
---|
394 | COMMON /D/ LCOWT !T=USE CODA WEIGHTS F=IGNORE CODA WEIGHTS (MAG) |
---|
395 | COMMON /D/ MLOGA0 !1=USE EATON LOGA0 RELATION 2=BAKUN & JOYNER(MAG |
---|
396 | C !3=RICHTER 4=BKY'S NORDQUIST (MAG) |
---|
397 | COMMON /D/ NLA0 !NUMBER OF COMPONENTS W/ UNIQUE LOGA0 REL.(LA0) |
---|
398 | COMMON /C/ CLA0(20) !COMPONENT CODES FOR UNIQUE LOGA0 (LA0) |
---|
399 | COMMON /D/ MLA0(20) !LOGA0 # (LIKE MLOGA0) FOR UNIQUE LOGA0S (LA0) |
---|
400 | COMMON /D/ DMA0 !CONSTANT IN TAU MAG (TAU) |
---|
401 | COMMON /D/ DMA1 !COEFF. OF LOG(TAU) IN TAU MAG (TAU) |
---|
402 | COMMON /D/ DMA2 !COEFF. OF LOG(TAU)**2 IN TAU MAG (TAU) |
---|
403 | COMMON /D/ DMLI !COEFF. OF TAU IN TAU MAG (TAU) |
---|
404 | COMMON /D/ DMZ !COEFF. OF DEPTH IN TAU MAG (TAU) |
---|
405 | COMMON /D/ DMGN !COEFF. OF GAIN/ATTEN IN TAU MAG (TAU) |
---|
406 | |
---|
407 | C--DURATION & AMPLITUDE MAGNITUDE CORRECTIONS |
---|
408 | CHARACTER CFCM*3,CXCM*3 |
---|
409 | COMMON /D/ ICDATE !INITIAL DATE FOR EXTRACTING ATTENS (ATE) |
---|
410 | COMMON /D/ IFDATE !INITIAL DATE FOR EXTRACTING FMCORS (FMC) |
---|
411 | COMMON /D/ IXDATE !INITIAL DATE FOR EXTRACTING XMCORS (XMC) |
---|
412 | COMMON /D/ LNOFMC !T TO USE STA W/ NO CORR, F TO OMIT (FMC) |
---|
413 | COMMON /D/ LNOXMC !T TO USE STA W/ NO CORR, F TO OMIT (XMC) |
---|
414 | COMMON /D/ CALSV(0:13) !CAL FACTORS FOR EACH ATTENUATION SETTING |
---|
415 | COMMON /D/ NFCM !# OF COMPONENTS TO CORRECT FMAGS (FCM) |
---|
416 | COMMON /C/ CFCM(10) !3-LETTER COMPONENTS TO CORRECT (FCM) |
---|
417 | COMMON /D/ AFCM(10) !FMAG CORRECTION FOR COMPONENT (FCM) |
---|
418 | COMMON /D/ NXCM !# OF COMPONENTS TO CORRECT XMAGS (XCM) |
---|
419 | COMMON /C/ CXCM(10) !3-LETTER COMPONENTS TO CORRECT (XCM) |
---|
420 | COMMON /D/ AXCM(10) !XMAG CORRECTION FOR COMPONENT (XCM) |
---|
421 | |
---|
422 | C--SELECTING WHICH COMPONENTS TO USE FOR THE 4 MAGNITUDES |
---|
423 | C--ALSO SELECT WHICH INSTRUMENT TYPES TO USE WITH THE 2 AMP MAGS |
---|
424 | CHARACTER*3 COMPF1,COMPF2,COMPX1,COMPX2 |
---|
425 | CHARACTER*1 LABF1,LABF2,LABX1,LABX2 |
---|
426 | COMMON /D/ NCPF1 !NUMBER OF COMPONENTS TO SELECT FOR FMAG1 (FC1) |
---|
427 | COMMON /C/ COMPF1(20) !COMPONENTS TO USE FOR FMAG1 (FC1) |
---|
428 | COMMON /C/ LABF1 !LABEL CODE TO IDENTIFY FMAG1 (FC1) |
---|
429 | COMMON /D/ NCPF2 !NUMBER OF COMPONENTS TO SELECT FOR FMAG2 (FC2) |
---|
430 | COMMON /C/ COMPF2(20) !COMPONENTS TO USE FOR FMAG2 (FC2) |
---|
431 | COMMON /C/ LABF2 !LABEL CODE TO IDENTIFY FMAG2 (FC2) |
---|
432 | COMMON /D/ NCPX1 !NUMBER OF COMPONENTS TO SELECT FOR XMAG1 (XC1) |
---|
433 | COMMON /C/ COMPX1(20) !COMPONENTS TO USE FOR XMAG1 (XC1) |
---|
434 | COMMON /C/ LABX1 !LABEL CODE TO IDENTIFY XMAG1 (XC1) |
---|
435 | COMMON /D/ NCPX2 !NUMBER OF COMPONENTS TO SELECT FOR XMAG2 (XC2) |
---|
436 | COMMON /C/ COMPX2(20) !COMPONENTS TO USE FOR XMAG2 (XC2) |
---|
437 | COMMON /C/ LABX2 !LABEL CODE TO IDENTIFY XMAG2 (XC2) |
---|
438 | COMMON /D/ MAG1TYPX !TYPE OF MAG FOR XMAG1 1=ML 2=MX (XMT) |
---|
439 | COMMON /D/ MAG2TYPX !TYPE OF MAG FOR XMAG2 1=ML 2=MX (XMT) |
---|
440 | |
---|
441 | C--SELECTING WHICH INSTRUMENT TYPES TO USE FOR THE 2 AMP MAGS |
---|
442 | LOGICAL LXCH |
---|
443 | COMMON /D/ LXCH !T FOR COMPONENT CHOICE, F INST TYPE CHOICE(XCH) |
---|
444 | COMMON /D/ NXTYP1 !NUMBER OF INST TYPES FOR XMAG1 (XTY) |
---|
445 | COMMON /D/ IXTYP1(3) !INST TYPES FOR XMAG1 (XTY) |
---|
446 | COMMON /D/ NXTYP2 !NUMBER OF INST TYPES FOR XMAG2 (XTY) |
---|
447 | COMMON /D/ IXTYP2(3) !INST TYPES FOR XMAG2 (XTY) |
---|
448 | |
---|
449 | C--PREFERRED EVENT MAGNITUDES |
---|
450 | C--MAG NUMBER TAGS: 1=FMAG 2=XMAG 3=PASS-THRU-MAG 4=XMAG2 5=FMAG2 |
---|
451 | CHARACTER LABPR*1 |
---|
452 | COMMON /D/ NMAGS !NUMBER OF PREFERRED MAGNITUDES (PRE) |
---|
453 | COMMON /D/ MPREF(10) !PREFERENCE ORDER OF MAG NUMBER TAGS (PRE) |
---|
454 | COMMON /D/ MNPREF(10) !MINIMUM NUMBER OF READINGS TO USE (PRE) |
---|
455 | COMMON /D/ AMPREF(10) !MIN MAG TO SELECT PREFERRED MAG (PRE) |
---|
456 | COMMON /D/ AXPREF(10) !MAX MAG TO SELECT PREFERRED MAG (PRE) |
---|
457 | COMMON /D/ PMAG !PREFERRED MAG |
---|
458 | C COMMON /D/ MPMAG !TOTAL WEIGHTS OF PREFERRED MAG *100 |
---|
459 | COMMON /D/ NPMAG !NUMBER OF PREFERRED MAGS W/ WEIGHTS>0 |
---|
460 | COMMON /C/ LABPR !LABEL CODE OF PREFERRED MAG |
---|
461 | COMMON /D/ PMMAD !MEDIAN-ABS-DIFF OF PREFERRED MAG |
---|
462 | |
---|
463 | C--MISCELLANEOUS PARAMETERS |
---|
464 | LOGICAL LCOMP1, LKEEP |
---|
465 | COMMON /D/ ITRDIS !ITERATION TO BEGIN DISTANCE WEIGHTING (DIS) |
---|
466 | COMMON /D/ DISCUT,DISW1,DISW2 !DIST. WEIGHT PARAMETERS (DIS) |
---|
467 | COMMON /D/ ITRDI1,DISCU1,DISW11,DISW21 !DIST. WEIGHT PARAMETERS (DI1) |
---|
468 | COMMON /D/ ITRRES !ITERATION TO BEGIN RESIDUAL WEIGHTING (RMS) |
---|
469 | COMMON /D/ RMSCUT,RMSW1,RMSW2 !RMS WEIGHTING PARAMS. (RMS) |
---|
470 | COMMON /D/ LJUNK !T TO KILL DIS & RES WEIGHTING WHEN M<4 (JUN) |
---|
471 | COMMON /D/ SWT !MULTIPLY ALL S WEIGHTS BY THIS (SWT) |
---|
472 | COMMON /D/ POS !VP/VS VELOCITY RATIO (POS) |
---|
473 | COMMON /D/ ZTR !TRIAL DEPTH (ZTR) |
---|
474 | COMMON /D/ RDERR !ESTIMATED READING ERROR IN SECS (ERR) |
---|
475 | COMMON /D/ ERCOF !WEIGHTING FACTOR OF RMS IN ERROR CALCS (ERC) |
---|
476 | COMMON /D/ NET !NET FOR ASSIGNING 3-LET. NAMES BASED ON LOCATION (NET) |
---|
477 | COMMON /D/ NSTLET !NUMBER OF STATION SITE LETTERS TO TEST (LET) |
---|
478 | COMMON /D/ NETLET !NUMBER OF STATION NET LETTERS TO TEST (LET) |
---|
479 | COMMON /D/ NCOMP !NUMBER OF COMPONENT LETTERS TO TEST (LET) |
---|
480 | COMMON /D/ NSLOC !# OF LOCATION CODE LETS TO TEST IN PHASE FILE (LET) |
---|
481 | COMMON /D/ NSLOC2 !# OF LOCATION CODE LETS TO TEST IN OTHER FILES (LET) |
---|
482 | COMMON /D/ LCOMP1 !T= 1 LETTER COMPONENT, F=OTHERWISE (LES) |
---|
483 | COMMON /D/ LKEEP !T=WRITE UNKNOWN STAS TO ARC FILE (KEP) |
---|
484 | COMMON /D/ WTVALS(4) !NUMERICAL WEIGHTS FOR CODES 0-3 (WET) |
---|
485 | |
---|
486 | C--PHASE & SHADOW DATA FOR UNKNOWN STATIONS. COMM BETWEEN HYPHS AND HYLST. |
---|
487 | PARAMETER (MAXUNK=40) |
---|
488 | CHARACTER PUNK*132,SUNK*132,LUNK*5 |
---|
489 | COMMON /UNK/ NLUNK !NUMBER OF UNK STAS IN "NO COMPLAIN" LIST (UNK) |
---|
490 | COMMON /UNC/ LUNK(10) !LIST OF STAS TO NOT COMPLAIN WHEN UNKNOWN (UNK) |
---|
491 | COMMON /UNK/ NUNK !NUMBER OF UNKNOWN STATIONS THIS EVENT |
---|
492 | COMMON /UNC/ PUNK(MAXUNK) !PHASE RECORDS OF UNKNOWN STATIONS |
---|
493 | COMMON /UNK/ NPUNK(MAXUNK) !LENGTH OF PHASE RECORDS |
---|
494 | COMMON /UNC/ SUNK(MAXUNK) !SHADOW RECORDS OF UNKNOWN STATIONS |
---|
495 | COMMON /UNK/ NSUNK(MAXUNK) !LENGTH OF SHADOW RECORDS |
---|
496 | |
---|
497 | C--INTERACTIVE PROCESSING OF SINGLE-EVENT FILES (PRO COMMAND) |
---|
498 | CHARACTER LSTFIL*40,LSTFOR*40 |
---|
499 | CHARACTER*8 EXTPHS,EXTARC,EXTSUM,EXTPRT |
---|
500 | COMMON /C/ LSTFIL !FILE FOR READING BASE EVENT NAMES (BAS) |
---|
501 | COMMON /D/ NCBASE !NUMBER OF CHARACTERS IN BASE NAME (BAS) |
---|
502 | COMMON /C/ LSTFOR !FORMAT FOR READING EVENT FILE (BAS) |
---|
503 | COMMON /C/ EXTPHS !PHASE FILENAME EXTENSION (BAS) |
---|
504 | COMMON /C/ EXTARC !ARCHIVE FILENAME EXTENSION (BAS) |
---|
505 | COMMON /C/ EXTSUM !SUMMARY FILENAME EXTENSION (BAS) |
---|
506 | COMMON /C/ EXTPRT !PRINTOUT FILENAME EXTENSION (BAS) |
---|
507 | COMMON /D/ IEDFLG !EDITOR FLAG 1=EDT 2=ANY ED (BAS) |
---|
508 | |
---|
509 | C--TRANSLATION ARRAYS FROM 3-LETTER DIGITIZERS TO 1-LETTER DATA SOURCE CODES |
---|
510 | CHARACTER DIG3*3, DIG1*1, DIGDEF*1 |
---|
511 | PARAMETER (MAXDIG=50) |
---|
512 | COMMON /C/ DIG3(MAXDIG) !3-LETTER DIGITIZER CODES IN MEM FILE (DIG) |
---|
513 | COMMON /C/ DIG1(MAXDIG) !CORRESPONDING 1-LETTER DATA SOURCE CODES (DIG) |
---|
514 | COMMON /C/ DIGDEF !DEFAULT DATA SOURCE CODE (DID) |
---|
515 | COMMON /D/ NDIG !NUMBER OF DATA SOURCE CODES (DID) |
---|
516 | |
---|
517 | C--P AMPLITUDE MAGNITUDE PARAMETERS |
---|
518 | LOGICAL LPMAG,LPPRT |
---|
519 | CHARACTER CPWM*3,PSOUR*1,PSOUR2*1, COMPP1*3,COMPP2*3, CCNTMM*1 |
---|
520 | CHARACTER LABP1*1,LABP2*1 |
---|
521 | COMMON /D/ LPMAG !T= COMPUTE PMAG, MUST HAVE SHADOW INPUT (PMA) |
---|
522 | COMMON /D/ LPPRT !T= ADD PMAG STA. LINES TO PRINT OUTPUT (PMA) |
---|
523 | COMMON /D/ CNT2MD !DEVEL. MM PER COUNT (P AMPS) RTP=.04 (PMA) |
---|
524 | COMMON /D/ CLPRAT !RATIO OF PMCLIP/PMUSED FOR CLIPPED MAG (PMA) |
---|
525 | COMMON /D/ LATYPP !LOGA0 FORMULA # FOR PMAGS (PMA) |
---|
526 | COMMON /D/ NPWM !NUMBER OF COMPONENTS WITH PAMAG WEIGHTS (PAC) |
---|
527 | COMMON /C/ CPWM(10) !COMPONENTS WITH DEFINED WEIGHTS (PAC) |
---|
528 | COMMON /D/ WPWM(10) !WEIGHTS FOR EACH COMPONENT (PAC) |
---|
529 | |
---|
530 | C--P AMP MAGNITUDE COMPONENT AND DATA SOURCE DEPENDENT PARAMETERS |
---|
531 | COMMON /C/ LABP1 !LABEL OF PRIMARY P AMP MAGNITUDE (PC1) |
---|
532 | COMMON /D/ PMA1 !COEFFICIENT IN PMAG = PMA1 +PMB1*LOGA (PC1) |
---|
533 | COMMON /D/ PMB1 !COEFFICIENT IN PMAG = PMA1 +PMB1*LOGA (PC1) |
---|
534 | COMMON /D/ NCPP1 !NUMBER OF COMPONENTS TO USE FOR PAMAG (PC1) |
---|
535 | COMMON /C/ COMPP1(10) !COMPONENT CODES FOR PRIMARY P MAG (PC1) |
---|
536 | COMMON /C/ LABP2 !LABEL OF SECONDARY P AMP MAGNITUDE (PC2) |
---|
537 | COMMON /D/ PMA2 !COEFFICIENT IN PMAG = PMA2 +PMB2*LOGA (PC1) |
---|
538 | COMMON /D/ PMB2 !COEFFICIENT IN PMAG = PMA2 +PMB2*LOGA (PC1) |
---|
539 | COMMON /D/ NCPP2 !NUMBER OF COMPONENTS TO USE FOR PAMAG2 (PC2) |
---|
540 | COMMON /C/ COMPP2(10) !COMPONENT CODES FOR SECONDARY P MAG (PC2) |
---|
541 | COMMON /D/ NCNTMM !NUMBER OF SOURCE SPECIFIC CONV FACTORS (PMC) |
---|
542 | COMMON /C/ CCNTMM(10) !DATA SOURCE CODES (PMC) |
---|
543 | COMMON /D/ CNT2MM(10) !COUNTS TO MM CONVERSION FOR SOURCE (PMC) |
---|
544 | |
---|
545 | C--PRIMARY P AMPLITUDE MAGNITUDE |
---|
546 | COMMON /D/ PAMAG !PRIMARY P AMPLITUDE MAGNITUDE FOR EVENT |
---|
547 | COMMON /D/ PAMAD !MEDIAN ABSOLUTE DIFFERENCE OF PMAG |
---|
548 | COMMON /D/ PMUSED !TOTAL OF WEIGHTS OF STATIONS USED IN PMAG |
---|
549 | COMMON /D/ PMCLIP !TOT OF WEIGHTS OF CLIPPED STAS USED IN PMAG |
---|
550 | COMMON /D/ MINPM !PMAG CLIP FLAG =1 OF PMAG IS A MIN, =0 OK |
---|
551 | COMMON /D/ PNORMN !L1 NORM FOR EVENT P MAG |
---|
552 | COMMON /C/ PSOUR !PRINCIPAL DATA SOURCE FOR P MAGS |
---|
553 | |
---|
554 | C--SECONDARY P AMPLITUDE MAGNITUDE (IE LOW GAIN) |
---|
555 | COMMON /D/ PAMAG2 !SECONDARY P AMPLITUDE MAGNITUDE FOR EVENT |
---|
556 | COMMON /D/ PAMAD2 !MEDIAN ABSOLUTE DIFFERENCE OF PMAG |
---|
557 | COMMON /D/ PMUSD2 !TOTAL OF WEIGHTS OF STATIONS USED IN PMAG |
---|
558 | COMMON /D/ PMCLP2 !TOT OF WEIGHTS OF CLIPPED STAS USED IN PMAG |
---|
559 | COMMON /D/ MINPM2 !PMAG CLIP FLAG =1 OF PMAG IS A MIN, =0 OK |
---|
560 | COMMON /D/ PNRMN2 !L1 NORM FOR EVENT P MAG |
---|
561 | COMMON /C/ PSOUR2 !PRINCIPAL DATA SOURCE FOR P MAGS |
---|