# data for legend
cat << EOT > _linp.txt
10 $ 10 $ tr $ 0.0 $ 0.0
10 $ 0.5 $ 6
Ground response (SR=100%) $ 1.5p            $ ---  $ --- $ ---
Ground response (SR=50%)  $ 1.5p,5_3:0      $ ---  $ --- $ ---
Tunnel lining             $ 1.5p            $ C0.2 $ 255 $ 1.0p
Limit supprt pressure     $ 1.0p,10_3_3_3:0 $ ---  $ --- $ ---
Plastic zone thickness    $ 1.5p,2_2:0      $ ---  $ --- $ ---
EOT
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _linp.txt > _lg.txt

gmt set FONT_ANNOT_PRIMARY 12p
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
xlabel='Displacement u@-a@- (mm)'
ylabel='Support pressure P@-i@- (MPa)'
scale=10/10
yga=a0.5f0.1
ylabelE='Plastic zone thickness (m)'
ygaE=a0.5f0.1

inp=out_A.txt
fig=fig_A.eps
range=0/5/0/2.5
rangeE=0/5/0/4
xga=a1f0.5
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
#Release rate 100%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$5}' $inp | gmt psxy -R -J -W1.5 -K -O >> $fig
#Release rate 50%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $1,$2}' $inp | gmt psxy -R -J -W1.5,5_3:0 -K -O >> $fig
#Sopport pressure
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -SC0.2 -W1 -G255 -N -K -O >> $fig
#Limit support pressure
awk 'NR==1{printf "0 %g 100 %g\n",$6,$6}' $inp | gmt psxy -R -JX  -Sv0.2+s -W1.0,10_3_3_3:0 -K -O  >> $fig
#Point a
awk '{if(NR==1){n=$1};if(n+3==NR)printf "%g %g 10p 0 LM a",$1,$2+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Point b
awk '{x=$1;y=$2}END{printf "%g %g 10p 0 LM b",x,y+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Plastic zone thickness
gmt psbasemap -R$rangeE -JX$scale -Bx$xga -By$ygaE+l"$ylabelE" -BE -K -O >> $fig
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$6/1000}' $inp | gmt psxy -R -J -W1.5,2_2:0 -K -O >> $fig
#
awk 'NR==1{x=$2;y=0.0;printf "0.55 0.72 10p 0 LM a(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$3;y=$4 ;printf "0.55 0.67 10p 0 LM b(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$6;y=$7 ;printf "0.55 0.62 10p 0 LM L.S.P=%.3fMPa(%.1fdays)",x,y}' $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
#
gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O << EOT >> $fig
0.95 1.02 10p 0 RB Support type: A
EOT
source _lg.txt


inp=out_B.txt
fig=fig_B.eps
range=0/10/0/2.5
rangeE=0/10/0/4
xga=a1f0.5
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
#Release rate 100%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$5}' $inp | gmt psxy -R -J -W1.5 -K -O >> $fig
#Release rate 50%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $1,$2}' $inp | gmt psxy -R -J -W1.5,5_3:0 -K -O >> $fig
#Sopport pressure
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -SC0.2 -W1 -G255 -N -K -O >> $fig
#Limit support pressure
awk 'NR==1{printf "0 %g 100 %g\n",$6,$6}' $inp | gmt psxy -R -JX  -Sv0.2+s -W1.0,10_3_3_3:0 -K -O  >> $fig
#Point a
awk '{if(NR==1){n=$1};if(n+3==NR)printf "%g %g 10p 0 LM a",$1,$2+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Point b
awk '{x=$1;y=$2}END{printf "%g %g 10p 0 LM b",x,y+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Plastic zone thickness
gmt psbasemap -R$rangeE -JX$scale -Bx$xga -By$ygaE+l"$ylabelE" -BE -K -O >> $fig
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$6/1000}' $inp | gmt psxy -R -J -W1.5,2_2:0 -K -O >> $fig
#
awk 'NR==1{x=$2;y=0.0;printf "0.55 0.72 10p 0 LM a(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$3;y=$4 ;printf "0.55 0.67 10p 0 LM b(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$6;y=$7 ;printf "0.55 0.62 10p 0 LM L.S.P=%.3fMPa(%.1fdays)",x,y}' $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
#
gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O << EOT >> $fig
0.95 1.02 10p 0 RB Support type: B
EOT
source _lg.txt


