| blob | 6a092b25b27529f10194bca1f8dc39dd570af278 |
1 Description of namelist variables
2 ---------------------------------
4 For WRF-NMM users, please see Chapter 5 of the WRF-NMM User's Guide for
5 information on NMM specific settings (http://www.dtcenter.org/wrf-nmm/users)
8 Note: variables followed by (max_dom) indicate that this variable needs to
9 be defined for the nests when max_dom > 1.
11 &time_control
12 run_days = 1, ; run time in days
13 run_hours = 0, ; run time in hours
14 Note: if it is more than 1 day, one may use both run_days and run_hours
15 or just run_hours. e.g. if the total run length is 36 hrs, you may
16 set run_days = 1, and run_hours = 12, or run_days = 0, and run_hours = 36
17 run_minutes = 0, ; run time in minutes
18 run_seconds = 0, ; run time in seconds
19 start_year (max_dom) = 2001, ; four digit year of starting time
20 start_month (max_dom) = 06, ; two digit month of starting time
21 start_day (max_dom) = 11, ; two digit day of starting time
22 start_hour (max_dom) = 12, ; two digit hour of starting time
23 start_minute (max_dom) = 00, ; two digit minute of starting time
24 start_second (max_dom) = 00, ; two digit second of starting time
25 Note: the start time is used to name the first wrfout file.
26 It also controls the start time for nest domains, and the time to restart
27 tstart (max_dom) = 00, ; FOR NMM: starting hour of the forecast
28 end_year (max_dom) = 2001, ; four digit year of ending time
29 end_month (max_dom) = 06, ; two digit month of ending time
30 end_day (max_dom) = 12, ; two digit day of ending time
31 end_hour (max_dom) = 12, ; two digit hour of ending time
32 end_minute (max_dom) = 00, ; two digit minute of ending time
33 end_second (max_dom) = 00, ; two digit second of ending time
34 It also controls when the nest domain integrations end
35 All start and end times are used by real.exe.
37 Note that one may use either run_days/run_hours etc. or
38 end_year/month/day/hour etc. to control the length of
39 model integration. But run_days/run_hours
40 takes precedence over the end times.
41 Program real.exe uses start and end times only.
43 interval_seconds = 10800, ; time interval between incoming real data, which will be the interval
44 between the lateral boundary condition file
45 input_from_file (max_dom) = T, ; whether nested run will have input files for domains other than 1
46 fine_input_stream (max_dom) = 0, ; field selection from nest input for its initialization
47 0: all fields are used; 2: only static and time-varying, masked land
48 surface fields are used. In V3.2, this requires the use of
49 io_form_auxinput2
50 history_interval (max_dom) = 60, ; history output file interval in minutes
51 frames_per_outfile (max_dom) = 1, ; number of output times per history output file,
52 used to split output into multiple files
53 into smaller pieces
54 restart = F, ; whether this run is a restart run
55 cycling = F, ; whether this run is a cycling run, if so, initializes look-up table for Thompson schemes only
56 restart_interval = 1440, ; restart output file interval in minutes
57 reset_simulation_start = F, ; whether to overwrite simulation_start_date with forecast start time
58 io_form_history = 2, ; 2 = netCDF
59 io_form_restart = 2, ; 2 = netCDF
60 io_form_input = 2, ; 2 = netCDF
61 io_form_boundary = 2, ; netCDF format
62 = 4, ; PHD5 format
63 = 5, ; GRIB1 format
64 = 10, ; GRIB2 format
65 = 11, ; pnetCDF format
66 frames_per_emissfile = 12, ; number of times in each chemistry emission file.
67 io_style_emiss = 1, ; style to use for the chemistry emission files.
68 ; 0 = Do not read emissions from files.
69 ; 1 = Cycle between two 12 hour files (set frames_per_emissfile=12)
70 ; 2 = Dated files with length set by frames_per_emissfile
71 debug_level = 0, ; 50,100,200,300 values give increasing prints
72 diag_print = 0, ; print out time series of model diagnostics
73 ; 0 = no print
74 ; 1 = no print
75 all_ic_times = .false., ; whether to write out wrfinput for all processing times
76 adjust_output_times = .false., ; adjust output times to the nearest hour
78 To choose between SI and WPS input to real for EM core:
79 auxinput1_inname = "met_em.d<domain>.<date>" ; Input to real from WPS (default since 3.0)
80 = "wrf_real_input_em.d<domain>.<date>" ; Input to real from SI
82 To choose between SI and WPS input to real for NMM core:
83 auxinput1_inname = "met_nm.d<domain>.<date>" ; Input to real from WPS
84 = "wrf_real_input_nm.d<domain>.<date>" ; Input to real from SI
86 Other output options:
88 auxhist2_outname = "rainfall" ; file name for extra output; if not specified,
89 auxhist2_d<domain>_<date> will be used
90 also note that to write variables in output other
91 than the history file requires Registry.EM file change
92 auxhist2_interval (max_dom) = 10, ; interval in minutes
93 io_form_auxhist2 = 2, ; output in netCDF
94 frames_per_auxhist2 = 1000, ; number of output times in this file
96 For SST updating (used only with sst_update=1):
98 auxinput4_inname = "wrflowinp_d<domain>"
99 auxinput4_interval = 360 ; minutes generally matches time given by interval_seconds
100 io_form_auxinput4 = 2 ; IO format, required in V3.2
102 For additional regional climate surface fields
104 output_diagnostics = 1 ; adds 36 surface diagnostic arrays (max/min/mean/std)
105 auxhist3_outname = 'wrfxtrm_d<domain>_<date>' ; file name for added diagnostics
106 io_form_auxhist3 = 2 ; netcdf
107 auxhist3_interval = 1440 ; minutes between outputs (1440 gives daily max/min)
108 frames_per_auxhist3 = 1 ; output times per file
110 For observation nudging:
111 auxinput11_interval = 10 ; interval in minutes for observation data. It should be
112 set as or more frequently as obs_ionf (with unit of
113 coarse domain time step).
114 auxinput11_end_h = 6 ; end of observation time in hours.
116 Options for run-time IO:
118 iofields_filename (max_dom) = "my_iofields_list.txt",
119 (example: +:h:21:rainc, rainnc, rthcuten)
120 ignore_iofields_warning = .true., ; what to do when encountering an error in the user-specified files
121 .false., : abort when encountering an error in iofields_filename file
123 Additional settings when running WRFVAR:
125 write_input = t, ; write input-formatted data as output
126 inputout_interval = 180, ; interval in minutes when writing input-formatted data
127 input_outname = 'wrfinput_d<domain>_<date>' ; you may change the output file name
128 inputout_begin_y = 0
129 inputout_begin_mo = 0
130 inputout_begin_d = 0
131 inputout_begin_h = 3
132 inputout_begin_m = 0
133 inputout_begin_s = 0
134 inputout_end_y = 0
135 inputout_end_mo = 0
136 inputout_end_d = 0
137 inputout_end_h = 12
138 inputout_end_m = 0
139 inputout_end_s = 0 ; the above shows that the input-formatted data are output
140 starting from hour 3 to hour 12 in 180 min interval.
142 &domains
143 time_step = 60, ; time step for integration in integer seconds
144 recommend 6*dx (in km) for typical real-data cases
145 time_step_fract_num = 0, ; numerator for fractional time step
146 time_step_fract_den = 1, ; denominator for fractional time step
147 Example, if you want to use 60.3 sec as your time step,
148 set time_step = 60, time_step_fract_num = 3, and
149 time_step_fract_den = 10
150 time_step_dfi = 60, ; time step for DFI, may be different from regular time_step
151 max_dom = 1, ; number of domains - set it to > 1 if it is a nested run
152 s_we (max_dom) = 1, ; start index in x (west-east) direction (leave as is)
153 e_we (max_dom) = 91, ; end index in x (west-east) direction (staggered dimension)
154 s_sn (max_dom) = 1, ; start index in y (south-north) direction (leave as is)
155 e_sn (max_dom) = 82, ; end index in y (south-north) direction (staggered dimension)
156 s_vert (max_dom) = 1, ; start index in z (vertical) direction (leave as is)
157 e_vert (max_dom) = 28, ; end index in z (vertical) direction (staggered dimension)
158 Note: this refers to full levels including surface and top
159 vertical dimensions need to be the same for all nests
160 Note: most variables are unstaggered (= staggered dim - 1)
161 dx (max_dom) = 10000, ; grid length in x direction; ARW: unit in meters, NMM: unit in degrees (e.g. 0.667)
162 dy (max_dom) = 10000, ; grid length in y direction; ARW: unit in meters, NMM: unit in degrees (e.g. 0.0658)
163 ztop (max_dom) = 19000. ; used in mass model for idealized cases
164 grid_id (max_dom) = 1, ; domain identifier
165 parent_id (max_dom) = 0, ; id of the parent domain
166 i_parent_start (max_dom) = 0, ; starting LLC I-indices from the parent domain
167 j_parent_start (max_dom) = 0, ; starting LLC J-indices from the parent domain
168 parent_grid_ratio (max_dom) = 1, ; parent-to-nest domain grid size ratio: for real-data cases
169 the ratio has to be odd; for idealized cases,
170 the ratio can be even if feedback is set to 0. (NMM: must be 3)
171 parent_time_step_ratio (max_dom) = 1, ; parent-to-nest time step ratio; it can be different
172 from the parent_grid_ratio (NMM: must be 3)
173 feedback = 1, ; feedback from nest to its parent domain; 0 = no feedback
174 smooth_option = 0 ; smoothing option for parent domain, used only with feedback
175 option on. 0: no smoothing; 1: 1-2-1 smoothing; 2: smoothing-desmoothing
177 Namelist variables specifically for the WPS input for real:
179 num_metgrid_soil_levels = 4 ; number of vertical soil levels or layers input
180 ; from WPS metgrid program
181 num_metgrid_levels = 27 ; number of vertical levels of 3d meteorological fields coming
182 ; from WPS metgrid program
183 interp_type = 2 ; vertical interpolation
184 ; 1 = linear in pressure
185 ; 2 = linear in log(pressure)
186 extrap_type = 2 ; vertical extrapolation of non-temperature fields
187 ; 1 = extrapolate using the two lowest levels
188 ; 2 = use lowest level as constant below ground
189 t_extrap_type = 2 ; vertical extrapolation for potential temperature
190 ; 1 = isothermal
191 ; 2 = -6.5 K/km lapse rate for temperature
192 ; 3 = constant theta
193 use_levels_below_ground = .true. ; in vertical interpolation, use levels below input surface level
194 ; T = use input isobaric levels below input surface
195 ; F = extrapolate when WRF location is below input surface value
196 use_surface = .true. ; use the input surface level data in the vertical interp and extrap
197 ; T = use the input surface data
198 ; F = do not use the input surface data
199 lagrange_order = 1 ; vertical interpolation order
200 ; 1 = linear
201 ; 2 = quadratic
202 zap_close_levels = 500 ; ignore isobaric level above surface if delta p (Pa) < zap_close_levels
203 lowest_lev_from_sfc = .false. ; place the surface value into the lowest eta location
204 ; T = use surface value as lowest eta (u,v,t,q)
205 ; F = use traditional interpolation
206 force_sfc_in_vinterp = 1 ; use the surface level as the lower boundary when interpolating
207 ; through this many eta levels
208 ; 0 = perform traditional trapping interpolation
209 ; n = first n eta levels directly use surface level
210 sfcp_to_sfcp = .false. ; Optional method to compute model's surface pressure when incoming
211 ; data only has surface pressure and terrain, but not SLP
212 smooth_cg_topo = .false. ; Smooth the outer rows and columns of domain 1's topography w.r.t.
213 ; the input data
214 use_tavg_for_tsk = .false. ; whether to use diurnally averaged surface temp as skin temp. The
215 diurnall averaged surface temp can be computed using WPS utility
216 avg_tsfc.exe. May use this option when SKINTEMP is not present.
217 aggregate_lu = .false. ; whetger to aggregate the grass, shrubs, trees in dominant landuse;
218 default is false.
219 rh2qv_wrt_liquid = .true., ; whether to compute RH with respect to water (true) or ice (false)
220 rh2qv_method = 1, ; which methed to use to computer mixing ratio from RH:
221 default is option 1, the old MM5 method; option 2 uses a WMO
222 recommended method (WMO-No. 49, corrigendum, August 2000) -
223 there is a difference between the two methods though small
224 p_top_requested = 5000 ; p_top (Pa) to use in the model
225 ptsgm = 42000. ; FOR NMM: defines the pressure interface dividing
226 ; the terrain following portion of the hybrid vertical
227 ; coordinate (p > ptsgm) and the purely
228 ; isobaric portion of the vertical coordinate (p < ptsgm)
229 vert_refine_fact = 1 ; vertical refinement factor for ndown
231 Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number
232 of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request
233 only the number of levels (with e_vert), and the real program will compute values. The computation assumes
234 a known first several layers, then generates equi-height spaced levels up to the top of the model.
236 eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946,
237 0.922, 0.894, 0.860, 0.817, 0.766,
238 0.707, 0.644, 0.576, 0.507, 0.444,
239 0.380, 0.324, 0.273, 0.228, 0.188,
240 0.152, 0.121, 0.093, 0.069, 0.048,
241 0.029, 0.014, 0.000,
242 eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
243 0.934, 0.909, 0.880, 0.845, 0.807,
244 0.765, 0.719, 0.672, 0.622, 0.571,
245 0.520, 0.468, 0.420, 0.376, 0.335,
246 0.298, 0.263, 0.231, 0.202, 0.175,
247 0.150, 0.127, 0.106, 0.088, 0.070,
248 0.055, 0.040, 0.026, 0.013, 0.000
250 Namelist variables for controling the specified moving nest:
251 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
252 to the ARCHFLAGS. The maximum number of moves, max_moves, is set to 50
253 but can be modified in source code file frame/module_driver_constants.F.
254 num_moves = 4 ; total number of moves
255 move_id(max_moves) = 2,2,2,2, ; a list of nest domain id's, one per move
256 move_interval(max_moves) = 60,120,150,180, ; time in minutes since the start of this domain
257 move_cd_x(max_moves) = 1,1,0,-1,; the number of parent domain grid cells to move in i direction
258 move_cd_y(max_moves) = 1,0,-1,1,; the number of parent domain grid cells to move in j direction
259 positive is to move in increasing i and j direction, and
260 negative is to move in decreasing i and j direction.
261 0 means no move. The limitation now is to move only 1 grid cell
262 at each move.
264 Namelist variables for controling the automatic moving nest:
265 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
266 and -DVORTEX_CENTER to the ARCHFLAGS. This option uses an mid-level vortex following algorthm to
267 determine the nest move. This option is experimental.
268 vortex_interval(max_dom) = 15 ; how often the new vortex position is computed
269 max_vortex_speed(max_dom) = 40 ; used to compute the search radius for the new vortex position
270 corral_dist(max_dom) = 8 ; how many coarse grid cells the moving nest is allowed to get
271 near the mother domain boundary
272 track_level = 50000 ; pressure value in Pa where the vortex is tracked
273 time_to_move(max_dom) = 0. ; time (in minutes) to start the moving nests
275 tile_sz_x = 0, ; number of points in tile x direction
276 tile_sz_y = 0, ; number of points in tile y direction
277 can be determined automatically
278 numtiles = 1, ; number of tiles per patch (alternative to above two items)
279 nproc_x = -1, ; number of processors in x for decomposition
280 nproc_y = -1, ; number of processors in y for decomposition
281 -1: code will do automatic decomposition
282 >1: for both: will be used for decomposition
284 Namelist variables for controlling the adaptive time step option:
285 These options are only valid for the ARW core.
286 use_adaptive_time_step = .false. ; T/F use adaptive time stepping, ARW only
287 step_to_output_time = .true. ; if adaptive time stepping, T/F modify the
288 time steps so that the exact history time is reached
289 target_cfl(max_dom) = 1.2,1.2 ; vertical and horizontal CFL <= to this value implies
290 no reason to reduce the time step, and to increase it
291 target_hcfl(max_dom) = .84,.84 ; horizontal CFL <= to this value implies
292 max_step_increase_pct(max_dom) = 5,51 ; percentage of previous time step to increase, if the
293 max(vert cfl, horiz cfl) <= target_cfl, then the time
294 will increase by max_step_increase_pct. Use something
295 large for nests (51% suggested)
296 starting_time_step(max_dom) = -1,-1 ; flag = -1 implies use 6 * dx (defined in start_em),
297 starting_time_step = 100 means the starting time step
298 for the coarse grid is 100 s
299 max_time_step(max_dom) = -1,-1 ; flag = -1 implies max time step is 3 * starting_time_step,
300 max_time_step = 100 means that the time step will not
301 exceed 100 s
302 min_time_step(max_dom) = -1,-1 ; flag = -1 implies max time step is 0.5 * starting_time_step,
303 min_time_step = 100 means that the time step will not
304 be less than 100 s
305 adaptation_domain = 1 ; default, all fine grid domains adaptive dt driven by coarse-grid
306 ; 2 = Fine grid domain #2 determines the fundamental adaptive dt.
308 &dfi_control
309 dfi_opt = 0 ; which DFI option to use (3 is recommended)
310 ; 0 = no digital filter initialization
311 ; 1 = digital filter launch (DFL)
312 ; 2 = diabatic DFI (DDFI)
313 ; 3 = twice DFI (TDFI)
314 dfi_nfilter = 7 ; digital filter type to use (7 is recommended)
315 ; 0 = uniform
316 ; 1 = Lanczos
317 ; 2 = Hamming
318 ; 3 = Blackman
319 ; 4 = Kaiser
320 ; 5 = Potter
321 ; 6 = Dolph window
322 ; 7 = Dolph
323 ; 8 = recursive high-order
324 dfi_write_filtered_input = .true. ; whether to write wrfinput file with filtered
325 ; model state before beginning forecast
326 dfi_write_dfi_history = .false. ; whether to write wrfout files during filtering integration
327 dfi_cutoff_seconds = 3600 ; cutoff period, in seconds, for the filter
328 dfi_time_dim = 1000 ; maximum number of time steps for filtering period
329 ; this value can be larger than necessary
330 dfi_bckstop_year = 2004 ; four-digit year of stop time for backward DFI integration
331 dfi_bckstop_month = 03 ; two-digit month of stop time for backward DFI integration
332 dfi_bckstop_day = 14 ; two-digit day of stop time for backward DFI integration
333 dfi_bckstop_hour = 12 ; two-digit hour of stop time for backward DFI integration
334 dfi_bckstop_minute = 00 ; two-digit minute of stop time for backward DFI integration
335 dfi_bckstop_second = 00 ; two-digit second of stop time for backward DFI integration
336 dfi_fwdstop_year = 2004 ; four-digit year of stop time for forward DFI integration
337 dfi_fwdstop_month = 03 ; two-digit month of stop time for forward DFI integration
338 dfi_fwdstop_day = 13 ; two-digit month of stop time for forward DFI integration
339 dfi_fwdstop_hour = 12 ; two-digit month of stop time for forward DFI integration
340 dfi_fwdstop_minute = 00 ; two-digit month of stop time for forward DFI integration
341 dfi_fwdstop_second = 00 ; two-digit month of stop time for forward DFI integration
342 dfi_radar = 0 ; DFI radar da switch
344 &physics
346 Note: even the physics options can be different in different nest domains,
347 caution must be used as what options are sensible to use
349 chem_opt = 0, ; chemistry option - use WRF-Chem
350 mp_physics (max_dom) microphysics option
351 = 0, no microphysics
352 = 1, Kessler scheme
353 = 2, Lin et al. scheme
354 = 3, WSM 3-class simple ice scheme
355 = 4, WSM 5-class scheme
356 = 5, Ferrier (new Eta) microphysics
357 = 6, WSM 6-class graupel scheme
358 = 7, Goddard GCE scheme (also uses gsfcgce_hail, gsfcgce_2ice)
359 = 8, Thompson scheme (new for V3.1)
360 = 9, Milbrandt-Yau 2-moment scheme (new for V3.2)
361 = 10, Morrison (2 moments)
362 = 13, SBU_YLIN scheme
363 = 14, WDM 5-class scheme
364 = 16, WDM 6-class scheme
366 For non-zero mp_physics options, to keep Qv .GE. 0, and to set the other moisture
367 fields .LT. a critcal value to zero
369 mp_zero_out = 0, ; no action taken, no adjustment to any moist field
370 = 1, ; except for Qv, all other moist arrays are set to zero
371 ; if they fall below a critical value
372 = 2, ; Qv is .GE. 0, all other moist arrays are set to zero
373 ; if they fall below a critical value
374 mp_zero_out_thresh = 1.e-8 ; critical value for moist array threshold, below which
375 ; moist arrays (except for Qv) are set to zero (kg/kg)
377 gsfcgce_hail = 0 ; for running gsfcgce microphysics with graupel
378 = 1 ; for running gsfcgce microphysics with hail
379 default value = 0
380 gsfcgce_2ice = 0 ; for running with snow, ice and graupel/hail
381 = 1 ; for running with only ice and snow
382 = 2 ; for running with only ice and graupel
383 (only used in very extreme situation)
384 default value = 0
385 gsfcgce_hail is ignored if gsfcgce_2ice is set to 1 or 2.
387 no_mp_heating = 0 ; normal
388 = 1 ; turn off latent heating from a microphysics scheme
390 ra_lw_physics (max_dom) longwave radiation option
391 = 0, no longwave radiation
392 = 1, rrtm scheme
393 = 3, cam scheme
394 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
395 = 4, rrtmg scheme
396 = 5, Goddard longwave scheme
397 = 31, Earth Held-Suarez forcing
398 = 99, GFDL (Eta) longwave (semi-supported)
399 also must use co2tf = 1 for ARW
401 ra_sw_physics (max_dom) shortwave radiation option
402 = 0, no shortwave radiation
403 = 1, Dudhia scheme
404 = 2, Goddard short wave
405 = 3, cam scheme
406 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
407 = 5, Goddard shortwave scheme
408 = 4, rrtmg scheme
409 = 99, GFDL (Eta) longwave (semi-supported)
410 also must use co2tf = 1 for ARW
412 radt (max_dom) = 30, ; minutes between radiation physics calls
413 recommend 1 min per km of dx (e.g. 10 for 10 km)
415 nrads (max_dom) = FOR NMM: number of fundamental timesteps between
416 calls to shortwave radiation; the value
417 is set in Registry.NMM but is overridden
418 by namelist value; radt will be computed
419 from this.
421 nradl (max_dom) = FOR NMM: number of fundamental timesteps between
422 calls to longwave radiation; the value
423 is set in Registry.NMM but is overridden
424 by namelist value.
426 co2tf CO2 transmission function flag only for GFDL radiation
427 = 0, read CO2 function data from pre-generated file
428 = 1, generate CO2 functions internally in the forecast
430 ra_call_offset radiation call offset
431 = 0 (no offset), =-1 (old offset)
433 cam_abs_freq_s = 21600 CAM clearsky longwave absorption calculation frequency
434 (recommended minimum value to speed scheme up)
435 levsiz = 59 for CAM radiation input ozone levels
436 paerlev = 29 for CAM radiation input aerosol levels
437 cam_abs_dim1 = 4 for CAM absorption save array
438 cam_abs_dim2 = value of e_vert for CAM 2nd absorption save array
440 sf_sfclay_physics (max_dom) surface-layer option (old bl_sfclay_physics option)
441 = 0, no surface-layer
442 = 1, Monin-Obukhov scheme
443 = 2, Monin-Obukhov (Janjic) scheme
444 = 3, NCEP Global Forecast System scheme (NMM only)
445 = 4, QNSE surface layer
446 = 5, MYNN surface layer
447 = 7, Pleim-Xiu surface layer (ARW only)
448 = 10, TEMF surface layer (ARW only)
450 sf_surface_physics (max_dom) land-surface option (old bl_surface_physics option)
451 = 0, no surface temp prediction
452 = 1, thermal diffusion scheme
453 = 2, Unified Noah land-surface model
454 = 3, RUC land-surface model
455 = 7, Pleim-Xiu LSM (ARW)
457 sf_urban_physics(max_dom) = 0, ; activate urban canopy model (in Noah LSM only)
458 = 0: no
459 = 1: Single-layer, UCM
460 = 2: Multi-layer, Building Environment Parameterization (BEP) scheme
461 (works only with MYJ and BouLac PBL)
462 = 3: Multi-layer, Building Environment Model (BEM) scheme
463 (works only with MYJ and BouLac PBL)
465 bl_pbl_physics (max_dom) boundary-layer option
466 = 0, no boundary-layer
467 = 1, YSU scheme
468 = 2, Mellor-Yamada-Janjic TKE scheme
469 = 3, NCEP Global Forecast System scheme (NMM only)
470 = 4, Quasi-Normal Scale Elimination PBL
471 = 5, MYNN 2.5 level TKE scheme, works with
472 sf_sfclay_physics=1 or 2 as well as 5
473 = 6, MYNN 3rd level TKE scheme, works only
474 MYNNSFC (sf_sfclay_physics = 5)
475 = 7, ACM2 (Pleim) PBL (ARW)
476 = 8, Bougeault and Lacarrere (BouLac) PBL
477 = 9, UW boundary layer scheme from CAM5 (CESM 1_0_1)
478 = 10, TEMF (Total Energy Mass Flux) scheme (ARW only)
479 = 99, MRF scheme
481 bldt (max_dom) = 0, ; minutes between boundary-layer physics calls
483 grav_settling = 0, ; MYNN PBL only; gravitational settling of fog/cloud droplets (1=yes)
484 nphs (max_dom) = FOR NMM: number of fundamental timesteps between
485 calls to turbulence and microphysics;
486 the value is set in Registry.NMM but is
487 overridden by namelist value; bldt will
488 be computed from this.
490 cu_physics (max_dom) cumulus option
491 = 0, no cumulus
492 = 1, Kain-Fritsch (new Eta) scheme
493 = 2, Betts-Miller-Janjic scheme
494 = 3, Grell-Devenyi ensemble scheme
495 = 4, Simplified Arakawa-Schubert scheme
496 = 5, Grell 3D ensemble scheme
497 = 6, Modifed Tiedtke scheme (ARW only)
498 = 7, Zhang-McFarlane scheme from CAM5 (CESM 1_0_1)
499 = 14, New GFS simplified Arakawa-Schubert scheme from YSU (ARW only)
500 = 99, previous Kain-Fritsch scheme
502 shcu_physics (max_dom) independent shallow cumulus option (not tied to deep convection)
503 = 0, no independent shallow cumulus
504 = 1, Grell 3D ensemble scheme (use with cu_physics=5) (PLACEHOLDER: SWITCH NOT YET IMPLEMENTED--use ishallow)
505 = 2, Park and Bretherton shallow cumulus from CAM5 (CESM 1_0_1)
507 ishallow = 1, Shallow convection used with Grell 3D ensemble scheme (cu_physics = 5)
508 clos_choice = 0, closure choice (place holder only)
510 cu_diag = 0, additional t-averaged stuff for cu physics (GD and G3 only)
511 convtrans_avglen_m = 30, averaging time for convective transport output variables (minutes) (GD and G3 only)
513 cudt = 0, ; minutes between cumulus physics calls
515 kfeta_trigger KF trigger option (cu_physics=1 only):
516 = 1, default option
517 = 2, moisture-advection based trigger (Ma and Tan [2009]) - ARW only
518 = 3, RH-dependent additional perturbation to option 1 (JMA)
520 cugd_avedx ; number of grid boxes over which subsidence is spread.
521 = 1, default, for large grid distances
522 = 3, for small grid distances (DX < 5 km)
524 ncnvc (max_dom) = FOR NMM: number of fundamental timesteps between
525 calls to convection; the value is set in Registry.NMM
526 but is overridden by namelist value; cudt will be
527 computed from this.
529 tprec (max_dom) = FOR NMM: number of hours in precipitation bucket
530 theat (max_dom) = FOR NMM: number of hours in latent heating bucket
531 tclod (max_dom) = FOR NMM: number of hours in cloud fraction average
532 trdsw (max_dom) = FOR NMM: number of hours in short wave buckets
533 trdlw (max_dom) = FOR NMM: number of hours in long wave buckets
534 tsrfc (max_dom) = FOR NMM: number of hours in surface flux buckets
535 pcpflg (max_dom) = FOR NMM: logical switch for precipitation assimilation
537 isfflx = 1, ; heat and moisture fluxes from the surface
538 (only works for sf_sfclay_physics = 1)
539 1 = with fluxes from the surface
540 0 = no flux from the surface
541 with bl_pbl_physics=0 this uses tke_drag_coefficient
542 and tke_heat_flux in vertical diffusion
543 2 = use drag from sf_sfclay_physics and heat flux from
544 tke_heat_flux with bl_pbl_physics=0
545 ifsnow = 0, ; snow-cover effects
546 (only works for sf_surface_physics = 1)
547 1 = with snow-cover effect
548 0 = without snow-cover effect
549 icloud = 1, ; cloud effect to the optical depth in radiation
550 (only works for ra_sw_physics = 1 and ra_lw_physics = 1)
551 1 = with cloud effect
552 0 = without cloud effect
553 swrad_scat = 1. ; scattering tuning parameter (default 1. is 1.e-5 m2/kg)
554 surface_input_source = 1, ; where landuse and soil category data come from:
555 1 = WPS/geogrid but with dominant categories recomputed
556 2 = GRIB data from another model (only possible
557 (VEGCAT/SOILCAT are in met_em files from WPS)
558 3 = use dominant land and soil categories from WPS/geogrid
560 num_soil_layers = 5, ; number of soil layers in land surface model
561 = 5: thermal diffusion scheme
562 = 4: Noah landsurface model
563 = 6: RUC landsurface model
564 = 2: Pleim-Xu landsurface model
565 num_land_cat = 24, ; number of land categories in input data.
566 24 - for USGS (default); 20 for MODIS
567 28 - for USGS if including lake category
568 21 - for MODIS if including lake category
569 num_soil_cat = 16, ; number of soil categories in input data
571 pxlsm_smois_init(max_dom) = 1 ; PXLSM Soil moisture initialization option
572 0 - From analysis, 1 - From MAVAIL
574 maxiens = 1, ; Grell-Devenyi only
575 maxens = 3, ; G-D only
576 maxens2 = 3, ; G-D only
577 maxens3 = 16 ; G-D only
578 ensdim = 144 ; G-D only
579 These are recommended numbers. If you would like to use
580 any other number, consult the code, know what you are doing.
581 seaice_threshold = 271 ; tsk < seaice_threshold, if water point and 5-layer slab
582 ; scheme, set to land point and permanent ice; if water point
583 ; and Noah scheme, set to land point, permanent ice, set temps
584 ; from 3 m to surface, and set smois and sh2o
585 sst_update = 0 ; time-varying sea-surface temp (0=no, 1=yes). If selected real
586 ; puts SST, XICE, ALBEDO and VEGFRA in wrflowinp_d01 file, and wrf updates
587 ; these from it at same interval as boundary file. Also requires
588 ; namelists in &time_control: auxinput4_interval, auxinput4_end_h,
589 ; auxinput4_inname = "wrflowinp_d<domain>",
590 ; and in V3.2 io_form_auxinput4
591 usemonalb = .true. ; use monthly albedo map instead of table value
592 ; (must be used for NMM and recommended for sst_update=1)
593 rdmaxalb = .true. ; use snow albedo from geogrid; false means using values from table
594 rdlai2d = .false. ; use LAI from input; false means using values from table
595 bucket_mm = -1. ; bucket reset value for water accumulations (value in mm, -1.=inactive)
596 bucket_J = -1. ; bucket reset value for energy accumulations (value in J, -1.=inactive)
597 tmn_update = 0 ; update deep soil temperature (1, yes; 0, no)
598 lagday = 150 ; days over which tmn is computed using skin temperature
599 sst_skin = 0 ; calculate skin SST
600 slope_rad (max_dom) = 0 ; slope effects for solar radiation (1=on, 0=off)
601 topo_shading (max_dom) = 0 ; neighboring-point shadow effects for solar radiation (1=on, 0=off)
602 shadlen = 25000. ; max shadow length in meters for topo_shading=1
603 omlcall = 0 ; activate simple ocean mixed layer model (0=no, 1=yes); works with
604 sf_surface_physics = 1 only
605 oml_hml0 = 50 ; oml model can be initialized with a constant depth everywhere (m)
606 oml_gamma = 0.14 ; oml deep water lapse rate (K m-1)
607 isftcflx = 0 ; alternative Ck, Cd formulation for tropical storm application (0=default, 1=new, 2=Garratt)
608 fractional_seaice = 0 ; treat sea-ice as fractional field (1) or ice/no-ice flag (0)
609 tice2tsk_if2cold = .false. ; set Tice to Tsk to avoid unrealistically low sea ice temperatures
610 iz0tlnd = 0 ; thermal roughness length for sfclay and myjsfc (0 - old, 1 - veg dependent Czil)
611 mp_tend_lim = 10., ; limit on temp tendency from mp latent heating from radar data assimilation
612 prec_acc_dt (max_dom) = 0., ; number of minutes in precipitation bucket (ARW only) - will add three
613 new 2d output fields: prec_acc_c, prec_acc_nc and snow_acc_nc
615 Options for wind turbine drag parameterization:
617 td_turbgridid = -1 ; which grid id has turbines in it
618 td_hubheight = 100. ; hub height (m)
619 td_diameter = 60. ; turbine diameter (m)
620 td_stdthrcoef = .158 ; standing thrust coefficient
621 td_cutinspeed = 4. ; cut-in speed (m/s)
622 td_cutoutspeed = 27. ; cut-out speed (m/s)
623 td_power = 2. ; turbine power (MW)
624 td_turbpercell = 1. ; number of turbines per cell
625 td_ewfx = 0 ; extent of wind farm in x-cells
626 td_ewfy = 0 ; extent of wind farm in y-cells
627 td_pwfx = 1 ; southwest corner of wind farm in x-cells
628 td_pwfy = 1 ; southwest corner of wind farm in y-cells
630 Options for stochastic kinetic-energy backscatter scheme:
632 stoch_force_opt (max_dom) = 0, : No stochastic parameterization
633 1, : Stochastic kinetic-energy backscatter scheme (SKEB)
634 stoch_vertstruc_opt (max_dom) = 0, : Constant vertical structure of random pattern generator
635 1, : Random phase vertical structure random pattern generator
636 tot_backscat_psi = 115200, ; Controls amplitude of rotational wind perturbations
637 tot_backscat_t = 2.E-6 ; Controls amplitude of potential temperature perturbations
638 nens = 1 ; an integer that controls the random number stream which will then
639 change the run. When running an ensemble, this can be
640 ensemble member number, so that each ensemble member gets a
641 different random number stream, hence a different perturbed run.
643 &fdda
644 grid_fdda (max_dom) = 1 ; grid-nudging fdda on (=0 off) for each domain
645 = 2 ; spectral nudging
646 gfdda_inname = "wrffdda_d<domain>" ; defined name in real
647 gfdda_interval_m (max_dom) = 360 ; time interval (in min) between analysis times (must use minutes)
648 gfdda_end_h (max_dom) = 6 ; time (in hours) to stop nudging after start of forecast
649 io_form_gfdda = 2 ; analysis data io format (2 = netCDF)
650 fgdt (max_dom) = 0 ; calculation frequency (minutes) for grid-nudging (0=every step)
651 if_no_pbl_nudging_uv (max_dom) = 0 ; 1= no nudging of u and v in the pbl, 0=nudging in the pbl
652 if_no_pbl_nudging_t (max_dom) = 0 ; 1= no nudging of temp in the pbl, 0=nudging in the pbl
653 if_no_pbl_nudging_q (max_dom) = 0 ; 1= no nudging of qvapor in the pbl, 0=nudging in the pbl
654 if_zfac_uv (max_dom) = 0 ; 0= nudge u and v in all layers, 1= limit nudging to levels above k_zfac_uv
655 k_zfac_uv (max_dom) = 10 ; 10=model level below which nudging is switched off for u and v
656 if_zfac_t (max_dom) = 0 ; 0= nudge temp in all layers, 1= limit nudging to levels above k_zfac_t
657 k_zfac_t (max_dom) = 10 ; 10=model level below which nudging is switched off for temp
658 if_zfac_q (max_dom) = 0 ; 0= nudge qvapor in all layers, 1= limit nudging to levels above k_zfac_q
659 k_zfac_q (max_dom) = 10 ; 10=model level below which nudging is switched off for qvapor
660 guv (max_dom) = 0.0003 ; nudging coefficient for u and v (sec-1)
661 gt (max_dom) = 0.0003 ; nudging coefficient for temp (sec-1)
662 gq (max_dom) = 0.0003 ; nudging coefficient for qvapor (sec-1)
663 if_ramping = 0 ; 0= nudging ends as a step function, 1= ramping nudging down at end of period
664 dtramp_min = 60.0 ; time (min) for ramping function, 60.0=ramping starts at last analysis time,
665 -60.0=ramping ends at last analysis time
666 grid_sfdda (max_dom) = 0 ; surface fdda switch (1, on; 0, off)
667 sgfdda_inname = "wrfsfdda_d<domain>" ; defined name for sfc nudgingi in input file (from program obsgrid)
668 sgfdda_end_h (max_dom) = 6 ; time (in hours) to stop sfc nudging after start of forecast
669 sgfdda_interval_m (max_dom) = 180 ; time interval (in min) between sfc analysis times (must use minutes)
670 io_form_sgfdda = 2 ; sfc analysis data io format (2 = netCDF)
671 guv_sfc (max_dom) = 0.0003 ; nudging coefficient for sfc u and v (sec-1)
672 gt_sfc (max_dom) = 0.0003 ; nudging coefficient for sfc temp (sec-1)
673 gq_sfc (max_dom) = 0.0003 ; nudging coefficient for sfc qvapor (sec-1)
674 rinblw = 250.0 ; radius of influence used to determine the confidence (or weights) for
675 the analysis, which is based on the distance between the grid point to the nearest
676 obs. The analysis without nearby observation is used at a reduced weight.
678 pxlsm_soil_nudge(max_dom) = 1 ; PXLSM Soil nudging option (requires wrfsfdda file)
680 The following are for spectral nudging:
681 fgdtzero (max_dom) = 0, ; 1= nudging tendencies are set to zero in between fdda calls
682 if_no_pbl_nudging_ph = 0, ; 1= no nudging of ph in the pbl, 0= nuding in the pbl
683 if_zfac_ph (max_dom) = 0, ; 0= nudge ph in all layers, 1= limit nudging to levels above k_zfac_ph
684 k_zfac_ph (max_dom) = 10, ; 10= model level below which nudging is switched off for ph
685 dk_zfac_uv (max_dom) = 1, ; depth in k between k_zfac_X to dk_zfac_X where nudging increases
686 linearly to full strength
687 dk_zfac_t (max_dom) = 1,
688 dk_zfac_ph (max_dom) = 1,
689 gph (max_dom) = 0.0003,
690 xwavenum (max_dom) = 3, ; top wave number to nudge in x direction
691 ywavenum (max_dom) = 3, ; top wave number to nudge in y direction
693 The following are for observation nudging:
694 obs_nudge_opt (max_dom) = 1 ; obs-nudging fdda on (=0 off) for each domain
695 also need to set auxinput11_interval and auxinput11_end_h
696 in time_control namelist
697 max_obs = 150000 ; max number of observations used on a domain during any
698 given time window
699 fdda_start = 0 ; obs nudging start time in minutes
700 fdda_end = 180 ; obs nudging end time in minutes
701 obs_nudge_wind (max_dom) = 1 ; whether to nudge wind: (=0 off)
702 obs_coef_wind = 6.E-4, ; nudging coefficient for wind, unit: s-1
703 obs_nudge_temp = 1 ; whether to nudge temperature: (=0 off)
704 obs_coef_temp = 6.E-4, ; nudging coefficient for temperature, unit: s-1
705 obs_nudge_mois = 1 ; whether to nudge water vapor mixing ratio: (=0 off)
706 obs_coef_mois = 6.E-4, ; nudging coefficient for water vapor mixing ratio, unit: s-1
707 obs_nudge_pstr = 0 ; whether to nudge surface pressure (not used)
708 obs_coef_pstr = 0. ; nudging coefficient for surface pressure, unit: s-1 (not used)
709 obs_rinxy = 200., ; horizonal radius of influence in km
710 obs_rinsig = 0.1, ; vertical radius of influence in eta
711 obs_twindo (max_dom) = 0.66667 ; half-period time window over which an observation
712 will be used for nudging (hours)
713 obs_npfi = 10, ; freq in coarse grid timesteps for diag prints
714 obs_ionf (max_dom) = 2 ; freq in coarse grid timesteps for obs input and err calc
715 obs_idynin = 0 ; for dynamic initialization using a ramp-down function to gradually
716 turn off the FDDA before the pure forecast (=1 on)
717 obs_dtramp = 40 ; time period in minutes over which the nudging is ramped down
718 from one to zero.
719 obs_prt_freq (max_dom) = 10, ; Frequency in obs index for diagnostic printout
720 obs_prt_max = 1000, ; Maximum allowed obs entries in diagnostic printout
721 obs_ipf_in4dob = .true. ; print obs input diagnostics (=.false. off)
722 obs_ipf_errob = .true. ; print obs error diagnostics (=.false. off)
723 obs_ipf_nudob = .true. ; print obs nudge diagnostics (=.false. off)
724 obs_ipf_init = .true. ; Enable obs init warning messages
726 obs_no_pbl_nudge_uv (max_dom) = 0 ; 1=no wind-nudging within pbl
727 obs_no_pbl_nudge_t (max_dom) = 0 ; 1=no temperature-nudging within pbl
728 obs_no_pbl_nudge_q (max_dom) = 0 ; 1=no moisture-nudging within pbl
729 obs_sfc_scheme_horiz = 0 ; horizontal spreading scheme for surf obs;
730 0=wrf scheme, 1=original mm5 scheme
731 obs_sfc_scheme_vert = 0 ; vertical spreading scheme for surf obs
732 0=regime vif scheme, 1=original simple scheme
733 obs_max_sndng_gap = 20 ; Max pressure gap between soundings, in cb
734 obs_nudgezfullr1_uv = 50 ; Vert infl full weight height for lowest model level (LML) obs, regime 1, winds
735 obs_nudgezrampr1_uv = 50 ; Vert infl ramp-to-zero height for LML obs, regime 1, winds
736 obs_nudgezfullr2_uv = 50 ; Vert infl full weight height for LML obs, regime 2, winds
737 obs_nudgezrampr2_uv = 50 ; Vert infl ramp-to-zero height for LML obs, regime 2, winds
738 obs_nudgezfullr4_uv = -5000 ; Vert infl full weight height for LML obs, regime 4, winds
739 obs_nudgezrampr4_uv = 50 ; Vert infl ramp-to-zero height for LML obs, regime 4, winds
740 obs_nudgezfullr1_t = 50 ; Vert infl full weight height for LML obs, regime 1, temperature
741 obs_nudgezrampr1_t = 50 ; Vert infl ramp-to-zero height for LML obs, regime 1, temperature
742 obs_nudgezfullr2_t = 50 ; Vert infl full weight height for LML obs, regime 2, temperature
743 obs_nudgezrampr2_t = 50 ; Vert infl ramp-to-zero height for LML obs, regime 2, temperature
744 obs_nudgezfullr4_t = -5000 ; Vert infl full weight height for LML obs, regime 4, temperature
745 obs_nudgezrampr4_t = 50 ; Vert infl ramp-to-zero height for LML obs, regime 4, temperature
746 obs_nudgezfullr1_q = 50 ; Vert infl full weight height for LML obs, regime 1, moisture
747 obs_nudgezrampr1_q = 50 ; Vert infl ramp-to-zero height for LML obs, regime 1, moisture
748 obs_nudgezfullr2_q = 50 ; Vert infl full weight height for LML obs, regime 2, moisture
749 obs_nudgezrampr2_q = 50 ; Vert infl ramp-to-zero height for LML obs, regime 2, moisture
750 obs_nudgezfullr4_q = -5000 ; Vert infl full weight height for LML obs, regime 4, moisture
751 obs_nudgezrampr4_q = 50 ; Vert infl ramp-to-zero height for LML obs, regime 4, moisture
752 obs_nudgezfullmin = 50 ; Min depth through which vertical infl fcn remains 1.0
753 obs_nudgezrampmin = 50 ; Min depth (m) through which vert infl fcn decreases from 1 to 0
754 obs_nudgezmax = 3000 ; Max depth (m) in which vert infl function is nonzero
755 obs_sfcfact = 1.0 ; Scale factor applied to time window for surface obs
756 obs_sfcfacr = 1.0 ; Scale factor applied to horiz radius of influence for surface obs
757 obs_dpsmx = 7.5 ; Max pressure change (cb) allowed within horiz radius of influence
758 /
760 &scm
761 scm_force = 1, ; switch for single column forcing (=0 off)
762 scm_force_dx = 4000. ; DX for SCM forcing (in meters)
763 num_force_layers = 8 ; number of SCM input forcing layers
764 scm_lu_index = 2 ; SCM landuse category (2 is dryland, cropland and pasture)
765 scm_isltyp = 4 ; SCM soil category (4 is silt loam)
766 scm_vegfra = 0.5 ; SCM vegetation fraction
767 scm_canwat = 0.0 ; SCM canopy water
768 scm_lat = 37.600 ; SCM latitude
769 scm_lon = -96.700 ; SCM longitude
770 scm_th_adv = .true. ; turn on theta advection in SCM
771 scm_wind_adv = .true. ; turn on wind advection in SCM
772 scm_qv_adv = .true. ; turn on moisture advection in SCM
773 scm_ql_adv = .true. ; turn on cloud liquid water advection in SCM
774 scm_vert_adv = .true. ; turn on vertical advection in SCM
775 num_force_soil_layers = 5, ; Number of SCM soil forcing layer
776 scm_soilT_force = .false. ; Turn on soil temp forcing in SCM
777 scm_soilq_force = .false. ; Turn on soil moisture forcing in SCM
778 scm_force_th_largescale = .false. ; Turn on large scale theta forcing in SCM
779 scm_force_qv_largescale = .false. ; Turn on large scale qv forcing in SCM
780 scm_force_ql_largescale = .false. ; Turn on large scale cloud water forcing in SCM
781 scm_force_wind_largescale = .false. ; Turn on large scale wind forcing in SCM
783 &dynamics
784 rk_ord = 3, ; time-integration scheme option:
785 2 = Runge-Kutta 2nd order
786 3 = Runge-Kutta 3rd order
787 diff_opt = 0, ; turbulence and mixing option:
788 0 = no turbulence or explicit
789 spatial numerical filters (km_opt IS IGNORED).
790 1 = evaluates 2nd order
791 diffusion term on coordinate surfaces.
792 uses kvdif for vertical diff unless PBL option
793 is used. may be used with km_opt = 1 and 4.
794 (= 1, recommended for real-data cases)
795 2 = evaluates mixing terms in
796 physical space (stress form) (x,y,z).
797 turbulence parameterization is chosen
798 by specifying km_opt.
799 km_opt = 1, ; eddy coefficient option
800 1 = constant (use khdif kvdif)
801 2 = 1.5 order TKE closure (3D)
802 3 = Smagorinsky first order closure (3D)
803 Note: option 2 and 3 are not recommended for DX > 2 km
804 4 = horizontal Smagorinsky first order closure
805 (recommended for real-data cases)
806 damp_opt = 0, ; upper level damping flag
807 0 = without damping
808 1 = with diffusive damping, maybe used for real-data cases
809 (dampcoef nondimensional ~0.01-0.1)
810 2 = with Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .003; idealized case only
811 not for real-data cases)
812 3 = with w-Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .05;
813 for real-data cases)
814 diff_6th_opt = 0, ; 6th-order numerical diffusion
815 0 = no 6th-order diffusion (default)
816 1 = 6th-order numerical diffusion (not recommended)
817 2 = 6th-order numerical diffusion but prohibit up-gradient diffusion
818 diff_6th_factor = 0.12, ; 6th-order numerical diffusion non-dimensional rate (max value 1.0
819 corresponds to complete removal of 2dx wave in one timestep)
820 dampcoef (max_dom) = 0., ; damping coefficient (see above)
821 zdamp (max_dom) = 5000., ; damping depth (m) from model top
822 w_damping = 0, ; vertical velocity damping flag (for operational use)
823 0 = without damping
824 1 = with damping
825 base_temp = 290., ; real-data, em ONLY, base sea-level temp (K)
826 base_pres = 10^5 ; real-data, em ONLY, base sea-level pres (Pa), DO NOT CHANGE
827 base_lapse = 50., ; real-data, em ONLY, lapse rate (K), DO NOT CHANGE
828 iso_temp = 0., ; real-data, em ONLY, reference temp in stratosphere
829 use_baseparam_fr_nml = .f., ; whether to use base state parameters from the namelist
830 khdif (max_dom) = 0, ; horizontal diffusion constant (m^2/s)
831 kvdif (max_dom) = 0, ; vertical diffusion constant (m^2/s)
832 smdiv (max_dom) = 0.1, ; divergence damping (0.1 is typical)
833 emdiv (max_dom) = 0.01, ; external-mode filter coef for mass coordinate model
834 (0.01 is typical for real-data cases)
835 epssm (max_dom) = .1, ; time off-centering for vertical sound waves
836 non_hydrostatic (max_dom) = .true., ; whether running the model in hydrostatic or non-hydro mode
837 pert_coriolis (max_dom) = .false., ; Coriolis only acts on wind perturbation (idealized)
838 top_lid (max_dom) = .false., ; Zero vertical motion at top of domain
839 mix_full_fields(max_dom) = .true., ; used with diff_opt = 2; value of ".true." is recommended, except for
840 highly idealized numerical tests; damp_opt must not be 1 if ".true."
841 is chosen. .false. means subtract 1-d base-state profile before mixing
842 mix_isotropic(max_dom) = 0 ; 0=anistropic vertical/horizontal diffusion coeffs, 1=isotropic
843 mix_upper_bound(max_dom) = 0.1 ; non-dimensional upper limit for diffusion coeffs
844 tke_drag_coefficient(max_dom) = 0., ; surface drag coefficient (Cd, dimensionless) for diff_opt=2 only
845 tke_heat_flux(max_dom) = 0., ; surface thermal flux (H/(rho*cp), K m/s) for diff_opt=2 only
846 h_mom_adv_order (max_dom) = 5, ; horizontal momentum advection order (5=5th, etc.)
847 v_mom_adv_order (max_dom) = 3, ; vertical momentum advection order
848 h_sca_adv_order (max_dom) = 5, ; horizontal scalar advection order
849 v_sca_adv_order (max_dom) = 3, ; vertical scalar advection order
851 ; advection options for scalar variables: 0=simple, 1=positive definite, 2=monotonic
852 moist_adv_opt (max_dom) = 1 ; for moisture
853 scalar_adv_opt (max_dom) = 1 ; for scalars
854 chem_adv_opt (max_dom) = 1 ; for chem variables
855 tracer_adv_opt (max_dom) = 1 ; for tracer variables (WRF-Chem activated)
856 tke_adv_opt (max_dom) = 1 ; for tke
858 time_step_sound (max_dom) = 4 / ; number of sound steps per time-step (0=set automatically)
859 (if using a time_step much larger than 6*dx (in km),
860 proportionally increase number of sound steps - also
861 best to use even numbers)
862 do_avgflx_em (max_dom) = 0, ; whether to output time-averaged mass-coupled advective velocities
863 0 = no (default)
864 1 = yes
865 do_avgflx_cugd (max_dom) = 0, ; whether to output time-averaged convective mass-fluxes from Grell-Devenyi ensemble scheme
866 0 = no (default)
867 1 = yes (only takes effect if do_avgflx_em=1 and cu_physics= 3
868 do_coriolis (max_dom) = .true., ; whether to do Coriolis calculations (idealized) (inactive)
869 do_curvature (max_dom) = .true., ; whether to do curvature calculations (idealized) (inactive)
870 do_gradp (max_dom) = .true., ; whether to do horizontal pressure gradient calculations (idealized) (inactive)
871 fft_filter_lat = 45. ; the latitude above which the polar filter is turned on
873 gwd_opt = 0 ; for running without gravity wave drag
874 = 1 ; for running the WRF-ARW with its gravity wave drag
875 = 2 ; for running the WRF-NMM with its gravity wave drag
876 sfs_opt (max_dom) = 0 ; nonlinear backscatter and anisotropy (NBA) off
877 = 1 ; NBA1 using diagnostic stress terms (km_opt=2,3 for scalars)
878 = 2 ; NBA2 using tke-based stress terms (km_opt=2 needed)
879 m_opt (max_dom) = 0 ; no added output
880 = 1 ; adds output of Mij stress terms when NBA is not used
881 tracer_opt(max_dom) = 0 ;
883 &bdy_control
884 spec_bdy_width = 5, ; total number of rows for specified boundary value nudging
885 spec_zone = 1, ; number of points in specified zone (spec b.c. option)
886 relax_zone = 4, ; number of points in relaxation zone (spec b.c. option)
887 specified (max_dom) = .false., ; specified boundary conditions (only can be used for domain 1)
888 the above 4 are used for real-data runs
889 spec_exp = 0. ; exponential multiplier for relaxation zone ramp for specified=.t.
890 (0.=linear ramp default, e.g. 0.33=~3*dx exp decay factor)
891 constant_bc = .false. ; constant boundary condition used with DFI
893 periodic_x (max_dom) = .false., ; periodic boundary conditions in x direction
894 symmetric_xs (max_dom) = .false., ; symmetric boundary conditions at x start (west)
895 symmetric_xe (max_dom) = .false., ; symmetric boundary conditions at x end (east)
896 open_xs (max_dom) = .false., ; open boundary conditions at x start (west)
897 open_xe (max_dom) = .false., ; open boundary conditions at x end (east)
898 periodic_y (max_dom) = .false., ; periodic boundary conditions in y direction
899 symmetric_ys (max_dom) = .false., ; symmetric boundary conditions at y start (south)
900 symmetric_ye (max_dom) = .false., ; symmetric boundary conditions at y end (north)
901 open_ys (max_dom) = .false., ; open boundary conditions at y start (south)
902 open_ye (max_dom) = .false., ; open boundary conditions at y end (north)
903 nested (max_dom) = .false., ; nested boundary conditions (must be used for nests)
904 polar = .false., ; polar boundary condition
905 (v=0 at polarward-most v-point)
906 euler_adv = .false., ; conservative Eulerian passive advection (NMM only)
907 idtadt = 1, ; fundamental timesteps between calls to Euler advection, dynamics (NMM only)
908 idtadc = 1 ; fundamental timesteps between calls to Euler advection, chemistry (NMM only)
912 &tc ; controls for tc_em.exe ONLY, no impact on real, ndown, or model
914 insert_bogus_storm = .false. ; T/F for inserting a bogus tropical storm (TC)
915 remove_storm = .false. ; T/F for only removing the original TC
916 num_storm = 1 ; Number of bogus TC
917 latc_loc = -999. ; center latitude of the bogus TC
918 lonc_loc = -999. ; center longitude of the bogus TC
919 vmax_meters_per_second(max_bogus) = -999. ; vmax of bogus storm in meters per second
920 rmax = -999. ; maximum radius outward from storm center
921 vmax_ratio(max_bogus) = -999. ; ratio for representative maximum winds, 0.75 for 45 km grid, and
922 0.9 for 15 km grid.
923 rankine_lid = -999. ; top pressure limit for the tc bogus scheme
925 &namelist_quilt This namelist record controls asynchronized I/O for MPI applications.
927 nio_tasks_per_group = 0, default value is 0: no quilting; > 0 quilting I/O
928 nio_groups = 1, default 1, don't change
931 &grib2:
932 background_proc_id = 255, ; Background generating process identifier, typically defined
933 by the originating center to identify the background data that
934 was used in creating the data. This is octet 13 of Section 4
935 in the grib2 message
936 forecast_proc_id = 255, ; Analysis or generating forecast process identifier, typically
937 defined by the originating center to identify the forecast process
938 that was used to generate the data. This is octet 14 of Section
939 4 in the grib2 message
940 production_status = 255, ; Production status of processed data in the grib2 message.
941 See Code Table 1.3 of the grib2 manual. This is octet 20 of
942 Section 1 in the grib2 record
943 compression = 40, ; The compression method to encode the output grib2 message.
944 Only 40 for jpeg2000 or 41 for PNG are supported