events_prune.f  examine_sac.f   read_summary.f  rotate_sac.f

events_prune.f
c  f77 -o /park/backup/bin/events_prune events_prune.f /park/backup/Ritz/jlib
c  cross-reference EVENT lines from Weed summary files
c  against data files
c  Goal is to remove all EVENT lines for which there is no data
c  (i.e. station was down)
c  look for agreement between event time and records time of two hours
c  must prune remainder by hand
c  events_rawlist is output of
c  "grep EVENT ???.summary" with leading string before the date removed.
c  infiles_raw is the list of datafiles, in threes
c  output of program is file 'events' with aligned strings

c  program rotate_sac
c  7/8/99 JJP  -- adapted for WEED output, after running event_prune
c  f77 -o /park/backup/bin/rotate_sac rotate_sac.f /park/backup/Ritz/jlib
c  to read filenames of data files (E,N,Z) for a single event
      open(7,file='infiles1',form='formatted')  ! sacfilenames
      open(18,file='events',form='formatted') ! source locations and times
      open(9,file='outfile',form='formatted')

c  program examine_sac
c  7/8/99 JJP - adapted for WEED output, after running event_prune & rotate_sac
c  f77 -o /park/backup/bin/examine_sac examine_sac.f /park/backup/Ritz/jlib
c  to read filenames of data files (bhr,t,z) for a single event
      open(7,file='infiles2',form='formatted')  ! sacfilenames  19*.bh[rtz]
      open(8,file='outfile1',form='formatted')  ! to be in_recfunk
      open(9,file='outfile',form='formatted') ! to be sac_examine
c

FOR SKS Splitting, there is a fourth routine to generate a SAC macro

c  program split_sac
c  xf77 -o /park/backup/bin/split_sac split_sac.f /park/backup/Ritz/jlib.a /park/backup/Plotxy/plotlib.a
c  5/26/03 -- adapted to make splitting estimates in sac macro,
c  calling sac_split from within sac
c  7/8/99 JJP - adapted for WEED output, after running event_prune & rotate_sac
     open(7,file='infiles3',form='formatted')  ! sacfiles  19*.BH[ENZ]*   20*.BH[ENZ]*
      open(9,file='outfile',form='formatted') ! to be sac_split

and the sac_split macro calls the FORTRAN program sac_splitting.f:
c  program sac_splitting
c  JJP 5/26/03
c  xf77 -o /park/backup/bin/sac_splitting sac_splitting.f /park/backup/Plotxy/plotlib.a /park/backup/Ritz/jlib.a
c  to read SAC-format data  HARDWIRED to read three filenames lowe lown lowz
c  do gridsearch over delay time and fast axis strike