rem ********************************************
rem Probability plotting of maximum daily rainfall
rem ********************************************
rem ===============================
rem making plotting data files
rem ===============================
gawk -f awk_prb_paper.awk 1 inp_RF_M.txt > dat_M_PRB_ND.txt
gawk -f awk_prb_paper.awk 2 inp_RF_M.txt > dat_M_PRB_LN.txt
gawk -f awk_prb_paper.awk 3 inp_RF_M.txt > dat_M_PRB_GD.txt
gawk -f awk_prb_paper.awk 4 inp_RF_M.txt > dat_M_PRB_WD.txt
gawk -f awk_prb_paper.awk 5 inp_RF_M.txt > dat_M_PRB_ED.txt

gawk -f awk_prb_paper.awk 1 inp_RF_K.txt > dat_K_PRB_ND.txt
gawk -f awk_prb_paper.awk 2 inp_RF_K.txt > dat_K_PRB_LN.txt
gawk -f awk_prb_paper.awk 3 inp_RF_K.txt > dat_K_PRB_GD.txt
gawk -f awk_prb_paper.awk 4 inp_RF_K.txt > dat_K_PRB_WD.txt
gawk -f awk_prb_paper.awk 5 inp_RF_K.txt > dat_K_PRB_ED.txt

gawk -f awk_prb_paper.awk 1 inp_RF_S.txt > dat_S_PRB_ND.txt
gawk -f awk_prb_paper.awk 2 inp_RF_S.txt > dat_S_PRB_LN.txt
gawk -f awk_prb_paper.awk 3 inp_RF_S.txt > dat_S_PRB_GD.txt
gawk -f awk_prb_paper.awk 4 inp_RF_S.txt > dat_S_PRB_WD.txt
gawk -f awk_prb_paper.awk 5 inp_RF_S.txt > dat_S_PRB_ED.txt



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



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="Converted variate [y(p)=z(p)]"
set pname=Normal distribution
set inp_1=dat_M_PRB_ND.txt
set fig=fig_M_PRB_ND.eps
call bat_plot
set inp_1=dat_K_PRB_ND.txt
set fig=fig_K_PRB_ND.eps
call bat_plot
set inp_1=dat_S_PRB_ND.txt
set fig=fig_S_PRB_ND.eps
call bat_plot



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="Converted variate [y(p)=z(p)]"
set pname=Log-Normal distribution
set inp_1=dat_M_PRB_LN.txt
set fig=fig_M_PRB_LN.eps
call bat_plot
set inp_1=dat_K_PRB_LN.txt
set fig=fig_K_PRB_LN.eps
call bat_plot
set inp_1=dat_S_PRB_LN.txt
set fig=fig_S_PRB_LN.eps
call bat_plot



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="Converted variate [y(p)=-ln{-ln(p)}]"
set pname=Gumbel distribution
set inp_1=dat_M_PRB_GD.txt
set fig=fig_M_PRB_GD.eps
call bat_plot
set inp_1=dat_K_PRB_GD.txt
set fig=fig_K_PRB_GD.eps
call bat_plot
set inp_1=dat_S_PRB_GD.txt
set fig=fig_S_PRB_GD.eps
call bat_plot



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="Converted variate [y(p)=ln{-ln(1-p)}]"
set pname=Weibull distribution
set inp_1=dat_M_PRB_WD.txt
set fig=fig_M_PRB_WD.eps
call bat_plot
set inp_1=dat_K_PRB_WD.txt
set fig=fig_K_PRB_WD.eps
call bat_plot
set inp_1=dat_S_PRB_WD.txt
set fig=fig_S_PRB_WD.eps
call bat_plot



rem ********************************************
rem Exponential distribution plotting
rem ********************************************
set range0=0/1/0.0010005/6.90776
set range1=0/500/0.0010005/6.90776
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="Converted variate [y(p)=-ln(1-p)]"
set pname=Exponential distribution
set inp_1=dat_M_PRB_ED.txt
set fig=fig_M_PRB_ED.eps
call bat_plot
set inp_1=dat_K_PRB_ED.txt
set fig=fig_K_PRB_ED.eps
call bat_plot
set inp_1=dat_S_PRB_ED.txt
set fig=fig_S_PRB_ED.eps
call bat_plot



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=
set pname=

del .gmt*
del _*