附录1 WRF模式参数配置说明

1、附录1 WRF模式参数配置说明由 wrfchina 于 星期五,2012-04-06 15:08 提交注意，参数选项名称后跟的(max_dom)是表示此参数需定义成嵌套形式。参数配置第一部分这部分参数仅用于由真实大气方案的预处理程序产生的输入数据。当输入数 据产生于理想大气试验方案时，这部分参数将会被忽略。对于大多数真实大气方 案来说，起止时间的分和秒都应该设为0。常用的小时和秒之间的换算关系有：3小时=10800秒；6小时= 21600秒；12小时= 43200秒。&time_controlrun_days运行的天数run_hours运行的小时数注意：如果模式积分时间大于1天，则可同时设置r

2、un_days和_run_hours，也 可设置run_hours 一个参数。比如：模式运行的总时间长度为36小时，则 可设置 run_days=1， 且 run_hours=12，或者设置 run_days=0， 且 run_hours=36。run_minutes运行的分钟数run_seconds运行的秒数start_year(max_dom) =2001四位数字表示的起始年份。start_month(max_dom) =04两位数字(01-12)表示的起始月份。start_day(max_dom) =20两位数字(01-31)表示的起始天数。start_hour(max_dom) =12

3、两位数字(00-23)表示的起始小时数。start_minute(max_dom) =00两位数字(00-59)表示的起始分钟数。start_second (max_dom) =00两位数字(00-59)表示的起始秒数。end_year(max_dom) =2001四位数字表示的终止年份。end_month(max_dom) = 04两位数字(01-12)表示的终止月份。end_day(max_dom) =21两位数字(01-31)表示的终止天数。end_hour =00两位数字(00-23)表示的终止小时数。end_minute =00两位数字(00-59)表示的终止分钟数。end_seco

4、nd =00两位数字(00-59)表示的终止秒数。说明：起止时间设置也可以用来控制模式的积分的起止。并且，real.exe的时 间控制信息是用起止时间参数来设定的。模式的积分时间可以用run_days、run_hours等来控制，也可以用end_year、 end_month等来控制。但前者run_days等优先与后者end_year等。而在 real.exe中只用end_year等来控制时间信息。intervalseconds =43200前处理程序的两次分析时间之间的时间间隔，以秒为单位。也即模式的实时输入 数据的时间间隔，一般为输入边界条件的文件的时间间隔。input_from_file

5、 (max_dom) = T嵌套初始场输入选项。嵌套时，指定嵌套网格是否用不同的初始场文件。fine_input_stream(max_dom) = 0选择从嵌套网格中的输入要素场，仅在嵌套网格时有用。0表示选择从子嵌 套网格中输入的所有要素场，2表示在子网格嵌套输入场中仅选择由通道2 (在注册表中定义)所指定的那些要素场。history_interval (max_dom) = 60此参数指定模式结果输出的时间间隔，以分钟为单位。frames_per_outfile (max_dom) = 1此参数指定每一个结果文件中保存输出结果的次数，因此可以将模式结果分成多 个文件保存，默认值为10。r

6、estart = F指定模式运行是否为断点重启方式。restart_interval = 1440此参数指定模式断点重启输出的时间间隔，以分钟为单位。io_form_history = 2指定模式结果输出的格式，2为netCDF格式io_form_restart = 2指定模式断点重启输出的格式，2为netCDF格式io_form_initial = 2指定模式初始场数据的格式，2为netCDF格式io_form_boundary = 2 指定模式边界条件数据的格式，2为netCDF格式，4为PHD5格式，5为GRIB1 格式(目前没有后处理程序)，1为二进制格式(目前没有后处理程序)。deb

7、ug_level =0此选项指定模式运行时的调试信息输出等级。取值可为0,50,100,200,300 ,数 值越大，调试信息输出就越多，默认值为0。auxhist2_outname = rainfall指定模式加密输出文件的文件名，缺省时取值为“auxhist2_d_”。另外，需要 指出的是，加密输出变量需要修改注册表文件Registry.EM。auxhist2_interval = 10此参数指定模式加密结果输出的时间间隔，以分钟为单位。io_form_auxhist2 = 2指定模式加密输出文件的格式，2为netCDF格式nocolons = .FALSE.在输出文件名中是否用下划线“_

