rem bat_bts_draw.bat
set range=1/10000/0/1000
set scale=15l/10
set xga=g3a1
set yga=g100a100
set xlabel="Return period (year)"
set ylabel="Maximum daily rainfall (mm/day)"
gmtset ANOT_FONT_SIZE 16
gmtset LABEL_FONT_SIZE 16
gmtset TICK_LENGTH 0c

echo 15 $ 10 $ tl $ 0.3 $ 0.3 > _inp_legend.txt
echo 14 $ 0.5 $ 6 >> _inp_legend.txt
echo Original X(T)        $ 5t15_5_5_5:0 $ C0.2 $ 255/255/255 $ 3 >> _inp_legend.txt
echo Bootstrap X@#B@#(T)  $ 5            $ ---  $ ---         $ --- >> _inp_legend.txt
echo Bootstrap 95%%%% limit $ 5t10_5:0   $ ---  $ ---         $ --- >> _inp_legend.txt
gawk -f \WANtaroHP_plot\gmt_basic\awk_klegend.awk _inp_legend.txt > legend.bat
rem ===============================
rem Plotting
rem ===============================
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$6}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -W5t15_5_5_5:0 -P -O -K >> %fig%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -SC0.2 -G255/255/255 -W3 -P -O -K >> %fig%
gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.1 16 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -JX12/10 -N -P -O -K >> %fig%
call legend.bat