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_leg.txt
12 $ 8 $ tr $ 0.5 $ 0.5
10 $ 0.5 $ 4
delta=1 (x,y)=(0,0.00) $ 1   $ C0.2  $ 255  $ 1
delta=1 (x,y)=(0,0.25) $ 1   $ S0.2  $ 255  $ 1
delta=1 (x,y)=(0,0.45) $ 1   $ D0.2  $ 255  $ 1
delta=1 (x,y)=(0,0.50) $ 1   $ T0.2  $ 255  $ 1
delta=5 (x,y)=(0,0.00) $ 1,. $ C0.2  $ 0  $ ---
delta=5 (x,y)=(0,0.25) $ 1,. $ S0.2  $ 0  $ ---
delta=5 (x,y)=(0,0.45) $ 1,. $ D0.2  $ 0  $ ---
delta=5 (x,y)=(0,0.50) $ 1,. $ T0.2  $ 0  $ ---
EOT
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _inp_leg.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_div20.eps
gmt psbasemap -R$range  -JX$scale -Bx$bax+l"$strx" -By$bay+l"$stry" -BSWne -P -K > $fig
# ********************************************
# (delta=1days) broken line+white circle
# ********************************************
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$13}' out_11_div20_10m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$18}' out_11_div20_10m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$22}' out_11_div20_10m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$23}' out_11_div20_10m.csv | gmt psxy -JX -R -W1 -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==951{print $2,$13}' out_11_div20_10m.csv | gmt psxy -JX -R -SC$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==951{print $2,$18}' out_11_div20_10m.csv | gmt psxy -JX -R -SS$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==951{print $2,$22}' out_11_div20_10m.csv | gmt psxy -JX -R -SD$psize -W1 -G255 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==951{print $2,$23}' out_11_div20_10m.csv | gmt psxy -JX -R -ST$psize -W1 -G255 -N -K -O >> $fig
# ********************************************
# (delta=5days) dotted line+black triangle
# ********************************************
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$13}' out_11_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$18}' out_11_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$22}' out_11_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==1851{print $2,$23}' out_11_div20_50m.csv | gmt psxy -JX -R -W1,. -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==871{print $2,$13}' out_11_div20_50m.csv | gmt psxy -JX -R -SC$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==871{print $2,$18}' out_11_div20_50m.csv | gmt psxy -JX -R -SS$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==871{print $2,$22}' out_11_div20_50m.csv | gmt psxy -JX -R -SD$psize -G0 -N -K -O >> $fig
awk 'BEGIN{FS=","}NR==851,NR==871{print $2,$23}' out_11_div20_50m.csv | gmt psxy -JX -R -ST$psize -G0 -N -K -O >> $fig

source _legend.txt