GMT色度条更改--附加代码

上一文章讲到GMT色度条的更改，接下来附加代码：

rem 2021.3.29/ Z.Q.WEN ECJTU

set R=60/96/30/50

set J=M20c

gmt gmtset FORMAT_GEO_MAP=ddd:mm:ssF

set FONT_ANNOT_PRIMARY  [20p,Helvetica, black]

set MAP_ORIGIN_X   = 0.5i

set MAP_ORIGIN_Y   = 0.5i

set PS="tarimu.ps"

REM gmt grdcut ETOPO1_Bed_g_gdal.grd -R%R% -GcutTopo.grd

REM gmt grdgradient cutTopo.grd -Ne0.7 -A50 -GcutTopo_i.grd

REM gmt grd2cpt cutTopo.grd  -Cglobe  -T1000/5000/10  -D -Z >t.cpt

gmt grdimage cutTopo.grd -IcutTopo_i.grd  -R%R% -J%J% -Xc -Yc  -Ctopo.cpt  -K >%PS%

gmt pscoast -R%R% -J%J% -Df -W0.1p -I1/0.8p,blue -O -K >>%PS%

gmt psxy tarimu.txt -J%J% -R%R% -Sc0.05c -W1,red -O -K>> %ps%

REM gmt xyz2grd  grace.txt -R%R% -I1  -Gtmp.grd

REM gmt grdsample tmp.grd -Gtmp.grd -I0.01

REM gmt makecpt -Crainbow -T-140/400/0.1 -D >tem.cpt

REM gmt psclip Amazon.txt  -J%J% -R%R% -Am  -O -K >>%PS%

REM gmt grdimage tmp.grd -R%R% -J%J% -Ctem.cpt -B5f1/5f1WenS -O -K>>%PS%

REM gmt Amazon.txt -J%J% -R%R% -Sc0.05c -W1,red  -O -K >> %ps%

REM gmt psclip -C -O -K  >>%PS%

gmt psbasemap -R%R% -J%J% -Bf10a10 -BWeSn -Bg10  --FONT_ANNOT_PRIMARY=20p,Helvetica,black -X0c -Y0c -O -K >>%PS%

gmt psscale -Dx0.5c/-01.2c+jTC+w8c/0.4c+h+e --FONT_ANNOT_PRIMARY=20p,Helvetica,black -X11c -Y3c  -O  -K -Ctopo.cpt -Ba4000 -By+lm>> %ps%

REM 绘制图中图

set Rg=-180/180/-90/90

set Jg=G70/40/6c

gmt grdimage ETOPO1_Bed_g_gdal.grd  -R%Rg% -J%Jg%  -Ct.cpt   -X-10c  -Y5c  -O -K >>%PS%

gmt pscoast -R%Rg% -J%Jg% -N1 -A1000 -B0  -Wthinnest  -O -K -X0c -Y0c  --MAP_FRAME_TYPE=plain>> %PS%

gmt psbasemap -R%Rg% -J%Jg% -D%R% -Fd+p2p,red+c0.1  -O>>%PS%

gmt psconvert  %PS% -A -Tg -P 

有任何问题，欢迎大家与我交流！！

QQ：226687758