gcc -o gccBTS.exe gccBTS.c
gccBTS 1 inp_PRB.txt out_B_GUM.txt inp_BH_GUM.txt
gccBTS 2 inp_PRB.txt out_B_GEV.txt inp_BH_GEV.txt
gccBTS 3 inp_PRB.txt out_B_SQR.txt inp_BH_SQR.txt
gccBTS 4 inp_PRB.txt out_B_LN3.txt inp_BH_LN3.txt
gccBTS 5 inp_PRB.txt out_B_LP3.txt inp_BH_LP3.txt


echo off
rem ********************************************
rem Drawing by GMT
rem ********************************************
set range=1/10000/0/1000
set scale=12l/10
set xga=g3a1
set yga=g100a100
set xlabel="Return period (year)"
set ylabel="Maximum daily rainfall (mm/day)"
gmtset ANOT_FONT_SIZE 12
gmtset LABEL_FONT_SIZE 12
gmtset TICK_LENGTH 0c

gawk -f awk_klegend.awk inp_legend_BTS.txt > legend.bat
rem ===============================
rem Plotting
rem ===============================
set inp_1=out_B_GUM.txt
set fig_out=fig_B_GUM.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$6}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -W5t15_5_5_5:0 -P -O -K >> %fig_out%
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_out%
gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -JX12/10 -N -P -O -K >> %fig_out%
call legend.bat



set inp_1=out_B_GEV.txt
set fig_out=fig_B_GEV.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$6}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -W5t15_5_5_5:0 -P -O -K >> %fig_out%
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_out%
gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -JX12/10 -N -P -O -K >> %fig_out%
call legend.bat



set inp_1=out_B_SQR.txt
set fig_out=fig_B_SQR.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$6}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -W5t15_5_5_5:0 -P -O -K >> %fig_out%
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_out%
gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -JX12/10 -N -P -O -K >> %fig_out%
call legend.bat



set inp_1=out_B_LN3.txt
set fig_out=fig_B_LN3.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$6}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -W5t15_5_5_5:0 -P -O -K >> %fig_out%
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_out%
gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -JX12/10 -N -P -O -K >> %fig_out%
call legend.bat



set inp_1=out_B_LP3.txt
set fig_out=fig_B_LP3.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$6}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$3}" %inp_1% | psxy -R -J -B -W5t15_5_5_5:0 -P -O -K >> %fig_out%
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_out%
gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -JX12/10 -N -P -O -K >> %fig_out%
call legend.bat



set range=
set scale=
set xga=
set yga=
set xlabel=
set ylabel=
set inp_1=
set fig_out=

del .gmt*
del _*