WRFnamelist说明mzm.docx
《WRFnamelist说明mzm.docx》由会员分享,可在线阅读,更多相关《WRFnamelist说明mzm.docx(30页珍藏版)》请在冰豆网上搜索。
WRFnamelist说明mzm
关于WRFnamelist的说明
ynkm_mzm(2014年5月)
1、
&share
wrf_core='ARW',(选内核,ARW用于科学研究,NMM用于业务预报)
max_dom=2,(模拟区域总数,包括母区域,1无嵌套,2-4表示嵌套区域及母区域的总数)
start_date='2006-08-16_12:
00:
00','2006-08-16_12:
00:
00',(输入资料的起始时间)
end_date='2006-08-16_18:
00:
00','2006-08-16_18:
00:
00',(输入资料的终止时间,同上)
注:
对于模拟起止时间,亦可用变量START_YEAR(四位数年)、START_MONTH(两位数月)、START_DAY(两位数日)、START_HOUR(两位数小时)以及END_YEAR、END_MONTH、END_DAY、END_HOUR来表示,若两种时间均指定,则start_date优先被采用,且均为UTC时间。
interval_seconds=21600(前处理程序的两次分析时间之间的时间间隔,以秒为单位,即模式的实时输入数据的时间间隔,一般为输入初始场文件的时间间隔。
)
、
io_form_geogrid=2(存数据的格式:
1forbinary;2forNetCDF;3forGRIB1。
缺省通常为2,ARWUserGuide上表述为geogrid产生的数据文件格式)
OPT_OUTPUT_FORM_GEODRID_PATH=’./’(geogrid输出的数据文件存储路径,通常不需指定,即用当前路径’./’)
/
&geogrid(建立“静态的”地面数据,定义模式域)
parent_id=1,1(嵌套区域的母区域的标号。
注意MOAD本身没有母区域,因此PARENT_ID的第一列总是设为1。
第二列必须等于1。
总列数必须等于NUM_DOMAINS)
parent_grid_ratio=1,3(嵌套时,母网格相对于嵌套网格的水平网格比例。
在真实大气方案中,此比例必须为奇数;在理想大气方案中,如果将反馈选项feedback设置为0的话,则此比例也可以为偶数)
i_parent_start=1,30(嵌套网格的左下角(LLC)在上一级网格(母网格)中x方向的起始位置)
j_parent_start=1,19(嵌套网格的左下角(LLC)在上一级网格(母网格)中y方向的起始位置)
s_we=1,1(x方向(西-东方向)的起始格点值,总是1)
s_sn=1,1(y方向(南-北方向)的起始格点值,总是1)
e_we=74,112(x方向(西-东方向)的终止格点值(通常为x方向的格点数),嵌套网格的格点数减去1要能被3整除)
e_sn=59,100(y方向(南-北方向)的终止格点值(通常为y方向的格点数),同上)
注:
关于嵌套格点计算如下图:
geog_data_res=‘10m’,‘2m’(区域对应选择的地表面静态数据,解像度分为10m,5m,2m,30s,也即大约分别为110km,55km,20km,1km)
dx=30000,10000(指定x方向的格距(单位为米))
dy=30000,10000(指定y方向的格距(单位为米))
注:
设置嵌套模拟时,只需设置首层区域的dx与dy值,嵌套区域的格点距将根据parent_grid_ratio自动计算,若设置嵌套网格的格点距则会出错。
对于ARW则‘polar’、’lambert’、’mercator’投影方式以米为单位,’lat-lon’投影以度为单位,对于NMM则以度为单位。
map_proj=‘lambert’(指定模拟区域地图投影方式)
注:
对于ARW可选择的方式有:
"polar"->极射赤面投影(以平面切割60°N,南极的光源将地球赤道以北的区域投影到地图上,适合高纬度地区);"lambert"->兰伯托等角投影(圆锥体与地球30°N和60°N相割,地心的光源投影北半球的区域,适合中维度地区);"mercator"->麦卡托投影(圆柱体与赤道正交从南北回归线切割,地心的光源投影南北回归线之间的区域(适合低纬度地区),"lat-lon"->等经纬度投影。
~
ref_lat=(模拟区域的中心纬度)
ref_lon=103(模拟区域的中心经度)
truelat1=(真实纬度1,lambert投影为第一个与地球相割的维度)
truelat2=(真实纬度2,lambert投影为第二个与地球相割的维度,其他投影则忽略truelat2而以truelat1指定唯一纬度)
stand_lon=103(标准经度,常设为与ref_lon一致)
geog_data_path=‘/…/geog‘(放置地表面静态数据路径)
/
&ungrib(解码GRIB气象数据,并归纳成一个中间格式文件)
out_format='WPS'(输出数据格式,通常为WPS,不需修改,若指定为SI或者MM5则生成相应的格式数据)
prefix='FILE'(指定生成的中间格式文件前缀,亦可加入路径信息,通常指定为FILE,不需修改)
—
/
&metgrid(将准备好的气象数据水平插入模式域)
fg_name='FILE'(ungrib中间格式文件的前缀或者路径,通常与ungrib一致为FILE,不需更改)
io_form_metgrid=2,(2=netCDF,megtgrid输出数据格式)
opt_output_form_metgrid_path=’./’(指定metgrid输出数据文件路径,缺省为当前路径)
/
&mod_levs
press_pa=201300,200100,100000,95000,90000,85000,80000,75000,70000,65000,60000,55000,50000,45000,40000,35000,30000,25000,20000,15000,10000,5000,1000
/
设置完成以上后,若不通过WRFDomainWizard进行预处理,则直接在终端执行如下步骤运行WPS准备模式边界场:
@
$./(建立静态的地面数据)
生成,如果WPS文件夹中没有Vtable,则需要进行下一步:
$ln–sfungrib/Variable_Tables/Vtable
(Variable_Tables中有各种原始资料运行的数据列表,可根据需要选择,Vtable是一个文本文件,内容以表的形式列出初始资料的各种参数配置情况,有此文件才能操作数据)
$.//…/fnl*(读入数据)
$./(解码GRIB气象数据,并归纳成一个intermediate文件。
)
生成FILE*的一系列数据文件
./(把气象数据水平插入模式领域内)
生成以met开头的一系列插值数据,即完成WPS的运行,且met开头的数据为WRF输入数据。
2、
|
&time_control(时间控制)
run_days=0(runtimeindays)
run_hours=12(runtimeinhours)
run_minutes=0(runtimeinminutes)
run_seconds=0(runtimeinseconds)
(注意:
这里的时间相加起来是等于用于运行的数据的时间总长度,若时间超过1天,则可用run_days、run_hours或者run_hours表示,如:
36小时可设置为run_days=1,run_hours=12或者run_hours=36)
start_year(max_dom)=2011(四位数的起始年份)
start_month(max_dom)=05(两位数的起始月份)
start_day(max_dom)=01(两位数的起始日数)
start_hour(max_dom)=00(两位数的起始小时)
start_minute(max_dom)=00(两位数的起始分钟)
[
start_second(max_dom)=00(两位数的起始秒数)
(Note:
thestarttimeisusedtonamethefirstwrfoutfile.Italsocontrolsthestarttimefornestdomains,andthetimetorestart)
end_year(max_dom)=2011(结束时间)
end_month(max_dom)=05
end_day(max_dom)=01,
end_hour(max_dom)=12
end_minute(max_dom)=00
end_second(max_dom)=00
注:
模式积分时间可用run_days、run_hours等来控制,也可用end_year、end_month等来控制,但前者优先于后者,而在中只用后者来控制。
interval_seconds=21600(预处理程序输入的连续两次分析资料之间的时间间隔,也是侧边界条件的时间间隔)
input_from_file(max_dom)=.true.(这是一个逻辑语句,指定嵌套网格是否使用不同的初始场文件,如三重嵌套是否用wrfinput_d02、wrfinput_d03)
~
fine_input_stream(max_dom)=0(选择从嵌套网格中的输入要素场,仅在嵌套网格时有用。
0表示选择从子嵌套网格中输入的所有要素场,2表示在子网格嵌套输入场中仅选择由通道2(在注册表中定义)所指定的那些要素场。
)
history_interval(max_dom)=180,60(指定模式结果输出的时间间隔,以分钟为单位,即每多长时间输出一次模拟结果)
frames_per_outfile(max_dom)=1000,1000(每个输出文件的输出时间,用来把输出文件分成几个小的文件)
restart=.false.(此次运行是否要重新运行,以下用于模式中断后的运行)
restart_interval=500(此参数指定模式断点重启输出的时间间隔,以分钟为单位。
)
io_form_history=2(指定模式结果输出的格式,2=netCDF,1=二进制格式,4=PHDF5格式,5=GRIB1,10=GRIB2,11=parallelnetCDF)
io_form_restart=2(指定模式断点重启输出的格式,2=netCDF,102=splitnetCDFfilesoneperprocessor)
io_form_input=2(2=netCDF,102=allowsprogramtoreadinsplitmet_em*files,andwritesplitwrfinputfiles.Nosplitfileforwrfbdy.)
io_form_boundary=2(指定模式边界条件数据的格式,2为netCDF格式,4为PHD5格式,5为GRIB1格式(目前没有后处理程序),1为二进制格式(目前没有后处理程序))
debug_level=0(指定模式运行时的调试信息输出等级。
取值可为0,50,100,200,300,数值越大,调试信息输出就越多,默认值为0)
auxhist2_outname="rainfall"(指定模式加密输出文件的文件名,缺省时取值为“auxhist2_d_”。
另外,需要指出的是,加密输出变量需要修改注册表文件。
)
<
auxhist2_interval=10(此参数指定模式加密结果输出的时间间隔,以分钟为单位。
)
io_form_auxhist2=2(指定模式加密输出文件的格式,2为netCDF格式)
nocolons=.FALSE.(在输出文件名中是否用下划线“_”代替冒号“:
”。
)
运行3DVAR时需要的额外参数:
write_input=T(指定模式是否输出用于3DVAR的输入数据格式)
inputout_interval=180(此参数指定模式结果输出用于3DVAR的输入数据的时间间隔,以分钟为单位。
)
input_outname=‘wrf_3dvar_input_d_’(指定模式输出用于3DVAR的输入数据文件名,缺省时取值为“wrf_3dvar_input_d_”。
)
inputout_begin_y=0(四位数字表示输出3DVAR数据开始年份。
)
inputout_begin_mo=0(两位数字表示输出3DVAR数据开始月份。
)
inputout_begin_d=0(两位数字表示输出3DVAR数据开始日期。
)
inputout_begin_h=3(两位数字表示输出3DVAR数据开始时次。
)
^
Inputout_begin_m=0(两位数字表示输出3DVAR数据开始分钟数。
)
inputout_begin_s=0(两位数字表示输出3DVAR数据开始秒数。
)
inputout_end_y=0(四位数字表示输出3DVAR数据终止年份。
)
inputout_end_mo=0(两位数字表示输出3DVAR数据终止月份。
)
inputout_end_d=0(两位数字表示输出3DVAR数据终止日期。
)
inputout_end_h=12(两位数字表示输出3DVAR数据终止时次。
)
Inputout_end_m=0(两位数字表示输出3DVAR数据终止分钟数。
)
inputout_end_s=0(两位数字表示输出3DVAR数据终止秒数。
)
说明:
输出用于3DVAR输入数据的时间控制以上面的默认设置为例,模式将从第3时次到第12时次每180分钟输出一次。
&domains(区域控制:
范围,嵌套,参数)
time_step=180(积分的时间步长,为整型数,单位为秒,在真实大气中推荐值为dx公里数的6倍)
#
time_step_fract_num=0(实数型时间步长的分子部分)
time_step_fract_den=1(实数型时间步长的分母部分)
说明:
如果想以秒作为积分时间步长,那么可以设置time_step=60,time_step_fract_num=3,并且设置time_step_fract_den=10。
其中time_step对应于时间步长的整数部分,time_step_fract_num/time_step_fract_den对应于时间步长的小数部分。
注:
在设置history_interval、time_step、parent_time_step_ratio
参数时,相互之间应配置合理。
若当history_interval=15minutes,而time_step=2minutes,则15/2为非整数,在第15minutes时无可用数据,结果将出错,即模拟输出只有第一个时次,故尽量设置history_interval与time_step为相同值。
max_dom=2(计算区域个数,使用嵌套时必须大于1)
s_we(max_dom)=1(x方向(西-东方向)的起始格点值(通常为1))
e_we(max_dom)=74,112(x方向(西-东方向)的终止格点值(通常为x方向的格点数))
s_sn(max_dom)=1(y方向(南-北方向)的起始格点值(通常为1))
e_sn(max_dom)=59,100(y方向(南-北方向)的终止格点值(通常为y方向的格点数))
s_vert(max_dom)=1(z方向(垂直方向)的起始格点值(通常为1))
.
e_vert(max_dom)=28(z方向(垂直方向)的终止格点值,即全垂直eta层的总层数,垂直层数在各嵌套网格中必须保持一致)
num_metgrid_levels=18(来自WPS的metgrid的输入数据的垂直层次数。
一般为WPS的三维变量的层数加上一层地面量,比如三维量是17层,那么总数应该是17+1=18层,可用ncdunp–h查询)
num_metgrid_soil_levels=4(WPS输出的插值的土壤层数)
eta_levels=,,…,(模式的eta层数值,仅用于来自WPS的输入数据。
此eta的数值个数要与模式的垂直层数(e_vert)相一致。
如果缺省,real程序会自动生成一套eta数值。
)
force_sfc_in_vinterp=1(在垂直插值时,在边界层低层,使用地面量作为模式面量的层数。
默认值时只有1层,即最低层使用地面量作为模式面量。
)
p_top_requested=5000(最高等压面,即模式顶气压,单位为帕,在WPS中有效)
interp_type=1(垂直插值的类型:
1,气压线性插值;2,对数气压线性插值)
lagrange_order=1(垂直插值的精度阶数:
1,线性;2,二次)
lowest_lev_form_sfc=.FALSE.(是否使用地面量作为模式最低层的值(u,v,t,q)。
.TRUE.:
使用;.FALSE.:
利用通常的插值方法插值)
dx(max_dom)=30000,1000(x方向网格距)
dy(max_dom)=30000,1000(y方向网格距)
]
ztop(max_dom)=10000(单位为米,此参数指定模式顶的高度。
通常取20000米。
在真实大气方案中,用于高度坐标动力框架模式,此高度值必须与WRFSI的数据或其他输入数据中的高度值相同。
在质量坐标动力框架中,此高度值仅用于理想实验方案)
grid_id(max_dom)=1,2(计算区域编号,一般是从1开始)
parent_id(max_dom)=0,1(嵌套网格的上一级网格(母网格)的编号,一般是从0开始)
i_parent_start(max_dom)=1(嵌套网格的左下角(LLC)在上一级网格(母网格)中x方向的起始位置)
j_parent_start(max_dom)=1,(嵌套网格的左下角(LLC)在上一级网格(母网格)中y方向的起始位置)
parent_grid_ratio=1,3(嵌套时,母网格相对于嵌套网格的水平网格比例。
在真实大气方案中,此比例必须为奇数;在理想大气方案中,如果将返馈选项feedback设置为0的话,则此比例也可以为偶数)
parent_time_step_ratio=1,3(嵌套时,母网格相对于嵌套网格的时间步长比例)
feedback=1(嵌套时,嵌套网格向母网格的反馈作用。
设置为0时,无反馈作用。
而反馈作用也只有在母网格和子网格的网格比例(parent_grid_ratio)为奇数时才起作用)
smooth_option=0(向上一级网格(母网格)反馈的平滑选项,只有设置了反馈选项为1时才起作用的,0:
不平滑;1:
1-2-1平滑;2:
smoothing-desmoothing)
移动嵌套参数:
用户指定移动:
编译时需要在ARCHFLAGS选项中添加“-DMOVE_NESTS”来激活。
允许的最大移动套网格移动次数为50,不过也可以在源程序frame/进行修改。
)
¥
num_moves=4(移动嵌套网格总移动次数。
)
move_id=2,2,2,2,(每一次移动嵌套网格区域编号列表。
)
move_interval=60,120,150,180,(每一次移动的启动时间列表,单位为分钟,自模式积分起始时刻算起。
)
move_cd_x=1,1,0,-1,(在i方向(即东西方向)每一次相对于父网格移动格点数。
)
move_cd_y=1,0,-1,1,(在j方向(即南北方向)每一次相对于父网格移动格点数。
)
注:
正整数表示顺着i/j值增大的方向,负值表示顺着i/j值减小的方向。
0表示不移动。
目前移动距离限制只能为一个网格单元。
自动移动:
编译时需要在ARCHFLAGS选项中添加“-DMOVE_NESTS”和“-DVORTEX_CENTER”来激活。
目前,这些参数是应用中等涡旋追随法(mid-levelvortexfollowingalgorithm)来确定嵌套网格的移动,还在测试阶段。
vortex_interval=15(经过多长时间计算一次涡旋的位置,单位为分钟)
max_vortex_speed=40(涡旋的最大移动速度,用于计算新涡旋位置的搜索半径)
corral_dist=8(移动嵌套网格靠近粗网格边界允许的最大网格单元数,此参数也就是规定了移动网格靠近粗网格允许的最大距离。
)
/
&physics(物理方案)
(
chem_opt(此选项指定是否使用化学过程方案,默认值为0)
mp_physics(max_dom)=2(微物理参数化方案)
0nomicrophysics(不采用微物理过程方案)
1Kessler方案(暖雨方案)
来自于COMMAS模式,是一个简单的暖云降水方案,考虑的微物理过程包括:
雨水的产生、降落以及蒸发,云水的增长,以及由凝结产生云水的过程,微物理过程中显式预报水汽、云水和雨水,无冰相过程。
2Lin等的方案(水汽、雨、雪、云水、冰、冰雹)
微物理过程中,包括了对水汽、云水、雨、云冰、雪和霰的预报,在结冰点以下,云水处理为云冰,雨水处理为雪。
所有的参数化项都是在Lin等人以及Rutledge和Hobbs的参数化方案的基础上得到的,某些地方稍有修改,饱和修正方案采用Tao的方法。
这个方案是WRF模式中相对比较成熟的方案,更适合于理论研究。
3WSM3类简单冰方案
该方案来自于旧的NCEP3方案的修正,包括冰的沉降和冰相的参数化。
和其它方案不同的是诊断关系所使用冰的数浓度是基于冰的质量含量而非温度。
方案包括三类水物质:
水汽、云水或云冰、雨水或雪。
在这种被称为是简单的冰方案里面,云水和云冰被作为同一类来计算。
它们的区别在于温度,也就是说当温度低于或等于凝结点时冰云存在,否则水云存在,雨水和雪也是这样考虑的。
该方案对于业务模式来说已足够有效。
4WSM5类方案
与WSM3类似的对NCEP5方案进行了修正,它代替了NCEP5版本。
@
5Ferrier(newEta)微物理方案(水汽、云水)
此方案预报模式平流项中水汽和总凝结降水的变化。
程序中,用一个局域数组变量来保存初始猜测场信息,然后从中分解出云水,雨水,云冰,以及降冰的变化的密度(冰的形式包括雪、霰或冰雹)。
降冰密度是根据存有冰的增长信息的局域数组来估计,其中,冰的增长与水汽凝结和液态水增长有关。
沉降过程的处理是将降水时间平均通量分离成格点单元的立体块。
这种处理方法,伴随对快速微物理过程处理方法的一些修改,使得方案在大时间步长时计算结果稳定。
根据Ryan的观测结果,冰的平均半径假定为温度函数。
冰水混合相仅在温度高于-10℃时考虑,而冰面饱和状态则假定在云体低于-10℃。
6WSM6类