NCEP资料简介

  • 格式:doc
  • 大小:48.50 KB
  • 文档页数:5

美国国家环境预报中心(NCEP)和国家大气研究中心(NCAR)联合执行的全球大气40年资料再分析计划通过CDC(Climate Dianogistic Center)利用磁带的形式向外发行。

南京大气资料服务中心通过NCEP朱跃建获得了磁带形式的40年再分析逐日资料。

现在把资料的基本情况作一简单的介绍。

1 资料分类该资料集分:等压面资料、地面资料、通量资料1.1等压面资料资料格距:2.50*2.50的经纬网格网格点数:144* 73个格点资料范围:900N~900S,O0E~357.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天第1~17层的值、…、该年最后一天第1~17层的值。

表1 等压面资料说明变量名缩写物理量单位air 温度 0.1Khgt 位势高度 mrhum 对湿度 %shum 比湿 0.00001 kg/kgomega 垂直速度 0.001 Pa/suwnd 纬向风速 0.1m/svwnd 经向风速 0.1m/s1.2地面资料资料格距:2.50*2.50的经纬网格网格点数:144*73个格点资料范围:900N~900S,00E~357.50E资料层:地面或近地层(0.995层)资料文件:文件名由变量名的缩写、层和年份组成,如:air.Sig995.83,表示1983年地面温度。

资料内容:各资料的说明见表2资料存放方式:每一个要素一年为一个数据文件;在文件中,先存放第1天的值、第2天的值、…、该年最后一天的值。

表2 地面资料说明变量名缩写和层物理量单位air.sig9 95 温度 0.1Klftx.sfc 地面抬升指数 0.1Klftx4.sfc 最佳(4层)地面抬升指数 0.1Komega.sig995 垂直速度 0.001Pa/spottmp.sig995 位温 0.1Kpr-wtr.eatm(整层气柱)可降水量 0.1kg/m2pres.sfc 地面气压 10Parhum.sig995 相对湿度 %slp 海平面气压 10Pauwnd.Sig995 纬向风速 0.1m/svwnd.sig995 经向风速 0.1m/shgt.sfc 地形高度 mland 海陆分布1.3通量资料资料网格:T62高斯格点,192*94个格点资料范围:88.5420N~88.5420S,00E~358.1250E资料文件:文件名由变量名的缩写和年份组成,如:pres.hcb.83,表示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 中云顶气压 10Patcdc.eatm 总云量 %ulwrf.ntat 向上长波辐射通量 W/m2uswrf.ntat 向上短波辐射通量 W/m22 资料使用以上资料用netCDF(Netware Communication Data Format)格式以压缩二进制的形式存储。

在机器上安装支持netCDF格式的GrADS绘图系统,我们以两个简单的例子来说明如何使用该资料集。

例1:显示1960年5月3日500 hPa高度场,假设该资料文件:hgt.60已在当前目录,进入GrADS:ga>sdfopen hgt.60 〈回车〉(扫开支件)ga>set z 6<回车>(第6层为500 hPa)ga>set t 123<回车>(第123日为5月3日)ga>d hgt<回车>(d表示显示水,hgt是变量缩写)此时屏幕显示的是1960年5月3日500hPa高度图。

例2:读出1960年5月3日~10日共8日的北半球500hPa位势高度场,进入GrADS:ga>sdfopen hgt.60<回车>(打开文件)ga>set x l l44<回车>(X方向整个纬圈)ga>set y 37 73<回车>(纬度0~900N)ga>set z 6<回车>(第6层为500hpa)ga>set t 123 130 〈回车〉(第123~130日为5月3~10日)ga>set fwrite h500.grd 〈回车〉(设置输出文件名为:h500.grd)ga>set gxout fwrite 〈回车〉(设置输出)ga>d hgt 〈回车〉(d表示显示,hgt是变量缩写)此时,文件:h500.grd中包含了l960年5月3日~10日共8日的北半球500hpa 位势高度场资料,相应的可以用如下的FORTRAN程序把它读出:DIMENSION H(144,37,8)OPEN(1,FILE=H500.GRD,FORM=UNFORMATTED,ACCESS=DIRECT,RECL=144*37*8*4)READ(1,REC=1)H . . .STOPEND继美国环境预报中心(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个降水正异常年份。

)NcEP/NCAR全球逐6h再分析资料一览表光盘编号变量名文件名文件容量(Mb)单位比率A 气温 air 521 K 10-1B 位势高度 hgt 521 m 100C 相对湿度 rhum 245 % 100C 比湿 shum 245 kg.kg-1 10e-5D 垂直速度 omega 367 Pa.s-1 10e-1E 纬向风 uwnd 521 m.s-1 10e-1F 经向风 vwnd 521 m.s-1 10e-1写论文,参考文献的修改很麻烦,删除一个,添加一个,就需要改一长串数字。

怎么办呢。

本人推荐一种简单方法:尾注法方法如下(以Word2003为例): 1.光标移到要插入参考文献的地方,菜单中“插入”——“引用”-“脚注和尾注”。

2.对话框中选择“尾注”,编号方式选“自动编号”,所在位置建议选“节的结尾”。

3.如“自动编号”后不是阿拉伯数字,选右下角的“选项”,在编号格式中选中阿拉伯数字。

4.确定后在该处就插入了一个上标“1”,而光标自动跳到文章最后,前面就是一个上标“1”,这就是输入第一个参考文献的地方。

5.将文章最后的上标“1”的格式改成正常(记住是改格式,而不是将它删掉重新输入,否则参考文献以后就是移动的位置,这个序号也不会变),再在它后面输入所插入的参考文献(格式按杂志要求来慢慢输,好像没有什么办法简化)。

6.对着参考文献前面的“1”双击,光标就回到了文章内容中插入参考文献的地方,可以继续写文章了。

7.在下一个要插入参考文献的地方再次按以上方法插入尾注,就会出现一个“2”(Word已经自动为你排序了),继续输入所要插入的参考文献。

8.所有文献都引用完后,你会发现在第一篇参考文献前面一条短横线(页面视图里才能看到),如果参考文献跨页了,在跨页的地方还有一条长横线,这些线无法选中,也无法删除。

这是尾注的标志,但一般科技论文格式中都不能有这样的线,所以一定要把它们删除。

9.切换到普通视图,菜单中“视图”——“脚注”,这时最下方出现了尾注的编辑栏。

10.在尾注右边的下拉菜单中选择“尾注分隔符”,这时那条短横线出现了,选中它,删除。

11.再在下拉菜单中选择“尾注延续分隔符”,这是那条长横线出现了,选中它,删除。

12.切换回到页面视图,参考文献插入已经完成了。

这时,无论文章如何改动,参考文献都会自动地排好序了。

如果删除了,后面的参考文献也会自动消失,绝不出错。

13.参考文献越多,这种方法的优势就体现的越大。

还没完,标号上的方括号如何加呢?很简单:在全文中,查找尾注标记^e,然后全部替换为[^&]即可;如果用了脚注就查找脚注标记^f,再全部替换为[^&]便可以了(注意查找时让“不限定格式”按钮为灰色)别看步骤多,操作一遍,就知道很简单了。