#
#    plate1.gmt on geomag in induced+remanence_model directory 
#
#   For manuscript 'The southern edge of cratonic North America:
#     Evidence from new satellite magnetometer observations'
#
blockmean tfsemm0+mf.allmap1a.xyz -: -I2.5/2.5 -R-180/180/-90/90 -L > tf1mean.xyz
surface tf1mean.xyz -Gtf1mean.grd -I2.5/2.5 -R-180/180/-90/90 -V 
grdsample tf1mean.grd -Gtf1mean_oneq.grd -I0.5/0.5 -R-180/180/-90/90 -Lg -V
grdimage tf1mean_oneq.grd -JG-90/30/3.2i -Cmag1.cpt -R-180/180/-90/90 \
 -X0.5 -Y7.5 -K -P > plate1.ps
pscoast -R-180/180/-90/90 -JG-90/30/3.2i -Bg30/g15 -Dc -O -K \
 -A2000 -W5 -P >> plate1.ps
#
grdmath Atp12_4.00_dec_av.grd 2 MUL = total.grd
grdimage total.grd -JG270/30/3.2i -Cmag1.cpt -R0/360/-90/90 -Bg30/g15 \
 -X3.3 -O -K -E300 -P >> plate1.ps
pscoast -R0/360/-90/90 -JG270/30/3.2i -Bg30/g15 -Dc -O -K \
 -A2000 -W5 -P >> plate1.ps
#
blockmean tfsemm0+mf.allmap1a_360.xyz -: -I2.5/2.5 -R0/360/-90/90 -L > tf1mean.xyz
surface tf1mean.xyz -Gtf1mean.grd -I2.5/2.5 -R0/360/-90/90 -V 
grdsample tf1mean.grd -Gtf1mean_oneq.grd -I0.5/0.5 -R0/360/-90/90 -Lg -V
grdmath Atp12_4.00_dec_av.grd 2 MUL tf1mean_oneq.grd ADD = total.grd
grdimage total.grd -JG270/30/3.2i -Cmag1.cpt -R0/360/-90/90 -Bg30/g15 \
 -P -X-3.3 -Y-3.3 -O -K -E300 >> plate1.ps
pscoast -R0/360/-90/90 -JG270/30/3.2i -Bg30/g15 -Dc -O \
 -A2000 -K -P -W5 >> plate1.ps
#
xyz2grd total_d1425.xyz -Gtotal_d1428.grd -R0/360/-90/90 -I0.5/0.5 -V
grdimage total_d1428.grd -JG270/30/3.2i -Cmag1.cpt -R0/360/-90/90 \
 -Bg30/g15 -X3.3 -O -K -E300 -P >> plate1.ps
pscoast -R0/360/-90/90 -JG270/30/3.2i -Bg30/g15 -Dc -O \
 -A2000 -K -W5 -P >> plate1.ps
#
blockmean magsat_d1526_tf400.xyz -: -L -I2.5/2.5 -R-180/180/-90/90 >tf1mean.xyz
surface tf1mean.xyz -Gtf1mean.grd  -I2.5/2.5 \
 -R-180/180/-90/90 -V 
grdsample tf1mean.grd -Gtf1mean_oneq.grd -I0.5/0.5 -R-180/180/-90/90 -Lg -V
grdimage tf1mean_oneq.grd -JG-90/30/3.2i -Cmag1.cpt -R-180/180/-90/90 \
 -X-3.3 -Y-3.3 -K -O -P >> plate1.ps
pscoast -R-180/180/-90/90 -JG-90/30/3.2i -Bg30/g15 -Dc -O \
 -A2000 -P -K -W5 >> plate1.ps
#
blockmean oersted_03c_01_d1526_tf400.xyz -: -L -I2.5/2.5 -R-180/180/-90/90 >tf1mean.xyz
surface tf1mean.xyz -Gtf1mean.grd -I2.5/2.5 -R-180/180/-90/90 -V
grdsample tf1mean.grd -Gtf1mean_oneq.grd -I0.5/0.5 -R-180/180/-90/90 -Lg -V
grdimage tf1mean_oneq.grd -JG-90/30/3.2i -Cmag1.cpt -R-180/180/-90/90 \
 -X3.3 -O -K -P >> plate1.ps
pscoast -R-180/180/-90/90 -JG-90/30/3.2i -Bg30/g15 -Dc -O \
 -A2000 -K -W5 -O >> plate1.ps
#
psscale -Cmag1.cpt -D3.5/4.5/5/0.3 -Xr -Yr -O -K -P -L -B::/:"nT":  >> plate1.ps
#
pstext -R-1/1/-1/1 -JX8i/11i -O -P -S0.5p -N <<EOF >> plate1.ps
-1.8 0.75 16 0 0 6 a)
-0.25  0.75 16 0 0 6 b)
-1.8 0.15 16 0 0 6 c)
-0.25  0.15 16 0 0 6 d)
-1.8 -0.45 16 0 0 6 e)
-0.25  -0.45 16 0 0 6 f)
EOF 
#
rm tf1mean.xyz
rm tf1mean.grd
rm tf1mean_oneq.grd
rm total.grd
rm total_d1428.grd