GrADS，NCL一些经验

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 再控制标记的分布

GRADS 风羽朝向 set hempref auto|nhem|shem

默认的是坐标在北半球，则风羽朝向风的左侧，在南半球，则风羽朝向风的右侧。一般情况下是没问题的。但是，如果画的是有特殊含义的图（坐标经纬度是相对的经纬度，那么就需要修改这一选项）