#!/bin/bash
#	$Id: EGM96.sh 11490 2013-05-16 06:26:21Z pwessel $
#
# Evaluate first 36 degrees of EGM96 geoid model

ps=EGM96.ps
gmt sph2grd EGM96_to_36.txt -I1 -Rg -Ng -Ggrid.nc
gmt grdmath grid.nc 1e5 MUL = grid.nc
gmt grd2cpt grid.nc -Crainbow -E16 -Z > t.cpt
gmt grdimage grid.nc -JH0/6i -P -B30 -Ct.cpt > $ps