gmt set FONT_LABEL 12p
gmt set MAP_TICK_LENGTH_PRIMARY 0

range=0/120/-400/400
scale=12/4
xga=g30a30
yga=g100a100
xlabel='time (sec)'
ylabel='Acceleration (gal)'

gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -Y21 -K > $fig
awk 'BEGIN{FS=",";i=-1}{if(3<NR){i=i+1;print i*0.01,$1}}' $inp_1 | gmt psxy -R -J -B -W0.25 -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(3<NR){i=i+1;print i*0.01,$2*100}}' $inp_1 | gmt psxy -R -J -B -W2,red -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(3<NR){i=i+1;print i*0.01,-$2*100}}' $inp_1 | gmt psxy -R -J -B -W2,red -O -K >> $fig
echo '0.01 420 10p 0 BL Original' | gmt pstext -R -J -F+f+a+j -N -O -K >> $fig

gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -Y-6 -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(3<NR){i=i+1;print i*0.01,$2}}' $inp_2 | gmt psxy -R -J -B -W0.25 -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(3<NR){i=i+1;print i*0.01,$2*300}}' $inp_1 | gmt psxy -R -J -B -W2,red -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(3<NR){i=i+1;print i*0.01,-$2*300}}' $inp_1 | gmt psxy -R -J -B -W2,red -O -K >> $fig
echo 0.01 420 10p 0 BL Simulated | gmt pstext -R -J -F+f+a+j -N -O -K >> $fig

range=0.01/1000/0.01/10000
scale=12l/10l
xga=g1a1
yga=g1a1
xlabel='Period (sec)'
ylabel='Response acceleration (gal)'
cat << EOT > _inp_legend.txt
12 $ 10 $ bl $ 0.3 $ 0.3
10 $ 0.7 $ 3
Target    $ 2,red     $ ---  $ ---  $ ---
Simulated $ 1         $ ---  $ ---  $ ---
Original  $ 2,blue    $ ---  $ ---  $ ---
EOT
awk -f awk_klegend.awk _inp_legend.txt > _legend.txt
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWSen -P -Y-12 -O -K >> $fig
awk 'BEGIN{FS=","}{if(3<NR)print $2,$5}' $inp_3 | gmt psxy -R -J -B -W2,red  -O -K >> $fig
awk 'BEGIN{FS=","}{if(3<NR)print $2,$4}' $inp_3 | gmt psxy -R -J -B -W2      -O -K >> $fig
awk 'BEGIN{FS=","}{if(3<NR)print $2,$3}' $inp_3 | gmt psxy -R -J -B -W2,blue -O -K >> $fig
echo '0.01 12000 10 0 BL Damping factor h=0.05' | gmt pstext -R -J -F+f+a+j -N -O -K >> $fig
source _legend.txt