GMT绘制重力异常等值线图

vfory

gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
gmtset FRAME_PEN 1.0p
gmtset FRAME_WIDTH 0.1

默认的边框宽度是0.2，感觉有点宽改作0.1
gmtset TICK_LENGTH 0.1
gmtset LABEL_FONT 5
gmtset LABEL_FONT_SIZE 10
gmtset ANOT_FONT_SIZE 8
set psf=cdm.ps
set inf=gra.cdm
surface -I2m -V %inf% -Gtmp.grd -R90/100/35/40
grd2cpt tmp.grd -Chaxby -Z -V > tmp.cpt
grdimage tmp.grd -R -JM15 -Ctmp.cpt -P -E135 -K -V -Y3.0 > %psf%
grdcontour tmp.grd -JM -C20 -A100+p+g255+s5 -S4 -B1a2 -W0.25p -P -V -O -K  >> %psf%
psscale -D8/-0.7/10/0.2h -Ctmp.cpt -B100/:"mgal": -O -P -I -V -K>> %psf%
del .gmt* tmp.*

网格化surface改为nearneighbor，再使用grdgradient产生个梯度文件，代码如下：
gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
gmtset FRAME_PEN 1.0p
gmtset FRAME_WIDTH 0.1
gmtset TICK_LENGTH 0.1
gmtset LABEL_FONT 5
gmtset LABEL_FONT_SIZE 10
gmtset ANOT_FONT_SIZE 8
set psf=cdm.ps
set inf=gra.cdm
nearneighbor -I2m -S5m -V %inf% -Gtmp.grd -R90/100/35/40
grd2cpt tmp.grd -Chaxby -Z -V > tmp.cpt
grdgradient tmp.grd -A45 -Nt -Gtmp.i
grdimage tmp.grd -R -Itmp.i -JM15 -Ctmp.cpt -P -E135 -K -V -Y3.0 > %psf%
grdcontour tmp.grd -JM -C20 -A100+p+g255+s5 -S4 -B1a2 -W0.25p,gray30 -P -V -O -K  >> %psf%
psscale -D8/-0.7/10/0.2h -Ctmp.cpt -B100/:"mgal": -O -P -I -V -K>> %psf%
del .gmt* tmp.*

所用数据请到http://topex.ucsd.edu/cgi-bin/get_data.cgi按经纬度范围下载。