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
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_div20.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==851,NR==1351{print $2,$3}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$4}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$8}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$11}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$13}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$15}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$18}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$22}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1351{print $2,$23}' out_22_div20_01m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
# ********************************************
# (delta=1days) broken line+white circle
# ********************************************
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$3}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$4}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$8}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$11}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$13}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$15}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$18}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$22}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$23}' out_22_div20_10m.csv | gmt psxy -JX -R -W1,- -K -O >> $fig
#
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$3}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$4}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$8}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$11}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$13}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$15}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$18}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$22}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==901{print $2,$23}' out_22_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
# ********************************************
# (delta=5days) dotted line+black triangle
# ********************************************
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$3}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$4}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$8}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$11}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$13}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$15}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$18}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$22}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$23}' out_22_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
#
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$3}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$4}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$8}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$11}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$13}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$15}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$18}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$22}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==861{print $2,$23}' out_22_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig
# ********************************************
# legend (pstext)
# ********************************************
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,-0.50)",$3}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,-0.45)",$4}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,-0.25)",$8}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,-0.10)",$11}' out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,0.00)",$13}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,0.10)",$15}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,0.25)",$18}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,0.45)",$22}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==901{printf "52 %.2f 10 0 ML (x,y)=(0,0.50)",$23}'  out_22_div20_10m.csv | gmt pstext -JX -R -F+f+a+j -N -K -O >> $fig

source _legend.txt