Dear Jim,

This is the whole 9 yards.  We use this
program for several things - primarily
the estimation of WAS amplitudes and
ground motions such as PGA, PGV, and PGD.
The main program is called redi_wave and
can be invoked with different options.  Of
relevance to you -w for WAS and -g for
the ground motions.

This codes does do all kinds of things,
so it may be difficult for you to work
through it.  Here is a general outline.

redi_wave is entirely command-line driven.
Various options tell it what to do.  The
final argument is the list of files to work
with.

Once the program figures out what it is supposed
to do (and it can do many different things with 
the same data files), it loops over the list of
data files.  It opens the files, reads in the
data and retrieves information about the station
and instrument.  It checks for clipping and 
optionally checks to make sure that it has the
correct amount of time.  Then it actually does
something.

est_wood is the subroutine that computes the WAS
amplitudes.  This routine calls coinst2, which is
the workhorse Fourier domain routine.  coinst2
does the deconvolution/reconvolution/filtering
operations.

Have fun!
