月平均的sst整理成季平均，输出后画图出错（已解决）

好猫

一个很简单的季平均处理，通过fwrite输出1951-2012年的季平均值，共62年，每年4个值，我存成了4个变量。但是在画图的时候不仅出现诡异的特大值，而且还跑到陆地上去了……
经river的帮助，我用fortran读了一下自己fwrite出来的东西，发现缺测值居然神奇的变了。缺测值并不是原nc文件中的“-9.9e~~blabla”，而是fortran读出来的那个哦~所以如果大家发现图画出来像我的第一张图这个样子，肯定就是缺测值又淘气了！

ctl只要把原来的
undef -1e30
改成
undef 32767
就好！（我的是这样，你的是不是一样呢？自己编个小程读一下吧~）

然后是图……
这张是对的！

river

目测是缺测值问题，改成-9.99e8试试吧

好猫

river 发表于 2013-12-12 20:24

                               
登录/注册后可看大图

目测是缺测值问题，改成-9.99e8试试吧

仍旧是不行
本来还想做一个季平均的时间经度剖面图，写出来了一样有问题
如果不是按这个思路，你觉得还有别的方法吗？

river

好猫 发表于 2013-12-13 18:54

                               
登录/注册后可看大图

仍旧是不行
本来还想做一个季平均的时间经度剖面图，写出来了一样有问题
如果不是按这个思路，你觉得还 ...

你用fortran把你提取出来的资料读出来，输出到屏幕上看看缺测值到底是哪个值？然后改一下ctl里的缺测值试试

好猫

river 发表于 2013-12-13 20:43

                               
登录/注册后可看大图

你用fortran把你提取出来的资料读出来，输出到屏幕上看看缺测值到底是哪个值？然后改一下ctl里的缺测值试 ...

做出来了！坑爹的缺测值啊！在nc文件的ctlinfo里明明是一个值，在fortran读出来的时候原来是32767！
太感谢了~又多了一份经验！

river

好猫 发表于 2013-12-13 21:46

                               
登录/注册后可看大图

做出来了！坑爹的缺测值啊！在nc文件的ctlinfo里明明是一个值，在fortran读出来的时候原来是32767！
太感 ...

这个主要是有运算过程了，所以可能就不一样了。正常只是提取出里面的变量的话，缺测值就应该是grads默认值-9.99e8

好猫

river 发表于 2013-12-13 22:52

                               
登录/注册后可看大图

这个主要是有运算过程了，所以可能就不一样了。正常只是提取出里面的变量的话，缺测值就应该是grads默认值 ...

学习了

lee309

http://bbs.06climate.com/forum.php?mod=viewthread&tid=31317
请教下楼主

supermann

你好楼主，我想问下你是怎么输出的季平均文件的呢，我的也是多年月平均海温，想输出每年的冬季数据