NCEP/NCAR再分析数据集是由美国⽓象环境预报中⼼(NCEP)和美国国家⼤⽓研究中⼼(NCAR)联合制作的,他们采⽤了当今最先进的全球资料同化系统和完善的数据库,对各种来源(地⾯、船舶、⽆线电探空、测风⽓球、飞机、卫星等)的观测资料进⾏质量控制和同化处理,获得了⼀套完整的再分析资料集,它不仅包含的要素多,范围⼴,⽽且延伸的时段长,是⼀个综合的数据集。依其统计时段不同,分以下⼏个类别。 1.NCEP/NCAR逐⽇再分析资料
提供者:NOAA-CIRES Climate Diagnostics Center
载体和容量:硬盘80GB
覆盖时段: 1948~2004年
覆盖地区:全球 (2.5º*2.5º)
内容简介:
NCEP/NCAR 逐⽇再分析资料有变量71个,全部资料数据按性质和变量归类,按年代排列,分为五⼤部分:
(1)等压⾯资料7要素
(2)地⾯资料11要素
(4)对流层顶资料2要素
(5)其它通量资料12要素
详见资料通讯第22期
2.NCEP/NCAR ⼀⽇四次再分析资料
提供者:NOAA-CIRES Climate Diagnostics Center
载体和容量:硬盘270GB河北大学人民武装学院
覆盖时段: 1948~2004年
覆盖地区:全球 (2.5º*2.5º)
内容简介:NCEP/NCAR ⼀⽇四次(00,06,12,18)再分析资料数据集的要素内容与NCEP/NCAR逐⽇再分析资料类同,包括
(1) 等压⾯资料7要素
(2)地⾯资料11要素
(3)地⾯通量资料,42要素
清贫思想详见资料通讯第23期
3. NCEP/NCAR⽉平均再分析资料
提供者:NOAA-CIRES Climate Diagnostics Center
载体和容量:硬盘5.5GB
覆盖时段: 1948~2003年
覆盖地区:全球 (2.5º*2.5º)
内容简介: NCEP/NCAR⽉平均再分析资料内容包括⽓候⽇、⽉平均,历年各⽉平均,标准差,四季昼夜距平,要素场正交展开的主分量等等
1.历年⽉平均(Monthly mean)是由1948年到2003年的每⽉⼀⽇4次资料计算所得。
2.长期⽉平均(Monthly Long-Term means)⼀般是由1968年到1996年的⽉平均资料计算所得。
3.长期⽇平均(Daily Long-Term mean)是由1968年到1996年的⽇平均资料计算所得,其中2⽉29⽇不予计算。
4.年际标准差(Interannual Standard Deviation)是由历年⽉平均对长期⽉平均求均⽅差,它表⽰历年平均相对于长期平均的平均偏离程度。
5.昼夜距平(Diurnal Cycle Anomaly)是由⼀⽇4次(00、06、12、18)资料计算所得。
美国国家环境预报中⼼(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.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,表⽰1983年地⾯温度。
资料内容:各资料的说明见表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 通量资料说明
变量名缩写和层物理量单位
pres.hcb ⾼云底⽓压 10Pa
pres.hct ⾼云顶⽓压 10Pa
pres.lcb 低云底⽓压 10Pa
pres.lct 底云顶⽓压 10Pa
wdc2010>物探与化探
< 中云顶⽓压 10Pa
tcdc.eatm 总云量 %
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位势⾼度场,进⼊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 〈回车〉(设置输出⽂件名为:d)
ga>set gxout fwrite 〈回车〉(设置输出)
ga>d hgt 〈回车〉(d表⽰显⽰,hgt是变量缩写)
此时,⽂件:d中包含了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
.
.
.
STOP
END