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.
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
50 history_interval (max_dom) = 60, ; history output file interval in minutes
51 frames_per_outfile (max_dom) = 1, ; output times per history output file, used to split output files
53 restart = F, ; whether this run is a restart run
54 cycling = F, ; whether this run is a cycling run, if so, initializes look-up table for Thompson schemes only
55 restart_interval = 1440, ; restart output file interval in minutes
56 reset_simulation_start = F, ; whether to overwrite simulation_start_date with forecast start time
57 io_form_history = 2, ; 2 = netCDF
58 io_form_restart = 2, ; 2 = netCDF
59 io_form_input = 2, ; 2 = netCDF
60 io_form_boundary = 2, ; netCDF format
64 = 11, ; pnetCDF format
65 frames_per_emissfile = 12, ; Number of times in each chemistry emission file.
66 io_style_emiss = 1, ; Style to use for the chemistry emission files.
67 ; 0 = Do not read emissions from files.
68 ; 1 = Cycle between two 12 hour files (set frames_per_emissfile=12)
69 ; 2 = Dated files with length set by frames_per_emissfile
70 debug_level = 0, ; 50,100,200,300 values give increasing prints
72 To choose between SI and WPS input to real for EM core:
73 auxinput1_inname = "met_em.d<domain>.<date>" ; Input to real from WPS
74 = "wrf_real_input_em.d<domain>.<date>" ; Input to real from SI
76 To choose between SI and WPS input to real for NMM core:
77 auxinput1_inname = "met_nm.d<domain>.<date>" ; Input to real from WPS
78 = "wrf_real_input_nm.d<domain>.<date>" ; Input to real from SI
82 auxhist2_outname = "rainfall" ; file name for extra output; if not specified,
83 auxhist2_d<domain>_<date> will be used
84 also note that to write variables in output other
85 than the history file requires Registry.EM file change
86 auxhist2_interval (max_dom) = 10, ; interval in minutes
87 io_form_auxhist2 = 2, ; output in netCDF
89 For SST updating (used only with sst_update=1):
91 auxinput4_inname = "wrflowinp_d<domain>"
92 auxinput4_interval = 360 ; minutes generally matches time given by interval_seconds
93 io_form_auxinput4 = 2 ; IO format, required in V3.2
95 Additional settings when running WRFVAR:
97 write_input = t, ; write input-formatted data as output
98 inputout_interval = 180, ; interval in minutes when writing input-formatted data
99 input_outname = 'wrfinput_d<domain>_<date>' ; you may change the output file name
101 inputout_begin_mo = 0
111 inputout_end_s = 0 ; the above shows that the input-formatted data are output
112 starting from hour 3 to hour 12 in 180 min interval.
115 time_step = 60, ; time step for integration in integer seconds
116 recommend 6*dx (in km) for typical real-data cases
117 time_step_fract_num = 0, ; numerator for fractional time step
118 time_step_fract_den = 1, ; denominator for fractional time step
119 Example, if you want to use 60.3 sec as your time step,
120 set time_step = 60, time_step_fract_num = 3, and
121 time_step_fract_den = 10
122 time_step_dfi = 60, ; time step for DFI, may be different from regular time_step
123 max_dom = 1, ; number of domains - set it to > 1 if it is a nested run
124 s_we (max_dom) = 1, ; start index in x (west-east) direction (leave as is)
125 e_we (max_dom) = 91, ; end index in x (west-east) direction (staggered dimension)
126 s_sn (max_dom) = 1, ; start index in y (south-north) direction (leave as is)
127 e_sn (max_dom) = 82, ; end index in y (south-north) direction (staggered dimension)
128 s_vert (max_dom) = 1, ; start index in z (vertical) direction (leave as is)
129 e_vert (max_dom) = 28, ; end index in z (vertical) direction (staggered dimension)
130 Note: this refers to full levels including surface and top
131 vertical dimensions need to be the same for all nests
132 Note: most variables are unstaggered (= staggered dim - 1)
133 dx (max_dom) = 10000, ; grid length in x direction; ARW: unit in meters, NMM: unit in degrees (e.g. 0.667)
134 dy (max_dom) = 10000, ; grid length in y direction; ARW: unit in meters, NMM: unit in degrees (e.g. 0.0658)
135 ztop (max_dom) = 19000. ; used in mass model for idealized cases
136 grid_id (max_dom) = 1, ; domain identifier
137 parent_id (max_dom) = 0, ; id of the parent domain
138 i_parent_start (max_dom) = 0, ; starting LLC I-indices from the parent domain
139 j_parent_start (max_dom) = 0, ; starting LLC J-indices from the parent domain
140 parent_grid_ratio (max_dom) = 1, ; parent-to-nest domain grid size ratio: for real-data cases
141 the ratio has to be odd; for idealized cases,
142 the ratio can be even if feedback is set to 0. (NMM: must be 3)
143 parent_time_step_ratio (max_dom) = 1, ; parent-to-nest time step ratio; it can be different
144 from the parent_grid_ratio (NMM: must be 3)
145 feedback = 1, ; feedback from nest to its parent domain; 0 = no feedback
146 smooth_option = 0 ; smoothing option for parent domain, used only with feedback
147 option on. 0: no smoothing; 1: 1-2-1 smoothing; 2: smoothing-desmoothing
149 Namelist variables specifically for the WPS input for real:
151 num_metgrid_soil_levels = 4 ; number of vertical soil levels or layers input
152 ; from WPS metgrid program
153 num_metgrid_levels = 27 ; number of vertical levels of 3d meteorological fields coming
154 ; from WPS metgrid program
155 interp_type = 2 ; vertical interpolation
156 ; 1 = linear in pressure
157 ; 2 = linear in log(pressure)
158 extrap_type = 2 ; vertical extrapolation of non-temperature fields
159 ; 1 = extrapolate using the two lowest levels
160 ; 2 = use lowest level as constant below ground
161 t_extrap_type = 2 ; vertical extrapolation for potential temperature
163 ; 2 = -6.5 K/km lapse rate for temperature
165 use_levels_below_ground = .true. ; in vertical interpolation, use levels below input surface level
166 ; T = use input isobaric levels below input surface
167 ; F = extrapolate when WRF location is below input surface value
168 use_surface = .true. ; use the input surface level data in the vertical interp and extrap
169 ; T = use the input surface data
170 ; F = do not use the input surface data
171 lagrange_order = 1 ; vertical interpolation order
174 zap_close_levels = 500 ; ignore isobaric level above surface if delta p (Pa) < zap_close_levels
175 lowest_lev_from_sfc = .false. ; place the surface value into the lowest eta location
176 ; T = use surface value as lowest eta (u,v,t,q)
177 ; F = use traditional interpolation
178 force_sfc_in_vinterp = 1 ; use the surface level as the lower boundary when interpolating
179 ; through this many eta levels
180 ; 0 = perform traditional trapping interpolation
181 ; n = first n eta levels directly use surface level
182 sfcp_to_sfcp = .false. ; Optional method to compute model's surface pressure when incoming
183 ; data only has surface pressure and terrain, but not SLP
184 smooth_cg_topo = .false. ; Smooth the outer rows and columns of domain 1's topography w.r.t.
186 use_tavg_for_tsk = .false. ; whether to use diurnally averaged surface temp as skin temp. The
187 diurnall averaged surface temp can be computed using WPS utility
188 avg_tsfc.exe. May use this option when SKINTEMP is not present.
189 p_top_requested = 5000 ; p_top (Pa) to use in the model
190 ptsgm = 42000. ; FOR NMM: defines the pressure interface dividing
191 ; the terrain following portion of the hybrid vertical
192 ; coordinate (p > ptsgm) and the purely
193 ; isobaric portion of the vertical coordinate (p < ptsgm)
194 vert_refine_fact = 1 ; vertical refinement factor for ndown
196 Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number
197 of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request
198 only the number of levels (with e_vert), and the real program will compute values. The computation assumes
199 a known first several layers, then generates equi-height spaced levels up to the top of the model.
201 eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946,
202 0.922, 0.894, 0.860, 0.817, 0.766,
203 0.707, 0.644, 0.576, 0.507, 0.444,
204 0.380, 0.324, 0.273, 0.228, 0.188,
205 0.152, 0.121, 0.093, 0.069, 0.048,
207 eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
208 0.934, 0.909, 0.880, 0.845, 0.807,
209 0.765, 0.719, 0.672, 0.622, 0.571,
210 0.520, 0.468, 0.420, 0.376, 0.335,
211 0.298, 0.263, 0.231, 0.202, 0.175,
212 0.150, 0.127, 0.106, 0.088, 0.070,
213 0.055, 0.040, 0.026, 0.013, 0.000
215 Namelist variables for controling the specified moving nest:
216 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
217 to the ARCHFLAGS. The maximum number of moves, max_moves, is set to 50
218 but can be modified in source code file frame/module_driver_constants.F.
219 num_moves = 4 ; total number of moves
220 move_id(max_moves) = 2,2,2,2, ; a list of nest domain id's, one per move
221 move_interval(max_moves) = 60,120,150,180, ; time in minutes since the start of this domain
222 move_cd_x(max_moves) = 1,1,0,-1,; the number of parent domain grid cells to move in i direction
223 move_cd_y(max_moves) = 1,0,-1,1,; the number of parent domain grid cells to move in j direction
224 positive is to move in increasing i and j direction, and
225 negative is to move in decreasing i and j direction.
226 0 means no move. The limitation now is to move only 1 grid cell
229 Namelist variables for controling the automatic moving nest:
230 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
231 and -DVORTEX_CENTER to the ARCHFLAGS. This option uses an mid-level vortex following algorthm to
232 determine the nest move. This option is experimental.
233 vortex_interval(max_dom) = 15 ; how often the new vortex position is computed
234 max_vortex_speed(max_dom) = 40 ; used to compute the search radius for the new vortex position
235 corral_dist(max_dom) = 8 ; how many coarse grid cells the moving nest is allowed to get
236 near the mother domain boundary
237 track_level = 50000 ; pressure value in Pa where the vortex is tracked
238 time_to_move(max_dom) = 0. ; time (in minutes) to start the moving nests
240 tile_sz_x = 0, ; number of points in tile x direction
241 tile_sz_y = 0, ; number of points in tile y direction
242 can be determined automatically
243 numtiles = 1, ; number of tiles per patch (alternative to above two items)
244 nproc_x = -1, ; number of processors in x for decomposition
245 nproc_y = -1, ; number of processors in y for decomposition
246 -1: code will do automatic decomposition
247 >1: for both: will be used for decomposition
249 Namelist variables for controlling the adaptive time step option:
250 These options are only valid for the ARW core.
251 use_adaptive_time_step = .false. ; T/F use adaptive time stepping, ARW only
252 step_to_output_time = .true. ; if adaptive time stepping, T/F modify the
253 time steps so that the exact history time is reached
254 target_cfl(max_dom) = 1.2,1.2 ; vertical and horizontal CFL <= to this value implies
255 no reason to reduce the time step, and to increase it
256 max_step_increase_pct(max_dom) = 5,51 ; percentage of previous time step to increase, if the
257 max(vert cfl, horiz cfl) <= target_cfl, then the time
258 will increase by max_step_increase_pct. Use something
259 large for nests (51% suggested)
260 starting_time_step(max_dom) = -1,-1 ; flag = -1 implies use 6 * dx (defined in start_em),
261 starting_time_step = 100 means the starting time step
262 for the coarse grid is 100 s
263 max_time_step(max_dom) = -1,-1 ; flag = -1 implies max time step is 3 * starting_time_step,
264 max_time_step = 100 means that the time step will not
266 min_time_step(max_dom) = -1,-1 ; flag = -1 implies max time step is 0.5 * starting_time_step,
267 min_time_step = 100 means that the time step will not
269 adaptation_domain = 1 ; default, all fine grid domains adaptive dt driven by coarse-grid
270 ; 2 = Fine grid domain #2 determines the fundamental adaptive dt.
273 dfi_opt = 0 ; which DFI option to use (3 is recommended)
274 ; 0 = no digital filter initialization
275 ; 1 = digital filter launch (DFL)
276 ; 2 = diabatic DFI (DDFI)
277 ; 3 = twice DFI (TDFI)
278 dfi_nfilter = 7 ; digital filter type to use (7 is recommended)
287 ; 8 = recursive high-order
288 dfi_write_filtered_input = .true. ; whether to write wrfinput file with filtered
289 ; model state before beginning forecast
290 dfi_write_dfi_history = .false. ; whether to write wrfout files during filtering integration
291 dfi_cutoff_seconds = 3600 ; cutoff period, in seconds, for the filter
292 dfi_time_dim = 1000 ; maximum number of time steps for filtering period
293 ; this value can be larger than necessary
294 dfi_bckstop_year = 2004 ; four-digit year of stop time for backward DFI integration
295 dfi_bckstop_month = 03 ; two-digit month of stop time for backward DFI integration
296 dfi_bckstop_day = 14 ; two-digit day of stop time for backward DFI integration
297 dfi_bckstop_hour = 12 ; two-digit hour of stop time for backward DFI integration
298 dfi_bckstop_minute = 00 ; two-digit minute of stop time for backward DFI integration
299 dfi_bckstop_second = 00 ; two-digit second of stop time for backward DFI integration
300 dfi_fwdstop_year = 2004 ; four-digit year of stop time for forward DFI integration
301 dfi_fwdstop_month = 03 ; two-digit month of stop time for forward DFI integration
302 dfi_fwdstop_day = 13 ; two-digit month of stop time for forward DFI integration
303 dfi_fwdstop_hour = 12 ; two-digit month of stop time for forward DFI integration
304 dfi_fwdstop_minute = 00 ; two-digit month of stop time for forward DFI integration
305 dfi_fwdstop_second = 00 ; two-digit month of stop time for forward DFI integration
306 dfi_radar = 0 ; DFI radar da switch
310 Note: even the physics options can be different in different nest domains,
311 caution must be used as what options are sensible to use
313 chem_opt = 0, ; chemistry option - use WRF-Chem
314 mp_physics (max_dom) microphysics option
317 = 2, Lin et al. scheme
318 = 3, WSM 3-class simple ice scheme
319 = 4, WSM 5-class scheme
320 = 5, Ferrier (new Eta) microphysics
321 = 6, WSM 6-class graupel scheme
322 = 7, Goddard GCE scheme (also uses gsfcgce_hail, gsfcgce_2ice)
323 = 8, Thompson scheme (new for V3.1)
324 = 9, Milbrandt-Yau 2-moment scheme (new for V3.2)
325 = 10, Morrison (2 moments)
326 = 13, SBU_YLIN scheme
327 = 14, WDM 5-class scheme
328 = 16, WDM 6-class scheme
329 = 98, Thompson scheme (version from V3.0)
331 For non-zero mp_physics options, to keep Qv .GE. 0, and to set the other moisture
332 fields .LT. a critcal value to zero
334 mp_zero_out = 0, ; no action taken, no adjustment to any moist field
335 = 1, ; except for Qv, all other moist arrays are set to zero
336 ; if they fall below a critical value
337 = 2, ; Qv is .GE. 0, all other moist arrays are set to zero
338 ; if they fall below a critical value
339 mp_zero_out_thresh = 1.e-8 ; critical value for moist array threshold, below which
340 ; moist arrays (except for Qv) are set to zero (kg/kg)
342 gsfcgce_hail = 0 ; for running gsfcgce microphysics with graupel
343 = 1 ; for running gsfcgce microphysics with hail
345 gsfcgce_2ice = 0 ; for running with snow, ice and graupel/hail
346 = 1 ; for running with only ice and snow
347 = 2 ; for running with only ice and graupel
348 (only used in very extreme situation)
350 gsfcgce_hail is ignored if gsfcgce_2ice is set to 1 or 2.
352 no_mp_heating = 0 ; normal
353 = 1 ; turn off latent heating from a microphysics scheme
355 ra_lw_physics (max_dom) longwave radiation option
356 = 0, no longwave radiation
359 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
361 = 5, Goddard longwave scheme
362 = 31, Earth Held-Suarez forcing
363 = 99, GFDL (Eta) longwave (semi-supported)
364 also must use co2tf = 1 for ARW
366 ra_sw_physics (max_dom) shortwave radiation option
367 = 0, no shortwave radiation
369 = 2, Goddard short wave
371 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
372 = 5, Goddard shortwave scheme
374 = 99, GFDL (Eta) longwave (semi-supported)
375 also must use co2tf = 1 for ARW
377 radt (max_dom) = 30, ; minutes between radiation physics calls
378 recommend 1 min per km of dx (e.g. 10 for 10 km)
380 nrads (max_dom) = FOR NMM: number of fundamental timesteps between
381 calls to shortwave radiation; the value
382 is set in Registry.NMM but is overridden
383 by namelist value; radt will be computed
386 nradl (max_dom) = FOR NMM: number of fundamental timesteps between
387 calls to longwave radiation; the value
388 is set in Registry.NMM but is overridden
391 co2tf CO2 transmission function flag only for GFDL radiation
392 = 0, read CO2 function data from pre-generated file
393 = 1, generate CO2 functions internally in the forecast
395 ra_call_offset radiation call offset
396 = 0 (no offset), =-1 (old offset)
398 cam_abs_freq_s = 21600 CAM clearsky longwave absorption calculation frequency
399 (recommended minimum value to speed scheme up)
400 levsiz = 59 for CAM radiation input ozone levels
401 paerlev = 29 for CAM radiation input aerosol levels
402 cam_abs_dim1 = 4 for CAM absorption save array
403 cam_abs_dim2 = value of e_vert for CAM 2nd absorption save array
405 sf_sfclay_physics (max_dom) surface-layer option (old bl_sfclay_physics option)
406 = 0, no surface-layer
407 = 1, Monin-Obukhov scheme
408 = 2, Monin-Obukhov (Janjic) scheme
409 = 3, NCEP Global Forecast System scheme (NMM only)
410 = 4, QNSE surface layer
411 = 5, MYNN surface layer
412 = 7, Pleim-Xiu surface layer (ARW only)
414 sf_surface_physics (max_dom) land-surface option (old bl_surface_physics option)
415 = 0, no surface temp prediction
416 = 1, thermal diffusion scheme
417 = 2, Unified Noah land-surface model
418 = 3, RUC land-surface model
419 = 7, Pleim-Xiu LSM (ARW)
421 sf_urban_physics(max_dom) = 0, ; activate urban canopy model (in Noah LSM only)
423 = 1: Single-layer, UCM
424 = 2: Multi-layer, Building Environment Parameterization (BEP) scheme
425 (works only with MYJ and BouLac PBL)
426 = 3: Multi-layer, Building Environment Model (BEM) scheme
427 (works only with MYJ and BouLac PBL)
429 bl_pbl_physics (max_dom) boundary-layer option
430 = 0, no boundary-layer
432 = 2, Mellor-Yamada-Janjic TKE scheme
433 = 3, NCEP Global Forecast System scheme (NMM only)
434 = 4, Quasi-Normal Scale Elimination PBL
435 = 5, MYNN 2.5 level TKE scheme, works with
436 sf_sfclay_physics=1 or 2 as well as 5
437 = 6, MYNN 3rd level TKE scheme, works only
438 MYNNSFC (sf_sfclay_physics = 5)
439 = 7, ACM2 (Pleim) PBL (ARW)
440 = 8, Bougeault and Lacarrere (BouLac) PBL
443 bldt (max_dom) = 0, ; minutes between boundary-layer physics calls
445 grav_settling = 0, ; MYNN PBL only; gravitational settling of fog/cloud droplets (1=yes)
446 nphs (max_dom) = FOR NMM: number of fundamental timesteps between
447 calls to turbulence and microphysics;
448 the value is set in Registry.NMM but is
449 overridden by namelist value; bldt will
450 be computed from this.
452 cu_physics (max_dom) cumulus option
454 = 1, Kain-Fritsch (new Eta) scheme
455 = 2, Betts-Miller-Janjic scheme
456 = 3, Grell-Devenyi ensemble scheme
457 = 4, Simplified Arakawa-Schubert scheme
458 = 5, Grell 3D ensemble scheme
459 = 6, Modifed Tiedtke scheme (ARW only, from Zhang and Wang)
460 = 14, New GFS simplified Arakawa-Schubert scheme from YSU (ARW only)
461 = 99, previous Kain-Fritsch scheme
463 ishallow = 1, Shallow convection used with Grell 3D ensemble scheme (cu_physics = 5)
465 cudt = 0, ; minutes between cumulus physics calls
467 kfeta_trigger KF trigger option (cu_physics=1 only):
469 = 2, moisture-advection based trigger (Ma and Tan [2009]) - ARW only
470 = 3, RH-dependent additional perturbation to option 1 (JMA)
472 cugd_avedx ; number of grid boxes over which subsidence is spread.
473 = 1, default, for large grid distances
474 = 3, for small grid distances (DX < 5 km)
476 ncnvc (max_dom) = FOR NMM: number of fundamental timesteps between
477 calls to convection; the value is set in Registry.NMM
478 but is overridden by namelist value; cudt will be
481 tprec (max_dom) = FOR NMM: number of hours in precipitation bucket
482 theat (max_dom) = FOR NMM: number of hours in latent heating bucket
483 tclod (max_dom) = FOR NMM: number of hours in cloud fraction average
484 trdsw (max_dom) = FOR NMM: number of hours in short wave buckets
485 trdlw (max_dom) = FOR NMM: number of hours in long wave buckets
486 tsrfc (max_dom) = FOR NMM: number of hours in surface flux buckets
487 pcpflg (max_dom) = FOR NMM: logical switch for precipitation assimilation
489 isfflx = 1, ; heat and moisture fluxes from the surface
490 (only works for sf_sfclay_physics = 1)
491 1 = with fluxes from the surface
492 0 = no flux from the surface
493 with bl_pbl_physics=0 this uses tke_drag_coefficient
494 and tke_heat_flux in vertical diffusion
495 2 = use drag from sf_sfclay_physics and heat flux from
496 tke_heat_flux with bl_pbl_physics=0
497 ifsnow = 0, ; snow-cover effects
498 (only works for sf_surface_physics = 1)
499 1 = with snow-cover effect
500 0 = without snow-cover effect
501 icloud = 1, ; cloud effect to the optical depth in radiation
502 (only works for ra_sw_physics = 1 and ra_lw_physics = 1)
503 1 = with cloud effect
504 0 = without cloud effect
505 swrad_scat = 1. ; scattering tuning parameter (default 1. is 1.e-5 m2/kg)
506 surface_input_source = 1, ; where landuse and soil category data come from:
507 1 = WPS/geogrid but with dominant categories recomputed
508 2 = GRIB data from another model (only possible
509 (VEGCAT/SOILCAT are in met_em files from WPS)
510 3 = use dominant land and soil categories from WPS/geogrid
512 num_soil_layers = 5, ; number of soil layers in land surface model
513 = 5: thermal diffusion scheme
514 = 4: Noah landsurface model
515 = 6: RUC landsurface model
516 = 2: Pleim-Xu landsurface model
517 num_land_cat = 24, ; number of land categories in input data
518 num_soil_cat = 16, ; number of soil categories in input data
520 pxlsm_smois_init(max_dom) = 1 ; PXLSM Soil moisture initialization option
521 0 - From analysis, 1 - From MAVAIL
523 maxiens = 1, ; Grell-Devenyi only
524 maxens = 3, ; G-D only
525 maxens2 = 3, ; G-D only
526 maxens3 = 16 ; G-D only
527 ensdim = 144 ; G-D only
528 These are recommended numbers. If you would like to use
529 any other number, consult the code, know what you are doing.
530 seaice_threshold = 271 ; tsk < seaice_threshold, if water point and 5-layer slab
531 ; scheme, set to land point and permanent ice; if water point
532 ; and Noah scheme, set to land point, permanent ice, set temps
533 ; from 3 m to surface, and set smois and sh2o
534 sst_update = 0 ; time-varying sea-surface temp (0=no, 1=yes). If selected real
535 ; puts SST, XICE, ALBEDO and VEGFRA in wrflowinp_d01 file, and wrf updates
536 ; these from it at same interval as boundary file. Also requires
537 ; namelists in &time_control: auxinput4_interval, auxinput4_end_h,
538 ; auxinput4_inname = "wrflowinp_d<domain>",
539 ; and in V3.2 io_form_auxinput4
540 usemonalb = .true. ; use monthly albedo map instead of table value
541 ; (must be used for NMM and recommended for sst_update=1)
542 rdmaxalb = .true. ; use snow albedo from geogrid; false means using values from table
543 rdlai2d = .false. ; use LAI from input; false means using values from table
544 bucket_mm = -1. ; bucket reset value for water accumulations (value in mm, -1.=inactive)
545 bucket_J = -1. ; bucket reset value for energy accumulations (value in J, -1.=inactive)
546 tmn_update = 0 ; update deep soil temperature (1, yes; 0, no)
547 lagday = 150 ; days over which tmn is computed using skin temperature
548 sst_skin = 0 ; calculate skin SST
549 slope_rad (max_dom) = 0 ; slope effects for solar radiation (1=on, 0=off)
550 topo_shading (max_dom) = 0 ; neighboring-point shadow effects for solar radiation (1=on, 0=off)
551 shadlen = 25000. ; max shadow length in meters for topo_shading=1
552 omlcall = 0 ; activate simple ocean mixed layer model (0=no, 1=yes); works with
553 sf_surface_physics = 1 only
554 oml_hml0 = 50 ; oml model can be initialized with a constant depth everywhere (m)
555 oml_gamma = 0.14 ; oml deep water lapse rate (K m-1)
556 isftcflx = 0 ; alternative Ck, Cd formulation for tropical storm application (0=default, 1=new, 2=Garratt)
557 fractional_seaice = 0 ; treat sea-ice as fractional field (1) or ice/no-ice flag (0)
558 tice2tsk_if2cold = .false. ; set Tice to Tsk to avoid unrealistically low sea ice temperatures
559 iz0tlnd = 0 ; thermal roughness length for sfclay and myjsfc (0 - old, 1 - veg dependent Czil)
560 mp_tend_lim = 10., ; limit on temp tendency from mp latent heating from radar data assimilation
561 prec_acc_dt (max_dom) = 0., ; number of minutes in precipitation bucket (ARW only) - will add three
562 new 2d output fields: prec_acc_c, prec_acc_nc and snow_acc_nc
565 grid_fdda (max_dom) = 1 ; grid-nudging fdda on (=0 off) for each domain
566 = 2 ; spectral nudging
567 gfdda_inname = "wrffdda_d<domain>" ; defined name in real
568 gfdda_interval_m (max_dom) = 360 ; time interval (in min) between analysis times (must use minutes)
569 gfdda_end_h (max_dom) = 6 ; time (in hours) to stop nudging after start of forecast
570 io_form_gfdda = 2 ; analysis data io format (2 = netCDF)
571 fgdt (max_dom) = 0 ; calculation frequency (minutes) for grid-nudging (0=every step)
572 if_no_pbl_nudging_uv (max_dom) = 0 ; 1= no nudging of u and v in the pbl, 0=nudging in the pbl
573 if_no_pbl_nudging_t (max_dom) = 0 ; 1= no nudging of temp in the pbl, 0=nudging in the pbl
574 if_no_pbl_nudging_q (max_dom) = 0 ; 1= no nudging of qvapor in the pbl, 0=nudging in the pbl
575 if_zfac_uv (max_dom) = 0 ; 0= nudge u and v in all layers, 1= limit nudging to levels above k_zfac_uv
576 k_zfac_uv (max_dom) = 10 ; 10=model level below which nudging is switched off for u and v
577 if_zfac_t (max_dom) = 0 ; 0= nudge temp in all layers, 1= limit nudging to levels above k_zfac_t
578 k_zfac_t (max_dom) = 10 ; 10=model level below which nudging is switched off for temp
579 if_zfac_q (max_dom) = 0 ; 0= nudge qvapor in all layers, 1= limit nudging to levels above k_zfac_q
580 k_zfac_q (max_dom) = 10 ; 10=model level below which nudging is switched off for qvapor
581 guv (max_dom) = 0.0003 ; nudging coefficient for u and v (sec-1)
582 gt (max_dom) = 0.0003 ; nudging coefficient for temp (sec-1)
583 gq (max_dom) = 0.0003 ; nudging coefficient for qvapor (sec-1)
584 if_ramping = 0 ; 0= nudging ends as a step function, 1= ramping nudging down at end of period
585 dtramp_min = 60.0 ; time (min) for ramping function, 60.0=ramping starts at last analysis time,
586 -60.0=ramping ends at last analysis time
587 grid_sfdda (max_dom) = 0 ; surface fdda switch (1, on; 0, off)
588 sgfdda_inname = "wrfsfdda_d<domain>" ; defined name for sfc nudgingi in input file (from program obsgrid)
589 sgfdda_end_h (max_dom) = 6 ; time (in hours) to stop sfc nudging after start of forecast
590 sgfdda_interval_m (max_dom) = 180 ; time interval (in min) between sfc analysis times (must use minutes)
591 io_form_sgfdda = 2 ; sfc analysis data io format (2 = netCDF)
592 guv_sfc (max_dom) = 0.0003 ; nudging coefficient for sfc u and v (sec-1)
593 gt_sfc (max_dom) = 0.0003 ; nudging coefficient for sfc temp (sec-1)
594 gq_sfc (max_dom) = 0.0003 ; nudging coefficient for sfc qvapor (sec-1)
595 rinblw = 250.0 ; radius of influence used to determine the confidence (or weights) for
596 the analysis, which is based on the distance between the grid point to the nearest
597 obs. The analysis without nearby observation is used at a reduced weight.
599 pxlsm_soil_nudge(max_dom) = 1 ; PXLSM Soil nudging option (requires wrfsfdda file)
601 The following are for spectral nudging:
602 fgdtzero (max_dom) = 0, ; 1= nudging tendencies are set to zero in between fdda calls
603 if_no_pbl_nudging_ph = 0, ; 1= no nudging of ph in the pbl, 0= nuding in the pbl
604 if_zfac_ph (max_dom) = 0, ; 0= nudge ph in all layers, 1= limit nudging to levels above k_zfac_ph
605 k_zfac_ph (max_dom) = 10, ; 10= model level below which nudging is switched off for ph
606 dk_zfac_uv (max_dom) = 1, ; depth in k between k_zfac_X to dk_zfac_X where nudging increases
607 linearly to full strength
608 dk_zfac_t (max_dom) = 1,
609 dk_zfac_ph (max_dom) = 1,
610 gph (max_dom) = 0.0003,
611 xwavenum (max_dom) = 3, ; top wave number to nudge in x direction
612 ywavenum (max_dom) = 3, ; top wave number to nudge in y direction
614 The following are for observation nudging:
615 obs_nudge_opt (max_dom) = 1 ; obs-nudging fdda on (=0 off) for each domain
616 also need to set auxinput11_interval and auxinput11_end_h
617 in time_control namelist
618 max_obs = 150000 ; max number of observations used on a domain during any
620 fdda_start = 0 ; obs nudging start time in minutes
621 fdda_end = 180 ; obs nudging end time in minutes
622 obs_nudge_wind (max_dom) = 1 ; whether to nudge wind: (=0 off)
623 obs_coef_wind = 6.E-4, ; nudging coefficient for wind, unit: s-1
624 obs_nudge_temp = 1 ; whether to nudge temperature: (=0 off)
625 obs_coef_temp = 6.E-4, ; nudging coefficient for temperature, unit: s-1
626 obs_nudge_mois = 1 ; whether to nudge water vapor mixing ratio: (=0 off)
627 obs_coef_mois = 6.E-4, ; nudging coefficient for water vapor mixing ratio, unit: s-1
628 obs_nudge_pstr = 0 ; whether to nudge surface pressure (not used)
629 obs_coef_pstr = 0. ; nudging coefficient for surface pressure, unit: s-1 (not used)
630 obs_rinxy = 200., ; horizonal radius of influence in km
631 obs_rinsig = 0.1, ; vertical radius of influence in eta
632 obs_twindo (max_dom) = 0.66667 ; half-period time window over which an observation
633 will be used for nudging (hours)
634 obs_npfi = 10, ; freq in coarse grid timesteps for diag prints
635 obs_ionf (max_dom) = 2 ; freq in coarse grid timesteps for obs input and err calc
636 obs_idynin = 0 ; for dynamic initialization using a ramp-down function to gradually
637 turn off the FDDA before the pure forecast (=1 on)
638 obs_dtramp = 40 ; time period in minutes over which the nudging is ramped down
640 obs_nobs_prt (max_dom) = 10, ; Number of current obs to print grid coord. info.
641 obs_ipf_in4dob = .true. ; print obs input diagnostics (=.false. off)
642 obs_ipf_errob = .true. ; print obs error diagnostics (=.false. off)
643 obs_ipf_nudob = .true. ; print obs nudge diagnostics (=.false. off)
644 obs_ipf_init = .true. ; Enable obs init warning messages
646 obs_no_pbl_nudge_uv (max_dom) = 0 ; 1=no wind-nudging within pbl
647 obs_no_pbl_nudge_t (max_dom) = 0 ; 1=no temperature-nudging within pbl
648 obs_no_pbl_nudge_q (max_dom) = 0 ; 1=no moisture-nudging within pbl
649 obs_nudgezfullr1_uv = 50 ; Vert infl full weight height for lowest model level (LML) ; obs, regime 1, winds
650 obs_nudgezrampr1_uv = 50 ; Vert infl ramp-to-zero height for LML obs, regime 1, winds
651 obs_nudgezfullr2_uv = 50 ; Vert infl full weight height for LML obs, regime 2, winds
652 obs_nudgezrampr2_uv = 50 ; Vert infl ramp-to-zero height for LML obs, regime 2, winds
653 obs_nudgezfullr4_uv = -5000 ; Vert infl full weight height for LML obs, regime 4, winds
654 obs_nudgezrampr4_uv = 50 ; Vert infl ramp-to-zero height for LML obs, regime 4, winds
655 obs_nudgezfullr1_t = 50 ; Vert infl full weight height for LML obs, regime 1, temperature
656 obs_nudgezrampr1_t = 50 ; Vert infl ramp-to-zero height for LML obs, regime 1, temperature
657 obs_nudgezfullr2_t = 50 ; Vert infl full weight height for LML obs, regime 2, temperature
658 obs_nudgezrampr2_t = 50 ; Vert infl ramp-to-zero height for LML obs, regime 2, temperature
659 obs_nudgezfullr4_t = -5000 ; Vert infl full weight height for LML obs, regime 4, temperature
660 obs_nudgezrampr4_t = 50 ; Vert infl ramp-to-zero height for LML obs, regime 4, temperature
661 obs_nudgezfullr1_q = 50 ; Vert infl full weight height for LML obs, regime 1, moisture
662 obs_nudgezrampr1_q = 50 ; Vert infl ramp-to-zero height for LML obs, regime 1, moisture
663 obs_nudgezfullr2_q = 50 ; Vert infl full weight height for LML obs, regime 2, moisture
664 obs_nudgezrampr2_q = 50 ; Vert infl ramp-to-zero height for LML obs, regime 2, moisture
665 obs_nudgezfullr4_q = -5000 ; Vert infl full weight height for LML obs, regime 4, moisture
666 obs_nudgezrampr4_q = 50 ; Vert infl ramp-to-zero height for LML obs, regime 4, moisture
667 obs_nudgezfullmin = 50 ; Min depth through which vertical infl fcn remains 1.0
668 obs_nudgezrampmin = 50 ; Min depth (m) through which vert infl fcn decreases from 1 to 0
669 obs_nudgezmax = 3000 ; Max depth (m) in which vert infl function is nonzero
670 obs_sfcfact = 1.0 ; Scale factor applied to time window for surface obs
671 obs_sfcfacr = 1.0 ; Scale factor applied to horiz radius of influence for surface obs
672 obs_dpsmx = 7.5 ; Max pressure change (cb) allowed within horiz radius of influence
676 scm_force = 1, ; switch for single column forcing (=0 off)
677 scm_force_dx = 4000. ; DX for SCM forcing (in meters)
678 num_force_layers = 8 ; number of SCM input forcing layers
679 scm_lu_index = 2 ; SCM landuse category (2 is dryland, cropland and pasture)
680 scm_isltyp = 4 ; SCM soil category (4 is silt loam)
681 scm_vegfra = 0.5 ; SCM vegetation fraction
682 scm_canwat = 0.0 ; SCM canopy water
683 scm_lat = 37.600 ; SCM latitude
684 scm_lon = -96.700 ; SCM longitude
685 scm_th_adv = .true. ; turn on theta advection in SCM
686 scm_wind_adv = .true. ; turn on wind advection in SCM
687 scm_qv_adv = .true. ; turn on moisture advection in SCM
688 scm_vert_adv = .true. ; turn on vertical advection in SCM
691 rk_ord = 3, ; time-integration scheme option:
692 2 = Runge-Kutta 2nd order
693 3 = Runge-Kutta 3rd order
694 diff_opt = 0, ; turbulence and mixing option:
695 0 = no turbulence or explicit
696 spatial numerical filters (km_opt IS IGNORED).
697 1 = evaluates 2nd order
698 diffusion term on coordinate surfaces.
699 uses kvdif for vertical diff unless PBL option
700 is used. may be used with km_opt = 1 and 4.
701 (= 1, recommended for real-data cases)
702 2 = evaluates mixing terms in
703 physical space (stress form) (x,y,z).
704 turbulence parameterization is chosen
705 by specifying km_opt.
706 km_opt = 1, ; eddy coefficient option
707 1 = constant (use khdif kvdif)
708 2 = 1.5 order TKE closure (3D)
709 3 = Smagorinsky first order closure (3D)
710 Note: option 2 and 3 are not recommended for DX > 2 km
711 4 = horizontal Smagorinsky first order closure
712 (recommended for real-data cases)
713 damp_opt = 0, ; upper level damping flag
715 1 = with diffusive damping, maybe used for real-data cases
716 (dampcoef nondimensional ~0.01-0.1)
717 2 = with Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .003; idealized case only
718 not for real-data cases)
719 3 = with w-Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .05;
721 diff_6th_opt = 0, ; 6th-order numerical diffusion
722 0 = no 6th-order diffusion (default)
723 1 = 6th-order numerical diffusion (not recommended)
724 2 = 6th-order numerical diffusion but prohibit up-gradient diffusion
725 diff_6th_factor = 0.12, ; 6th-order numerical diffusion non-dimensional rate (max value 1.0
726 corresponds to complete removal of 2dx wave in one timestep)
727 dampcoef (max_dom) = 0., ; damping coefficient (see above)
728 zdamp (max_dom) = 5000., ; damping depth (m) from model top
729 w_damping = 0, ; vertical velocity damping flag (for operational use)
732 base_temp = 290., ; real-data, em ONLY, base sea-level temp (K)
733 base_pres = 10^5 ; real-data, em ONLY, base sea-level pres (Pa), DO NOT CHANGE
734 base_lapse = 50., ; real-data, em ONLY, lapse rate (K), DO NOT CHANGE
735 iso_temp = 0., ; real-data, em ONLY, reference temp in stratosphere
736 khdif (max_dom) = 0, ; horizontal diffusion constant (m^2/s)
737 kvdif (max_dom) = 0, ; vertical diffusion constant (m^2/s)
738 smdiv (max_dom) = 0.1, ; divergence damping (0.1 is typical)
739 emdiv (max_dom) = 0.01, ; external-mode filter coef for mass coordinate model
740 (0.01 is typical for real-data cases)
741 epssm (max_dom) = .1, ; time off-centering for vertical sound waves
742 non_hydrostatic (max_dom) = .true., ; whether running the model in hydrostatic or non-hydro mode
743 pert_coriolis (max_dom) = .false., ; Coriolis only acts on wind perturbation (idealized)
744 top_lid (max_dom) = .false., ; Zero vertical motion at top of domain
745 mix_full_fields(max_dom) = .true., ; used with diff_opt = 2; value of ".true." is recommended, except for
746 highly idealized numerical tests; damp_opt must not be 1 if ".true."
747 is chosen. .false. means subtract 1-d base-state profile before mixing
748 mix_isotropic(max_dom) = 0 ; 0=anistropic vertical/horizontal diffusion coeffs, 1=isotropic
749 mix_upper_bound(max_dom) = 0.1 ; non-dimensional upper limit for diffusion coeffs
750 tke_drag_coefficient(max_dom) = 0., ; surface drag coefficient (Cd, dimensionless) for diff_opt=2 only
751 tke_heat_flux(max_dom) = 0., ; surface thermal flux (H/(rho*cp), K m/s) for diff_opt=2 only
752 h_mom_adv_order (max_dom) = 5, ; horizontal momentum advection order (5=5th, etc.)
753 v_mom_adv_order (max_dom) = 3, ; vertical momentum advection order
754 h_sca_adv_order (max_dom) = 5, ; horizontal scalar advection order
755 v_sca_adv_order (max_dom) = 3, ; vertical scalar advection order
757 ; advection options for scalar variables: 0=simple, 1=positive definite, 2=monotonic
758 moist_adv_opt (max_dom) = 1 ; for moisture
759 scalar_adv_opt (max_dom) = 1 ; for scalars
760 chem_adv_opt (max_dom) = 1 ; for chem variables
761 tracer_adv_opt (max_dom) = 1 ; for tracer variables (WRF-Chem activated)
762 tke_adv_opt (max_dom) = 1 ; for tke
764 time_step_sound (max_dom) = 4 / ; number of sound steps per time-step (0=set automatically)
765 (if using a time_step much larger than 6*dx (in km),
766 proportionally increase number of sound steps - also
767 best to use even numbers)
768 do_avgflx_em (max_dom) = 0, ; whether to output time-averaged mass-coupled advective velocities
771 do_avgflx_cugd (max_dom) = 0, ; whether to output time-averaged convective mass-fluxes from Grell-Devenyi ensemble scheme
773 1 = yes (only takes effect if do_avgflx_em=1 and cu_physics= 3
774 do_coriolis (max_dom) = .true., ; whether to do Coriolis calculations (idealized) (inactive)
775 do_curvature (max_dom) = .true., ; whether to do curvature calculations (idealized) (inactive)
776 do_gradp (max_dom) = .true., ; whether to do horizontal pressure gradient calculations (idealized) (inactive)
777 fft_filter_lat = 45. ; the latitude above which the polar filter is turned on
779 gwd_opt = 0 ; for running without gravity wave drag
780 = 1 ; for running the WRF-ARW with its gravity wave drag
781 = 2 ; for running the WRF-NMM with its gravity wave drag
782 sfs_opt (max_dom) = 0 ; nonlinear backscatter and anisotropy (NBA) off
783 = 1 ; NBA1 using diagnostic stress terms (km_opt=2,3 for scalars)
784 = 2 ; NBA2 using tke-based stress terms (km_opt=2 needed)
785 m_opt (max_dom) = 0 ; no added output
786 = 1 ; adds output of Mij stress terms when NBA is not used
787 tracer_opt(mac_dom) = 0 ;
790 spec_bdy_width = 5, ; total number of rows for specified boundary value nudging
791 spec_zone = 1, ; number of points in specified zone (spec b.c. option)
792 relax_zone = 4, ; number of points in relaxation zone (spec b.c. option)
793 specified (max_dom) = .false., ; specified boundary conditions (only can be used for domain 1)
794 the above 4 are used for real-data runs
795 spec_exp = 0. ; exponential multiplier for relaxation zone ramp for specified=.t.
796 (0.=linear ramp default, e.g. 0.33=~3*dx exp decay factor)
797 constant_bc = .false. ; constant boundary condition used with DFI
799 periodic_x (max_dom) = .false., ; periodic boundary conditions in x direction
800 symmetric_xs (max_dom) = .false., ; symmetric boundary conditions at x start (west)
801 symmetric_xe (max_dom) = .false., ; symmetric boundary conditions at x end (east)
802 open_xs (max_dom) = .false., ; open boundary conditions at x start (west)
803 open_xe (max_dom) = .false., ; open boundary conditions at x end (east)
804 periodic_y (max_dom) = .false., ; periodic boundary conditions in y direction
805 symmetric_ys (max_dom) = .false., ; symmetric boundary conditions at y start (south)
806 symmetric_ye (max_dom) = .false., ; symmetric boundary conditions at y end (north)
807 open_ys (max_dom) = .false., ; open boundary conditions at y start (south)
808 open_ye (max_dom) = .false., ; open boundary conditions at y end (north)
809 nested (max_dom) = .false., ; nested boundary conditions (must be used for nests)
810 polar = .false., ; polar boundary condition
811 (v=0 at polarward-most v-point)
812 euler_adv = .false., ; conservative Eulerian passive advection (NMM only)
813 idtadt = 1, ; fundamental timesteps between calls to Euler advection, dynamics (NMM only)
814 idtadc = 1 ; fundamental timesteps between calls to Euler advection, chemistry (NMM only)
818 &tc ; controls for tc_em.exe ONLY, no impact on real, ndown, or model
820 insert_bogus_storm = .false. ; T/F for inserting a bogus tropical storm (TC)
821 remove_storm = .false. ; T/F for only removing the original TC
822 num_storm = 1 ; Number of bogus TC
823 latc_loc = -999. ; center latitude of the bogus TC
824 lonc_loc = -999. ; center longitude of the bogus TC
825 vmax_meters_per_second = -999. ; vmax of bogus storm in meters per second
826 rmax = -999. ; maximum radius outward from storm center
827 vmax_ratio = -999. ; ratio for representative maximum winds, 0.75 for 45 km grid, and
830 &namelist_quilt This namelist record controls asynchronized I/O for MPI applications.
832 nio_tasks_per_group = 0, default value is 0: no quilting; > 0 quilting I/O
833 nio_groups = 1, default 1, don't change
837 background_proc_id = 255, ; Background generating process identifier, typically defined
838 by the originating center to identify the background data that
839 was used in creating the data. This is octet 13 of Section 4
841 forecast_proc_id = 255, ; Analysis or generating forecast process identifier, typically
842 defined by the originating center to identify the forecast process
843 that was used to generate the data. This is octet 14 of Section
844 4 in the grib2 message
845 production_status = 255, ; Production status of processed data in the grib2 message.
846 See Code Table 1.3 of the grib2 manual. This is octet 20 of
847 Section 1 in the grib2 record
848 compression = 40, ; The compression method to encode the output grib2 message.
849 Only 40 for jpeg2000 or 41 for PNG are supported