gmt set FONT_ANNOT_PRIMARY 12
gmt set MAP_ANNOT_OFFSET_PRIMARY 0.3c
gmt set FONT_LABEL 12p
gmt set MAP_LABEL_OFFSET 0.3c
gmt set MAP_TICK_LENGTH_PRIMARY -0.2c

cat << EOT > _inp.txt
12 $ 8 $ br $ 0.5 $ 0.5
10 $ 0.5 $ 4
Theory         $ 1     $ ---   $ ---  $ ---
FEM(delta=1.0) $ ---   $ C0.2  $ 255  $ 1
EOT
awk 'BEGIN{for(t=0;t<=50;t=t+0.1){print t,20+40*(1-exp(-0.2*t))}}' > z_theory.txt
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _inp.txt > _legend.txt
range=0/50/20/70
scale=12/8
bax=a10f5
bay=a10f5
strx="time (hour)"
stry="Temperature (@%12%\260@%%C)"
psize=0.150
fig=fig_gmt_00_div4.eps
gmt psbasemap -R$range  -JX$scale -Bx$bax+l"$strx" -By$bay+l"$stry" -BSWne -P -K > $fig
gmt psxy z_theory.txt -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==101{print $2,$5}' out_test00_div4.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
source _legend.txt


cat << EOT > _inp.txt
12 $ 8 $ tr $ 0.5 $ 0.5
10 $ 0.5 $ 4
delta=1 (x,y)=(0,0.00) $ ---   $ C0.2  $ 255  $ 1
delta=1 (x,y)=(0,0.25) $ ---   $ S0.2  $ 255  $ 1
delta=1 (x,y)=(0,0.50) $ ---   $ T0.2  $ 255  $ 1
EOT
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _inp.txt > _legend.txt
range=0/100/0/60
scale=12/8
bax=a10f5
bay=a10f5
strx="time (hour)"
stry="Temperature (@%12%\260@%%C)"
psize=0.150
fig=fig_gmt_11_div4.eps
gmt psbasemap -R$range  -JX$scale -Bx$bax+l"$strx" -By$bay+l"$stry" -BSWne -P -K > $fig
awk 'BEGIN{FS=","}NR==51,NR==151{print $2,$5}' out_test11_div4.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==151{print $2,$4}' out_test11_div4.csv | gmt psxy -JX -R -SS$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==151{print $2,$3}' out_test11_div4.csv | gmt psxy -JX -R -ST$psize -W1 -G255 -N -K -O >> $fig
source _legend.txt


cat << EOT > _inp.txt
10 $ 10 $ ot $ 0.5 $ 0.5
10 $ 0.5 $ 8
delta=0.1 $ 1   $ ---   $ ---  $ ---
delta=1.0 $ 1,- $ C0.2  $ 255  $ 1
delta=5.0 $ 1,. $ T0.2  $ 0    $ ---
EOT
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _inp.txt > _legend.txt
range=0/50/-5/45
scale=10/10
bax=a10f5
bay=a5f1
strx="time (hour)"
stry="Temperature (@%12%\260@%%C)"
psize=0.150
fig=fig_gmt_22_div4.eps
gmt psbasemap -R$range  -JX$scale -Bx$bax+l"$strx" -By$bay+l"$stry" -BSWne -P -K > $fig
# ********************************************
# (delta=0.1days) solid line
# ********************************************
awk 'BEGIN{FS=","}NR==51,NR==5511{print $2,$3}' out_test22_div4_01.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==5511{print $2,$4}' out_test22_div4_01.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==5511{print $2,$5}' out_test22_div4_01.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==5511{print $2,$6}' out_test22_div4_01.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==5511{print $2,$7}' out_test22_div4_01.csv | gmt psxy -JX -R -W1 -K -O >> $fig
# ********************************************
# (delta=1days) white circle
# ********************************************
awk 'BEGIN{FS=","}NR==51,NR==101{print $2,$3}' out_test22_div4_10.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==101{print $2,$4}' out_test22_div4_10.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==101{print $2,$5}' out_test22_div4_10.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==101{print $2,$6}' out_test22_div4_10.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==101{print $2,$7}' out_test22_div4_10.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
# ********************************************
# (delta=5days) black triangle
# ********************************************
awk 'BEGIN{FS=","}NR==51,NR==61{print $2,$3}' out_test22_div4_50.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==61{print $2,$4}' out_test22_div4_50.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==61{print $2,$5}' out_test22_div4_50.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==61{print $2,$6}' out_test22_div4_50.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==51,NR==61{print $2,$7}' out_test22_div4_50.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
# ********************************************
# legend
# ********************************************
awk 'BEGIN{FS=","}NR==101{printf "52 %.2f 10 0 ML (x,y)=(0,-0.50)",$3}' out_test22_div4_10.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==101{printf "52 %.2f 10 0 ML (x,y)=(0,-0.25)",$4}' out_test22_div4_10.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==101{printf "52 %.2f 10 0 ML (x,y)=(0,0.00)",$5}'  out_test22_div4_10.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==101{printf "52 %.2f 10 0 ML (x,y)=(0,0.25)",$6}'  out_test22_div4_10.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==101{printf "52 %.2f 10 0 ML (x,y)=(0,0.50)",$7}'  out_test22_div4_10.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
source _legend.txt