NCEP资料简介.doc

美国国家环境预报中心（NCEP）和国家大气研究中心（NCAR）联合执行的全球大气40年资料再分析计划通过CDC（Climate Dianogistic Center）利用磁带的形式向外发行。南京大气资料服务中心通过NCEP朱跃建获得了磁带形式的40年再分析逐日资料。现在把资料的基本情况作一简单的介绍。1 资料分类 该资料集分：等压面资料、地面资料、通量资料11等压面资料 资料格距：2.50*2.50的经纬网格 网格点数：144* 73个格点 资料范围：900N900S，O0E357.50E 等压面层：共17层（hPa），1000，925，850，700，600，500，400，300，250，200，150，100，70，50，30，20，10 资料文件：文件名由变量名的缩写和年份组成，如：air.83，表示1983年各等压面温度。资料内容:各资料的说明见表1 资料存放方式：每一个要素一年为一个数据文件；在文件中，先存放第1天第1层（l000）的值、第17层（10）的值；第2天第117层的值、该年最后一天第117层的值。表1 等压面资料说明变量名缩写 物理量 单位air 温度0.1Khgt 位势高度 mrhum 对湿度 %shum 比湿 0.00001 kg/kgomega 垂直速度 0.001 Pa/suwnd 纬向风速 0.1m/svwnd 经向风速0.1m/s 12地面资料 资料格距：250*250的经纬网格 网格点数：144*73个格点 资料范围：900N900S,00E357.50E 资料层：地面或近地层（0995层）资料文件：文件名由变量名的缩写、层和年份组成，如：airSig995.83，表示1983年地面温度。 资料内容：各资料的说明见表2 资料存放方式：每一个要素一年为一个数据文件；在文件中，先存放第1天的值、第2天的值、该年最后一天的值。表2 地面资料说明变量名缩写和层 物理量 单位airsig9 95 温度 0.1Klftx.sfc 地面抬升指数0.1Klftx4.sfc 最佳（4层） 地面抬升指数 0.1Komega.sig995 垂直速度 0.001Pa/spottmp.sig995 位温 0.1Kpr-wtreatm（整层气柱） 可降水量 0.1kg/m2pres.sfc地面气压 10Parhum.sig995 相对湿度 %slp 海平面气压 10PauwndSig995 纬向风速 0.1m/svwnd.sig995 经向风速 0.1m/shgt.sfc 地形高度 m land 海陆分布13通量资料 资料网格：T62高斯格点，192*94个格点 资料范围：88.5420N88.5420S，00E358.1250E 资料文件：文件名由变量名的缩写和年份组成，如：preshcb83，表示1983年高云底的气压。 资料内容：各资料的说明见表3 资料存放方式：每一个要素一年为一个数据文件；在文件中，先存放第1天的值、第2大的值、该年最后一天的值。表3 通量资料说明变量名缩写和层物理量 单位csulf.ntat 净空向上长波通量 W/m2csusf.ntat 净空向上短波通量 W/m2dswrf.ntat 向下短波辐射通量 W/m2pres.hcb 高云底气压 10Papres.hct 高云顶气压 10Papres.lcb 低云底气压 10Papres.lct 底云顶气压10Papres.mcb 中云底气压10Papres.mct中云顶气压 10Patcdceatm 总云量 %ulwrf.ntat 向上长波辐射通量 W/m2uswrf.ntat 向上短波辐射通量 W/m22 资料使用 以上资料用netCDF(Netware Communication Data Format)格式以压缩二进制的形式存储。在机器上安装支持netCDF格式的GrADS绘图系统，我们以两个简单的例子来说明如何使用该资料集。 例1：显示1960年5月3日500 hPa高度场，假设该资料文件：hgt.60已在当前目录，进入GrADS: gasdfopen hgt.60 回车 （扫开支件） ga>set z 6<回车>（第6层为500 hPa） gaset t 123<回车>（第123日为5月3日） ga>d hgt<回车>（d表示显示水，hgt是变量缩写） 此时屏幕显示的是1960年5月3日500hPa高度图。 例2：读出1960年5月3日10日共8日的北半球500hPa位势高度场，进入GrADS: gasdfopen hgt60<回车>（打开文件） gaset x l l44<回车>（X方向整个纬圈） gaset y 37 73<回车>（纬度0900N） gaset z 6<回车>（第6层为500hpa） gaset t 123 130 回车 （第123130日为5月310日） gaset fwrite h500grd 回车 （设置输出文件名为：h500.grd） gaset gxout fwrite 回车 （设置输出） gad hgt 回车 （d表示显示，hgt是变量缩写） 此时，文件：h500.grd中包含了l960年5月3日10日共8日的北半球500hpa位势高度场资料，相应的可以用如下的FORTRAN程序把它读出: DIMENSION H(144,37,8) OPEN(1，FILE=H500GRD，FORMUNFORMATTED，ACCESSDIRECT，RECL=144*37*8*4） READ（1，REC=1）H . . . STOP END 继美国环境预报中心(NCEP)和国家大气研究中心(NCAR)联合推出月平均再分析资料（1949-01）和日平均再分析资料（1958-01-01）后，最近又正式对外发布逐6h再分析资料集（1958-01-01T00），这为研究天气尺度和中尺度系统变化过程提供了良好的条件，并能为中尺度数值模式及区域气候模式提供初始场和侧边界条件。 该资料以netCDF(netware Common Data Form)的形式存储，通常可用COLA发展的自由软件GrADS或PMEL发展的可视化自由软件FERRET和RSI提供的商用软件包IDL进行操作处理，逐6h再分析资料集包括7个基本气象要素，根据变量名逐年生成独立文件。气温、位势高度、纬向风和经向风文件容量均为521Mb/a（17层）；相对湿度和比湿文件容量均为245Mb/a(8层）；垂自速度文件容量为367Mb/a（12层）。因此全部42年7要素逐6h再分析资料需要124Gb的介质存储，资料量极其庞大。中心近期已获得1969、1980、1991和1998年完整的资料集（注：这些年份是在近30年来长江流域4个降水正异常年份。）NcEPNCAR全球逐6h再分析资料一览表 光盘编号 变量名 文件名 文件容量（Mb） 单位比率 A气温 air 521K10-1 B 位势高度 hgt 521 m 100 C 相对湿度 rhum 245 % 100 C 比湿 shum 245 kg.kg-1 10e-5 D 垂直速度 omega 367Pa.s-110e-1 E纬向风 uwnd 521m.s-1 10e-1 F 经向风 vwnd 521m.s-1 10e-1 写论文，参考文献的修改很麻烦，删除一个，添加一个，就需要改一长串数字。怎么办呢。本人推荐一种简单方法：尾注法方法如下（以Word2003为例）：1光标移到要插入参考文献的地方，菜单中“插入”“引用”-“脚注和尾注”。 2对话框中选择“尾注”，编号方式选“自动编号”，所在位置建议选“节的结尾”。 3如“自动编号”后不是阿拉伯数字，选右下角的“选项”，在编号格式中选中阿拉伯数字。 4确定后在该处就插入了一个上标“1”，而光标自动跳到文章最后，前面就是一个上标“1”，这就是输入第一个参考文献的地方。 5将文章最后的上标“1”的格式改成正常（记住是改格式，而不是将它删掉重新输入，否则参考文献以后就是移动的位置，这个序号也不会变），再在它后面输入所插入的参考文献（格式按杂志要求来慢慢输，好像没有什么办法简化）。 6对着参考文献前面的“1”双击，光标就回到了文章内容中插入参考文献的地方，可以继续写文章了。 7在下一个要插入参考文献的地方再次按以上方法插入尾注，就会出现一个“2”（Word已经自动为你排序了），继续输入所要插入的参考文献。 8所有文献都引用完后，你会发现在第一篇参考文献前面一条短横线（页面视图里才能看到），如果参考文献跨页了，在跨页的地方还有一条长横线，这些线无法选中，也无法删除。这是尾注的标志，但一般科技论文格式中都不能有这样的线，所以一定要把它们删除。 9切换到普通视图，菜单中“视图”“脚注”，这时最下方出现了尾注的编辑栏。 10在尾注右边的下拉菜单中选择“尾注分隔符”，这时那条短横线出现了，选中它，删除。 11再在下拉菜单中选择“尾注延续分隔符”，这是那条长横线出现了，选中它，删除。 12切换回到页面