gcc -o gccWEI.exe gccWEI.c
gccWEI inp_PRB.txt out_WEI.txt



rem ********************************************
rem Drawing by GMT
rem ********************************************
set range=0/400/0/400
set scale=10/10
set xga=f50a100
set yga=f50a100
set xlabel="Observed data (mm/day)"
set ylabel="Calculated data (mm/day)"
set inp_1=out_WEI.txt
set fig_out=fig_WEI.eps
gmtset ANOT_FONT_SIZE 12
gmtset LABEL_FONT_SIZE 12
gmtset TICK_LENGTH -0.2c



gawk -f 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_out%
gawk "BEGIN{FS=\",\"}{if(6<=NR)print $2,$3}" %inp_1% | psxy -R -J -B -SC0.2 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(6<=NR)print $2,$4}" %inp_1% | psxy -R -J -B -ST0.2 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(6<=NR)print $2,$5}" %inp_1% | psxy -R -J -B -SS0.2 -P -O -K >> %fig_out%
echo 0 0 > _temp.txt
echo 400 400 >> _temp.txt
psxy _temp.txt -R -J -B -W1 -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==2)printf \"0.5 9.5 11 0 0 ML %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==3)printf \"0.5 9.0 11 0 0 ML %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==4)printf \"0.5 8.5 11 0 0 ML %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "END{printf \"0.5 10.1 11 0 0 BL %%s\n\",\"Weibull distribution\"}" %inp_1% | pstext -R0/10/0/10 -J -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 _*