ncep资料介绍 2010-01-21 21:19 美国国家环境预报中心(NCEP)和国家大气研究中心(NCAR)联合执行的全球大气40年资料再分析计划通过CDC(Climate Dianogistic Center)利用磁带的形式向外发行。南京大气资料服务中心通过NCEP朱跃建获得了磁带形式的40年再分析逐日资料。现在把资料的基本情况作一简单的介绍。 表1 等压面资料说明 变量名缩写 物理量 单位 air 温度 0.1K 资料文件:文件名由变量名的缩写、层和年份组成,如:air.Sig995.83,表示1983年地面温度。 表2 地面资料说明 变量名缩写和层 物理量 单位 air.sig995 温度 0.1K 1.3通量资料 表3 通量资料说明 变量名缩写和层 物理量 单位 csulf.ntat 净空向上长波通量 W/m2 2 资料使用 nc格式文件的读取 下面给出使用Grads读取nc格式文件的一个范例脚本文件(.gs文件): 以位势高度场为例(红色字体为解释说明语句;蓝色字体为输入的命令,也就是gs文件的内容) ‘reinit’ ‘sdfopen d:\xxxx\hgt.mon.mean.nc ‘ (打开需要的nc文件,注意路径) ‘set gxout fwrite’ ‘set fwrite d:\xxxxxx\hgt.bin’ (指定要输出的二进制文件名和路径) ‘set lon 0 ‘set lat -90 90 ‘ (设定写出区域范围,这是全球范围,第一句是必需的,如果是用默认,则会在x方向输出为145个点,要注意) ‘set lev ‘set t 1 ‘d hgt’ (输出) ‘reinit’ 这样就生成了需要的二进制文件,其对应的ctl文件需要另外给出。 以上是输出单层的二进制文件,把多层的数据放在同一个二进制文件中,就要使用循环,根据grads规定的数据存放次序(依次为经度、纬度、高度、变量、时次),如下: ‘reinit’ ‘sdfopen d:\xxxx\hgt.mon.mean.nc’ (打开需要的nc文件,注意路径) ‘set gxout fwrite’ ‘set fwrite d:\xxxxxx\hgt.bin’ (指定要输出的二进制文件名和路径) ‘set lon 0 ‘set lat -90 t=a (a为起始时次) while(t<=b) (开始循环,b为结束时次) ‘set t’ t (设定时次) ‘set lev ‘d hgt ‘ (输出) ‘set lev ‘d hgt ‘ (输出) ‘set lev ‘d hgt ‘ (输出) ‘set lev ‘d hgt ‘ (输出) ‘set lev ‘d hgt ‘ (输出) ‘set lev ‘d hgt ‘ (输出) 。 。 。 (可以添加任意的层次) 。 (由于高度分布不连续而且数量有限,所以用这种方法) 。 t=t 1 endwhile (结束循环) ‘reinit’ 继美国环境预报中心(NCEP)和国家大气研究中心(NCAR)联合推出月平均再分析资料(1949-01)和日平均再分析资料( NcEP/NCAR全球逐6h再分析资料一览表 光盘编号 变量名 文件名 文件容量(Mb) 单位 比率 A 气温 air 521 K 10-1 B 位势高度 hgt 521 m 100 C 相对湿度 rhum 245 % 100 C 比湿 shum 245 kg.kg-1 10-5 D 垂直速度 omega 367 Pa.s-1 10-1 E 纬向风 uwnd 521 m.s-1 10-1 F 经向风 vwnd 521 m.s-1 10-1 |
|