python3 py_fig_data.py out_torus.txt # Displacement mode 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="Radial distance r (m)" ylabel="Elevation z (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 _mesh.txt -R -J -W1.0,#aaaaaa,10_5:0 -O -K >> $fig gmt pstext -R -J -F+f+a+j -N -K -O << EOT >> $fig #x y size ang MC text 6 5 10p 0 ML 0 4 7 10p 90 MR 90 2 5 10p 0 MR 180 4 3 10p 90 ML 270 EOT gmt psxy _disp.txt -R -J -W1.0,#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 # Stress distribution 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 _p1.txt -R -J -W2,#000000 -K -O >> $fig gmt psxy _st.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