# ********************************************
# Drawing by GMT (response spectrum)
# ********************************************
gmt set FONT_LABEL 12p
gmt set MAP_TICK_LENGTH_PRIMARY 0

range=0/50/-180/180
scale=12/4
xga=g10a10
yga=g90a90
xlabel='Frequency (Hz)'
ylabel='Phase angle (degree)'

fig=fig_fps1.eps
inp=out_fps1.csv
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
awk 'BEGIN{FS=","}{if(2<=NR){print $2,$3}}' $inp | gmt psxy -R -J -B -W0.25 -O -K >> $fig
echo '1 190 12p 0 BL Fourier Phase Spectrum (input wave)' | gmt pstext -R -J -F+f+a+j -N -O >> $fig

fig=fig_fps2.eps
inp=out_fps2.csv
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
awk 'BEGIN{FS=","}{if(2<=NR){i=i+1;print $2,$3}}' $inp | gmt psxy -R -J -B -W0.25 -O -K >> $fig
echo '1 190 12p 0 BL Fourier Phase Spectrum (simulated wave)' | gmt pstext -R -J -F+f+a+j -N -O >> $fig

range=0/30/-150/150
scale=12/4
xga=g10a10
yga=g50a50
xlabel='Time (sec)'
ylabel='Acceleration (gal)'

fig=fig_fps3.eps
inp=inp_ran1.csv
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
awk 'BEGIN{FS=",";i=-1}{if(2<=NR){i=i+1;print 0.01*i,$1}}' $inp | gmt psxy -R -J -B -W0.25 -O -K >> $fig
echo '1 160 12p 0 BL Acceleration data (input wave)' | gmt pstext -R -J -F+f+a+j -N -O >> $fig

fig=fig_fps4.eps
inp=out_pw2.csv
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
awk 'BEGIN{FS=",";i=-1}{if(2<=NR){i=i+1;print 0.01*i,$1}}' $inp | gmt psxy -R -J -B -W0.25 -O -K >> $fig
echo '1 160 12p 0 BL Phase wave Acceleration data (simulated wave)' | gmt pstext -R -J -F+f+a+j -N -O >> $fig