#!/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 sts2.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 10 amin 0.001
quit
EOF
plotnps -F7 -W10 -EPS -K < P001.PLT > sts2.eps
convert -trim sts2.eps sts2.png

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

gsac > /dev/null << EOF
qdp off
r impulse.sac
transfer from none to polezero subtype ss1.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 10 amin 0.001
EOF
plotnps -F7 -W10 -EPS -K < P001.PLT > ss1.eps
convert -trim ss1.eps ss1.png

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

gsac > /dev/null << EOF
qdp off
r impulse.sac
transfer from none to polezero subtype epi.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 10 amin 0.001
quit
EOF
plotnps -F7 -W10 -EPS -K < P001.PLT > epi.eps
convert -trim epi.eps epi.png



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