# ********************************************
# Drawing by GMT
# ********************************************
xlabel="Radius  R ( mm )"
ylabel="Load  P ( MPa )"
# ===============================
# Making of batch file for legend
# ===============================
gmt set FONT_ANNOT_PRIMARY 12
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.0c
# ===============================
# Plotting (no hinges)
# ===============================
range=1300/1700/0/10
scale=10/12
xga=g100a100
yga=g1a1

inp_1=_dat_circ060in.csv
fig=fig_circ060in.eps
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWSen -P -X5 -Y6 -K > $fig
# *** (Line) ***
awk 'BEGIN{FS=","}{if(1<=NR)print $1, $2+1650}' $inp_1 | gmt psxy -: -R -J -B -W1 -P -O -K >> $fig
awk 'BEGIN{FS=","}{if(1<=NR)print $1, $3+1350}' $inp_1 | gmt psxy -: -R -J -B -W1 -P -O -K >> $fig
# *** (Symbol) ***
awk 'BEGIN{FS=","}{if(1<=NR&&$1<10)print $1,$2+1650}' $inp_1 | gmt psxy -: -R -J -SC0.15 -G255 -W1 -N -P -O -K >>$fig
awk 'BEGIN{FS=","}{if(1<=NR&&$1<10)print $1,$3+1350}' $inp_1 | gmt psxy -: -R -J -ST0.15 -G255 -W1 -N -P -O -K >>$fig
#
awk 'BEGIN{FS=","}{if(1<=NR&&$1<10){yy=$1;xx=$2+1650}}END{printf "%g %g 12 0 TR u+R@-0R@-",xx+60,yy}' $inp_1 | gmt pstext -R -J -F+f+a+j -N -P -O -K >> $fig
awk 'BEGIN{FS=","}{if(1<=NR&&$1<10){yy=$1;xx=$3+1350}}END{printf "%g %g 12 0 TL v+R@-0T@-",xx-60,yy}' $inp_1 | gmt pstext -R -J -F+f+a+j -N -P -O -K >> $fig
echo 1650 0.2 12 90 ML R@-0R@-=1650 | gmt pstext -R -J -F+f+a+j -N -P -O -K >> $fig
echo 1350 0.2 12 90 ML R@-0T@-=1350 | gmt pstext -R -J -F+f+a+j -N -P -O -K >> $fig
# *** marginal note ***
echo 1300 10.1 12 0 BL Ring with large initial deflection | gmt pstext -R -J -F+f+a+j -N -P -O >> $fig