#
#    plate2a.gmt on snooks in induced+remanence_model directory 
#     Alternate version of Plate 2
#
#   For manuscript 'Global lithospheric magnetization models 
#     validated and refined using new satellite observations"
#
/bin/awk '{if ($1>0 && $1<80) print $2, $1, $11; }' icosfid2.f2c >t2
/bin/awk '{if ($1>-160 && $1<-30) print $1, $2, $3; }' t2 >t3
surface t3 -Gtf1mean.grd -I2.5/2.5 -R-160/-30/0/80 
grdsample tf1mean.grd -Gtf1mean_oneq.grd -I0.5 -R-160/-30/0/80
grdcut tf1mean_oneq.grd -Gtemp.grd -R-130/-60/20/60
pscoast -R-130/-60/20/60 -JQ-95/3i -Ba20WSen -Dc -X.8 -Y7.5 \
-K -W5 -P -A2000 > plate2a.ps
grdcontour temp.grd -JQ-95/3i -C5 -P -O -K -A5 \
 -R-130/-60/20/60 >> plate2a.ps
#
xyz2grd total_d1425.xyz -Gtotal_d1428.grd -R0/360/-90/90 -I0.5/0.5 -V
grdcut total_d1428.grd -Gtotal_d1428_na.grd -R-130/-60/20/60
grdimage total_d1428_na.grd -JQ-95/3i -Cmag2.cpt -R-130/-60/20/60 \
 -P -X3.3 -O -K -E300 >> plate2a.ps
pscoast -R-130/-60/20/60 -JQ-95/3i -Ba20WSen -Dl -O -K -P \
 -N1 -N2 -W2 -A2000 >> plate2a.ps
#
blockmean magsat_d1526_tf400.xyz -: -I2.5/2.5 -L -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
grdcut tf1mean_oneq.grd -Gtf1mean_oneh.grd -R-130/-60/20/60
grdinfo tf1mean_oneh.grd
grdimage tf1mean_oneh.grd -JQ-95/3i -Cmag2.cpt -R-130/-60/20/60 \
 -X-3.3 -Y-3.3 -K -O -P >> plate2a.ps
pscoast -R-130/-60/20/60 -JQ-95/3i -Ba20WSen -Dl -O -P \
 -N1 -N2 -K -W2 -A2000 >> plate2a.ps
#
blockmean oersted_03c_01_d1526_tf400.xyz -: -I2.5/2.5 -L -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
grdcut tf1mean_oneq.grd -Gtf1mean_oneh.grd -R-130/-60/20/60
grdinfo tf1mean_oneh.grd
grdimage tf1mean_oneh.grd -JQ-95/3i -Cmag2.cpt -R-130/-60/20/60 \
 -X3.3 -O -K -P >> plate2a.ps
pscoast -R-130/-60/20/60 -JQ-95/3i -Ba20wsen -Dl -O -K \
 -N1 -N2 -W2 -A2000 >> plate2a.ps
#
/bin/awk '{print $2, $1, $12; }' icosfid2.f2c_namod_wrem >t5
blockmean t5 -I2.5/2.5 -L -R0/360/-90/90 >tf1mean.xyz
surface tf1mean.xyz -Gtf1mean.grd -I2.5 -R0/360/-90/90
grdsample tf1mean.grd -Gtf1mean_oneq.grd -I0.5 -R0/360/-90/90
blockmean orsted_semm_delintmag.xyz -: -I5/5 -L -R0/360/-90/90 >tf2mean.xyz
surface tf2mean.xyz -Gtf2mean.grd -I5 -R0/360/-90/90
grdsample tf2mean.grd -Gtf2mean_oneq.grd -I0.5 -R0/360/-90/90
grdmath tf1mean_oneq.grd tf2mean_oneq.grd ADD = total_mag.grd
grdcut total_mag.grd -Gtemp.grd -R-130/-60/20/60
pscoast -R-130/-60/20/60 -JQ-95/3i -Ba20WSen -Dc -X-3.3 -Y-3.3 \
-K -W5 -P -A2000 -O >> plate2a.ps
grdcontour temp.grd -JQ-95/3i -C5 -P -O -K -A5 \
 -R-130/-60/20/60 >> plate2a.ps
#
xyz2grd total_iter1_d1425.xyz -Gtotal_d1428.grd -R0/360/-90/90 -I0.5 -V
grdcut total_d1428.grd -Gtotal_d1428_na.grd -R-130/-60/20/60
grdimage total_d1428_na.grd -JQ-95/3i -Cmag2.cpt -R-130/-60/20/60 \
 -P -X3.3 -O -K -E300 >> plate2a.ps
pscoast -R-130/-60/20/60 -JQ-95/3i -Ba20wsen -Dl -O -K \
 -N1 -N2 -W2 -P -A2000 >> plate2a.ps
#
psscale -Cmag2.cpt -D3.5/4.5/5/0.3 -Xr -Yr -O -K -P -L -B::/:"nT": >> plate2a.ps
#
rm tf1mean.xyz
rm tf1mean.grd
rm tf1mean_oneq.grd
rm tf1mean_oneh.grd
rm tf2mean.xyz
rm tf2mean.grd
rm tf2mean_oneq.grd
rm total_d1428.grd
rm total_d1428_na.grd
rm total_mag.grd
rm t2
rm t3
rm t5
rm temp.grd
#
pstext -R-1/1/-1/1 -JX8i/11i -O -P -S0.5p -N <<EOF >> plate2a.ps
-1.8 0.6 16 0 0 6 a)
-0.25 0.6 16 0 0 6 b)
-1.8 0.0 16 0 0 6 c)
-0.25 0.0 16 0 0 6 d)
-1.8 -0.6 16 0 0 6 e)
-0.25 -0.6 16 0 0 6 f)
EOF