# ********************************************
# Drawing by GMT
# ********************************************
fig=fig_func.eps
range=0/150/0/1.3
scale=12/6
xga=a30f5
yga=a1f0.1
xlabel='time  t (sec)'
ylabel='e(t)'
# ===============================
# Making of batch file for legend
# ===============================
gmt set FONT_LABEL 12p
gmt set MAP_LABEL_OFFSET 0.0c
gmt set MAP_TICK_LENGTH_PRIMARY 0.2c
# ===============================
# Plotting
# ===============================
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
awk 'BEGIN{FS=",";i=-1}{if(4<=NR){i=i+1;print i*0.01,$2}}' inp_ran1.csv | gmt psxy -R -J -B -W2,5_5:0 -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(4<=NR){i=i+1;print i*0.01,$2}}' inp_ran2.csv | gmt psxy -R -J -B -W2,5_5:0 -O -K >> $fig
awk 'BEGIN{FS=",";i=-1}{if(4<=NR){i=i+1;print i*0.01,$2}}' inp_ran3.csv | gmt psxy -R -J -B -W2 -O -K >> $fig
gmt psxy -R -J -B -W1 -O -K << EOT >> $fig
  0 0.1
125 0.1
EOT
gmt pstext -R -JX -F+f+a+j -N -O << EOT >> $fig
126 0.10 11p  0 ML e=0.1
  5 1.03 11p 90 ML a=5
 15 1.03 11p 90 ML b=15
 25 1.03 11p 90 ML b=25
 35 1.03 11p 90 ML b=35
 30 0.13 11p 90 ML c=30
 60 0.13 11p 90 ML c=60
120 0.13 11p 90 ML c=120
  3 1.30 12p  0 BL Envelop function of acceleration wave due to earthquake 
EOT