GrADS画特定经线方法
set clevs 23.5; d lat
set clevs 120;d lon
GrADS 查看多个打开文件的ctl
q ctlinfo 1
q ctlinfo 2
...
2个不同的缺测值 fwrite更改默认缺测值
GrADS用fwrite输出时,默认的缺测值是-9.9900000E+08,如果输入源的缺测值与之不一致,则会出现两种缺测值的情况。解决的办法,可以在fwrite 声明后面设置默认缺测值设置。比如
‘reinit‘
‘open J:/ec79-15/ansnablameansst79_15.ctl‘
‘set t 1 12‘
‘save=ave(absnablamsst,t+0,t=444,12)‘
‘modify save seasonal‘
‘set gxout fwrite‘
‘set undef 9.999E+20‘
‘set x 1 480‘
‘set y 1 89‘
‘set fwrite J:/ec79-15/seasonal.grd‘
‘d save‘
‘disable fwrite‘
‘reinit‘
;
(注:modify save seasonal命令意为将save变量修改成季节变量,方便与其它不同大小的变量计算(比如求季节距平),此外modify还有另外一个选项diurnal日变化)
NCL使用自定义色标
清风给的调色盘确实好用,生成的rgb格式也可以。但是最近再次遇到了ncl无法正常读取colormap错误的问题,参考了一下官网。最后发现是这样的。官网给出的文件有“ncolors=8”这一行,用NCL调色盘生成的脚本也有这一句话。但是在我的ncl里面无法使用。最后删了这一句话之后就好了。另外,注意前两行是用来定义底色和填色的。需要手动添加。比如加的是“255 255 255”和“0 0 0”这两行。最后能正常生成了。(一句话,用GrADS调色盘生成的rgb文件需要删掉nclolors这一行,并在文件最前面加上“255 255 255”和“0 0 0”这两行(定义底色))
(注:今天发现rgb文件最后一行上也要加,貌似是这样的:第一行定义,第二行定义填充线色,最后一行表示背景色)
GrADS使用NCL色标
(20171102注:GrADS使用colormap命令使用NCL色标时候,出现了一个问题,发现colormap打开的相关rgb文件在ncolors那一行后面还要有 # r g b 这一行,否则可能无法正常使用。
Cressman插值(oacre函数)
oacres函数实际上就是对数据进行多次cressman插值,其插值的影响半径是逐渐减少的,默认的影响半径(单位是格点)是10,7,4,2,1。如果某格点影响半径范围内的站点数目小于3个,那么该格点的值设置为缺测值。(参考http://cola.gmu.edu/grads/gadoc/gradfuncoacres.html)
oacres函数源代码 ftp://cola.gmu.edu/grads/2.1/grads-2.1.1.b0-src.tar.gz,解压后src目录下gafunc.c文件,第5313-5602行
GrADS 坐标
调字体
set xlopts 4 1 0.1 设置X轴标记数字的颜色、粗细和大小(英寸)
set ylopts 4 1 0.1 设置Y轴标记数字的颜色、粗细和大小(英寸)
调坐标间隔
set xlint interval
set ylint interval
set xlevs/ylevs 再控制标记的分布