#!/bin/bash
#
#       $Id: gspline_1.sh 12114 2013-09-03 19:19:00Z fwobbe $

ps=gspline_1.ps

# Figure 2 in Wessel, P. (2009), A general-purpose Green's function-based
#	interpolator, Computers & Geosciences, 35, 1247–1254.


R=-2000/25000/0/11
T=Table_4.2.txt
gmt psbasemap -R-2/25/0/11 -JX6i/3i -P -K -Bx5f1+l"Distance (km)" -By2f1+l"Mg (ppm)" -BWSne -X1.25i -Y2i --FONT_LABEL=18p > $ps
gmt psxy -R$R -J -O -K $T -Sc0.075i -Gblack >> $ps
gmt greenspline -R-2000/25000 -I100 $T -Sl -D0 | gmt psxy -R$R -JX -O -K -Wthin,. >> $ps
gmt greenspline -R-2000/25000 -I100 $T -Sc -D0 | gmt psxy -R$R -JX -O -K -Wthin,- >> $ps
gmt greenspline -R-2000/25000 -I100 $T -St0.25 -D0 | gmt psxy -R$R -J -O -K -Wthin >> $ps
gmt pslegend -R$R -JX -O -K -F+p -Dx5.9i/2.9i/2.05i/TR --FONT_ANNOT_PRIMARY=12p << EOF >> $ps
S 0.2i - 0.35i - 0.5p 0.5i Tension (@%6%t@%% = 0.25)
S 0.2i - 0.35i - 0.5p,- 0.5i No tension
S 0.2i - 0.35i - 0.5p,. 0.5i Linear interpolation
EOF
gmt psxy -R$R -JX -O -T >> $ps