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

river

提米苏拉 发表于 2013-4-17 21:25

                               
登录/注册后可看大图

第一步就是想把60年逐日资料整合到一起，但是ncep资料是一年一个，不知道怎么用grads整合成一个grd文件。 ...

文件名连续可以使用循环来打开或者编写一个批量描述的ctl就可以直接在gs里合并了，提取冬半年的资料也可以在gs脚本完成

提米苏拉

river 发表于 2013-4-17 21:36

                               
登录/注册后可看大图

文件名连续可以使用循环来打开或者编写一个批量描述的ctl就可以直接在gs里合并了，提取冬半年的资料也可以 ...

还是不太清楚啊，有没有类似的例子借鉴一下呢？

river

提米苏拉 发表于 2013-4-17 21:58

                               
登录/注册后可看大图

还是不太清楚啊，有没有类似的例子借鉴一下呢？

批量描述nc文件，自己在论坛里搜一下，很多例子

提米苏拉

river 发表于 2013-4-17 22:00

                               
登录/注册后可看大图

批量描述nc文件，自己在论坛里搜一下，很多例子

'reinit'
year=1948
while(year<=2012)
'sdfopen e:\coldair\%y4.nc'

'set gxout fwrite'
'set fwrite e:\coldair\1948-

2012.grd'
'set lat 30 160'
'set lon 10 80'
'set t 1 23742'
'd air'
year=year+1
endwhile
'disable fwrite'
'reinit'
;
这是我按照你的提示，参考一些帖子编的，但不知道哪里错了，生成不了文件啊？

/猪头哼哼

提米苏拉 发表于 2013-4-18 19:59

                               
登录/注册后可看大图

'reinit'
year=1948
while(year

{:3_41:}上下QQ呗

river

提米苏拉 发表于 2013-4-18 19:59

                               
登录/注册后可看大图

'reinit'
year=1948
while(year

完全没有学到方法嘛，估计你都没怎么看懂这个循环啊，grads基础亟需加强。你的nc资料名字是什么样的，只有年份？应该不是吧，我的都是air.1948.nc   air.1949.nc   air.1950.nc·······以此为例吧：
'reinit'
year=1948
while(year<=2012)
'sdfopen e:\coldair\air.'year'.nc'

'set gxout fwrite'
'set fwrite e:\coldair\1948-2012.grd

'set lat 30 160'
'set lon 10 80'
tt=1
while(tt<=365)
'set t 'tt''
'd air'
tt=tt+1
year=year+1
endwhile
endwhile
'disable fwrite'
;
还有不建议使用'set lat 30 160' 'set lon 10 80'设置经纬度，最好改成set x 和set y这样的。但是上面的gs并不完整，遇到闰年的话就会少一天的资料，所以要在其中再加上一个闰年的判断才算圆满。就交给楼主了，如果楼主还是弄不明白的话，那就去加固一下grads的基础吧

提米苏拉

river 发表于 2013-4-18 20:28

                               
登录/注册后可看大图

完全没有学到方法嘛，估计你都没怎么看懂这个循环啊，grads基础亟需加强。你的nc资料名字是什么样的，只有 ...

我想我把批处理ctl文件与循环合并nc文件弄混了，%y4应该在描述文件中用的。哎，自己确实基础不牢，概念混淆，多谢大神指点，还有，那个合并好的65年的逐日温度资料存放顺序就是从1到23742吧？

提米苏拉

river 发表于 2013-4-18 20:28

                               
登录/注册后可看大图

完全没有学到方法嘛，估计你都没怎么看懂这个循环啊，grads基础亟需加强。你的nc资料名字是什么样的，只有 ...

万分感谢。。。。

river

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

                               
登录/注册后可看大图

我想我把批处理ctl文件与循环合并nc文件弄混了，%y4应该在描述文件中用的。哎，自己确实基础不牢，概念混 ...

是啊，从打开的第一个nc文件的第一天资料为1，然后往后推

﹏末歌尽微凉

看评论  来学习