2024年2月22日发(作者：)
摘要： 查了一顿也没找到fnl数据到底包含多少气象数据，这是用wgrib从fnl数据里面读出来的，自己整理的。希望哪位大侠给指导 下，fnl数据。 1 PRES=Pressure
气压 2 HGT=Geopotential height [gpm] 位势高度 3 TMP=Temp. [K] 温度 4
VVEL=Pressure vertical velocit
-
查了一顿也没找到fnl数据到底包含多少气象数据，这是用wgrib从fnl数据里面读出来的，自己整理的。希望哪位大侠给指导 下，fnl数据。
1 PRES=Pressure 气压
2 HGT=Geopotential height [gpm] 位势高度
3 TMP=Temp. [K] 温度
4 VVEL=Pressure vertical velocity [Pa/s] 气压垂直速度
5 RH=Relative humidity [%] 相对湿度
6 ABSV=Absolute vorticity [/s] 绝对涡度
7 O3MR=Ozone mixing ratio [kg/kg] 臭氧混合比
8 CLWMR=Cloud water [kg/kg] 云水
9 5WAVH=5-wave geopotential height [gpm]
10 UGRD=u wind [m/s] 纬向风
11 VGRD=v wind [m/s] 经向风
12 SPFH=Specific humidity [kg/kg] 比湿
13 PWAT=Precipitable water [kg/m^2] 降雨量
14 VWSH=Vertical speed shear [1/s] 垂向风切变
15 LFTX=Surface lifted index [K] 地面举升指数
16 CAPE=Convective Avail. Pot. Energy [J/kg] 对流可用位能
17 CIN=Convective inhibition [J/kg] 对流抑制能
18 4LFTX=Best (4-layer) lifted index [K]
19 PRMSL=Pressure reduced to MSL [Pa]
20 POT=Potential temp. [K 位温
21 TOZNE=Total ozone [Dobson] 总臭氧
22 CWAT=Cloud water [kg/m^2] 云水
23 SOILW=Volumetric soil moisture [fraction]
24 WEASD=Accum. snow [kg/m^2] 雪深
25 LAND=Land cover (land=1;sea=0) [fraction]土地标示
26 ICEC=Ice concentration (ice=1;no ice=0) [fraction] 冰密集度
27 HPBL=Planetary boundary layer height
行星边界层高度
28TCDC=Total cloud cover [%] 总云水
29 GPA=Geopotential height anomaly [gpm] 位势高度异常
30 5WAVA=5-wave geopot. height anomaly [gpm]
ｇｒａｄｓ处理多个ｃｔｌ
(2010-07-29 15:02:41)
标签：
杂谈
先生成一个ctl
grib1:
C:Program FilesPCGrADSwin32>grib2ctl grib_file > grib_
or
gs>!grib2ctl grib_file > grib_
grib2:
For analyses:
$ g2ctl -O grib2_file >grib2_
$ gribmap -O -i grib2_
$ grads
Landscape mode? (no for portrait):
ga-> open grib2_
For forecasts:
$ g2ctl grib2_file >grib2_
$ gribmap -i grib2_
$ grads
Landscape mode? (no for portrait):
ga-> open grib2_
用grads处理多个相同格式的数据时若单个单个处理非常麻烦，当文件非常多的时候是单个处理是不实际的。下面介绍一种方法;
第一步，在这种情况下可以重新写一个ctl描述文件，其文件变量都和已知的ctl相同，若原来的n文件只是时间不同，那么新描述文件的时间维数是所有原文件的时间的和。同样，若其他维数不同时也用同样的方法处理。
第二步，在第一行之后添加一行： options template 表示多个时间序列原始数据文件想用一个描述文件统一地描述。这些原数据的原文件名由dset定义的形势命名文件名。
第三步，修改dset 的文件名。原路径不变，把文件名用％表示。其中：
％y2 代表两位数年
％y4 代表四位数年
％m1 代表一位或者两位数的月
％m2 代表两位数月（用0补齐1位数）
％mc 3个字符月份的缩写
％d1 1或2位天
％d2 两位天
％h1 1或者2位时
％h2 2位时
例如：
原文件其中之一的文件名为gdas2006050812f00，且所有文件只有天和时的变化
那么新描述文件的文件名为：gdas200605%d2%h2f00
另外如果源文件里有index项的话，需要修改其idx的文件名，假设改成。并用在dos下用gribmap函数生成一个新的idx文件。gribmap -e-i （加绝对路径）
open 就可以打开所有文件。
安装WRF及其组件
2010-09-16 17:53:54 来源:NoBodY 【大 中 小】 评论: 条
摘要： 网上有很多指导安装WRF及其组件的文章，也有很多帮助解决安装过程中各种困难的帖子，但大多分散且不系统，下面我就以过来人的身份把我在PC上安装WRF的步骤及遇到的各种问题和答案放在模式联盟论坛上，希望能帮助那些正在挣扎于安装WRF泥淖中的各位，同时也纪
网上有很多指导安装WRF及其组件的文章，也有很多帮助解决安装过程中各种困难的帖子，但大多分散且不系统，下面我就以过来人的身份把我在PC上安装WRF的步骤及遇到的各种问题和答案放在模式联盟论坛上，希望能帮助那些正在
挣扎于安装WRF泥淖中的各位，同时也纪念下我之前所经历的这种痛苦
1，安装平台 fedora 8
2，所需的各种组件：
(1) PGI 7.15 netcdf4.0 NCL(== ncarg) # 其实ifort的性能比PGI更好，但是其他的都过了，只在安装WRFDA
# 的过程中一直编译不过去，因此暂不介绍用ifort编译WRF的过程
(2) WRFV3.1 WPSV3.0.1 WRFDAV3.0.1(if needed)
(3) RIP4
在安装前需要说明的是，我在linux下自己建了两个目录，一个是Model，这是放源程序的（上述所有组件），另一个是Program，这个是放安装好的程序的
3, 安装PGI
（1）解压缩PGI压缩包
$cd /home/name/Model
$mkdir PGI
$tar -xvf -C PGI
$cd PGI
$./install 出现提示信息
Do you accept these terms?[accept,decline]
accept
1,single system install
2,Network install
1
Install the ACML?[y/n]
n
Installation directory?[/opt/pgi]
/home/name/Program/pgi
Do you wish to install MPICH1?[y/n]
n
Do you want the files in the install directory to be read-only?[y/n]
n
install complete
把与安装版本相适应的拷贝到你安装的目录下
配置.bashrc里的环境变量
export PGI=/home/name/Program/pgi/linux86/7.1-5/
export MANPATH=$MANPATH: $PGI/man
export LM_LICENSE_FILE=/home/name/Program/pgi/
export PATH=$PATH: $PGI/bin
$source .bashrc
为了试验pgi是否安装成功，可以打如下命令
$pgf90
如果有如下提示信息，说明已安装成功
$pgf90-Warning-No files to process
如果大家是在更高版本的fedora下装pgi，估计pgi的版本也会相应的升高，这里是一个在fedora11下安装pgi_9.0-1的步骤
/
4，安装netcdf
（1）解压netcdf-4.0压缩包
$mkdir /home/name/Program/netcdf
$cd /home/name/Model
$tar -xvf netcdf-4.0
$cd netcdf-4.0
$./configure --prefix=/home/name/Program/netcdf FC=pgf90
$make
$make install
配置.bashrc里的环境变量
export NETCDF=/home/name/Program/netcdf
export NETCDF_LIB=$NETCDF/lib
export NETCDF_INC=$NETCDF/include
export PATH=$PATH: $NETCDF/bin
$source .bashrc
试验netcdf安装是否成功
$ncdump
如果出现使用说明则说明安装成功
5，安装ncl（目前的ncl已经集合了ncarg，所以大家没必要装ncarg）
在安装ncl之前，需要装三个外部库
下载地址：/?tid=31&extra=page=1&frombbs=1
$su
password: 如果以root身份登陆可跳过此步
#rpm -i *.rpm *代表库名
#exit
安装好之后，开始安装ncl，如果从论坛ftp下载的ncl,解压出来的bin，include，lib可以直接放到/home/name/Program/ncl下即可
配置.bashrc中关于ncl和ncarg的环境变量
export NCARG_ROOT=/home/name/Program/ncl
export NCARG_LIB=/home/name/Program/ncl/lib
export NCARG_INC=/home/name/Program/ncl/include
export PATH=$PATH:/home/name/Program/ncl/bin
$source .bashrc
如想试验安装是否成功
$ncl 如进入程序就说明已经成功
6，安装WRF
在安装之前，请在.bashrc中配置环境变量
export WRFIO_NCD_LARGE_FILE_SUPPORT=1
$source .bashrc
然后按照/里关于的解决办法（2）修改（务必）
解压WRFV3.1压缩包
$ tar -xvf -C /home/name/Program/
$cd /home/name/Program/WRFV3
$./configure
7
$./compile em_real (这时要用到csh，如果PC上还没有安装，可以从这里下载/)
如果/WRFV3/test/em_real里有和说明WRF安装成功了
7，安装WPS
安装WPS前要先安装三个外部库，在以后识别GRIB2格式的数据时会用得到
这三个库的下载地址：/?tid=31&pid=49&page=1&extra=page=1
分别解压缩这三个压缩包
$tar -xvf
$tar -xvf
$tar -xvf
然后开始安装
$su
password: 输入密码，如一开始就以root进入可跳过这步
#cd jasper-1.701.0
#./configure
#make
#make install
#cd ..
#cd libpng-1.2.12
#./configure
#make
#make install
#cd ..
#cd zlib-1.2.3
#./configure
#make
#make install
#exit
现在可以安装WPS了
$cd /home/name/Model
$tar -xvf -C /home/name/Program/
$cd /home/name/Program/WPS
$cd arch
$vi ts
找到跟PGI compiler有关的那几项，分别是serial NO GRIB2, serial, DM
parallel NO GRIB2, DM parallel
修改其中的NCARG_LIBS2项，修改后的结果是：
NCARG_LIBS2 = -L$(PGI)/lib -lpgftnrtl -lpgc
-L/usr/lib/gcc/i386-redhat-linux/4.1.2 -lg2c
需要说明的是，并不是所有的系统都要这么改，修改的原则是与用户系统的路径一致即可，具体可通过进入/usr/lib后找关于gcc及子目录
如果安装前已经打算以后只编译上述四项中的一项或几项，可相应修改，或者一次性全部改过来，一劳永逸
$vi preamble
修改NCARG_LIBS选项，修改后结果如下：
NCARG_LIBS = -L$(NCARG_ROOT)/lib -lncarg -lncarg_gks
-lncarg_c
-L/usr/lib -lX11
如/usr目录下有X11R6则无须做任何改动
这样安装WPS前的准备工作就做完了，下面开始正式安装
$cd ..
$./configure
1
$./compile
当目录下出现, , and 时，说明安装成功了
8，安装WRFDA
安装WRFDAV3前要先装三个外部库bufr（如不用bufr格式的数据，可以不用装，这是3.1版本之前的WRFDA所不允许的），lapack，和blas
前一个的可在论坛ftp中下载，后两个的下载地址：
/?tid=31&pid=50&page=1&extra=page=1
$cd /home/name/Model
$tar -xvf -C /home/name/Program
$cd /home/name/Program
$mkdir BUFR
$tar -xvf /home/name/Model/ -C BUFR
$cd BUFR
$cpp -traditional -I. -C -DRECL1 wrdlen.F>wrdlen.f
$cpp -traditional -I. -C -DRECL1 irev.F>irev.f
$ksh
/home/name/Program/WRFDA/var/scripts/compile_external/pgi_i386/setup_bufr
(这时要用到ksh，如果PC上还没有安装，可以从这里下载/)
$cd ..
此时 bufr安装成功
$tar -xvf /home/name/Model/ ./
$cd BLAS
$ksh
/home/name/Program/WRFDA/var/scripts/compile_external/pgi_i386/setup_blas
$cd ..
blas安装成功
$tar -xvf /home/name/Model/ ./
$cd lapack-3.1.1
$cp INSTALL/ ./
修改，把其中g77或关于fortran的文字全部改成pgf90，把OPT改成-O2 -r8
$ksh
/home/name/Program/WRFDA/var/scripts/compile_external/pgi_i386/setup_lapack
lapack安装成功
配置.bashrc环境变量
export BLAS=/home/name/Program/BLAS
export LAPACK=/home/name/Program/lapack-3.1.1
export BUFR=/home/name/Program/BUFR
export WRF_DA_CORE=1
$source .bashrc
现在开始安装WRFDA
$cd /home/name/Program/WRFDA
$./configure wrfda
$./compile all_wrfvar
如果/home/name/Program/WRFDA/var/da 里面有da_,
da_update_及它共26个.exe文件，那就说明WRFDA安装成功了
9， 安装RIP4
解压缩RIP4压缩包
$tar -xvf /home/name/Model/RIP4 -C /home/name/Program/
$cd /home/name/Program/RIP4
$vi Makefile
把linux（第94行）项中的NETCDFLIB和NETCDFINC和LIBS改成相应目录
$make linux
如果生成rip, ripdp_mm5, ripdp_wrfarw等11个.exe文件则说明安装成功了。
配置.bashrc环境变量
export RIP_ROOT=/home/name/Program/RIP4
export PATH=$PATH: $RIP_ROOT
$source .bashrc
到此为止，WRF组件全部安装完成，剩下的就是设置各个组件的namelist并运行可执行程序了
需要说明的是，不同的Linux版本可能略有不同，所以步骤上也可能有差别，如果有如此问题可以发帖，共同讨论
WRF物理过程参数化方案简介(WRF V2)
2010-07-12 22:03:49 来源:未知 【大 中 小】 评论: 条
摘要： 作者：胡向军, 陶健红 ,郑飞 ,王娜,张铁军,刘世祥,尚大成 1 辐射过程参数化 1.1 RRTM长波辐射方案 来自于MM5模式,采用了Mlawer等人的方法。它是利用一个预先处理的对照表来表示由于水汽、臭氧、二化碳和其他气体,以及云的光学厚度引起的长波过程。 1.2 Dudhi
作者：胡向军, 陶健红 ,郑飞 ,王娜,张铁军,刘世祥,尚大成
1 辐射过程参数化
1.1 RRTM长波辐射方案
来自于MM5模式,采用了Mlawer等人的方法。它是利用一个预先处理的对照表来表示由于水汽、臭氧、二化碳和其他气体,以及云的光学厚度引起的长波过程。
1.2 Dudhia 短波辐射方案
来自于MM5模式,采用Dudhia的方法,它是简单地累加由于干净空气散射、水汽吸收、云反射和吸收所引起的太阳辐射通量。采用了Stephens的云对照表。
1.3 Goddard短波辐射方案
它是由Chou和Suarez发展的一个复杂光学方案。包括了霰的影响,适用于云分辨模式。
1.4 Eta Geophysical Fluid Dynamics Laboratory(GFDL)长波辐射方案
这个辐射方案来自于GFDL。它将Fels和Schwarzkopf的两个方案简单的结合起来了,计算了二氧化碳、水汽、臭氧的光谱波段。
1. 5 Eta Geophysical Fluid Dynamics Laboratory(GFDL) 短波辐射方案
这个短波辐射方案是Lacis和Hansen参数化的GFDL版本。用Lacis和Hansen的方案计算大气水汽、臭氧的作用。用Sasamori等人的方案计算二氧化碳的作用。云是随机重叠考虑的。短波计算用到时间间隔太阳高度角余弦的日平均。
2 微物理过程参数化
2.1 Kessler暖云方案
来自于COMMAS模式,是一个简单的暖云降水方案,考虑的微物理过程包括:雨水的产生、降落以及蒸发,云水的增长,以及由凝结产生云水的过程,微物理过程中显式预报水汽、云水和雨水,无冰相过程。
2.2 Purdue Lin方案
微物理过程中,包括了对水汽、云水、雨、云冰、雪和霰的预报,在结冰点以下,云水处理为云冰,雨水处理为雪。所有的参数化项都是在L in等人以及Rutledge和Hobbs的参数化方案的基础上得到的,某些地方稍有修改,饱和修正方案采用Tao的方法。这个方案是WRF模式中相对比较成熟的方案,更适合于理论研究。
2.3 Eta Ferrier方案
此方案预报模式平流项中水汽和总凝结降水的变化。程序中,用一个局域数组变量来保存初始猜测场信息,然后从中分解出云水,雨水,云冰,以及降冰的变化的密度(冰的形式包括雪、霰或冰雹) 。降冰密度是根据存有冰的增长信息的局域数组来估计,其中,冰的增长与水汽凝结和液态水增长有关。沉降过程的处理是将降水时间平均通量分离成格点单元的立体块。这种处理方法,伴随对快速微物理过程处理方法的一些修改,使得方案在大时间步长时计算结果稳定。根据Ryan的观测结果,冰的平均半径假定为温度函数。冰水混合相仅在温度高于-10℃时考虑,而冰面饱和状态则假定在云体低于-10℃。
2.4 WRF Single_Moment_3_class (WSM3) 方案
该方案来自于旧的NCEP3方案的修正,包括冰的沉降和冰相的参数化。和其它方案不同的是诊断关系所使用冰的数浓度是基于冰的质量含量而非温度。方案包括三类水物质:水汽、云水或云冰、雨水或雪。在这种被称为是简单的冰方案里面,云水和云冰被作为同一类来计算。它们的区别在于温度,也就是说当温度低于或等于凝结点时冰云存在,否则水云存在,雨水和雪也是这样考虑的。该方案对于业务模式来说已足够有效。
2.5 WSM5 方案
与WSM3类似的对NCEP5方案进行了修正,它代替了NCEP5版本。
2.6 WSM6 方案
该方案扩充了WSM5方案,它还包括有霰和与它关联的一些过程。这些过程的参数化大多数和Lin等人的方案相似,在计算增长和其它参数上有些差别。为了增加垂直廓线的精度,在下降过程中会考虑凝结/融化过程。过程的顺序会最优化选择,是为了减少方案对模式时间步长的敏感性。和WSM3、WSM5一样,饱和度调节按照Dudhia和Hong等人的方案分开处理冰和水的饱和过程。
2.7 Thompson方案
该方案改进了较早的Reisner方案,还作了广泛的测试。该方案还被用来做理想试验研究和中纬度冬季观测资料的比较。这个方案被设计用来提高冻雨天气情况下航天安全保障的预报。
3 边界层参数化方案
3.1 MM5相似理论近地面层方案
这个方案用了Paulson、Dyer和Webb稳定性函数来计算地面热量、湿度、动力的交换系数。用Beljaars提出的对流速度来加强地面热量和湿度通量。常与MRF或YSU边界层方案联合使用。
3.2 ETA相似理论近地面层方案
基于Monin-Obukhov理论,在水面上,粘性下层显式参数化,在陆地近地面层上,粘性下层则考虑了变化的位势高度对温度和湿度的作用,近地面通量通过迭代途径进行计算,并用Beljaars修正法来避免在不稳定表面层和无风时出现的奇异性。常与ETAM-Y-J TKE边界层方案联合使用。
3.3 Eta Mellor-Yamada-Janjic TKE边界层方案
此方案用边界层和自由大气中的湍流参数化过程代替Mellor-Yamada的2. 5阶湍流闭合模型。这是将用于Eta模式中的Mellor-Yamada-Janjic方案引入该模式的一种边界层方案,它预报湍流动能,并有局地垂直混合。该方案调用SLAB (薄层)模式来计算地面的温度;在SLAB之前,用相似理论计算交换系数,在SLAB之后,用隐式扩散方案计算垂直通量。
3.4 Medium Range ForecastModel (MRF) 边界层方案
该方案在不稳定状态下使用反梯度通量来处理热量和水汽。在行星边界层中使用增强的垂直通量系数,行星边界层高度由临界Richardson 数决定。它利用一个基于局地自由大气Ri的隐式局地方案来处理垂直扩散项。
3.5 YonseiUniversity (YSU) 边界层方案
YSU是MRF边界层方案的第二代。对于MRF增加了处理边界层顶部夹卷层的方法。
4 积云对流参数化
4.1 浅对流EtaKain-Fritsch方案
在Eta模式中对Kain-Fritsch 方案进行了调整,利用了一个简单的云模式伴随水汽的上升和下沉,同时包括了卷入和卷出,以及相对粗糙的微物理过程的作用。
4.2 Betts-Miller-Janjic方案
对Betts-Miller方案进行了调整和改进,在一给定的时段,对热力廓线进行张弛调整,在张弛时间内,对流的质量通量可消耗一定的有效浮力。此方案为对流调整方案,浅对流调整是该方案的重要部分。
4.3 Kain-Frisch方案
此方案是KF方案的修正方案。与老的KF方案一样,此方案也用了一个简单的包含水汽抬升和下沉运动的云模式,包括卷出、卷吸、气流上升和气流下沉现象。
4.4 Grell-Devenyi集合方案
该方案是质量通量类型,用不同的上升、下沉、卷入、卷出的参数和降水率。静态控制的不同结合了动态控制的不同,这是决定云质量通量的方法。
5 陆面过程参数化
5.1 热量扩散方案
基于MM5的5层土壤温度模式,分别是1、2、4、8和16cm,在这些层下温度固定为日平均值。能量计算包括辐射、感热和潜热通量,同时也允许雪盖效应。
5.2 Noah 方案
Noah陆面过程参数化是OSU的后继版,与原先的相比,可以预报土壤结冰、积雪影响,提高了处理城市地面的能力,考虑了地面发射体的性质,这些是OSU所没有的。
5.3 Rapid Update Cycle (RUC) 方案
这个方案有六个土壤层和两个雪层。它考虑了土壤结冰过程、不均匀雪地、雪的温度和密度差异,以及植被效应和冠层水。
6 讨论
1)参数化的选取与模式的分辨率有关,应根据模式网格设计情况选取相适应的参数化方案。如在高分辨率情况下,对流已不再完全是次网格尺度现象,这时应考虑选择合理的纯显式云物理方案 。对于格距小于5km的情况,一般建议不采用积云参数化方案。
2)由于各种参数化方案在设计原理、复杂程度、计算耗费机时和成熟程度等方面存在差异,研究者应根据研究目的和计算条件等情况来综合判断、对比选择。如对中尺度系统的考察,积云参数化需包括湿下沉气流、中上层的云卷出和非降水性浅对
流,显式云物理方案则需同时加入含有水相和冰相的预报方程, 以计入水负荷、凝结蒸发、冻结融化和凝华升华的影响。
3)由于模式在区域预报的效果与参数化方案的适应性至关重要,目前很多参数化方案对中小尺度系统描述能力不足,所以开发适合本地域特点的参数化方案意义重大。
发布评论