inp=out_C.txt
fig=fig_C.eps
range=0/20/0/2.5
rangeE=0/20/0/4
xga=a5f1
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
#Release rate 100%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$5}' $inp | gmt psxy -R -J -W1.5 -K -O >> $fig
#Release rate 50%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $1,$2}' $inp | gmt psxy -R -J -W1.5,5_3:0 -K -O >> $fig
#Sopport pressure
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -SC0.2 -W1 -G255 -N -K -O >> $fig
#Limit support pressure
awk 'NR==1{printf "0 %g 100 %g\n",$6,$6}' $inp | gmt psxy -R -JX  -Sv0.2+s -W1.0,10_3_3_3:0 -K -O  >> $fig
#Point a
awk '{if(NR==1){n=$1};if(n+3==NR)printf "%g %g 10p 0 LM a",$1,$2+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Point b
awk '{x=$1;y=$2}END{printf "%g %g 10p 0 LM b",x,y+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Plastic zone thickness
gmt psbasemap -R$rangeE -JX$scale -Bx$xga -By$ygaE+l"$ylabelE" -BE -K -O >> $fig
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$6/1000}' $inp | gmt psxy -R -J -W1.5,2_2:0 -K -O >> $fig
#
awk 'NR==1{x=$2;y=0.0;printf "0.55 0.72 10p 0 LM a(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$3;y=$4 ;printf "0.55 0.67 10p 0 LM b(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$6;y=$7 ;printf "0.55 0.62 10p 0 LM L.S.P=%.3fMPa(%.1fdays)",x,y}' $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
#
gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O << EOT >> $fig
0.95 1.02 10p 0 RB Support type: C
EOT
source _lg.txt


inp=out_D.txt
fig=fig_D.eps
range=0/30/0/2.5
rangeE=0/30/0/4
xga=a5f1
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
#Release rate 100%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$5}' $inp | gmt psxy -R -J -W1.5 -K -O >> $fig
#Release rate 50%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $1,$2}' $inp | gmt psxy -R -J -W1.5,5_3:0 -K -O >> $fig
#Sopport pressure
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -SC0.2 -W1 -G255 -N -K -O >> $fig
#Limit support pressure
awk 'NR==1{printf "0 %g 100 %g\n",$6,$6}' $inp | gmt psxy -R -JX  -Sv0.2+s -W1.0,10_3_3_3:0 -K -O  >> $fig
#Point a
awk '{if(NR==1){n=$1};if(n+3==NR)printf "%g %g 10p 0 LM a",$1,$2+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Point b
awk '{x=$1;y=$2}END{printf "%g %g 10p 0 LM b",x,y+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Plastic zone thickness
gmt psbasemap -R$rangeE -JX$scale -Bx$xga -By$ygaE+l"$ylabelE" -BE -K -O >> $fig
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$6/1000}' $inp | gmt psxy -R -J -W1.5,2_2:0 -K -O >> $fig
#
awk 'NR==1{x=$2;y=0.0;printf "0.55 0.72 10p 0 LM a(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$3;y=$4 ;printf "0.55 0.67 10p 0 LM b(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$6;y=$7 ;printf "0.55 0.62 10p 0 LM L.S.P=%.3fMPa(%.1fdays)",x,y}' $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
#
gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O << EOT >> $fig
0.95 1.02 10p 0 RB Support type: D
EOT
source _lg.txt


inp=out_E.txt
fig=fig_E.eps
range=0/100/0/2.5
rangeE=0/100/0/4
xga=a20f5
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWS -P -K > $fig
#Release rate 100%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$5}' $inp | gmt psxy -R -J -W1.5 -K -O >> $fig
#Release rate 50%
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $1,$2}' $inp | gmt psxy -R -J -W1.5,5_3:0 -K -O >> $fig
#Sopport pressure
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{if(NR==1){n=$1};if(n+3<=NR)print $1,$2}' $inp        | gmt psxy -R -J -SC0.2 -W1 -G255 -N -K -O >> $fig
#Limit support pressure
awk 'NR==1{printf "0 %g 100 %g\n",$6,$6}' $inp | gmt psxy -R -JX  -Sv0.2+s -W1.0,10_3_3_3:0 -K -O  >> $fig
#Point a
awk '{if(NR==1){n=$1};if(n+3==NR)printf "%g %g 10p 0 LM a",$1,$2+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Point b
awk '{x=$1;y=$2}END{printf "%g %g 10p 0 LM b",x,y+0.07}' $inp | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
#Plastic zone thickness
gmt psbasemap -R$rangeE -JX$scale -Bx$xga -By$ygaE+l"$ylabelE" -BE -K -O >> $fig
awk '{if(NR==1){n=$1};if(2<=NR&&NR<=n+2)print $4,$6/1000}' $inp | gmt psxy -R -J -W1.5,2_2:0 -K -O >> $fig
#
awk 'NR==1{x=$2;y=0.0;printf "0.55 0.72 10p 0 LM a(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$3;y=$4 ;printf "0.55 0.67 10p 0 LM b(ua,Pi)=(%.3f,%.3f)",x,y}'    $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
awk 'NR==1{x=$6;y=$7 ;printf "0.55 0.62 10p 0 LM L.S.P=%.3fMPa(%.1fdays)",x,y}' $inp | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
#
gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O << EOT >> $fig
0.95 1.02 10p 0 RB Support type: E
EOT
source _lg.txt

rm _*