Sunday, October 16, 2016

New Versions of rfmig_mcboot.f, rfmig_boot.f, rfmig_mboot.f, rfmig_cboot.f, rfmig_mcbp.f, rfmig_mcbf.f, all with improved handling of larger data sets (500+ events).

Avoided some NaNs in the harmonic-regression calculations.  Earlier versions of the codes included some single-precision computations that led the algorithm into darkness.  Tested with a 1000-record dataset.  The link takes you to the Google Drive folder where the codes reside.

Wednesday, June 15, 2016

Picking P waves, PKP waves, PP waves, etc for Recfunk codes

c program recfunk_pick_through
c version of recfunk_pick that reads seismic record start times from
c the sac header, and writes to in_recpick the filenames that have been picked
c using EITHER the A-marker or the T1, T2, or T3 markers in the SAC header, placed by running plotpk
c within SAC. Set by pressing "T" and 1,2 or 3 in plotpk sequentially
c 10/13/05 JJP
c input file: in_recpick_raw
c output file: in_recpick