8、”代替冒号“：”。运行3DVAR时需要的额外参数：write input = T指定模式是否输出用于3DVAR的输入数据格式inputout_interval = 180此参数指定模式结果输出用于3DVAR的输入数据的时间间隔，以分钟为单位。input_outname = wrf_3dvar_input_d_指定模式出用于3DVAR的输入数据文件名，缺省时取值为 “wrf_3dvar_input_d_ ”。inputout_begin_y = 0inputout_begin_mo = 0 两位数字表示输出3DVAR数据开始月份。inputout_begin_d = 0两位数字表示输出3DVAR

9、数据开始日期。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 = 12Inputout_end_m = 0 两位数字表示输出3DVAR数据终止分钟

10、数。inputout_end_s = 0两位数字表示输出3DVAR数据终止秒数。说明：输出用于3DVAR输入数据的时间控制以上面的默认设置为例，模式将从第3时次到 第12时次每180分钟输出一次。参数配置第二部分&domainstime_step = 60积分的时间步长，为整型数，单位为秒。time_step_fract_num = 0实数型时间步长的分子部分。time_step_fract_den = 1实数型时间步长的分母部分。说明：如果想以60.3秒作为积分时间步长，那么可以设置time_step=60, time_step_fract_num=3，并且设置 time_step_frac

11、t_den=10。其中 time_step对应与时间步长的整数部分， time_step_fract_num/time_step_fract_den 对应于时间步长的小数部分。max_dom =1计算区域个数。计算区域默认值为1，如果使用嵌套功能，则max_dom大于 1。s_we（max_dom） =1x方向（西-东方向）的起始格点值（通常为1）.e_we（max_dom） =32x方向（西-东方向）的终止格点值（通常为x方向的格点范围）。s_sn （max_dom） =1y方向（南-北方向）的起始格点值（通常为1）.e_sn （max_dom） =32y方向（南-北方向）的终止格点值（通常

12、为y方向的格点范围）。s_vert （max_dom） =1z方向（垂直方向）的起始格点值。e_vert （max_dom） =31z方向（垂直方向）的终止格点值，即全垂直eta层的总层数。垂直层数在各 嵌套网格中必须保持一致。num_metgrid_levels =18来自WPS的metgrid的输入数据的垂直层次数。一般为WPS的三维变量的 层数加上一层地面量，比如三维量是17层，那么总数应该是17+1 = 18层。eta_levels =1.0, 0.997,，0.0模式的eta层数值，仅用于来自WPS的输入数据。此eta的数值个数要与 模式的垂直层数（e_vert）相一致。如果缺省，r

13、eal程序会自动生成一套 eta数值。force_sfc_in_vinterp =1在垂直插值时，在边界层低层，使用地面量作为模式面量的层数。默认值 时只有1层，即最低层使用地面量作为模式面量。p_top_requested =5000模式的顶部气压，单位为帕。interp_type =1垂直插值的类型：1，气压线性插值；2,对数气压线性插值lagrange_order =1垂直插值的精度阶数：1，线性；2，二次 lowest_lev_form_sfc =.FALSE.是否使用地面量作为模式最低层的值(u,v,t,q)。.TRUE.:使用；.FALSE.： 利用通常的插值方法插值。dx (ma

14、x_dom) =10000 (单位为米)指定x方向的格距。在真实大气方案中，此参数值必须与输入数据中的x 方向格距一致。dy (max_dom) = 10000 (单位为米)指定y方向的格距。通常与x方向格距相同。ztop (max_dom) = 10000 (单位为米)此参数指定模式顶的高度。通常取20000米。在真实大气方案中，用于高 度坐标动力框架模式，此高度值必须与WRF SI的数据或其他输入数据中的 高度值相同。在质量坐标动力框架中，此高度值仅用于理想实验方案。grid_id (max_dom) =1计算区域的编号。一般是从1开始。parent_id (max_dom) = 0嵌套网

