rem ============================
rem Time-sequence diagram (Temperature & Rainfall)
rem ============================
set fig=fig_gmt_tsd_TR.eps
gmtset PLOT_DATE_FORMAT o 
gmtset TIME_FORMAT_PRIMARY abbreviated
gmtset INPUT_DATE_FORMAT yyyy/mm/dd
gmtset ANNOT_FONT_SIZE_PRIMARY 7
gmtset ANNOT_FONT_SIZE_SECONDARY 14
gmtset LABEL_FONT_SIZE 14
gmtset HEADER_FONT_SIZE 18
gmtset HEADER_OFFSET -0.3c
gmtset TICK_LENGTH -0.1c
rem ---------------------------------------
rem Create batch file for legend
echo 20 $ 10 $ tr $ 0.5 $ 0.3 > inp_legend_tsd_TR.txt
echo 10 $ 0.5 $ 6 >> inp_legend_tsd_TR.txt
echo Temperature  $ 5/255/0/0    $ ---  $ ---  $ --- >> inp_legend_tsd_TR.txt
echo Rainfall     $ 20/0/255/255 $ ---  $ ---  $ --- >> inp_legend_tsd_TR.txt
gawk -f \WANtaroHP_plot\gmt_basic\awk_klegend.awk inp_legend_tsd_TR.txt > legend1.bat
rem ---------------------------------------
psbasemap -R1999-01-01T/2010-01-01T/-50/50 -JX20/10 -U"%fig%" -Bsa1Yg1Y/S -Bpa3Of1o/Sn:."Temperature and Rainfall in Maebashi": -K > %fig%
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $1,$2}" inp_TR_maebashi.txt | psxy -R -J -Bs/a20f20g20:"Temperature (@%%12%%\260@%%%%C)":W -W5,red -O -K >> %fig%
gawk "BEGIN{FS=\",\"}{if(2<=NR)print $1,$5}" inp_TR_maebashi.txt | psxy -R1999-01-01T/2010-01-01T/0/1000 -J -Gcyan -W2,magenta -Sb30u -Bs/a100f100:"Rainfall (mm)":E -O -K >> %fig%
call legend1.bat
set fig=
del _*
del .gmt*