wq_ocean_的博客

set plot=coast.psgmt pscoast -JM10c -R70/140/2/60 -Bafg -W1/thin -Ggray -P -K > %plot%rem -JM10c 表示使用墨卡托投影，整个地图的宽度是10厘米，地图的高度由-R和-J自动确定rem -R70/140/2/60表示X方向范围是70到140，Y方向范围是2到60°rem -Bafg会根据当前...

set ps=psvelo_ex2.psrem Long. Lat. Evel Nvel Esig Nsig CorEN SITErem (deg) (deg) (mm/yr) (mm/yr)echo -8. 0. 5.0 0.0 4.0 6.0 0.500 4x6 > tmp.txtecho ...

psvelo模块时用来绘制速度矢量图，但把其他类型的矢量当作速度就好，代码如下：set ps=psvelo_ex2.psrem Long. Lat. Evel Nvel rem (deg) (deg) (mm/yr) echo -8. 0. 5.0 0.0 > tmp.txtecho -5. 0. ...

set plot=ave.psgmt gmtinfo 2017050600.txtrem 查看文件信息，方便制作cpt，文件格式为 x y zgmt makecpt -Cjet -T20/60/5 > std.cptpscontour 2017050600.txt -R110/118/20/26 -JM6i -Cstd.cpt -Bafg -P -I >>%plot...