15、格的上一级网格(母网格)的编号。一般是从0开始。i_parent_start (max_dom) = 0嵌套网格的左下角(LLC)在上一级网格(母网格)中x方向的起始位置。j_parent_start (max_dom) = 0嵌套网格的左下角(LLC)在上一级网格(母网格)中y方向的起始位置。parent_grid_ratio(max_dom) = 1嵌套时，母网格相对于嵌套网格的水平网格比例。在真实大气方案中，此 比例必须为奇数；在理想大气方案中，如果将返馈选项feedback设置为0 的话，则此比例也可以为偶数。parent_time_step_ratio (max_dom) = 1嵌套

16、时，母网格相对于嵌套网格的时间步长比例。feedback = 嵌套时，嵌套网格向母网格得反馈作用。设置为0时，无反馈作用。而反 馈作用也只有在母网格和子网格的网格比例（parent_grid_ratio）为奇数 时才起作用。smooth_option = 0向上一级网格（母网格）反馈的平滑选项，只有设置了反馈选项为1时才 起作用的。0：不平滑；1： 1-2-1 平滑；2: smoothing-desmoothing移动网格控制参数说明：控制移动嵌套网格的方式有两种：1.用户指定移动网 格，此方式是指嵌套网格的每一次移动都是由用户通过参数来指定；2.自 动移动网格，此方式是指移动嵌套网格在启动后

17、，根据模式的计算状态， 自动判断下一次的移动参数。用户指定移动：编译时需要在ARCHFLAGS选项中添加“一D MOVE_NESTS”来激活。 允许的最大移动套网格移动次数为50，不过也可以在源程序 frame/module_driver_constants.F进行修改。）num_moves = 4移动嵌套网格总移动次数。move_id = 2,2,2,2,每一次移动嵌套网格区域编号列表。move_interval = 60,120,150,180,每一次移动的启动时间列表，单位为分钟，自模式积分起始时刻算起。move_cd_x = 1,1,0,-1,在i方向（即东西方向）每一次相对于父网格移

