日志
GrADS中差点没注意的细节-转换数据时缺测值的变化
热度 6 ||
今天用GrADS求了一下季节平均，源数据中的缺测值是-9999.0，求完之后画图一看，不对了，不是这个数了，还好先用的命令行测试的，下面直接输出了当前的数据范围和间隔，最小值是-9.99e+08，然后又固定到一个格点输出数值看了下，果然是这个值，修改ctl中的缺测值为这个，出图就正常了！
我不确定是不是所有文件都是这样，以前转换缺测值为-9.99e+33的时候貌似没这个问题，改天再试一下~
发表评论 评论 (12 个评论)
- 回复 刘振
- 今天也碰到这个问题，而且诡异的发现是用我之前老师给的grads的版本没有出现这个问题，但用的气象家园的清风斑竹改善的grads出现了这个问题，不知道这个问题是否具有普遍性。因为之前本科做论文也有同学因为用的是气象家园的grads出现了这个bug，所以提出来问问版主
- 回复 mofangbao
- 可能是你老是给的版本比较旧 新版的都是这样的 并不是打包才造成这样的 打包只是加了些gs和常用文件 没有对现有文件做修改
刘振: 而且我用的是ncep的海温数据，缺测值是-9.99e+33额，建议版主用其它版本的grads重做你的上述问题，谢谢
- 回复 mofangbao
- 不是bug 这样是很方便的 比如很多种资料批处理的时候就不用一个个的考虑处理后的缺测值了 便于统一 老版本的功能上当然不够新版本
刘振: 哦哦，这样啊，谢谢啊。那老版本做打包是不是更好啊，新版本这种结果我觉得本身就算是bug的吧
- 回复 afterrain_29
- 我最近也发现这个问题，hadisst_sst.nc本来陆地上的缺省值是-1e+30，但我在grads（气象家园整合版）上把冬季12，1，2这3个月的数据取出来写成grd后再用fortran读出来，那些陆地上的缺省值就变成-9.99+e08，以前用老版本的grads就没有这个问题，缺省值是不会变的。