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

river

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

                               
登录/注册后可看大图

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

批量处理就是在一个ctl里描述所有的资料,时间的设置就是从第一个资料的起始时间开始到最后一个资料的最后一个时次,所以不需要判断什么闰年

提米苏拉

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

                               
登录/注册后可看大图

批量处理就是在一个ctl里描述所有的资料,时间的设置就是从第一个资料的起始时间开始到最后一个资料的最后 ...

那我批量描述后，就可以直接打开该ctl，再直接写成grd的文件啦？

提米苏拉

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

                               
登录/注册后可看大图

批量处理就是在一个ctl里描述所有的资料,时间的设置就是从第一个资料的起始时间开始到最后一个资料的最后 ...

dset  e:\coldair\daily\%y4.nc
undef 1e+15

options template
title nc data
dtype netcdf
xdef 53 linear 30 2.5
ydef 29 linear 10 2.5
zdef 1 linear  1 1
tdef 23742 linear 01jan1948 1dy
vars 1
air 1 99 temperature            

                  
endvars

'reinit'
'open e:\coldair\air.ctl'
'set gxout fwrite'
'set fwrite e:\coldair\1948-

2012.grd'
'set lat 10 80'
'set lon 30 160'
'set z 1'
'set t 1 23742'
'd air'
'disable fwrite'
'reinit'
;
为什么提示错误呢

river

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

                               
登录/注册后可看大图

dset  e:\coldair\daily\%y4.nc
undef 1e+15

你写完ctl可不可以先在grads里面打开，看看能不能出图，出来的图是不是正确的，然后再进行下一步啊。而且你的air的资料是有一层吗？
air 1 99 temperature 注意红色的部分,他代表的意思是该变量有多少个层次,你写1,难道air的nc资料只有一层吗？一般的都是17层。还有这个资料的缺测值你是怎么判断出来是1e+15，有验证过吗？还有你的这个nc资料只是从东经30度开始到162.5度？纬度只是从北纬10度到北纬41.5吗？一般不都是全球的吗？

上面这些问题不经过确认，那你的ctl很可能完全就是错误的，更别说下一步的东西了。难道没有一个完整的思路吗，一步一步都得经过验证，才能进行下一步。不然的话就是一步错，步步错

提米苏拉

river 发表于 2013-4-22 20:16

                               
登录/注册后可看大图

你写完ctl可不可以先在grads里面打开，看看能不能出图，出来的图是不是正确的，然后再进行下一步啊。而且 ...

好吧，一直觉得自己逻辑感不强，被程序弄的我头昏脑胀的。我再改改吧。还有上次的grd经人指点终于整出来了，原因在于应该打开一个文件关一个文件，否则会一直读取第一年的数据。在最后一个endwhile前面加一个‘close 1’就好了。还要特别感谢river，不厌其烦的帮助我。我以后会尽量改掉这些毛病，否则我未来的两年会很痛苦。。。

提米苏拉

river 发表于 2013-4-22 20:16

                               
登录/注册后可看大图

你写完ctl可不可以先在grads里面打开，看看能不能出图，出来的图是不是正确的，然后再进行下一步啊。而且 ...

还有我这个是sig995层温度资料，只有一层，nc资料的缺测值我是看别的帖子得到的。我也不知道怎么得到nc资料的缺测值。

提米苏拉

river 发表于 2013-4-22 20:16

                               
登录/注册后可看大图

你写完ctl可不可以先在grads里面打开，看看能不能出图，出来的图是不是正确的，然后再进行下一步啊。而且 ...

我知道了，应该 q ctlinfo 来得知nc资料的一些具体信息

river

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

                               
登录/注册后可看大图

我知道了，应该 q ctlinfo 来得知nc资料的一些具体信息

只有一层啊，那挺好的。自己试试吧，我今天没有午休，暴躁的很，我要疯啦快。赶紧睡觉去了

river

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

                               
登录/注册后可看大图

我知道了，应该 q ctlinfo 来得知nc资料的一些具体信息

好像搞定了哦，写成下面那样就能提取成一个grd了。注意改路径啊

1. 'reinit'
   
2. 'set gxout fwrite'
   
3. 'set fwrite  1948-2012.grd'
   
4. 
   
5. yy=1948
   
6. while(yy<=2000)
   
7. 'sdfopen   'yy'.nc'
   
8. 'set lat 30 160'
   
9. 'set lon 10 80'
   
10. if (math_mod(yy,4)=0|math_mod(yy,400)=0)
    
11. tt=1
    
12. while(tt<=366)
    
13. 'set t 'tt''
    
14. 'set z 1'
    
15. 'd air'
    
16. tt=tt+1
    
17. endwhile
    
18. else
    
19. tt=1
    
20. while(tt<=365)
    
21. 'set t 'tt''
    
22. 'set z 1'
    
23. 'd air'
    
24. tt=tt+1
    
25. endwhile
    
26. endif

复制代码

提米苏拉

river 发表于 2013-4-23 19:19

                               
登录/注册后可看大图

好像搞定了哦，写成下面那样就能提取成一个grd了。注意改路径啊

每打开一个文件还要关一个文件，即在最后一个endwhile前加一个‘close 1’才行，否则会一直读取第一年的资料哈，不过还是谢谢你给了这么好的思路。