博文
GMT绘制grd格网图
||
REM ANOT_FONT_SIZE 表示在图上标记的数字的大小
REM GRID_CROSS_SIZE 表示在图上标记的边框的大小
REM pscoast -R100/140/0/40 -JM -N1 -Di -A10000 -W -V -K -O >>china.ps 画边界线
REM psscale -Ct.cpt -D2.05i/-0.40i/3.75i/0.15ih -I1.5 -V -O -B/:nT: >>china.ps 画色棒
REM GRID_CROSS_SIZE 表示在图上标记的边框的大小
REM pscoast -R100/140/0/40 -JM -N1 -Di -A10000 -W -V -K -O >>china.ps 画边界线
REM psscale -Ct.cpt -D2.05i/-0.40i/3.75i/0.15ih -I1.5 -V -O -B/:nT: >>china.ps 画色棒
gmtset GRID_CROSS_SIZE 0 ANOT_FONT_SIZE 10 FRAME_WIDTH 0.09c
gmtset GRID_PEN 0.25pta
grdreformat kms02g.grd temp0.grd=nf -V
makecpt -Cwysiwyg -T-350/400/100 -Z >t.cpt
grdimage temp0.grd -R100/140/0/40 -E200 -JM4.0i -B5f5 -X3.5i -Y2.5i -K -V -Ct.cpt >china.ps
pscoast -R -JM -N1 -Di -A10000 -W -V -K -O >>china.ps
psscale -Ct.cpt -D2.05i/-0.40i/3.75i/0.15ih -I1.5 -V -O -B/:nT: >>china.ps
gmtset GRID_PEN 0.25pta
grdreformat kms02g.grd temp0.grd=nf -V
makecpt -Cwysiwyg -T-350/400/100 -Z >t.cpt
grdimage temp0.grd -R100/140/0/40 -E200 -JM4.0i -B5f5 -X3.5i -Y2.5i -K -V -Ct.cpt >china.ps
pscoast -R -JM -N1 -Di -A10000 -W -V -K -O >>china.ps
psscale -Ct.cpt -D2.05i/-0.40i/3.75i/0.15ih -I1.5 -V -O -B/:nT: >>china.ps
pause
西太平洋海域地磁异常图
gmtset GRID_CROSS_SIZE 0 ANOT_FONT_SIZE 15 FRAME_WIDTH 0.2c
gmtset GRID_PEN 0.25pta
grdreformat kms02_grv_bin.grd temp0.grd=nf -V
makecpt -Cwysiwyg -T-300/300/1 -D >t.cpt
grdimage temp0.grd -R100/140/0/40 -E200 -JM6.0i -B5f5/5f5WeSn -X1.5i -Y2.5i -K -P -V -Ct.cpt >china.ps
pscoast -R -JM -N1 -Di -A1000 -Wthin -V -K -O >>china.ps
psxy region1.txt -R -JM -Wthickest,black -K -O >>china.ps
psxy region2.txt -R -JM -Wthickest,black -K -O >>china.ps
rem psscale -Ct.cpt -D3.0i/-0.45i/5.75i/0.18ih -I1.5 -V -O -B150/:nT: >>china.ps
psscale -Ct.cpt -D3.0i/-0.45i/5.5i/0.18ih -I1.5 -V -O -B50/:mgal: >>china.ps
pause
西太平洋海域重力异常图
region1.txt
134 11
134 6
139 6
139 11
134 11
region2.txt
130 28130 23
135 23
135 28
130 28
https://blog.sciencenet.cn/blog-51026-566643.html
上一篇：【GAMIT-BLOCK】sh_upd_stnfo 用法
下一篇：【GMTSAR】GMTSAR软件处理实例终端输出参考
扫一扫，分享此博文