labelS="Time (hour)"
labelW="Discharge (m@+3@+/s)"
labelE="Rainfall (mm/hr)"
labelN=''
# ===============================
# Making of batch file for legend
# ===============================
cat << EOT > _leg.txt
12 $ 8 $ tr $ 0.3 $ 0.3
10 $ 0.5 $ 4
Rainfall (Observed) $ ---     $ S0.3 $ 229  $ 0.5
Given n1 and TL     $ 1,3_3:0 $ C0.1 $ 255  $ 1
n1=1 and TL=0       $ 1       $ C0.1 $ 0    $ ---
EOT
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _leg.txt > _legend.txt
# **************************************************************************
# * Execute GMT
# **************************************************************************
gmt set FONT_ANNOT_PRIMARY 12
gmt set MAP_ANNOT_OFFSET_PRIMARY 0.3c
gmt set FONT_LABEL 12p
gmt set MAP_LABEL_OFFSET 0.5c
gmt set MAP_TICK_LENGTH_PRIMARY -0.2c
gmt psbasemap -R$rangeNE -JX$scaleNE -Bx$afgN -By$afgE+l"$labelE" -BnE -P -K > $fig
awk 'BEGIN{FS=","}5<=NR{print $2+0.5,$3}' $inp | gmt psxy -J -R -Sb1ub0 -W0.5 -G#eeeeee -O -K >> $fig
gmt psbasemap -R$rangeSW -JX$scaleSW -Bx$afgS+l"$labelS" -By$afgW+l"$labelW" -BSW -O -K >> $fig
awk 'BEGIN{FS=","}5<=NR{print $2+1,$4}' $inp | gmt psxy -J -R -B -W1,3_3:0 -O -K >> $fig
awk 'BEGIN{FS=","}5<=NR{print $2+1,$4}' $inp | gmt psxy -J -R -B -SC0.1 -G255 -W1 -O -K >> $fig
awk 'BEGIN{FS=","}5<=NR{print $2+1,$5}' $inp | gmt psxy -J -R -B -W1 -O -K >> $fig
awk 'BEGIN{FS=","}5<=NR{print $2+1,$5}' $inp | gmt psxy -J -R -B -SC0.1 -G0 -O -K >> $fig
awk 'BEGIN{FS=","}NR==1{printf "0.02 1.02 12p 0 BL %s\n",$0}' $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -O -K >> $fig
source _legend.txt