#
#  Global version of Plate 1 in GRL paper, Satellite magnetic
#  anomalies related to seafloor spreading in the South Atlantic
#   Ocean
#          fig1x2a_d1540a_global.script
#
psbasemap -R0/7/0/9.2 -Jx1 -B0 -P -K  > fig1x2_d1540a_global.ps
#
xyz2grd m102389_d1540_tf400.data -Gtemp.grd -I2/2 -R-180/180/-88/88
grdsample temp.grd -Gm102389_d1540_tf400_global.grd -I0.5/0.5 -R-180/180/-70/70
grdimage m102389_d1540_tf400_global.grd -Clangel.cpt -Jm0.015d -P -O \
 -K -R-180/180/-70/70 -X0.7 -Y0.33 >> fig1x2_d1540a_global.ps 
pscoast -R-180/180/-70/70 -Ba60f60/a35f35WSen -Dl -O -K -A1 -Jm0.015d \
 -P -G0 >> fig1x2_d1540a_global.ps
psxy isochrons_gmt.latlon -R-180/180/-70/70 -: -H -Jm0.015d -M -O -P \
 -K >>fig1x2_d1540a_global.ps
psxy platebound_gmt.latlon -R-180/180/-70/70 -: -H  -Jm0.015d -M -O -K \
 -P -W4 >>fig1x2_d1540a_global.ps
echo "200. 0. 14 90.0 1 MC Observations" | pstext -O -K -P -R-180/180/-70/70 \
 -Jm0.015d >>fig1x2_d1540a_global.ps
#
#xyz2grd tf_d1439.model2_all.xyz -Gtf_d1439.model2_all.grd -I0.5/0.5 \
 -R-180/180/-90/90
grdmath tf_d1439.model2_all.grd 5.0 DIV = tf_d1439.model2_all.div10.grd
grdimage tf_d1439.model2_all.div10.grd -Clangel.cpt -Jm0.015 -O -K \
 -R-180/180/-70/70 -Y5.1 >> fig1x2_d1540a_global.ps 
pscoast -R-180/180/-70/70 -Jm0.015 -Ba30f30/a35f35WSen -Dl -O -K -A1 \
 -G0 >> fig1x2_d1540a_global.ps
psxy isochrons_gmt.latlon -R-180/180/-70/70 -: -H -M -O -K \
 -Jm0.015  >>fig1x2_d1540a_global.ps
psxy platebound_gmt.latlon -R-180/180/-70/70 -: -H -M -O -K -Jm0.015 \
 -W4 >>fig1x2_d1540a_global.ps
echo "200. 0. 14 90.0 1 MC Model" | pstext -O -K -P -R-180/180/-70/70 \
 -Jm0.015d >>fig1x2_d1540a_global.ps
# 
psscale -Clangel.cpt -D2.4/-0.8/4/0.25h -O \
 -B:nT::.FIELD: >>fig1x2_d1540a_global.ps