gawk 'BEGIN{FS=","}{if(1090<=NR&&NR<=1449)print $3/1000,$2/1000}' out_toroidal.csv > _inp1.txt
gawk 'BEGIN{FS=","}{if(NR==1090)print $3/1000,$2/1000}' out_toroidal.csv >> _inp1.txt
gawk 'BEGIN{FS=","}{if(1090<=NR&&NR<=1449)print ($3+200*$5)/1000,($2+200*$4)/1000}' out_toroidal.csv > _inp2.txt
gawk 'BEGIN{FS=","}{if(NR==1090)print ($3+200*$5)/1000,($2+200*$4)/1000}' out_toroidal.csv >> _inp2.txt
gmtset ANOT_FONT_SIZE 10
gmtset LABEL_FONT_SIZE 10
gmtset TICK_LENGTH -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
psbasemap -R$range -JX$scale -B$xga:"$xlabel":/$yga:"$ylabel":ES -P -K > $fig
psxy _inp1.txt -R -J -W5t10_5:0/#777777 -O -K >> $fig
psxy _inp2.txt -R -J -W5/#000000 -O -K >> $fig
echo "0.05 0.98 10 0 0 TR @~d@~ x 200" | pstext -R0/1/0/1 -J -N -K -O >> $fig
echo "0 0" | psxy -R0/1/0/1 -J -Sp -O >> $fig


gawk 'BEGIN{FS=","}{if(1812<=NR&&NR<=2171)print $1,$7}' out_toroidal.csv > _inp1.txt
gawk 'BEGIN{FS=","}{if(1812<=NR&&NR<=2171)print $1,$5}' out_toroidal.csv > _inp2.txt
gmtset ANOT_FONT_SIZE 10
gmtset LABEL_FONT_SIZE 10
gmtset TICK_LENGTH -0.1c
range=0/360/0/350
scale=10/6
xga=a90f10
yga=a50f10
xlabel="@~f@~ (degree)"
ylabel="Stress (N/mm@+2@+)"

fig=fig_graph2.eps
psbasemap -R$range -JX$scale -B$xga:"$xlabel":/$yga:"$ylabel":SW -P -K > $fig
psxy _inp1.txt -R -J -W5/#000000 -O -K >> $fig
psxy _inp2.txt -R -J -W5t10_5:0/#777777 -O -K >> $fig
echo "30 190 12 0 0 MC @~s@~@-@~j@~@-" | pstext -R -J -N -K -O >> $fig
echo "30  80 12 0 0 MC @~s@~@-@~q@~@-" | pstext -R -J -N -K -O >> $fig
psbasemap -R -JX -B/EN -K -O >> $fig
echo "0 0" | psxy -R0/1/0/1 -J -Sp -O >> $fig