# ********************************************
# Drawing by GMT
# ********************************************
xlabel="Displacement  ( u/L or v/L )"
ylabel="Load  ( PL@+2@+/EI )"
# ===============================
# Making of batch file for legend
# ===============================
gmt set FONT_ANNOT_PRIMARY 10
gmt set MAP_ANNOT_OFFSET_PRIMARY 0.3c
gmt set FONT_LABEL 10p
gmt set MAP_LABEL_OFFSET 0.3c
gmt set MAP_TICK_LENGTH_PRIMARY 0.0c
# ===============================
# Plotting (no hinges)
# ===============================
range=0/1/-20/40
scale=10/12
xga=g0.2a0.2
yga=g10a10
inp=_dat_lee.csv
fig=fig_lee.eps
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BWSen -P -X5 -Y6 -K > $fig
# *** Elastica (Line) ***
awk 'BEGIN{FS=","}{if(1<=NR)print -$1*1e6/2e5/833, $2/1000}' $inp | gmt psxy -: -R -J -B -W1 -P -O -K >> $fig
awk 'BEGIN{FS=","}{if(1<=NR)print -$1*1e6/2e5/833,-$3/1000}' $inp | gmt psxy -: -R -J -B -W1 -P -O -K >> $fig
# *** Calculation (Symbol) ***
awk 'BEGIN{FS=","}{if(1<=NR&&-$1*1e6/2e5/833<40)print -$1*1e6/2e5/833, $2/1000}' $inp | gmt psxy -: -R -J -SC0.15 -G255 -W1 -N -P -O -K >>$fig
awk 'BEGIN{FS=","}{if(1<=NR&&-$1*1e6/2e5/833<40)print -$1*1e6/2e5/833,-$3/1000}' $inp | gmt psxy -: -R -J -ST0.15 -G255 -W1 -N -P -O -K >>$fig
#
awk 'BEGIN{FS=","}{if(1<=NR&&-$1*1e6/2e5/833<40){yy=-$1*1e6/2e5/833;xx= $2/1000}}END{printf "%g %g 12 0 TR u/L",xx-0.02,yy}' $inp | gmt pstext -R -J -F+f+a+j -N -P -O -K >> $fig
awk 'BEGIN{FS=","}{if(1<=NR&&-$1*1e6/2e5/833<40){yy=-$1*1e6/2e5/833;xx=-$3/1000}}END{printf "%g %g 12 0 TL v/L",xx+0.02,yy}' $inp | gmt pstext -R -J -F+f+a+j -N -P -O -K >> $fig
# *** marginal note ***
echo "0.02 41 12 0 BL Lee's frame" | gmt pstext -R -J -F+f+a+j -N -P -O >> $fig