18、动格点数。move_cd_y = 1,0,-1,1,在j方向（即南北方向）每一次相对于父网格移动格点数。正整数表示顺着i/j值增大的方向，负值表示顺着i/j值减小的方向。0 表示不移动。目前移动距离限制只能为一个网格单元。自动移动：编译时需要在ARCHFLAGS选项中添加“一D MOVE_NESTS”和 “-DVORTEX_CENTER”来激活。目前，这些参数是应用中等涡旋追随法）来确定嵌套网格的移(mid-level vortex following algorithm 动，还在测试阶段。vortex_interval = 15经过多长时间计算一次涡旋的位置，单位为分钟max_vortex_

19、speed = 40涡旋的最大移动速度，用于计算新涡旋位置的搜索半径corral_dist = 8移动嵌套网格靠近粗网格边界允许的最大网格单元数，此参数也就是规定了移动 网格靠近粗网格允许的最大距离。参数配置第三部分虽然不同的嵌套网格可以使用不同的物理方案，但必须注意没中方案的使用条件 和范围。&physicschem_opt此选项指定是否使用化学过程方案，默认值为0。mp_physics （max_dom）此选项设置微物理过程方案，默认值为0。目前的有效选择值为：=0,不采用微物理过程方案=1, Kessler方案（暖雨方案）=2, Lin等的方案（水汽、雨、雪、云水、冰、冰雹）=3, WS

20、M 3类简单冰方案=4, WSM 5类方案=5, Ferrier（new Eta）微物理方案（水汽、云水）=6, WSM 6类冰雹方案=8, Thompson 等方案=98, NCEP 3类简单冰方案(水汽、云/冰和雨/雪)(将放弃) =99, NCEP 5类方案(水汽、雨、雪、云水和冰)(将放弃)新添参数：mp_zero_out = 0,选用微物理过程时，保证Qv .GE. 0,以及当其他一些水汽变量小于临界值 时，将其设置为0。=0,表示不控制，=1,除了Qv外，所有的其他水汽变量当其小于临界值时，则设置为0=2,确保Qv .GE. 0,并且所有的其他水汽变量当其小于临界值时，则 设置为0

21、。mp_zero_out_thresh = 1.e-8水汽变量(Qv除外)的临界值，低于此值时，则设置为0 (kg/kg)。ra_lw_physics(max_dom)此选项指定长波辐射方案，默认值为0。有效选择值如下：=0,不采用长波辐射方案=1, rrtm 方案=99, GFDL (Eta)长波方案(semi-supported)ra_sw_physics (max_dom)此选项指定短波辐射方案，默认值为0。有效选择值如下：=0,不采用短波辐射方案=1, Dudhia 方案=2, Goddard短波方案=99, GFDL (Eta)短波方案(semi-supported)radt (ma

22、x_dom)此参数指定调用辐散物理方案的时间间隔，默认值为0,单位为分钟。通 常比较合理的间隔值为30分钟。当网格水平分辨率提高时，则需将间隔时 间相应地缩短。建议为水平分辨率的1倍，如dx=10km，则取10分钟。nrads (max_dom)=用于NMM版本WRF。指定调用短波辐散过程的时间间隔，单位为粗网格的时 间步数。默认值在Registry.NMM进行设置，但也可在namelist.input文 件中重新赋值。模式会根据此数值计算radt的值。nradl (max_dom)=用于NMM版本WRF。指定调用长波辐散过程的时间间隔，单位为粗网格的时 间步数。默认值在Registry.NM

23、M进行设置，但也可在namelist.input文 件中重新赋值。模式会根据此数值计算radt的值。co2tf = 0CO2的传输函数参数，用于GFDL辐射方案=0,从预生成的数据文件中读取CO2函数数据=1,模式自己生成CO2函数sf_sfclay_physics (max_dom)此选项指定近地面层(surface-layer)方案，默认值为0。旧参数表中的 bl_sfclay_physics，有效选择值有：=0,不采用近地面层方案=1, Monin-Obukhov 方案=2, MYJ Monin-Obukhov方案(仅用于MYJ边界层方案)=3, NCEP Global Forecast

24、 System schemesf_surface_physics(max_dom)此选项指定陆面过程方案，默认值为0。旧参数表中的 bl_surface_physics，有效选择值有：=0,不采用陆面过程方案=1,热量扩散方案=2, Noah陆面过程方案=3, RUC陆面过程方案bl_pbl_physics (max_dom)此选项指定边界层方案，默认值为0。有效选择值有：=0,不采用边界层方案=1, YSU方案=2, Eta Mellor-Yamada-Janjic TKE(湍流动能)方案=3, NCEP Global Forecast System scheme=99, MRF方案(将放弃

25、)bldt (max_dom)此参数指定调用边界层物理方案的时间间隔，默认值为0，单位为分钟。0 (推荐值)表示每一个时间步长都调用边界层物理方案。cu_physics (max_dom)此选项指定积云参数化方案，默认值为0。有效选择值有：=0,不采用积云参数化方案=1, 浅对流 Kain-Fritsch (new Eta)方案=2, Betts-Miller-Janjic 方案=3, Grell-Devenyi 集合方案=4,简化 Arakawa-Schubert 方案=99,老 Kain-Fritsch 方案cudt(max_dom)此参数设定积云参数化方案的调用时间间隔，默认值为0,单位

26、为分钟。一般的积云参数化方案是每一步都要调用，但如果是用Kain-Fritsch方案 (cu_physics=1)，则可以设 cudt=5。用于NMM版本WRF。指定调用积云参数化过程的时间间隔，单位为粗网格的 时间步数。默认值在Registry.NMM进行设置，但也可在namelist.input 文件中重新赋值。模式会根据此数值计算cudt的值。ISFFLX -在使用扰动边界层时有效，即sf_sfclay_physics = 1此选项指定在选用扰动边界层和陆面物理过程时，是否考虑地面热量和水 汽通量，默认值为1：=考虑地面通量0 =不考虑地面通量IFSNOW此选项指定是否考虑雪盖效应。考虑

27、雪盖效应时，必须要有雪盖输入场。默认值为0，只有在利用扰动边界层PBL预报土壤温度是才有效，即 sf_surface_physics = 1。=考虑雪盖效应0 =不考虑雪盖效应ICLOUD此参数指定辐射光学厚度中是否考虑云的影响，默认值为1。仅当 ra_sw_physics = 1 和 ra_lw_physics = 1 时有效。=考虑云的影响0 =不考虑云的影响surface_input_source此参数知道土地利用类型和土壤类型数据的来源格式，默认值为1。= SI/gridgen （由 SI 的 gridgen_model.exe 程序产生）=其他模式产生的GRIB码数据（VEGCAT/

28、SOILCAT 数据都在由 SI 产生的 wrf_real_input_em 文件中）num_soil_layers指定陆面模式中的土壤层数，默认值为5=5:热量扩散方案 =4: Noah陆面过程方案=6: RUC陆面过程方案maxiens默认值为1,仅用于积云参数化方案中的Grell-Devenyi集合方案maxens默认值为3,仅用于积云参数化方案中的Grell-Devenyi集合方案maxens2默认值为3,仅用于积云参数化方案中的Grell-Devenyi集合方案maxens3默认值为16,仅用于积云参数化方案中的Grell-Devenyi集合方案ensdim默认值为144,仅用于积云

29、参数化方案中的Grell-Devenyi集合方案说明：以上这些用于Grell-Devenyi方案的默认值，都是一下推荐使用的数值。 如果要改变数值，请谨慎修改。seaice_threshold = 271海冰温度临界值。当TSK小于此临界值时，如果模式格点是水体，陆面过 程选用5层的SLAB方案，则将此模式格点设置为陆地，且为永久性冰体； 如果模式格点是水体，陆面过程选用Noah方案，则将此模式格点设置为陆 地，且为永久性冰体，并将设置03米的TEMPS,以及设置SMOIS和SH2O。sst_update = 0时变海温控制参数。0表示不用，1表示使用。如果选择使用时变海温，则 real.ex

30、e 会从 wrflowinp_d01 文件中读取 SST 和 VEGFRA 数据，wrf.exe 则 会以更新边条件数据相同的时间间隔来更新这些数据。要使用此功能，则 在参数列表文件namelist.input的时间控制区还必须包含 auxinput5_interval, auxinput5_end_h, 和 auxinput5_inname = wrflowinp_d。参数配置第四部分dynamicsdyn_opt模式框架配置选项，默认值为2：=欧拉高度坐标（已经放弃）=欧拉质量坐标=半拉格朗日（目前还没完成）rk_ord本参数定义Runge-Kutta时间积分方案阶数，默认值为3：= Ru

31、nge-Kutta 二阶= Runge-Kutta 三阶（推荐）diff_opt湍流和混合作用选项，默认值为0:0 =没有湍流或者显式空间数值滤波（km_opt将被忽略）（explicit spatial numerical filters）1=老扩散方案，计算坐标面上二阶扩散项。如果没有指定PBL选项，则用kvdif 选项当作垂直扩散系数。通常用于km_opt=1或者4。（在真实大 气方案的水平分辨率格距小于10km时，推荐使用1）2 =新扩散方案，计算物理空间（x,y,z）中的混合作用项（应力形式）。用km_opt 来指明湍流参数化过程。km_opt湍涡系数选项，默认值为1：1 =固定不变（用参数配置第三部分的khdif, kvdif参数值）与diff_opt=1 的区别在于km_opt的水平扩散作用不在模式的zeta面上。因此， 只有在没有地形的情况下，这两个选项的作用才是相同的。1.5阶TKE