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.5c
gmt set MAP_TICK_LENGTH_PRIMARY -0.2c

cat << EOT > _inp.txt
12 $ 12 $ tr $ 0.5 $ 0.5
12 $ 0.5 $ 8
Structure Line $ 2,#000000       $ --- $ ---
Water Level    $ 2,#0000ff       $ --- $ ---
Critical Depth $ 2,#0000ff,3_3:0 $ --- $ ---
EOT
awk -f /Users/kk/DATA/GMT/awk_klegend.awk _inp.txt > _leg.txt

rangeSW=-1/7/802/810
scale=12/12
afgS=a1f0.5
afgW=a1f0.5
labelS='Horizontal distance (m)'
labelW='Elevation (EL.m)'
fig=fig_gra_1.eps
gmt psbasemap -R$rangeSW -JX$scale -Bx$afgS+l"$labelS" -By$afgW+l"$labelW" -BWSen -P -K > $fig
gmt psxy -R -J -Gp300/62:F000000B- -W2 -K -O << EOT >> $fig
-0.5 803.5
-0.5 804.0
 0.0 804.0
 0.0 805.0
 0.5 805.0
 4.5 804.0
 7.0 804.0
 7.0 803.5
EOT
gmt psxy -R -J -Gp300/62:F000000B- -W2 -K -O << EOT >> $fig
5.0 808.0
5.0 806.5
7.0 806.5
7.0 807.0
5.5 807.0
5.5 808.0
EOT
gmt psxy -R -J -W2 -K -O << EOT >> $fig
-1.0 804.0
-0.5 804.0
EOT
gmt psxy -R -J -W2 -K -O << EOT >> $fig
0.0 805.0
0.0 808.0
5.0 808.0
EOT
awk '{print $1,$3}' out.txt | gmt psxy -R -J -W0.5,#000000 -K -O >> $fig
awk '{print $1,$3+$4}' out.txt | gmt psxy -R -J -W2,#0000ff -K -O >> $fig
awk '{print $1,$3+$5}' out.txt | gmt psxy -R -J -W2,#0000ff,3_3:0 -K -O >> $fig
# Streight line
gmt psxy -R -JX -Sv0.2+s -W0.5 -N -K -O << EOT >> $fig
0.0 803.5 0.0 802.4
0.5 803.5 0.5 802.9
4.5 803.5 4.5 802.9
5.5 803.5 5.5 802.4
0.5 805.0 2.5 805.0
7.0 807.0 7.4 807.0
7.0 806.5 7.4 806.5
7.0 804.0 7.4 804.0
7.0 803.5 7.4 803.5
EOT
# One side arrow
gmt set MAP_VECTOR_SHAPE 0
gmt psxy -R -JX -Sv0.2+s+e+a60 -W1 -G0 -N -O -K << EOT >> $fig
5.5 803.0 6.5 803.0
EOT
# Both side arrow
gmt set MAP_VECTOR_SHAPE 1
gmt psxy -R -JX -Sv0.2+s+b+e+a45 -W0.5 -G0 -N -O -K << EOT >> $fig
0.0 803.0 0.5 803.0
0.5 803.0 4.5 803.0
4.5 803.0 5.5 803.0
0.0 802.5 5.5 802.5
0.0 807.6 5.0 807.6
7.3 807.0 7.3 806.5
7.3 806.5 7.3 804.0
7.3 804.0 7.3 803.5
EOT
gmt psxy -R -JX -Sv0.2+s+e+a60 -W1 -G0 -N -O -K << EOT >> $fig
6.0 807.6 5.5 807.6
EOT
# Text
gmt pstext -R -J -F+f+a+j -N -K -O << EOT >> $fig
0.25 803.1 12p,Helvetica 0 BC 0.500
2.50 803.1 12p,Helvetica 0 BC 4.000 (i=1:4.0)
5.00 803.1 12p,Helvetica 0 BC 1.000
6.25 803.1 12p,Helvetica 0 BC i=0.010
2.75 802.6 12p,Helvetica 0 BC Inlet (L=5.5m)
2.50 805.1 12p,Helvetica 0 BR EL.805.0
2.75 808.1 12p,Helvetica 0 BC EL.808.0
5.00 804.1 12p,Helvetica 0 BC EL.804.0
2.50 807.5 12p,Helvetica 0 TC 5.000
5.70 807.5 12p,Helvetica 0 TL 0.500
7.50 806.75 12p,Helvetica 90 TC 0.500
7.50 805.25 12p,Helvetica 90 TC 2.500
7.50 803.75 12p,Helvetica 90 TC 0.500
EOT
source _leg.txt
#echo '-1 802' | gmt psxy -R -J -Sp -O >> $fig


cat << EOT > _inp.txt
17 $ 10 $ tr $ 0.5 $ 0.5
12 $ 0.5 $ 8
Structure Line $ 1,#000000       $ --- $ ---
Water Level    $ 1,#0000ff       $ --- $ ---
Critical Depth $ 1,#0000ff,3_3:0 $ --- $ ---
EOT
awk -f /Users/kk/DATA/GMT/awk_klegend.awk _inp.txt > _leg.txt
rangeSW=0/1700/790/810
scale=17/10
afgS=a100f10g100
afgW=a5f1g5
labelS='Horizontal distance (m)'
labelW='Elevation (EL.m)'
fig=fig_gra_2.eps
gmt psbasemap -R$rangeSW -JX$scale -Bx$afgS+l"$labelS" -By$afgW+l"$labelW" -BWSen -P -K > $fig
gmt psxy -R -J -W1 -K -O << EOT >> $fig
0.0 808.0
5.5 808.0
5.5 806.5
391.0 802.505
1607.0 795.503
EOT
awk '{print $1,$3}' out.txt | gmt psxy -R -J -W1,#000000 -K -O >> $fig
awk '{print $1,$3+$4}' out.txt | gmt psxy -R -J -W1,#0000ff -K -O >> $fig
awk '{print $1,$3+$5}' out.txt | gmt psxy -R -J -W1,#0000ff,3_3:0 -K -O >> $fig
# Streight line
gmt psxy -R -JX -Sv0.2+s -W0.5 -N -K -O << EOT >> $fig
5.5 804.0 5.5 791.5
391.0 802.505 391.0 791.5
1607.0 795.503 1607.0 791.5
EOT
# Both side arrow
gmt set MAP_VECTOR_SHAPE 1
gmt psxy -R -JX -Sv0.2+s+b+e+a45 -W0.5 -G0 -N -O -K << EOT >> $fig
5.5 792.0 391.0 792.0
391.0 792.0 1607.0 792.0
EOT
# Text
gmt pstext -R -J -F+f+a+j -N -K -O << EOT >> $fig
203.75 792.3 12p,Helvetica 0 BC 386.000 (i=0.010)
203.75 791.8 12p,Helvetica 0 TC (except Inlet)
999.00 792.3 12p,Helvetica 0 BC 1216.000 (i=0.006)
391.0 803.0 12p,Helvetica 90 ML Upper Edge
431.0 804.0 12p,Helvetica 90 ML of Reservoir
1607.0 796.0 12p,Helvetica 90 ML Stert Point of Spoil Bank
EOT
source _leg.txt