awk '{if($1==2.7)print $1,0,$1,$3}' out_sB.txt > _line_d.txt
awk '{if($1==2.7)print $1,$3,0,$3}' out_sB.txt >> _line_d.txt
awk '{if($1==3.3)print $1,0,$1,$3}' out_sC.txt >> _line_d.txt
awk '{if($1==3.3)print $1,$3,0,$3}' out_sC.txt >> _line_d.txt
awk '{if($1==4.5)print $1,0,$1,$3}' out_sD.txt >> _line_d.txt
awk '{if($1==4.5)print $1,$3,0,$3}' out_sD.txt >> _line_d.txt
awk '{if($1==6.8)print $1,0,$1,$3}' out_sE.txt >> _line_d.txt
awk '{if($1==6.8)print $1,$3,0,$3}' out_sE.txt >> _line_d.txt

awk '{if($1==2.7)print $2,0,$2,$3}' out_sB.txt > _line_u.txt
awk '{if($1==2.7)print $2,$3,0,$3}' out_sB.txt >> _line_u.txt
awk '{if($1==3.3)print $2,0,$2,$3}' out_sC.txt >> _line_u.txt
awk '{if($1==3.3)print $2,$3,0,$3}' out_sC.txt >> _line_u.txt
awk '{if($1==4.5)print $2,0,$2,$3}' out_sD.txt >> _line_u.txt
awk '{if($1==4.5)print $2,$3,0,$3}' out_sD.txt >> _line_u.txt
awk '{if($1==6.8)print $2,0,$2,$3}' out_sE.txt >> _line_u.txt
awk '{if($1==6.8)print $2,$3,0,$3}' out_sE.txt >> _line_u.txt

awk '{if($1==2.7)printf "%g %g 10p 90 LB %.1f\n",$1,0.05,$1}' out_sB.txt > _text_d.txt
awk '{if($1==3.3)printf "%g %g 10p 90 LB %.1f\n",$1,0.05,$1}' out_sC.txt >> _text_d.txt
awk '{if($1==4.5)printf "%g %g 10p 90 LB %.1f\n",$1,0.05,$1}' out_sD.txt >> _text_d.txt
awk '{if($1==6.8)printf "%g %g 10p 90 LB %.1f\n",$1,0.05,$1}' out_sE.txt >> _text_d.txt

awk '{if($1==2.7)printf "%g %g 10p 90 LB %.1f\n",$2,0.05,$2}' out_sB.txt > _text_u.txt
awk '{if($1==3.3)printf "%g %g 10p 90 LB %.1f\n",$2,0.05,$2}' out_sC.txt >> _text_u.txt
awk '{if($1==4.5)printf "%g %g 10p 90 LB %.1f\n",$2,0.05,$2}' out_sD.txt >> _text_u.txt
awk '{if($1==6.8)printf "%g %g 10p 90 LB %.1f\n",$2,0.05,$2}' out_sE.txt >> _text_u.txt

awk '{if($1==7.0)printf "%g %g 12p 0 LT B\n",7.1,$3}' out_sB.txt >> _text_d.txt
awk '{if($1==7.0)printf "%g %g 12p 0 LM C\n",7.1,$3}' out_sC.txt >> _text_d.txt
awk '{if($1==7.0)printf "%g %g 12p 0 LM D\n",7.1,$3}' out_sD.txt >> _text_d.txt
awk '{if($1==7.0)printf "%g %g 12p 0 LM E\n",7.1,$3}' out_sE.txt >> _text_d.txt

awk '{if($2<20.0)y=$3}END{printf "%g %g 12p 0 RT B\n",20.1,y}' out_sB.txt >> _text_u.txt
awk '{if($2<20.0)y=$3}END{printf "%g %g 12p 0 RM C\n",20.1,y}' out_sC.txt >> _text_u.txt
awk '{if($2<20.0)y=$3}END{printf "%g %g 12p 0 RM D\n",20.1,y}' out_sD.txt >> _text_u.txt
awk '{if($2<20.0)y=$3}END{printf "%g %g 12p 0 RM E\n",20.1,y}' out_sE.txt >> _text_u.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

fig=fig_s.eps
xlabel='Roof settlement (mm)'
#ylabel='Limit Support pressure (MPa)'
range=0/20/0/2.0
scale=-8/10
xga=a5f1
yga=a0.2f0.1
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -BS -P -K > $fig
inpB=out_sB.txt
inpC=out_sC.txt
inpD=out_sD.txt
inpE=out_sE.txt
awk '{print $2,$3}' $inpB | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{print $2,$3}' $inpC | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{print $2,$3}' $inpD | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{print $2,$3}' $inpE | gmt psxy -R -J -W1.5 -K -O >> $fig
gmt psxy _line_u.txt -R -JX  -Sv0.2+s -W0.5 -K -O >> $fig
gmt pstext _text_u.txt -R -J -F+f+a+j -N -K -O >> $fig

xlabel='Time (day)'
ylabel='Support pressure (MPa)'
range=0/7/0/2.0
scale=8/10
xga=a1
yga=a0.2f0.1
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga -BWS -P -X8 -K -O >> $fig
inpB=out_sB.txt
inpC=out_sC.txt
inpD=out_sD.txt
inpE=out_sE.txt
awk '{print $1,$3}' $inpB | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{print $1,$3}' $inpC | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{print $1,$3}' $inpD | gmt psxy -R -J -W1.5 -K -O >> $fig
awk '{print $1,$3}' $inpE | gmt psxy -R -J -W1.5 -K -O >> $fig
gmt psxy _line_d.txt -R -JX  -Sv0.2+s -W0.5 -K -O >> $fig
gmt pstext _text_d.txt -R -J -F+f+a+j -N -K -O >> $fig
echo '0 2.1 12p 0 CM Limit Support pressure (MPa)' | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
echo '0 0' | gmt psxy -R -J -Sp -O >> $fig