GMT绘制大地水准面图

葡萄

img2mercgrd geoid.img.9.2 -Gg.grd -R-180/180/-72/72 -T1 -N3 -V
grdgradient g.grd -Gintg.grd -A0/270 -Ne0.6
grd2cpt g.grd -Z -Crainbow > grav.cpt
grdimage g.grd -Iintg.grd -Cgrav.cpt -Jx0.07 -X2 -Y4 -K -V>geoid.ps
psbasemap -R-180/180/-72/72 -Jm0.07 -Ba30 -O -K -U"xxq">> geoid.ps

psscale -D12.5/-1/20/0.6h -Cgrav.cpt -B2000 -I -O -V -K >>geoid.ps
echo 0 77 25 0 7 CM GLOBAL GEOID | pstext -G125 -R -Jm -O -N -K>>geoid.ps
echo 0.2 77.2 25 0 7 CM GLOBAL GEOID | pstext -R -Jm -O -N>>geoid.ps
del .gmt* *.grd *.cpt

                     数据在此！