Assignment: WKBJ synthetics

All files for this assignment are on the Macs in
~shearer/PROG/CLASS/227C/WKBJ/

The files are as follows:

sacwkbj.f  =  WKBJ Fortran source code
              This version outputs into SAC format
sacwkbj    =  Executable version
Makefile   =  Makefile to compile sacwkbj (if you ever want to make changes!)
diskio.o   =  Compiled C i/o library necessary for writing SAC files
vz.socal227  =  velocity depth model for southern California
                (depth  alpha  beta  den) 
                linear velocity gradients assumed between model points
do.sacwkbj  =  Example command file to run sacwkbj for socal model

Last year I ran this program on the Suns.  It compiled with a bunch
of warning messages on the Mac but seems to work.  Apparently, there 
is also something irregular about the SAC headers, which gives a warning 
when the files are read into SAC.  (The warning is not very helpful 
because it does not say which field in the header has the problem)


Your assignment:

(1) Run sacwkbj for the vz.socal227 model.  Output seismograms every
20 km from 0 to 300 km.  Use a sample interval of 0.02 s.  Plot your 
results using the SAC program and print out hard copies of the synthetics 
both with and without a reduction velocity of 8 km/s.  Sketch the travel 
time curves on top of your plots.  Label the Pg, PmP and Pn phases.  Which
phases are Hilbert transformed?

(2) What is the approximate Pg/Pn amplitude ratio at 200 km?  Now suppose
that you are trying to match data in which the Pg/Pn amplitude ratio is 2.
Experiment with modifying the velocity model so that the synthetics
better match this constraint.  Make a plot of your v(z) model compared
to the original v(z) model.  Also include a comparison plot of the seismograms
at 200 km.