rem ********************************************
rem Probability plotting of maximum daily rainfall
rem ********************************************
rem ===============================
rem making plotting data files
rem ===============================
gawk -f awk_Gumbel.awk inp_PRB_DM.txt > dat_PRB_G_DM.txt
gawk -f awk_Normal.awk inp_PRB_DM.txt > dat_PRB_N_DM.txt
gawk -f awk_Lognrm.awk inp_PRB_DM.txt > dat_PRB_L_DM.txt
gawk -f awk_Weibull.awk inp_PRB_DM.txt > dat_PRB_W_DM.txt



rem ********************************************
rem Drawing by GMT
rem ********************************************
gmtset ANOT_FONT_SIZE 14
gmtset LABEL_FONT_SIZE 14
gmtset TICK_LENGTH -0.2c



rem ********************************************
rem Gumbel distribution plotting
rem ********************************************
set range0=0/1/-1.93265/6.90726
set range1=0/500/-1.93265/6.90726
set range2=0/12/0/15
set scale0=12/15
set scale1=12/15
set xga=g100a100
set yga=a1
set xlabel="Maximun daily rainfall [x] (mm/day)"
set ylabel="Canonical variate [y=-ln{-ln(p)}]"
set inp_1=dat_PRB_G_DM.txt
set fig_out=fig_PRB_G_DM.eps
rem ---------------------------
rem y-axis text and grid
rem ---------------------------
psbasemap -R%range0% -JX%scale0% -G255 -P -X5 -Y6 -K > %fig_out%
gawk "{if(12<=NR&&NR<=22){sub(/#/,\"\",$0);print $0}}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "{if(24<=NR&&NR<=56){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W0.25p -M -N -P -O -K >> %fig_out%
rem ---------------------------
rem plotting datas
rem ---------------------------
psxy %inp_1% -R%range1% -JX%scale1% -B%xga%:%xlabel%:/%yga%:%ylabel%:wSEn -SC0.2 -N -P -O -K >> %fig_out%
gawk "{if(9<=NR&&NR<=10){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W1 -N -P -O -K >> %fig_out%
rem ---------------------------
rem comment and y-axis label
rem ---------------------------
gawk "{if(NR==1){sub(/#/,\"\",$0);printf \"0.3 14.5 12 0 0 ML %%s\n\",$0}}" %inp_1% > _val.txt
gawk "{if(NR==2){sub(/#/,\"\",$0);printf \"0.3 14.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==3){sub(/#/,\"\",$0);printf \"0.3 13.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==4){sub(/#/,\"\",$0);printf \"0.3 13.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==5){sub(/#/,\"\",$0);printf \"0.3 12.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
pstext _val.txt -R%range2% -JX%scale0% -N -P -O -K >> %fig_out%
echo 0.2 15.2 14 0 0 BL Gumbel distribution | pstext -R -J -N -P -O -K >> %fig_out%
echo -2 7.5 14 90 0 MC Non-exceeding probability [p] | pstext -R -J -N -P -O >> %fig_out%



rem ********************************************
rem Normal distribution plotting
rem ********************************************
set range0=0/1/-3.09023/3.09023
set range1=0/500/-3.09023/3.09023
set range2=0/12/0/15
set scale0=12/15
set scale1=12/15
set xga=g100a100
set yga=a1
set xlabel="Maximun daily rainfall [x] (mm/day)"
set ylabel="Canonical variate [y]"
set inp_1=dat_PRB_N_DM.txt
set fig_out=fig_PRB_N_DM.eps
rem ---------------------------
rem y-axis text and grid
rem ---------------------------
psbasemap -R%range0% -JX%scale0% -G255 -P -X5 -Y6 -K > %fig_out%
gawk "{if(12<=NR&&NR<=22){sub(/#/,\"\",$0);print $0}}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "{if(24<=NR&&NR<=56){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W0.25p -M -N -P -O -K >> %fig_out%
rem ---------------------------
rem plotting datas
rem ---------------------------
psxy %inp_1% -R%range1% -JX%scale1% -B%xga%:%xlabel%:/%yga%:%ylabel%:wSEn -SC0.2 -N -P -O -K >> %fig_out%
gawk "{if(9<=NR&&NR<=10){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W1 -N -P -O -K >> %fig_out%
rem ---------------------------
rem comment and y-axis label
rem ---------------------------
gawk "{if(NR==1){sub(/#/,\"\",$0);printf \"0.3 14.5 12 0 0 ML %%s\n\",$0}}" %inp_1% > _val.txt
gawk "{if(NR==2){sub(/#/,\"\",$0);printf \"0.3 14.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==3){sub(/#/,\"\",$0);printf \"0.3 13.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==4){sub(/#/,\"\",$0);printf \"0.3 13.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==5){sub(/#/,\"\",$0);printf \"0.3 12.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
pstext _val.txt -R%range2% -JX%scale0% -N -P -O -K >> %fig_out%
echo 0.2 15.2 14 0 0 BL Normal distribution | pstext -R -J -N -P -O -K >> %fig_out%
echo -2 7.5 14 90 0 MC Non-exceeding probability [p] | pstext -R -J -N -P -O >> %fig_out%



rem ********************************************
rem Logarithmic normal distribution plotting
rem ********************************************
set range0=0/1/-3.09023/3.09023
set range1=10/1000/-3.09023/3.09023
set range2=0/12/0/15
set scale0=12/15
set scale1=12l/15
set xga=g3a2
set yga=a1
set xlabel="Maximun daily rainfall [x] (mm/day)"
set ylabel="Canonical variate [y]"
set inp_1=dat_PRB_L_DM.txt
set fig_out=fig_PRB_L_DM.eps
rem ---------------------------
rem y-axis text and grid
rem ---------------------------
psbasemap -R%range0% -JX%scale0% -G255 -P -X5 -Y6 -K > %fig_out%
gawk "{if(12<=NR&&NR<=22){sub(/#/,\"\",$0);print $0}}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "{if(24<=NR&&NR<=56){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W0.25p -M -N -P -O -K >> %fig_out%
rem ---------------------------
rem plotting datas
rem ---------------------------
psxy %inp_1% -R%range1% -JX%scale1% -B%xga%:%xlabel%:/%yga%:%ylabel%:wSEn -SC0.2 -N -P -O -K >> %fig_out%
gawk "{if(9<=NR&&NR<=10){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W1 -N -P -O -K >> %fig_out%
rem ---------------------------
rem comment and y-axis label
rem ---------------------------
gawk "{if(NR==1){sub(/#/,\"\",$0);printf \"0.3 14.5 12 0 0 ML %%s\n\",$0}}" %inp_1% > _val.txt
gawk "{if(NR==2){sub(/#/,\"\",$0);printf \"0.3 14.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==3){sub(/#/,\"\",$0);printf \"0.3 13.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==4){sub(/#/,\"\",$0);printf \"0.3 13.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==5){sub(/#/,\"\",$0);printf \"0.3 12.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
pstext _val.txt -R%range2% -JX%scale0% -N -P -O -K >> %fig_out%
echo 0.2 15.2 14 0 0 BL Logarithmic normal distribution | pstext -R -J -N -P -O -K >> %fig_out%
echo -2 7.5 14 90 0 MC Non-exceeding probability [p] | pstext -R -J -N -P -O >> %fig_out%



rem ********************************************
rem Weibull distribution plotting
rem ********************************************
set range0=0/1/-6.90726/1.93264
set range1=10/1000/-6.90726/1.93264
set range2=0/12/0/15
set scale0=12/15
set scale1=12l/15
set xga=g3a2
set yga=a1
set xlabel="Maximun daily rainfall [x] (mm/day)"
set ylabel="Canonical variate [y=ln{-ln(1-p)}]"
set inp_1=dat_PRB_W_DM.txt
set fig_out=fig_PRB_W_DM.eps
rem ---------------------------
rem y-axis text and grid
rem ---------------------------
psbasemap -R%range0% -JX%scale0% -G255 -P -X5 -Y6 -K > %fig_out%
gawk "{if(12<=NR&&NR<=22){sub(/#/,\"\",$0);print $0}}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "{if(24<=NR&&NR<=56){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W0.25p -M -N -P -O -K >> %fig_out%
rem ---------------------------
rem plotting datas
rem ---------------------------
psxy %inp_1% -R%range1% -JX%scale1% -B%xga%:%xlabel%:/%yga%:%ylabel%:wSEn -SC0.2 -N -P -O -K >> %fig_out%
gawk "{if(9<=NR&&NR<=10){sub(/#/,\"\",$0);print $0}}" %inp_1% | psxy -R -J -W1 -N -P -O -K >> %fig_out%
rem ---------------------------
rem comment and y-axis label
rem ---------------------------
gawk "{if(NR==1){sub(/#/,\"\",$0);printf \"0.3 14.5 12 0 0 ML %%s\n\",$0}}" %inp_1% > _val.txt
gawk "{if(NR==2){sub(/#/,\"\",$0);printf \"0.3 14.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==3){sub(/#/,\"\",$0);printf \"0.3 13.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==4){sub(/#/,\"\",$0);printf \"0.3 13.0 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
gawk "{if(NR==5){sub(/#/,\"\",$0);printf \"0.3 12.5 12 0 0 ML %%s\n\",$0}}" %inp_1% >> _val.txt
pstext _val.txt -R%range2% -JX%scale0% -N -P -O -K >> %fig_out%
echo 0.2 15.2 14 0 0 BL Weibull distribution | pstext -R -J -N -P -O -K >> %fig_out%
echo -2 7.5 14 90 0 MC Non-exceeding probability [p] | pstext -R -J -N -P -O >> %fig_out%



set range0=
set range1=
set range2=
set scale0=
set scale1=
set xga=
set yga=
set xlabel=
set ylabel=
set inpl=
set inp_1=
set inp_2=
set fig_out=

del .gmt*
del _*