rem bat_hist_bootstrap.bat set wd=20 set xlabel="Maximum daily rainfall (mm/day)" set range=0/1500/0/3000 set xga=f50a250 set yga=f100a500 rem ===Command========================== minmax %inp_org% -H1 > %out_text% pshistogram %inp_org% -W%wd% -IO -H1 >> %out_text% rem ==================================== rem ===Command========================== gawk "{if(NR==1){nd=$4;str0=$5}}END{sub(//,\"\",str0);split(str0,a,\"/\");print nd,a[1],a[2]}" %out_text% > _tempdat.txt rem ==================================== rem ===Command========================== gawk "{printf \"9.5 9.5 16 0 0 MR n=%%g\n\",$1}" _tempdat.txt > _txtdat.txt gawk "{printf \"9.5 9.0 16 0 0 MR min=%%g\n\",$2}" _tempdat.txt >> _txtdat.txt gawk "{printf \"9.5 8.5 16 0 0 MR max=%%g\n\",$3}" _tempdat.txt >> _txtdat.txt rem ==================================== rem ===Command========================== gmtset ANOT_FONT_SIZE 16 gmtset LABEL_FONT_SIZE 16 gmtset TICK_LENGTH -0.2c pshistogram %inp_org% -R%range% -JX15/10 -B%xga%:%xlabel%:/%yga%:"Frequency":WSen -H1 -W%wd% -L0.5p -G230 -Z0 -P -K > %fig% pstext _txtdat.txt -R0/10/0/10 -J -N -P -O -K >> %fig% gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.2 16 0 0 BL %%s\n\",$0}" %inp_org% | pstext -R -J -N -P -O >> %fig% rem ====================================