.gz ... 变化 如果 无法使用 更改 2 个 open sea
GrADS 画特定经线方法
GrADS 查看多个打开文件的 ctl
set clevs 23.5; d lat
set clevs 120;d lon
2 个不同的缺测值 fwrite 更改默认缺测值
q ctlinfo 1
q ctlinfo 2
...
GrADS 用 fwrite 输出时,默认的缺测值是 - 9.9900000E+08,如果输入源的缺测值与之不一致,则会出现两种缺测值的情况.解决的办法,可以在 fwrite 声明后面设置默认缺测值设置.比如
(注:modify save seasonal 命令意为将 save 变量修改成季节变量,方便与其它不同大小的变量计算 (比如求季节距平),此外 modify 还有另外一个选项 diurnal 日变化)
'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'
;
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 xlevs/ylevs 再控制标记的分布
set xlint interval
set ylint interval
GrADS,NCL 一些经验
来源: http://www.bubuko.com/infodetail-2461505.html