rfmigrate.f -- recfunk version with migration
Click the title of the post to access the code rfmigrate.f on Google Drive directory JParkCodes.
c program rfmigrate
c 9/20/00 JJP -- adapted from recfunk --- UPDATED 08/19/07 JJP
c updated again 07/22/15 JJP
c migrates MTC receiver functions in the frequency domain.
c requires a stacking model in the anirec format
c such a model may have anisotropy parameters,
c but migration code only uses the isotropic velocities.
c
c has kluge to cheat the pre-event noise for synthetic records 3/12/00 JJP
c check to see if the kluge is commented out, because anirec_synth was updated
c in 2015 to have a proper-length pre-event time window.
c
c this version of the RF code reads a file of data filenames
c you have two choices: either read the time intervals in the filename file
c or read them in the sac header
c the data must be binary SAC format
c horizontals must be rotated to radial and transverse
c
c for start times in the file:
c the file defaults to "in_recfunk" and has lines of the form:
c
c 1997.070.19.33.bh? <-- code="" div="" r="" replaces="" t="" with="" z="">
c 57 52 <-- analysis="" div="" duration="" of="" sec="" start="" time="" window="">
-->
c 1997.076.08.15.bh?
c 62 62
c ...
c ...
c ...
c stop <-- 799="" code="" data="" div="" events="" finished="" is="" max="" tells="" that="">
-->
c
c
c for start times in the SAC header, the default file is "in_recpick"
c reads seismic record start times from the sac header
c will search the SAC header for specific markers of P phases
c T1 - P, Pdiff ahead(12)
c T2 - PKP,PKIKP ahead(13)
c T3 - PP ahead(14)
c T1=T2=T3=0 ==> use original A-marker ahead(9)
c
c legacy data has timing in the
c A-marker (ahead(9)) which doesnt indicate which phase was marked,
c For the A-marker, the phase is identified as P/Pdiff for DELTA<120 div="" nbsp="">
120>
-->
c PKP/PKIKP for DELTA>120
c
c code does NOT combine data with different sample rates
c it is possible to spline interpolate the data files for a uniform sample rate
c see the code in rfmig_mcboot.f, which should be spliceable into this code.
c data files are limited to 99K pnts. To increase, see common block /datastuff/
c
c many intermediate quantities are plotted with PLOTIT as the code proceeds.
c other intermediate quantities can be plotted by uncommenting calls to PLOTIT
c
c the code writes the BAZ- and EPICEN-dependent RFs to files
c in a format easily digested by GMT (traces are separated by ">" lines)
c
c filenames: out[rt]_baz.grid oum[rt]_epi.grid oum[rt]_baz.grid
c
c out --> no moveout correction; oum --> moveout correction
c
c these files are overwritten everytime you run the program
c so rename them to save them
c
c xf77 -o /park/backup/bin/rfmigrate rfmigrate.f /park/backup/Plotxy/plotlib.a /park/backup/Ritz/eislib.a /park/backup/Ritz/jlib.a
c xf77 -o /Users/jjpark/bin/rfmigrate rfmigrate.f /Users/jjpark/Plotxy/plotlib.a /Users/jjpark/Ritz/eislib.a /Users/jjpark/Ritz/jlib.a
c
No comments:
Post a Comment