NCEP资料的介绍

【引用】NCEP资料的介绍来源：张岳军的日志

美国国家环境预报中心（NCEP）和国家大气研究中心（NCAR）联合执行的全球大气4 0年资料再分析计划通过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，25 0，200，150，100，70，50，30，20，10

资料文件：文件名由变量名的缩写和年份组成，如：air.83，表示1983年各等压面温度。

资料内容:各资料的说明见表1

资料存放方式：每一个要素一年为一个数据文件；在文件中，先存放第1天第1层（l000）的值、…、第17层（10）的值；第2天第1～17层的值、…、该年最后一天第1～17层的值。

1．2地面资料

资料格距：2．50 * 2．50的经纬网格

网格点数：144 * 73个格点

资料范围：900N～900S,00 E～357.50E

资料层：地面或近地层（0．995层）

表1 等压面资料说明

变量名缩写物理量单位

air 温度0.1K

hgt 位势高度m

rhum 相对湿度%

shum 比湿0.00001 kg/kg

omega 垂直速度0.001 Pa/s

uwnd 纬向风速0.1m/s

vwnd 经向风速0.1m/s

资料文件：文件名由变量名的缩写、层和年份组成，如：air．Sig995.83，表示19 83年地面温度。

资料内容：各资料的说明见表2

资料存放方式：每一个要素一年为一个数据文件；在文件中，先存放第1天的值、第2天的值、…、该年最后一天的值。

表2 地面资料说明

变量名缩写和层物理量单位

air．sig995 温度0.1K

lftx.sfc 地面抬升指数0.1K

lftx4.sfc 最佳（4层）地面抬升指数0.1K

omega.sig995 垂直速度0.001Pa/s

pottmp.sig995 位温 0.1K

pr-wtr．eatm 可降水量（整层气柱）0.1kg/m2

pres.sfc 地面气压10Pa

rhum.sig995 相对湿度%

slp 海平面气压10Pa

uwnd．Sig995 纬向风速 0.1m/s

vwnd.sig995 经向风速0.1m/s

hgt.sfc 地形高度m

land 海陆分布

1．3通量资料

资料网格：T62高斯格点，192 * 94个格点

资料范围：88.5420N～88.5420S，00E～358.1250E

资料文件：文件名由变量名的缩写和年份组成，如：pres．hcb．83，表示1983年高云底的气压。

资料内容：各资料的说明见表3

资料存放方式：每一个要素一年为一个数据文件；在文件中，先存放第1天的值、第2大的值、…、该年最后一天的值。

表3 通量资料说明

变量名缩写和层物理量单位

csulf.ntat 净空向上长波通量W/m2

csusf.ntat 净空向上短波通量W/m2

dswrf.ntat 向下短波辐射通量W/m2

pres.hcb 高云底气压10Pa

pres.hct 高云顶气压10Pa

pres.lcb 低云底气压10Pa

pres.lct 底云顶气压10Pa

pres.mcb 中云底气压10Pa

pres.mct 中云顶气压10Pa

tcdc．eatm 总云量%

ulwrf.ntat 向上长波辐射通量W/m2

uswrf.ntat 向上短波辐射通量W/m2

2 资料使用

以上资料用netCDF(Netware Communication Data Format)格式以压缩二进制的形式存储。

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

例1：显示1960年5月3日500 hPa高度场，假设该资料文件：hgt.60已在当前目录，进入CrADS:

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位势高度场，进入Gr ADS:

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’，R ECL=144*37*8*4）

READ（1，REC=1）H