【已解决】grads出图的时候数值特别离谱

suifeng1106 · 发表于 2015-11-14 15:45:36

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 suifeng1106 于 2015-11-15 15:18 编辑

用的SST月平均资料，求1948-2007年的海温距平。已经成功得到了1948-2007年年平均SST，也得到了多年平均SST。但在求距平的时候遇到了问题。我可以确定的是自己的fortran程序没有错误，大致确定为缺测值的问题。现在的问题就是能出图，不过特别离谱。

这是出的图

river · 发表于 2015-11-14 16:10:33

既然确定问题了，那楼主就改吧

suifeng1106 · 发表于 2015-11-14 17:17:17

river 发表于 2015-11-14 16:10
既然确定问题了，那楼主就改吧

我换过了几个缺测值，可是结果还是不理想

river · 发表于 2015-11-14 17:47:53

suifeng1106 发表于 2015-11-14 17:17
我换过了几个缺测值，可是结果还是不理想

你的计算过程使用了fortran，那缺测值是哪个难道不应该用fortran来看吗，靠猜要猜到什么时候去

拖把 · 发表于 2015-11-14 19:41:43

找到对的缺测值，重新定义

叩容 · 发表于 2015-11-15 10:35:40

在fortran计算里，先定义缺测值，然后在数据里检查，遇到缺测值不计算，直接赋缺测。
具体程序结构如下：

undef=...
do i=1,N
if(x(i)...) then
...=undef
end if
然后是你自己要做的计算
end do

出现这种错误首先是你写程序计算的时候不够严谨，没有考虑到数据缺测时候的处理。

suifeng1106 · 发表于 2015-11-15 11:24:21

叩容发表于 2015-11-15 10:35
在fortran计算里，先定义缺测值，然后在数据里检查，遇到缺测值不计算，直接赋缺测。
具体程序结构如下：
...

谢谢指导我现在试试看！

洗雨 · 发表于 2015-11-17 23:59:44

缺测：32767，-9.99e+08,-9.99e+33,-999，好像写描述时用-9.99e+08

suifeng1106 · 发表于 2015-11-18 22:03:21

洗雨发表于 2015-11-17 23:59
缺测：32767，-9.99e+08,-9.99e+33,-999，好像写描述时用-9.99e+08

没有没有，我一个一个点定位发现了缺测值是-832，虽然已经解决了，还是谢谢你的帮忙

suifeng1106 · 发表于 2015-11-18 22:03:38

洗雨发表于 2015-11-17 23:59
缺测：32767，-9.99e+08,-9.99e+33,-999，好像写描述时用-9.99e+08

没有没有，我一个一个点定位发现了缺测值是-832，虽然已经解决了，还是谢谢你的帮忙

		自动登录	找回密码
密码			立即注册