Plane wave synthetics

This assignment requires that you run the Fortran program
"superplane" which is on the Macs in ~shearer/PROG/CLASS/227C/PLANE/
and on the Suns in ~shearer/CLASS/227C/d.plane/.  It is also
available via anonymous FTP to mahi.ucsd.edu in the /pub/227C/PLANE directory.
The other required files are also in these directories.  

SUPERPLANE computes the plane-wave response of a stack of layers.  
An example of a command file that runs this program is:

do.superplane

This operates on the velocity model file vz.mod1, which looks like:

air       0.0       0.0
2.000     1.150     1.5       99999.    99999.    0.5
5.000     2.890     2.5       99999.    99999.    1.0


The first line is for a free surface with air above.  The following
lines are for the layers in the model, specifying

pvel      svel      density   Qp        Qs        h (layer thickness, km)

The first line is in a fixed format (things must be in
the correct columns); the following lines are read with
a free-format read.

The command file do.superplane will generate the output file out.mod1
The first 6 lines of this file are header information.  The following
lines give the frequency and amplitude points.

This can be plotted using plotxy using the command file do.plotxy
or with your favorite alternative plotting program.


(1) Copy do.superplane and vz.mod1 into your own 
    directory and verify that you can run the program and plot the output.

    (non Sun or Mac users) See if you can compile the superplane program
    on your machine and get it to run.  You will need to include (or
    link with) the subs.f and cinvert.f subroutines.  
    
(2) The city of Santa Tembloria is planning to build a new hospital.
    The city commissions your geotechnical consulting firm to
    evaluate three possible sites.  From borehole sampling, you
    find the following velocity and attenuation profiles:

SITE 1
1.600     0.500     1.3       99999.    30        0.1
3.000     1.73      2.0       99999.    100       1.0
5.000     2.890     2.5       99999.    99999.    1.0

SITE 2
2.2       1.2       1.8       9999.     50        0.4
3.0       1.73      2.0       9999.     100       1.0
5.000     2.890     2.5       99999.    99999.    1.0

SITE 3
1.8       0.7       1.7       9999.     30        0.05
3.0       1.73      2.0       9999.     100       1.0
5.000     2.890     2.5       99999.    99999.    1.0

(a) Plot the S-velocity vs. depth profiles for the three sites.

(b) Run superplane on these models, using an input S-wave of
    ray parameter p=0.1.  Plot the response spectra for the
    Ux (horizontal component) of motion.

(c) Buildings are typically most sensitive to damage at frequencies
    around 1 Hz.  Which site do you recommend to the city?