gcc -o gccSTA.exe gccSTA.c
gccSTA inp_PRB.txt out_PRB.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_PRB.txt
gmtset ANOT_FONT_SIZE 12
gmtset LABEL_FONT_SIZE 12
gmtset TICK_LENGTH -0.2c

rem ===============================
rem Plotting 1
rem ===============================
set fig_out=fig_PRB_GUM.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(9<=NR)print $2,$3}" %inp_1% | psxy -R -J -B -SC0.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 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.5 12 0 0 ML n=%%s\n\",$2}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.0 12 0 0 ML r=%%s\n\",$3}" %inp_1% | pstext -R -J -N -P -O >> %fig_out%
rem ===============================
rem Plotting 2
rem ===============================
set fig_out=fig_PRB_GEV.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(9<=NR)print $2,$4}" %inp_1% | psxy -R -J -B -SC0.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==3)printf \"0.5 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.5 12 0 0 ML n=%%s\n\",$2}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.0 12 0 0 ML r=%%s\n\",$4}" %inp_1% | pstext -R -J -N -P -O >> %fig_out%
rem ===============================
rem Plotting 3
rem ===============================
set fig_out=fig_PRB_SQR.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(9<=NR)print $2,$5}" %inp_1% | psxy -R -J -B -SC0.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==4)printf \"0.5 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.5 12 0 0 ML n=%%s\n\",$2}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.0 12 0 0 ML r=%%s\n\",$5}" %inp_1% | pstext -R -J -N -P -O >> %fig_out%
rem ===============================
rem Plotting 4
rem ===============================
set fig_out=fig_PRB_LN3.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(9<=NR)print $2,$6}" %inp_1% | psxy -R -J -B -SC0.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==5)printf \"0.5 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.5 12 0 0 ML n=%%s\n\",$2}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.0 12 0 0 ML r=%%s\n\",$6}" %inp_1% | pstext -R -J -N -P -O >> %fig_out%
rem ===============================
rem Plotting 5
rem ===============================
set fig_out=fig_PRB_LP3.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
gawk "BEGIN{FS=\",\"}{if(9<=NR)print $2,$7}" %inp_1% | psxy -R -J -B -SC0.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==6)printf \"0.5 10.1 12 0 0 BL %%s\n\",$0}" %inp_1% | pstext -R0/10/0/10 -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.5 12 0 0 ML n=%%s\n\",$2}" %inp_1% | pstext -R -J -N -P -O -K >> %fig_out%
gawk "BEGIN{FS=\",\"}{if(NR==8)printf \"0.5 9.0 12 0 0 ML r=%%s\n\",$7}" %inp_1% | pstext -R -J -N -P -O >> %fig_out%

del .gmt*
del _*