gcc -o gccJKK.exe gccJKK.c
gccJKK 1 inp_PRB.txt out_J_GUM.txt inp_JH_GUM.txt
gccJKK 2 inp_PRB.txt out_J_GEV.txt inp_JH_GEV.txt
gccJKK 3 inp_PRB.txt out_J_SQR.txt inp_JH_SQR.txt
gccJKK 4 inp_PRB.txt out_J_LN3.txt inp_JH_LN3.txt
gccJKK 5 inp_PRB.txt out_J_LP3.txt inp_JH_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_JKK.txt > legend.bat
rem ===============================
rem Plotting 1
rem ===============================
set inp_1=out_J_GUM.txt
set fig_out=fig_J_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,$4+$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4-$5}" %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
rem ===============================
rem Plotting 2
rem ===============================
set inp_1=out_J_GEV.txt
set fig_out=fig_J_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,$4+$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4-$5}" %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
rem ===============================
rem Plotting 3
rem ===============================
set inp_1=out_J_SQR.txt
set fig_out=fig_J_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,$4+$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4-$5}" %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
rem ===============================
rem Plotting 4
rem ===============================
set inp_1=out_J_LN3.txt
set fig_out=fig_J_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,$4+$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4-$5}" %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
rem ===============================
rem Plotting 5
rem ===============================
set inp_1=out_J_LP3.txt
set fig_out=fig_J_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,$4+$5}" %inp_1% | psxy -R -J -B -W5t10_5:0 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(3<=NR)print $1,$4-$5}" %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 _*