爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10640|回复: 22

批量描述itpcas的nc数据

[复制链接]

新浪微博达人勋

发表于 2014-1-8 20:07:50 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
青藏高原所自己生成的nc数据,1981~2010年,每三小时数据,0.1*0.1,以气温为例:temp_ITPCAS-CMFD_V0105_B-01_198101.nc~temp_ITPCAS-CMFD_V0105_B-01_201012.批量描述nc数据主要是缺测值的问题:
第一步,随便打开一个nc数据,sdfopen ^.nc,
第二步,q ctlinfo, 把ctl文件导出
'reinit'
'sdfopen ^temp_ITPCAS-CMFD_V0105_B-01_199010.nc'
'q ctlinfo'
rc=write('^temp.ctl',result)
;
利用这个命令可以直接导出ctl文件,在在此基础上进行改动
捕获.PNG
第三步,q attr,得到具体信息
捕获.PNG
其中的missing_value -32767是实际缺测值,在进行后期运算时还要注意scale_factor 0.01,以及add_offset 273.15
第四步,对导出的ctl文件进行修改,包括dset,undef,options和tdef
之后进行检验。直接打开nc文件temp是开尔文,打开批量ctl文件时气温是摄氏度/0.01,因此要注意gs文件的编写。

我给出了批量的ctl文件以及用ctl文件编写的求日最高气温的一个gs
temp.gs (415 Bytes, 下载次数: 17)

评分

参与人数 1金钱 +15 贡献 +5 收起 理由
mofangbao + 15 + 5

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-1-8 20:09:38 | 显示全部楼层
支持自己,小沙发,还有小板凳
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-8 20:14:25 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-1-8 21:38:19 | 显示全部楼层
andrewsoong 发表于 2014-1-8 20:14
青高所阳坤老师的高分辨率的风压湿温降水辐射?

是滴,就是阳坤老师的那个资料
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-8 21:57:43 | 显示全部楼层
哎哟,顶一个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-8 22:49:28 | 显示全部楼层
谢谢分享,很不错的例子
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-9 21:24:07 | 显示全部楼层
学习了!!!谢谢!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-3 11:39:03 | 显示全部楼层
感谢楼主,好东西
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-3 13:31:35 | 显示全部楼层
多谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-2 21:15:31 | 显示全部楼层
{:eb513:}{:eb513:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表