#!/bin/sh

#####
#     create a temporary directory to store the seismograms
#####

if [ ! -d P_1 ]
then
	mkdir P_1
else
	rm P_1/*
fi

#####
# create a P -wave velocity model such that
# Layer 1:  VP = 6.0 at the top of the layer
#           VP = 6.8 at the bottom of the layer 
#           Layer 1 is 40 km thick
# Layer 2:  VP = 8.0 at the top of the layer
#           VP = 9.0 at the bottom of the layer
#           Layer 2 is 80 km thick
#####

#####
#    for the model format look at
#    PROGRAMS.330/DOC/OVERVIER.pdf/CPS330o.pdf
#####
#    since we will not use S waves, the entries
#    for VS are jsut 1/2 VP
#####

cat > prob1.mod << EOF
MODEL01.1
TEST MODEL
ISOTROPIC
KGS
FLAT EARTH
1-D
VARIABLE VELOCITY
LINE08
LINE09
LINE10
LINE11
 H   VP VS  RHO  QP  QS  ETAP ETAS FREFP FREFS
40. 6.0 3.0 2.2  0.0 0.0 0.0  0.0  1.0   1.0
40. 6.8 3.4 2.5  0.0 0.0 0.0  0.0  1.0   1.0
80. 8.0 4.0 3.8  0.0 0.0 0.0  0.0  1.0   1.0
80. 9.0 4.5 3.3  0.0 0.0 0.0  0.0  1.0   1.0
EOF

#####
#    this si to show the ray diagram at 10 and 200 km
#####
cat > dfile << EOF
 10.0 0.1 30 -10 6
200.0 0.1 30 -10 6
EOF
cprep96 -M prob1.mod -DOP  -N 10 -d dfile
cseis96 -R > cseism96.out
cray96 -ZMAX 50 -XMAX 300 -XMIN -20
plotnps -F7 -W10 -BGFILL -K -EPS < RAY_1.PLT > teps
convert -trim teps RAY_1.png ; rm teps


cat > dfile.profile << EOF
 10.0 0.05 400 -10 6
 20.0 0.05 400 -10 6
 30.0 0.05 400 -10 6
 40.0 0.05 400 -10 6
 50.0 0.05 400 -10 6
 60.0 0.05 400 -10 6
 70.0 0.05 400 -10 6
 80.0 0.05 400 -10 6
 90.0 0.05 400 -10 6
100.0 0.05 400 -10 6
110.0 0.05 400 -10 6
120.0 0.05 400 -10 6
130.0 0.05 400 -10 6
140.0 0.05 400 -10 6
150.0 0.05 400 -10 6
160.0 0.05 400 -10 6
170.0 0.05 400 -10 6
180.0 0.05 400 -10 6
190.0 0.05 400 -10 6
200.0 0.05 400 -10 6
210.0 0.05 400 -10 6
220.0 0.05 400 -10 6
230.0 0.05 400 -10 6
240.0 0.05 400 -10 6
250.0 0.05 400 -10 6
260.0 0.05 400 -10 6
270.0 0.05 400 -10 6
280.0 0.05 400 -10 6
290.0 0.05 400 -10 6
300.0 0.05 400 -10 6
310.0 0.05 400 -10 6
320.0 0.05 400 -10 6
330.0 0.05 400 -10 6
340.0 0.05 400 -10 6
350.0 0.05 400 -10 6
360.0 0.05 400 -10 6
370.0 0.05 400 -10 6
380.0 0.05 400 -10 6
390.0 0.05 400 -10 6
EOF
cprep96 -M prob1.mod -DOP  -N 8 -d dfile.profile
cseis96  > cseism96.out
cpulse96 -EXF -p -V -l 2 | ( cd P_1 ; f96tosac -G )

gsac << EOF
r P_1/*ZEX
bg plt
prs p 0.16667 vl 0 500 tl -10 10 sa 0. amp 0.2 shd pos
q
EOF
mv PRS001.PLT PRS_1.PLT
plotnps -BGFILL -F7 -W10 -EPS -K < PRS_1.PLT > t.eps ; convert -trim t.eps 03_1.png