Introduction

The purpose of this problem is to make some wholespace synthetic seismograms

1. Read the manual that is located on crust.eas.slu.edu  at /home/rbh/PROGRAMS.330/DOC/OVERVIEW.pdf/cps330.pdf, e.g.,
  
     acroread /home/rbh/PROGRAMS.330/DOC/OVERVIEW.pdf/cps330.pdf

2.  Using the program mkmod96 to create the following velocity model (NOTE give this the name 2half.mod):

    H Vp Vs Density Qp Qs Etap Etas FrefP FrefS
    -4 1 0 1.0 0 0 0 0 1 1
      2 2 0 1.5 0 0 0 0 1 1

3. Plot the model using shwmod96

4. Create a distance file, which you will name dfile, with entries "Distance DT NPTS T0 VRED" for the distances used for figure 3.2.2 in the text, e.g., you will require entries such as
    1.2 0.05 256 0.0 0.0
    1.6 0.05 256 0.0 0.0
   ......
   6.8 0.05 256 0.0 0.0

5. Use the program hprep96 to generate the ray paths between the source and receiver:
     hprep96 -M 2half.mod -HS -1.0 -HR -1.5 -ALL -TH -BH -d dfile

6. Then run the program hwhole96

     hwhole96

      and convolve with the source time function
   
     hpulse96 -V -p -l 4 > file96

     and them plot the figure using fprof96

     fprof96 -S 2.0 < file96

     You can view the plot with the following commands:

     plotxvig < FPROF96.PLT

     or you can get a PostScript file that you can plot using the command

     plotnps -F7 -W10 -K < FPROF96.PLT > my.ps

  Of course you can also look at the traces using gsac by creating the SAC waveform files using the command f96tosac -B < file96

7. The plots that will be created will be the vertical and radial component of the particle velocity, and since this is a fluid problem, you will also see the pressure field in the fluid.

Questions

What is the difference between the shape of the velocity and pressure synthetic time histories?
How dows amplitude change with distance>