#!/bin/sh
set -x


#####
#	create an impulse with unit area
#####
NPTS=16384
DELTA=0.0125

gsac > /dev/null << EOF
funcgen impulse npts ${NPTS} delta ${DELTA}
# gsac creates a pulse with unit amplitude spectrum
w impulse.sac
EOF

#####
#	create a plot of the sts2.rsp
#####

gsac > /dev/null << EOF
qdp off
r impulse.sac
transfer from none to polezero subtype asts2.rsp
ch KSTNM STS2
ch KCMPNM VEL_SENS
w sts2.sac
pctl ylab "STS2 Vel. Sens. (cts/m/s)"
fft
bg plt
grid on minor on
psp fmin 0.01 fmax 50 amax 1.0e+10 amin 1.0e+7
quit
EOF
CAL << EOF
NEWPEN 
1
CENTER
5.3 7.2 0.15 'Velocity Sensitivity (CTS/M/S)' 0.0
PEND
EOF
cat P001.PLT CAL.PLT | plotnps -F7 -W10 -EPS -K > asts2.eps
convert -trim asts2.eps asts2.png

#####
#	create a plot of the ss1.rsp
#####

gsac > /dev/null << EOF
qdp off
r impulse.sac
transfer from none to polezero subtype ass1.rsp
ch KSTNM SS1
ch KCMPNM VEL_SENS
w ss1.sac
fft
pctl ylab "SS1 Vel. Sens. (cts/m/s)"
bg plt
grid on minor on
psp fmin 0.01 fmax 50 amax 1.0e+9 amin 1.0e+6
EOF
CAL << EOF
NEWPEN 
1
CENTER
5.3 7.2 0.15 'Velocity Sensitivity (CTS/M/S)' 0.0
PEND
EOF
cat P001.PLT CAL.PLT | plotnps -F7 -W10 -EPS -K > ass1.eps
convert -trim ass1.eps ass1.png

#####
#	create a plot of the epi.rsp
#####

gsac > /dev/null << EOF
qdp off
r impulse.sac
transfer from none to polezero subtype aepi.rsp
ch KSTNM EPI
ch KCMPNM ACC_SENS
w epi.sac
fft
pctl ylab "Epi Accel. Sens. (cts/m/s/s)"
bg plt
grid on minor on
psp fmin 0.01 fmax 50 amax 1.0e+7 amin 1.0e+4
quit
EOF
CAL << EOF
NEWPEN 
1
CENTER
5.3 7.2 0.15 'Acceleration Sensitivity (CTS/M/S/S)' 0.0
PEND
EOF
cat P001.PLT CAL.PLT | plotnps -F7 -W10 -EPS -K > aepi.eps
convert -trim aepi.eps aepi.png



#####
#	clean up
#####
rm -f impulse.sac
rm -f *.PLT *.eps CAL.cmd
rm -f *.sac