rem gcc -o gccTES1_LP3.exe gccTES1_LP3.c
rem gccTES1_LP3 inp_PRB.txt out_TES1_5.csv


rem ********************************************
rem Drawing by GMT
rem ********************************************
set range=1/10000000/0.98/1.10
set scale=12l/10
set xga=f3a1p
set yga=f0.01a0.02
set xlabel="Coefficient (b)"
set ylabel="Ratio of quantile (WH/Gamma)"
gmtset ANOT_FONT_SIZE 12
gmtset LABEL_FONT_SIZE 12
gmtset TICK_LENGTH -0.2c

rem ===============================
rem Plotting
rem ===============================
set fig_out=fig_TES1_LP3.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out%
set inp_1=out_TES1_1.csv
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $9,$10}" %inp_1% | psxy -R -J -B -SC0.1 -P -O -K >> %fig_out%
set inp_1=out_TES1_2.csv
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $9,$10}" %inp_1% | psxy -R -J -B -SC0.1 -P -O -K >> %fig_out%
set inp_1=out_TES1_3.csv
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $9,$10}" %inp_1% | psxy -R -J -B -SC0.1 -P -O -K >> %fig_out%
set inp_1=out_TES1_4.csv
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $9,$10}" %inp_1% | psxy -R -J -B -SC0.1 -P -O -K >> %fig_out%
set inp_1=out_TES1_5.csv
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $9,$10}" %inp_1% | psxy -R -J -B -SC0.1 -P -O >> %fig_out%