#
# Plate 1 in GRL paper, A magnetic map of Mars
#        in colored form with an accompanying scale
#
#   /home/purucker/manuscripts/mgs_mag_esmap/plate1_mgs_mag.gmt 
#       on snooks.gsfc.nasa.gov
#       Postscript output is 7 Mb in size
#
gmtset GRID_CROSS_SIZE 0.1i COLOR_MODEL RGB PAPER_MEDIA A4+
xyz2grd br_spoabgt3_it87_mars_rdipole_200km_pm180_clip.xyz -Gtemp.grd \
 -L -I1/1 -R0/360/-90/90
grdsample temp.grd -Gtemp_pixel.grd -I1/1 -R0/360/-90/90 -T -Lg 
grdsample temp_pixel.grd -Gbr_spoabgt3_it87_mars_rdipole_200km.grd \
 -I0.25/0.25 -R0/360/-90/90 -Lg
grdimage br_spoabgt3_it87_mars_rdipole_200km.grd -Cmag.cpt -JX9/4.5d \
 -Iquartgrad.grd -Bf10a60WeSn -R0/360/-90/90 -K > plate1_mgs_mag.eps
pstext -R0/360/-90/90 -JX9/4.5d -O -K -S0.75p << END >> plate1_mgs_mag.eps
285 -20 12 0 1 6 V
335 -12 12 0 1 6 G
195  35 12 0 1 6 A
335  38 12 0 1 6 C
END
psxy -R0/360/-90/90 -JX9/4.5d -G0 -O -K -SV <<EOF >> plate1_mgs_mag.eps
288 -20 70 0.41 
332 -12 280 0.50
EOF
psscale -Cmag.cpt -D4.5/6.5/5/0.25h -B:nT: -O -K -L >> plate1_mgs_mag.eps
cat << EOF | pstext  -R0/11/0/8.5 -Jx1i -O -N >> plate1_mgs_mag.eps
9 4.65 12 0.0 1 2 0\312
6.75 4.65 12 0.0 1 2 90\312
4.5 4.65 12 0.0 1 2 180\312 
2.25 4.65 12 0.0 1 2 270\312
0 4.65 12 0.0 1 2 360\312
4.5 5 16 0.0 1 2 West Longitude
4.5 -0.5 16 0.0 1 2 East Longitude
EOF