GMT in Arctic

Northern terrestrial


pscoast -r0/360/-90/90 -jg150/90/3.5i -b30g30f -dc -a1000 -g0 -p > arctic.ps

Click the fig to get postscript.

Northern Hemisphere

First, prepare the shaded grid.


#! /bin/sh
grd=/somewhere/ETOPO2.pacific.grd
grdgradient $grd  -A315 -Gtmp315.grad -V -Ne0.8
grdhisteq tmp315.grad -Gtmp315.hst -N -V
grdmath tmp315.hst 3  / = tmp315.hst3

Then run grdimage.


#! /bin/sh
gmtset PAPER_MEDIA a1+
gmtset DEGREE_FORMAT 3
output=arctic_img.ps
cpt=cptjapangeogra
grd=/somewhere/ETOPO2.pacific.grd
range=0/360/-90/90 
jpara=150/90/15i
grdimage $grd -JG$jpara -R$range -B10g10f -C$cpt -Itmp315.hst3 -P > $output

I used cptjapangeogra.
Jpeg file (1.1Mb)

Southern terrestrial


pscoast -R0/360/-90/90 -JG150/-90/3.5i -B30g30f -Dc -A1000 -G0 -P > antiarctic.ps

Click the fig to get postscript.