- 积分
- 7945
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-11-13
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2013-10-8 18:55:12
|
显示全部楼层
本帖最后由 lxlxllx89 于 2013-10-8 19:10 编辑
给大家简单举个例子:
假设:提取若干年的逐日冬季资料(这个对大家来说是个比较常见的问题)
困难:冬季跨年且存在闰年的问题
情景:1948-2013年逐日高度场资料:hgt.1948.nc, hgt.1949.nc , .... , hgt.2013.nc
其他解决办法:NCL、Fortran、GrADS可以批量读取并且去掉闰年,或者有条件判断,脚本较长
cdo解决办法如下:
1、将所有年份合并到一起:(参见UsersGuide P25)
cdo cat hgt.*.nc hgt.1948-2013.nc
生成的hgt.1948-2013.nc为1948-2013年所有时次的整合nc文件。
2、提取所有年份的12,1,2月:(参见UsersGuide P36)
cdo selmon,1,2,12 hgt.1948-2013.nc hgt.1948-2013.1212.nc
生成的hgt.1948-2013.1212.nc为1948-2013年所有12-2月的整合nc文件。
3、由于冬季的第一年是从1948年12月开始,那么1948年的1,2月以及2013年的12月份(假设数据到2013年12月份)需要剔除:(参见UsersGuide P36)
cdo seldate,1948-12-01,2013-02-28 hgt.1948-2013.1212.nc hgt.1948-2012.DJF.nc
生成的hgt.1948-2012.DJF.nc 即为1948/1949至2012/2013的逐日冬季资料。
简单三句话,完爆以上所有软件!
有了冬季资料之后我们无论是想做季平均也好,算任何东西也好也只需简简单单一句话,详细方法参见附件,以上例子只是这款软件功能的冰山一角,他可以通过bash脚本达到任何想要达到的效果,甚至包括模式的后处理等等,功能十分强大
|
评分
-
查看全部评分
|