怎样处理逐日ncep温度数据来计算冷空气频次

river

qxtlyf 发表于 2013-4-19 23:06

                               
登录/注册后可看大图

while(yy

那还用说，肯定是2012嘛，真要到了20012咱们连骨灰都找不到了

提米苏拉

river 发表于 2013-4-19 22:56

                               
登录/注册后可看大图

之前那个有些问题，我又改了一下，应该可以了。
'reinit'
yy=1948

if (math_mod(yy,4)=0&math_mod(yy,400)=0)
yy不应该写成‘yy’吗？

river

提米苏拉 发表于 2013-4-20 20:09

                               
登录/注册后可看大图

if (math_mod(yy,4)=0&math_mod(yy,400)=0)
yy不应该写成‘yy’吗？

不是的，必须写成yy ，要不然会报错的。这是grads中的全局变量，随循环变化，加了引号就是动态数组形式了

river

提米苏拉 发表于 2013-4-20 20:09

                               
登录/注册后可看大图

if (math_mod(yy,4)=0&math_mod(yy,400)=0)
yy不应该写成‘yy’吗？

还是有些问题啊if (math_mod(yy,4)=0&math_mod(yy,400)=0)里面是并列的关系，成立一个就可以。应该是if (math_mod(yy,4)=0|math_mod(yy,400)=0)
其实你这个最好的办法还是使用一个ctl来描述多个文件，来得快也来得好。这样循环很容易出问题

wc619lf

提米苏拉 发表于 2013-4-19 11:23

                               
登录/注册后可看大图

按照你的方法，终于整合了一个文件，但看文件的大小感觉只整合了19年的资料，这是为什么呢？

二进制的文件都可以自己算下大小，来判断数据是否是对的，x格点数*y格点数*t*var*4，最后是总的文件大小，单位是字节

提米苏拉

river 发表于 2013-4-21 16:32

                               
登录/注册后可看大图

还是有些问题啊if (math_mod(yy,4)=0&math_mod(yy,400)=0)里面是并列的关系，成立一个就可以。应该是if ( ...

按我之前的那种笨方法，整合出来一个文件，而且大小按照wc619lf的方法算了也是合适的，但是我用整合的grd文件画图与直接用nc画的图不一致，不知道这是什么原因。我觉得fortran处理nc文件比较麻烦，所以想先整合成grd文件再用fortran计算，你说的ctl描述多个文件，是先批量描述，再open ctl文件吗？可是我最终是想生成grd便于fortran计算的。不好意思，不太懂你的意思啊。

提米苏拉

wc619lf 发表于 2013-4-22 09:35

                               
登录/注册后可看大图

二进制的文件都可以自己算下大小，来判断数据是否是对的，x格点数*y格点数*t*var*4，最后是总的文件大小， ...

万分谢谢，好方法，我记下了。

river

提米苏拉 发表于 2013-4-22 11:20

                               
登录/注册后可看大图

按我之前的那种笨方法，整合出来一个文件，而且大小按照wc619lf的方法算了也是合适的，但是我用整合的grd ...

批量描述就相当于所有的资料在一个文件里了，所以一般的计算什么的处理起来就比较方便了，就不需要再合并了。如果不明白就论坛里搜索一下，很多人都分享过或者实用手册里都有

提米苏拉

river 发表于 2013-4-22 11:58

                               
登录/注册后可看大图

批量描述就相当于所有的资料在一个文件里了，所以一般的计算什么的处理起来就比较方便了，就不需要再合并 ...

刚才找了一下批量描述nc的帖子，发现没有考虑闰年的情况，不太适合啊

提米苏拉

提米苏拉 发表于 2013-4-22 19:06

                               
登录/注册后可看大图

刚才找了一下批量描述nc的帖子，发现没有考虑闰年的情况，不太适合啊

我看错了，应该是考虑了的