anirec_synth.f --- a code to compute synthetic P coda from a collection of SAC-format data

click title to see Google Drive directory JParkCodes

c anirec_synth - anirec to compute synthetics for a specified model and
c a specified collection of seismic records (read from in_recpick-format file of filenames)
c modified to accept SACfiles with start time info in the header, in A-marker (ahead(9))
c because the timing information is not used, the value of A-marker is not important.
c However, in the synthetics, the A-marker and the T1-marker (or T2) variables in the header are set to 5.0 seconds
c for easy analysis by recfunk codes that use the header information for timing.
c code divides between T1 (P) or T2 (PKP) markers based on epicentral distance, with 120 degrees the divide
c the synthetic files are s_filename.bh[zrt], where filename is the name of the datafile.
c the code writes in_recpick_synth, a list of the filenames for recfunk_pick and other recfunk codes that
c expect the start time in the SAC header
c

The goal of this code is to generate P-coda form a specified model to replicate the earthquake distribution of a real data set. This allows the analyst to test whether the oddities in an RF sweep are likely caused by an imperfect data distribution, or by a shortcoming of his/her model.