rem ============================================ rem Comparison of Climate rem ============================================ rem ***** Temperature data ***** gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;T[i]=$2}END{for(i=1;i<=12;i++){print M[i]-0.5,T[i]}}" inp_ASAHIKAWA.txt >_inpT_ASAHIKAWA.txt gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;T[i]=$2}END{for(i=1;i<=12;i++){print M[i]-0.5,T[i]}}" inp_AKITA.txt >_inpT_AKITA.txt gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;T[i]=$2}END{for(i=1;i<=12;i++){print M[i]-0.5,T[i]}}" inp_TOKYO.txt >_inpT_TOKYO.txt gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;T[i]=$2}END{for(i=1;i<=12;i++){print M[i]-0.5,T[i]}}" inp_NAHA.txt >_inpT_NAHA.txt rem ***** Rainfall data ***** gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;R[i]=$3}END{for(i=1;i<=12;i++){print M[i]-0.5,R[i]}}" inp_ASAHIKAWA.txt >_inpR_ASAHIKAWA.txt gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;R[i]=$3}END{for(i=1;i<=12;i++){print M[i]-0.5,R[i]}}" inp_AKITA.txt >_inpR_AKITA.txt gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;R[i]=$3}END{for(i=1;i<=12;i++){print M[i]-0.5,R[i]}}" inp_TOKYO.txt >_inpR_TOKYO.txt gawk "BEGIN{i=0}3<=NR{i=i+1;M[i]=$1;R[i]=$3}END{for(i=1;i<=12;i++){print M[i]-0.5,R[i]}}" inp_NAHA.txt >_inpR_NAHA.txt rem ***** Execute GMT command ***** set fig_out=fig_gmt_climate.eps gmtset ANOT_FONT_SIZE_PRIMARY 14 gmtset ANOT_OFFSET_PRIMARY 0.2c gmtset LABEL_FONT_SIZE 14 gmtset LABEL_OFFSET 0.2c gmtset HEADER_FONT_SIZE 16 gmtset HEADER_OFFSET 0.3c gmtset TICK_LENGTH 0c gmtset TICK_PEN 0.5p rem ***** Title ***** echo 14 19.5 16 0 0 MC Comparison of Climate | pstext -R0/28/0/20 -JX28/20 -X0 -Y0 -N -K > %fig_out% echo 0.5 0.5 > _inp1.txt echo 1.5 0.5 >> _inp1.txt echo 0.5 1.5 > _inp2.txt echo 1.5 1.5 >> _inp2.txt echo 1.0 1.5 > _inp3.txt rem ***** Legend ***** psbasemap -R0/5/0/2 -JX5/1.2 -B/:: -X11.5 -Y17.5 -K -O>> %fig_out% psxy _inp1.txt -R -J -W40,cyan -K -O >> %fig_out% psxy _inp2.txt -R -J -W5,red -K -O >> %fig_out% psxy _inp3.txt -R -J -SC0.2,red -W1,red -Gred -K -O >> %fig_out% echo 2.0 0.5 14 0 0 ML Rainfall | pstext -R -J -N -K -O >> %fig_out% echo 2.0 1.5 14 0 0 ML Temperature | pstext -R -J -N -K -O >> %fig_out% rem ***** Retuen to origin ***** echo -1 -1 | psxy -R0/1/0/1 -JX1/1 -Sp -W255 -X-11.5 -Y-17.5 -K -O >> %fig_out% rem ***** Plot-1 ***** set datinpR=_inpR_TOKYO.txt set datinpT=_inpT_TOKYO.txt set xx=2.5 set yy=1.5 set locname=Tokyo call bat_climate_imp rem ***** Plot-2 ***** set datinpR=_inpR_NAHA.txt set datinpT=_inpT_NAHA.txt set xx=14 set yy=0 set locname=Naha call bat_climate_imp rem ***** Plot-3 ***** set datinpR=_inpR_ASAHIKAWA.txt set datinpT=_inpT_ASAHIKAWA.txt set xx=-14 set yy=9 set locname=Asahikawa call bat_climate_imp rem ***** Plot-4 ***** set datinpR=_inpR_AKITA.txt set datinpT=_inpT_AKITA.txt set xx=14 set yy=0 set locname=Akita call bat_climate_imp rem ***** End of plot ***** echo 0 0 | psxy -R -J -Sp -O >> %fig_out% rem ***** Delete work files ***** set fig_out= del _* del .gmt*