GMT绘制等值线图
set range=-R110/150/10/40set ofn=LUPIT.psrem Plot SST
surfacesst_1020.dat %range% -T0.01 -I0.1/0.1 -C0.001 -Gsst.grd
makecpt -Crainbow -T16/30/0.5 -Z > sst_colors.cptgrdimagesst.grd %range% -B5WSne -JM12 -Csst_colors.cpt -X5 -Y15 -K -P> %ofn%
pscoast -R -JM -Dh -Glightgray -Wthinnest -K -O >> %ofn%
gawk "{print $8,$6}" LUPIT_track.dat | psxy -W6/255lta -R -JM -K -O >> %ofn%
gawk "/ACT/{print $8,$6}" LUPIT_track.dat | psxy -W6/125 -R -JM -K -O >> %ofn%
gawk "NR%5==0{print $8-0.3,$6-0.1,"12 40 7 RM",$2"/"$3," ",substr($4,1,2),":00"}" LUPIT_track.dat | pstext -R -JM -K -O >> %ofn%
gawk "NR%5==0{print $8,$6}" LUPIT_track.dat | psxy -Sc0.2 -G0/0/255 -R -JM -K -O >> %ofn%
echo 134 15 12 -10 1 LM Actual Track of LUPIT | pstext -R -JM -K -O -G0/0/255>> %ofn%
echo 111 18 12 10 1 LM Predicted Track of LUPIT | pstext -R -JM -K -O -G255>> %ofn%
echo 130 42 25 0 7 CM Sea Surface Temperature in 10/20 | pstext -R -JM -N -K -O >> %ofn%
psscale -Csst_colors.cpt -D13/5/10/0.25 -Ba2g1/:(@+o@+C): -K -O >> %ofn%rem Plot SLA
surfacesla_0921.dat %range% -T0.01 -I0.1/0.1 -C0.001 -Gsla.grd
makecpt -Cpolar -T-40/40/1 -Z > sla_colors.cptgrdimagesla.grd -R -B5WSne -JM12 -Csla_colors.cpt -Y-13 -K -O >> %ofn%
grdcontour sla.grd -R -JM -A10f6 -C5 -G8 -Wa3/255 -Wc1/255 -O -K >> %ofn%
pscoast -R -JM -Dh -Glightgray -Wthinnest -K -O >> %ofn%
echo 130.85 19.85 19 0 1 LM --Cold Eddy | pstext -R -JM -K -O -G0>> %ofn%
echo 131 20 19 0 1 LM --Cold Eddy | pstext -R -JM -K -O -G0/0/255>> %ofn%
echo 122.85 21.85 19 0 1 LM --Warm Eddy | pstext -R -JM -K -O -G0>> %ofn%
echo 123 22 19 0 1 LM --Warm Eddy" | pstext -R -JM -K -O -G255/0/0>> %ofn%echo 130 42 25 0 7 CM Sea Level Anomaly in 09/21 | pstext -R -JM -N -K -O >> %ofn%
psscale -Csla_colors.cpt -D13/5/10/0.25 -Ba10g5/:(cm): -O >> %ofn%del .gmt* *.grd *.cpt------------------------------------------------------------------------------------------------------------讲解：surfacesst_1020.dat -R$lon1/$lon2/$lat1/$lat2 -T0.01 -I0.1/0.1 -C0.001 -Gsst.grd
--将SST资料(sst_1020.dat)网格化并存为grd文件(sst.grd)
-T: 设定tension factor
-I: 设定网格大小
-C: 设定收敛限制值
makecpt -Crainbow -T16/30/0.5 -Z > sst_colors.cpt
--产生等值图所需的颜色渐变文件(*.cpt)
-C: 选择参考色-T: 定义新cpt文件的范围,例如-Tz0/z1/dz为,起始值(z0),终止值(z1),间距(dz)
-Z: 让颜色具有连续性(smooth)
* 在此范例中使用了rainbow(预设)与polar两种色表
grdimagesst.grd -R$lon1/$lon2/$lat1/$lat2 -B5WSne -JM12 -Csst_colors.cpt -X5 -Y15 -K -P> $ofn
--使用sst.grd的资料绘出以颜色表示的等值图,色表使用sst_colors.cpt
grdcontour sla.grd -R -JM -A10f6 -C5 -G8 -Wa3/255 -Wc1/255 -O -K >> $ofn
--使用sst.grd的资料会出等值线
-A10f6 每隔10画出等值线并标出其值,f指定字型
-C5 每隔5画出等值线但不标出其值
-G8 設定沿著等值线标出值的间隔
-Wa3/255 设定-A所标出的等值线宽度及颜色
-Wc1/255 设定-C所标出的等值线宽度及颜色
pscoast -R -JM -Dh -Glightgray -Wthinnest -K -O >> $ofn
--海岸线
psscale -Csst_colors.cpt -D13/5/10/0.25 -Ba2g1/:'(@+o@+C)': -K -O >> $ofn
--绘制sst_colors.cpt的sebiao-D: 标出位置,基准点为上承panel左下点,-D13/5/10/0.25意思为基准点右移13厘米,上移5厘米,长度10,宽度0.25
-Ba2g1: 每隔单位2标出数值,每隔单位1标线,標出温度单位为(℃)
* 两个@+之间的字为"上标"字
------------------------------------------------------------------------------------------------------------------------------------------------
页:
[1]