inp=out_fem_torus.csv
awk 'BEGIN{FS=","}{if(1090<=NR&&NR<=1449)print $3/1000,$2/1000}' $inp > _inp1.txt
awk 'BEGIN{FS=","}{if(NR==1090)print $3/1000,$2/1000}' $inp >> _inp1.txt
awk 'BEGIN{FS=","}{if(1090<=NR&&NR<=1449)print ($3+200*$5)/1000,($2+200*$4)/1000}' $inp > _inp2.txt
awk 'BEGIN{FS=","}{if(NR==1090)print ($3+200*$5)/1000,($2+200*$4)/1000}' $inp >> _inp2.txt
gmt set FONT_LABEL 10p
gmt set MAP_TICK_LENGTH_PRIMARY -0.1c
range=0/7/1/8
scale=-7/7
xga=a1f1
yga=a1f1
xlabel="Distance (m)"
ylabel="Elevation (EL.m)"

fig=fig_graph1.eps
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BES -P -K > $fig
gmt psxy _inp1.txt -R -J -W2,#777777,10_5:0 -O -K >> $fig
gmt psxy _inp2.txt -R -J -W2,#000000 -O -K >> $fig
echo "0.05 0.98 10p,Helvetica 0 TR @~d@~ x 200" | gmt pstext -R0/1/0/1 -J -F+f+a+j -N -K -O >> $fig
echo "0 0" | gmt psxy -R0/1/0/1 -J -Sp -O >> $fig


awk 'BEGIN{FS=","}{if(1812<=NR&&NR<=2171)print $1,$7}' $inp > _inp1.txt
awk 'BEGIN{FS=","}{if(1812<=NR&&NR<=2171)print $1,$5}' $inp > _inp2.txt
range=0/360/0/350
scale=10/6
xga=a90f10g90
yga=a50f10g50
xlabel="@~f@~ (degree)"
ylabel="Stress (N/mm@+2@+)"

fig=fig_graph2.eps
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BSWen -P -K > $fig
gmt psxy _inp1.txt -R -J -W2,#000000 -K -O >> $fig
gmt psxy _inp2.txt -R -J -W2,#777777,10_5:0 -O -K >> $fig
echo "30 190 12p,Helvetica 0 MC @~s@~@-@~j@~@-" | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
echo "30  80 12p,Helvetica 0 MC @~s@~@-@~q@~@-" | gmt pstext -R -J -F+f+a+j -N -K -O >> $fig
echo "0 0" | gmt psxy -R0/1/0/1 -J -Sp -O >> $fig