rem ************************** rem ヒストグラム描画 rem ************************** set wd=20 set xlabel="Maximum daily rainfall (mm/day)" set range=0/1500/0/3000 set xga=f50a250 set yga=f100a500 set inp_org=inp_BH_GUM.txt set out_text=out_BH_GUM.txt set fig_out=fig_BH_GUM.eps 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 12 0 0 MR n=%%g\n\",$1}" _tempdat.txt > _txtdat.txt gawk "{printf \"9.5 9.0 12 0 0 MR min=%%g\n\",$2}" _tempdat.txt >> _txtdat.txt gawk "{printf \"9.5 8.5 12 0 0 MR max=%%g\n\",$3}" _tempdat.txt >> _txtdat.txt rem ==================================== rem ===Command========================== gmtset ANOT_FONT_SIZE 14 gmtset LABEL_FONT_SIZE 14 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_out% pstext _txtdat.txt -R0/10/0/10 -J -N -P -O -K >> %fig_out% gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.2 12 0 0 BL %%s\n\",$0}" %inp_org% | pstext -R -J -N -P -O >> %fig_out% rem ==================================== set inp_org=inp_BH_GEV.txt set out_text=out_BH_GEV.txt set fig_out=fig_BH_GEV.eps 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 12 0 0 MR n=%%g\n\",$1}" _tempdat.txt > _txtdat.txt gawk "{printf \"9.5 9.0 12 0 0 MR min=%%g\n\",$2}" _tempdat.txt >> _txtdat.txt gawk "{printf \"9.5 8.5 12 0 0 MR max=%%g\n\",$3}" _tempdat.txt >> _txtdat.txt rem ==================================== rem ===Command========================== gmtset ANOT_FONT_SIZE 14 gmtset LABEL_FONT_SIZE 14 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_out% pstext _txtdat.txt -R0/10/0/10 -J -N -P -O -K >> %fig_out% gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.2 12 0 0 BL %%s\n\",$0}" %inp_org% | pstext -R -J -N -P -O >> %fig_out% rem ==================================== set inp_org=inp_BH_SQR.txt set out_text=out_BH_SQR.txt set fig_out=fig_BH_SQR.eps 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 12 0 0 MR n=%%g\n\",$1}" _tempdat.txt > _txtdat.txt gawk "{printf \"9.5 9.0 12 0 0 MR min=%%g\n\",$2}" _tempdat.txt >> _txtdat.txt gawk "{printf \"9.5 8.5 12 0 0 MR max=%%g\n\",$3}" _tempdat.txt >> _txtdat.txt rem ==================================== rem ===Command========================== gmtset ANOT_FONT_SIZE 14 gmtset LABEL_FONT_SIZE 14 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_out% pstext _txtdat.txt -R0/10/0/10 -J -N -P -O -K >> %fig_out% gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.2 12 0 0 BL %%s\n\",$0}" %inp_org% | pstext -R -J -N -P -O >> %fig_out% rem ==================================== set inp_org=inp_BH_LN3.txt set out_text=out_BH_LN3.txt set fig_out=fig_BH_LN3.eps 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 12 0 0 MR n=%%g\n\",$1}" _tempdat.txt > _txtdat.txt gawk "{printf \"9.5 9.0 12 0 0 MR min=%%g\n\",$2}" _tempdat.txt >> _txtdat.txt gawk "{printf \"9.5 8.5 12 0 0 MR max=%%g\n\",$3}" _tempdat.txt >> _txtdat.txt rem ==================================== rem ===Command========================== gmtset ANOT_FONT_SIZE 14 gmtset LABEL_FONT_SIZE 14 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_out% pstext _txtdat.txt -R0/10/0/10 -J -N -P -O -K >> %fig_out% gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.2 12 0 0 BL %%s\n\",$0}" %inp_org% | pstext -R -J -N -P -O >> %fig_out% rem ==================================== set inp_org=inp_BH_LP3.txt set out_text=out_BH_LP3.txt set fig_out=fig_BH_LP3.eps 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 12 0 0 MR n=%%g\n\",$1}" _tempdat.txt > _txtdat.txt gawk "{printf \"9.5 9.0 12 0 0 MR min=%%g\n\",$2}" _tempdat.txt >> _txtdat.txt gawk "{printf \"9.5 8.5 12 0 0 MR max=%%g\n\",$3}" _tempdat.txt >> _txtdat.txt rem ==================================== rem ===Command========================== gmtset ANOT_FONT_SIZE 14 gmtset LABEL_FONT_SIZE 14 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_out% pstext _txtdat.txt -R0/10/0/10 -J -N -P -O -K >> %fig_out% gawk "BEGIN{FS=\",\"}{if(NR==1)printf \"0.2 10.2 12 0 0 BL %%s\n\",$0}" %inp_org% | pstext -R -J -N -P -O >> %fig_out% rem ==================================== set inp_org= set out_text= set wd= set xlabel= set fig_out= set range= se xga= set yga= del .gmt* del _*