Changeset 3172
- Timestamp:
- 12/14/07 14:17:22 (13 years ago)
- Location:
- trunk/src/seismic_processing/hyp2000
- Files:
-
- 1 added
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/seismic_processing/hyp2000/common.inc
r2768 r3172 2 2 LOGICAL LSUM,LARC,LPRT,LERR,LEJCT,LREP,DONE,LCUSP,LPRALL 3 3 LOGICAL LAPP,LJUNK,LMULT,LTBIG,LMAG,LBSTA,LBCRU,SUBMOD,LP153 4 LOGICAL L2000,HYPOFIX,ALLFIX 4 LOGICAL L2000,HYPOFIX,ALLFIX, LASK 5 5 6 6 C--SHORT INTEGER AND LOGICAL DECLARATIONS … … 43 43 CHARACTER*60 ATNFIL,FMCFIL,XMCFIL,MAGFIL,BSTAFL,BCRUFL,TERMIN 44 44 CHARACTER*60 CALFIL 45 CHARACTER FORID*30 45 CHARACTER FORID*30, GREETING*70 46 46 COMMON /C/ SUMFIL !HYPOINVERSE SUMMARY OUTPUT FILENAME (SUM) 47 47 COMMON /C/ ARCFIL !ARCHIVE OUTPUT FILENAME (ARC) … … 88 88 COMMON /D/ ICENT !DEFAULT CENTURY OLD PHASE FORMAT 89 89 COMMON /D/ IAMPU !DEFAULT AMP UNITS CODE 90 COMMON /C/ GREETING !GREETING MESSAGE FOR TERMINAL & PRT FILE 90 91 91 92 C--HYPOINVERSE CRUSTAL MODELS ******************************************** … … 225 226 COMMON /C/ KRMK (MAXPHS) !ANY 1 LETTER STATION REMARK 226 227 COMMON /C/ KRMK6 (MAXPHS) !ORIGINAL 6-LETTER REMARK 227 C COMMON /C/ KPSWT (MAXPHS) !STATION WEIGHT-OUT CODE NO LONGER USED228 228 COMMON /C/ KSHAD (MAXPHS) !OPTIONAL SHADOW PHASE RECORD 229 229 COMMON /D/ KLSHA (MAXPHS) !LENGTH OF KSHAD … … 299 299 300 300 C--MISCELLANEOUS DATA FOR THIS EVENT 301 CHARACTER REMK*3,RMK1*1,RMK2*1,FULNAM*25,SHADO*9 2,SHAD1*92301 CHARACTER REMK*3,RMK1*1,RMK2*1,FULNAM*25,SHADO*95,SHAD1*95 302 302 CHARACTER*1 SOUCOD,FMSOU,XMSOU,FMSOU2,XMSOU2 303 303 COMMON /D/ RMS !RMS RESIDUAL AFTER RESIDUAL WEIGHTING -
trunk/src/seismic_processing/hyp2000/hybda.for
r2487 r3172 4 4 5 5 C--GENERAL PARAMETERS 6 DATA GREETING 7 2 /'12/2007 VERSION 1.2 (Can fix origin time; g77 compiler)'/ 6 8 DATA PI,RDEG,LJUNK,LMULT /3.14159,57.2958,2*.FALSE./ 7 9 DATA NSTLET,NETLET,NCOMP,NSLOC,NSLOC2 /4,4*0/ -
trunk/src/seismic_processing/hyp2000/hycmd.for
r2768 r3172 3 3 INCLUDE 'common.inc' 4 4 C--LASK IS A LOGICAL FUNCTION. THE OS2 COMPILER COMPLAINS WITHOUT THESE LINES 5 LOGICAL LASK 6 EXTERNAL LASK 5 C LOGICAL LASK 6 C EXTERNAL LASK 7 7 8 CHARACTER TEMPSTR*80 8 9 PARAMETER (NCMD=92) -
trunk/src/seismic_processing/hyp2000/hyp.for
r2768 r3172 1 C--HYP IS THE VAX VERSION OF THE LOCATION PROGRAM HYPOINVERSE. THE PROGRAM1 C--HYP IS THE VAX/SUN VERSION OF THE LOCATION PROGRAM HYPOINVERSE. THE PROGRAM 2 2 C IS DESIGNED TO BE FAST, FILE ORIENTED AND COMMAND DRIVEN. 3 3 4 C--WRITTEN BY FRED KLEIN , 10/2006 VERSION4 C--WRITTEN BY FRED KLEIN. SEE GREETING IN HYBDA.FOR FOR VERSION # 5 5 6 6 C++++++++++++++++ LIST OF REQUIRED SUBROUTINES +++++++++++++++++ … … 135 135 C MAXSTA > MAXPHS > MMAX. THESE CAN BE MADE SMALLER FOR SMALL NETWORKS. 136 136 137 C--VERY PARTIAL VERSION HISTORY 138 C 139 C 1978 VERSION 0.x (Eclipse computer, simple & compact) 140 C 1985 VERSION 0.x (VAX and pro-350 computer, compact) 141 C 1989 VERSION 0.x (multiple crust models) 142 C 2002 VERSION 1.0 (Full documentation, completely Y2000 capable) 143 C 2/2007 VERSION 1.1 (Can fix origin time, many other changes) 144 C 5/2007 VERSION 1.11 (Can fix origin time, greeting) 145 C 12/2007 VERSION 1.2 (Can fix origin time; g77 compiler) 146 137 147 C++++++++++++++++ I/O DEVICE NUMBERS USED ++++++++++++++++++++++ 138 148 C 5 TERMINAL INPUT. … … 162 172 C--SEND A MESSAGE TO THE TERMINAL. THIS ALSO ASSIGNS UNIT 5 TO TERMINAL 163 173 WRITE (6,1000) 164 1000 FORMAT (' HYPOINVERSE 2000 STARTING' /165 2' 2/2007 VERSION 1.1 (Can fix origin time)')174 1000 FORMAT (' HYPOINVERSE 2000 STARTING') 175 WRITE (6,'(A)') GREETING 166 176 167 177 C--INITIALIZE VARIABLES NOT INITIALIZED IN BLOCK DATA -
trunk/src/seismic_processing/hyp2000/hypro.for
r2020 r3172 6 6 CHARACTER PRVSTA*5, PRVNET*2, PRVCMP*3, PRVLOC*2 7 7 INCLUDE 'common.inc' 8 LOGICAL FOUNDIT 8 9 C--LASK IS A LOGICAL FUNCTION. THE OS2 COMPILER COMPLAINS WITHOUT THESE LINES 9 LOGICAL LASK, FOUNDIT 10 EXTERNAL LASK10 c LOGICAL LASK 11 c EXTERNAL LASK 11 12 12 13 BASE=' ' -
trunk/src/seismic_processing/hyp2000/hystl.for
r2768 r3172 12 12 C--PRINT CURRENT DATE & TIME FIRST 13 13 CALL HYTIME (CURTIM) 14 WRITE (15,1000) CURTIM,RUNLAB, CDOMAN,CPVERS15 1000 FORMAT (' HYPOINVERSE 2000 (10/2006 VERSION)RUN ON ',A28,14 WRITE (15,1000) GREETING, CURTIM,RUNLAB, CDOMAN,CPVERS 15 1000 FORMAT (' HYPOINVERSE 2000'/1X,A/' RUN ON ',A28, 16 16 2 ' RUN LABEL=',A1/' PROCESSING DOMAIN=',A2, 17 17 3 ' PROCESSING VERSION=',A2/) -
trunk/src/seismic_processing/hyp2000/makefile.ux
r2902 r3172 1 1 2 # 3 # THIS FILE IS UNDER RCS - DO NOT MODIFY UNLESS YOU HAVE 4 # CHECKED IT OUT USING THE COMMAND CHECKOUT. 5 # 6 # $Id$ 7 # 8 # Revision history: 9 # $Log$ 10 # Revision 1.5 2007/03/29 18:09:58 paulf 11 # patched makefile.ux to use hyedit.f and hydelt.f instead of .for 12 # 13 # Revision 1.4 2007/03/29 18:01:50 paulf 14 # removed .fpp files as per update from Fred Klein 15 # 16 # Revision 1.3 2007/03/02 18:39:43 dietz 17 # removed make entries relating to test versions of hypoinverse 18 # 19 # Revision 1.2 2006/06/06 18:37:46 paulf 20 # clean compile on Linux, but still problems with run 21 # 22 # Revision 1.1 2005/07/20 00:53:48 friberg 23 # added in makefile.ux 24 # 25 # Revision 1.2 2000/09/12 22:15:59 dietz 26 # *** empty log message *** 27 # 28 # Revision 1.1 2000/09/12 20:53:03 dietz 29 # Initial revision 30 # 31 # Revision 1.1 2000/02/14 18:41:39 lucky 32 # Initial revision 33 # 34 # 35 # 36 37 # HYP HYPOINVERSE EARTHQUAKE LOCATION PROGRAM. 2 3 # some earthworm vars 38 4 39 5 B = $(EW_HOME)/$(EW_VERSION)/bin … … 41 7 42 8 43 hyp2000 : hyp.o hybda.o hycmd.o hysta.o hydel.o hyate.o hyfmc.o hycal.o \ 9 # HYP HYPOINVERSE EARTHQUAKE LOCATION PROGRAM. 10 hypg : hyp.o hybda.o hycmd.o hysta.o hydel.o hyate.o hyfmc.o hycal.o \ 44 11 hyxmc.o hycrh.o hycrt.o hystl.o hyopen.o hyinit.o hyphs.o hysou.o \ 45 12 hycin.o hycout.o hytrl.o hyloc.o hysol.o hysvd.o hytra.o hytrh.o hytrt.o \ 46 13 hymag.o hyrep.o hylst.o hysum.o hyinp.o hypro.o medwt.o utmcal.o \ 47 init_event.o hydelt.o hyedit.o hytime.o hybeg.o hypref.o hymagp.o hyfile.o \ 48 readq.o jdate.o spawn.o openr.o askr.o openw.o errset.o leng.o dayjl.o \ 49 lask.o jask.o klas.o upstr.o askc.o box2.o box3.o klass.o ksic.o 50 f77 hyp.o hybda.o hycmd.o hysta.o hydel.o hyate.o hyfmc.o hypref.o \ 14 init_event.o hydelt.o hyedit.o hytime.o hybeg.o hypref.o hymagp.o \ 15 hyfile.o allsubs.o 16 g77 -g hyp.o hybda.o hycmd.o hysta.o hydel.o hyate.o hyfmc.o hypref.o \ 51 17 hyxmc.o hycrh.o hycrt.o hystl.o hyopen.o hyinit.o hyphs.o hycal.o \ 52 18 hycin.o hytrl.o hyloc.o hysol.o hysvd.o hytra.o hytrh.o hytrt.o \ 53 19 hymag.o hyrep.o hylst.o hysum.o hyinp.o hypro.o medwt.o utmcal.o \ 54 init_event.o hydelt.o hyedit.o hytime.o hybeg.o hycout.o hymagp.o hyfile.o \ 55 hysou.o readq.o jdate.o spawn.o openr.o askr.o openw.o errset.o \ 56 leng.o dayjl.o lask.o jask.o klas.o upstr.o askc.o box2.o box3.o \ 57 klass.o ksic.o \ 58 -Bstatic -o $B/hyp2000 20 init_event.o hydelt.o hyedit.o hytime.o hybeg.o hycout.o \ 21 hymagp.o hysou.o hyfile.o \ 22 -Bstatic allsubs.o -o $B/hyp2000 23 24 # line for loading all "common" subroutines instead of library: 25 # -Bstatic allsubs.o -o hypg 26 # 27 # line for loading all "common" subroutines in a library using f77: 28 # -Bstatic -L/home/calnet/klein/subs -lfs -o hyp2000 29 30 # allsubs - collection of general subroutines into one file 31 # This is the g77 version of readq subroutine (no q format) 32 # Sun version (readq.f) built into 33 allsubs.o : allsubs.f 34 g77 -c allsubs.f 59 35 60 36 # HYP * MAIN HYPOINVERSE PROGRAM. 61 37 hyp.o : hyp.for common.inc 62 f77 -c hyp.for 63 64 # HYPOINV * MAIN HYPOINVERSE SUBROUTINE. 65 hypoinv.o : hypoinv.for common.inc 66 f77 -c hypoinv.for 38 g77 -c hyp.for 67 39 68 40 # HYBDA * BLOCK DATA INITIALIZATION OF COMMON. 69 41 hybda.o : hybda.for common.inc 70 f77 -c hybda.for42 g77 -c hybda.for 71 43 72 44 # HYBEG * INITIALIZATION OF OTHER VARIABLES. 73 45 hybeg.o : hybeg.f common.inc 74 f77 -c hybeg.f46 g77 -c hybeg.f 75 47 76 48 # HYCMD * GETS AND PROCESSES COMMANDS. 77 49 hycmd.o : hycmd.for common.inc 78 f77 -c hycmd.for 79 80 # HYFILE * DETERMINE THE FILE TYPE AND FORMAT. 81 hyfile.o : hyfile.for common.inc 82 f77 -c hyfile.for 50 g77 -c hycmd.for 83 51 84 52 # HYSTA * READS IN STATIONS. 85 53 hysta.o : hysta.for common.inc 86 f77 -c hysta.for54 g77 -c hysta.for 87 55 88 56 # HYDEL * READS IN STATION DELAYS (FOR MULTIPLE MODELS). 89 57 hydel.o : hydel.for common.inc 90 f77 -c hydel.for58 g77 -c hydel.for 91 59 92 60 # HYATE * READS IN STATION ATTENUATION HISTORY. 93 61 hyate.o : hyate.for common.inc 94 f77 -c hyate.for62 g77 -c hyate.for 95 63 96 64 # HYCAL * READS IN STATION CAL FACTOR HISTORY. 97 65 hycal.o : hycal.for common.inc 98 f77 -c hycal.for66 g77 -c hycal.for 99 67 100 68 # HYFMC * READS IN STATION FMAG CORRECTIONS. 101 69 hyfmc.o : hyfmc.for common.inc 102 f77 -c hyfmc.for 70 g77 -c hyfmc.for 71 72 # HYFILE * DETERMINES INPUT FILE TYPE. 73 hyfile.o : hyfile.for 74 g77 -c hyfile.for 103 75 104 76 # HYXMC * READS IN STATION XMAG CORRECTIONS. 105 77 hyxmc.o : hyxmc.for common.inc 106 f77 -c hyxmc.for78 g77 -c hyxmc.for 107 79 108 80 # HYCRH * READS IN HOMOGENOUS LAYER CRUSTAL MODELS. 109 81 hycrh.o : hycrh.for common.inc 110 f77 -c hycrh.for82 g77 -c hycrh.for 111 83 112 84 # HYCRT * READS IN TRAVEL-TIME-TABLE CRUSTAL MODELS. 113 85 hycrt.o : hycrt.for common.inc 114 f77 -c hycrt.for86 g77 -c hycrt.for 115 87 116 88 # HYSTL * OUTPUTS STATIONS, CRUST & PARAMETERS TO PRINT FILE. 117 89 hystl.o : hystl.for common.inc 118 f77 -c hystl.for90 g77 -c hystl.for 119 91 120 92 # HYOPEN * OPENS FILES FOR LOCATION RUN. 121 93 hyopen.o : hyopen.for common.inc 122 f77 -c hyopen.for94 g77 -c hyopen.for 123 95 124 96 # HYINIT * INITIALIZES SOME VARIABLES FOR LOCATION RUN. 125 97 hyinit.o : hyinit.for common.inc 126 f77 -c hyinit.for98 g77 -c hyinit.for 127 99 128 100 # HYPHS * READS IN PHASE DATA FOR ONE EVENT. 129 101 hyphs.o : hyphs.for common.inc 130 f77 -c hyphs.for102 g77 -c hyphs.for 131 103 132 104 # HYCIN * INPUTS PHASE DATA FROM CUSP MEM FILES (ALTERNATE TO HYPHS) 133 105 hycin.o : hycin.f 134 f77 -c hycin.f106 g77 -c hycin.f 135 107 136 108 # HYCOUT * OUTPUTS PHASE DATA TO CUSP MEM FILES (ALTERNATE TO HYARC) 137 109 hycout.o : hycout.f 138 f77 -c hycout.f110 g77 -c hycout.f 139 111 140 112 # HYTRL * SETS TRIAL HYPOCENTER. 141 113 hytrl.o : hytrl.for common.inc 142 f77 -c hytrl.for114 g77 -c hytrl.for 143 115 144 116 # HYLOC * LOCATES ONE EVENT. 145 117 hyloc.o : hyloc.for common.inc 146 f77 -c hyloc.for118 g77 -c hyloc.for 147 119 148 120 # HYSOL * PERFORMS INVERSION FOR ONE ITERATION. 149 121 hysol.o : hysol.for common.inc 150 f77 -c hysol.for122 g77 -c hysol.for 151 123 152 124 # HYSVD CANNED SINGLE-VALUE-DECOMPOSITION ROUTINE. 153 125 hysvd.o : hysvd.f 154 f77 -c hysvd.f126 g77 -c hysvd.f 155 127 156 128 # HYTRA * MANAGE CRUST MODEL CHOICE & AVERAGING. 157 129 hytra.o : hytra.for common.inc 158 f77 -c hytra.for130 g77 -c hytra.for 159 131 160 132 # HYTRH * CALC TRAVEL TIMES AND DERIVS FOR HOMO LAYER MODEL. 161 133 hytrh.o : hytrh.for common.inc 162 f77 -c hytrh.for134 g77 -c hytrh.for 163 135 164 136 # HYTRT * CALC TRAVEL TIMES AND DERIVS FROM TRAV-TIME TABLE. 165 137 hytrt.o : hytrt.for common.inc 166 f77 -c hytrt.for138 g77 -c hytrt.for 167 139 168 140 # HYMAG * COMPUTES DUR & MAX AMP MAGNITUDE. 169 141 hymag.o : hymag.for common.inc 170 f77 -c hymag.for142 g77 -c hymag.for 171 143 172 144 # HYMAGP * COMPUTES P AMP MAGNITUDE. 173 145 hymagp.o : hymagp.for common.inc 174 f77 -c hymagp.for146 g77 -c hymagp.for 175 147 176 148 # HYPREF * SELECTS PREFERRED MAGNITUDE. 177 149 hypref.o : hypref.for common.inc 178 f77 -c hypref.for150 g77 -c hypref.for 179 151 180 152 # HYREP * REPORTS A LOCATION ON THE TERMINAL. 181 153 hyrep.o : hyrep.for common.inc 182 f77 -c hyrep.for154 g77 -c hyrep.for 183 155 184 156 # HYSOU * TABULATES MOST COMMON DATA SOURCES. 185 157 hysou.o : hysou.for common.inc 186 f77 -c hysou.for158 g77 -c hysou.for 187 159 188 160 # HYLST * OUTPUTS DATA BY STATION TO PRINT & ARCHIVE FILES. 189 161 hylst.o : hylst.for common.inc 190 f77 -c hylst.for162 g77 -c hylst.for 191 163 192 164 # HYSUM * OUTPUTS SUMMARY RECORD (FINAL LOCATION). 193 165 hysum.o : hysum.for common.inc 194 f77 -c hysum.for166 g77 -c hysum.for 195 167 196 168 # HYINP FOR INTERACTIVE ENTRY OF PHASE DATA. 197 169 hyinp.o : hyinp.for 198 f77 -c hyinp.for170 g77 -c hyinp.for 199 171 200 172 # HYPRO * INTERACTIVELY PROCESSES A SERIES OF EVENTS. 201 173 hypro.o : hypro.for common.inc 202 f77 -c hypro.for174 g77 -c hypro.for 203 175 204 176 # MEDWT COMPUTES THE WEIGHTED MEDIAN OF A SERIES OF MAGNITUDES. 205 177 medwt.o : medwt.for 206 f77 -c medwt.for178 g77 -c medwt.for 207 179 208 180 # UTMCAL CALCULATE STATION DISTANCES ON A UTM GRID 209 181 utmcal.o : utmcal.for 210 f77 -c utmcal.for182 g77 -c utmcal.for 211 183 212 184 # INIT_EVENT INITIALIZES CUSP FOR READING OF MEM DATA FILES (VAX ONLY) 213 185 init_event.o : init_event.f 214 f77 -c init_event.f186 g77 -c init_event.f 215 187 216 188 # HYDELETE DELETES FILES IN INTERACTIVE PROCESSING 217 189 hydelt.o : hydelt.f 218 f77 -c hydelt.f190 g77 -c hydelt.f 219 191 220 192 # HYEDTIOR RUNS AN EDTIOR WITHIN A PROCESS 221 193 hyedit.o : hyedit.f 222 f77 -c hyedit.f194 g77 -c hyedit.f 223 195 224 196 # HYTIME GETS CURRENT SYSTEM TIME FOR LABELING PRINT FILE 225 197 hytime.o : hytime.f 226 f77 -c hytime.f 227 228 readq.o : subs/readq.for 229 f77 -c subs/readq.for 230 231 jdate.o : subs/jdate.for 232 f77 -c subs/jdate.for 233 234 spawn.o : subs/spawn.f 235 f77 -c subs/spawn.f 236 237 openr.o : subs/openr.f 238 f77 -c subs/openr.f 239 240 askr.o : subs/askr.for 241 f77 -c subs/askr.for 242 243 openw.o : subs/openw.f 244 f77 -c subs/openw.f 245 246 errset.o : subs/errset.f 247 f77 -c subs/errset.f 248 249 leng.o : subs/leng.for 250 f77 -c subs/leng.for 251 252 dayjl.o : subs/dayjl.for 253 f77 -c subs/dayjl.for 254 255 lask.o : subs/lask.for 256 f77 -c subs/lask.for 257 258 jask.o : subs/jask.for 259 f77 -c subs/jask.for 260 261 klas.o : subs/klas.for 262 f77 -c subs/klas.for 263 264 upstr.o : subs/upstr.for 265 f77 -c subs/upstr.for 266 267 askc.o : subs/askc.for 268 f77 -c subs/askc.for 269 270 box2.o : subs/box2.for 271 f77 -c subs/box2.for 272 273 box3.o : subs/box3.for 274 f77 -c subs/box3.for 275 276 klass.o : subs/klass.for 277 f77 -c subs/klass.for 278 279 ksic.o : subs/ksic.for 280 f77 -c subs/ksic.for 198 g77 -c hytime.f 199 200 # ALSO USES THESE GENERAL PURPOSE SUBROUTINES FROM THE SUBROUTINE LIBRARY 201 # KLAS ASSIGNS A NAME AND NUMBER TO AN EVENT BASED ON LOCATION. 202 # KLASS (NET 1), BOX2 (NET2) - USED BY KLAS. 203 # ISITIN DETERMINES IF EVENT IS IN A PARTICULAR REGION. 204 # UPSTR CONVERTS A STRING TO UPPER CASE. 205 # JASK INTERACTIVE PROMPT & ENTRY OF AN INTEGER. 206 # ASKC INTERACTIVE PROMPT AND ENTRY OF A STRING. 207 # ASKR INTERACTIVE PROMPT AND ENTRY OF A REAL VALUE. 208 # LASK INTERACTIVE PROMPT AND ENTRY OF A LOGICAL VALUE. 209 # LENG DETERMINES THE NON-BLANK LENGTH OF A STRING. 281 210 282 211 # … … 285 214 # 286 215 clean: 287 rm -f a.out core *.o *.obj *% *~ subs/*.o subs/*.obj216 rm -f a.out core *.o *.obj 288 217 289 218 clean_bin: 290 rm -f $B/hyp2000 *219 rm -f $B/hyp2000
Note: See TracChangeset
for help on using the changeset viewer.