立即注册 登录
气象家园 返回首页

zaozedi的个人空间 http://bbs.06climate.com/?9149 [收藏] [复制] [分享] [RSS]

日志

NCEP/NCAR再分析计划提供的资料介绍

已有 685 次阅读2015-4-29 11:32 | 计划, 资料

美国国家环境预报中心(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层)

 等压面资料说明

变量名缩写               物理量                   单位

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天的值、…、该年最后一天的值。

 地面资料说明

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

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                 地形高度                        
land                    海陆分布

1.3通量资料
    资料网格:T62高斯格点,192 * 94个格点
    资料范围:88.5420N~88.5420S,00E~358.1250E
    资料文件:文件名由变量名的缩写和年份组成,如:pres.hcb.83,表示1983年高云底的气压。
    资料内容:各资料的说明见表3
    资料存放方式:每一个要素一年为一个数据文件;在文件中,先存放第1天的值、第2大的值、…、该年最后一天的值。

 通量资料说明

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

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

 资料使用
    以上资料用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   〈回车〉  (设置输出文件名为: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
      .
    .
    .
    STOP
    END

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 立即注册

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部