gmtset FRAME_PEN 3p BASEMAP_TYPE PLAIN DEGREE_FORMAT 2

#img2grd data/bathy/SmithSandwell/topo_8.2.img -Gdata/bathy/SmithSandwell/topo5km44-52.grd -R-132/-122/44/52 -T1 -V

psbasemap -Jm1 -R-132/-122/44/52 -K -X1 -Y1 -P -Ba3f1/a1f1wESn -V > topo5km44-52.ps

#grdinfo data/bathy/SmithSandwell/topo5km44-52.grd
makecpt -T-4500/0/100 -Z > grey_nepac.cpt

grdimage data/bathy/SmithSandwell/topo5km44-52.grd -Jm -R -Cgrey_nepac.cpt -B -V -O -K >> topo5km44-52.ps

grdcontour data/bathy/SmithSandwell/topo5km44-52.grd -Jm -R -C500 -G3m -A1000 -B -V -O -K >> topo5km44-52.ps

# plot the coast line
pscoast -Jm -R -Dh -A9 -O -W -N1 -N2 -I1 -V >> topo5km44-52.ps
#pscoast -Jm -R -Dh -A9 -O -W -G0/200/0 -N1 -N2 -I1 -V >> topo5km44-52.ps

grd2xyz data/bathy/SmithSandwell/topo5km44-52.grd -V > data/bathy/SmithSandwell/topo5km44-52.xyz