| blob | 6e7ebc1a3ddf167651afd80cb3aae66e451ce72f |
1 Description of namelist variables
3 ---------------------------------
5 Note: variables followed by (max_dom) indicate that this variable needs to
6 be defined for the nests when max_dom > 1.
8 &time_control
9 run_days = 0, ! run time in days
10 run_hours = 0, ! run time in hours
11 Note: if it is more than 1 day, one may use both run_days and run_hours
12 or just run_hours. e.g. if the total run length is 36 hrs, you may
13 set run_days = 1, and run_hours = 12, or run_days = 0, and run_hours = 36.
14 run_minutes = 0, ! run time in minutes
15 run_seconds = 0, ! run time in seconds
16 start_year (max_dom) = 2001, ! four digit year of starting time
17 start_month (max_dom) = 06, ! two digit month of starting time
18 start_day (max_dom) = 11, ! two digit day of starting time
19 start_hour (max_dom) = 12, ! two digit hour of starting time
20 start_minute (max_dom) = 00, ! two digit minute of starting time
21 start_second (max_dom) = 00, ! two digit second of starting time
22 Note: the start time is used to name the first wrfout file.
23 It also controls the start time for nest domains, and the time to restart
24 end_year (max_dom) = 2001, ! four digit year of ending time
25 end_month (max_dom) = 06, ! two digit month of ending time
26 end_day (max_dom) = 12, ! two digit day of ending time
27 end_hour (max_dom) = 12, ! two digit hour of ending time
28 end_minute (max_dom) = 00, ! two digit minute of ending time
29 end_second (max_dom) = 00, ! two digit second of ending time
30 It also controls when the nest domain integrations end
31 All start and end times are used by real.exe.
33 Note that one may use either run_days/run_hours etc. or
34 end_year/month/day/hour etc. to control the length of
35 model integration. But run_days/run_hours
36 takes precedence over the end times.
37 Program real.exe uses start and end times only.
39 interval_seconds = 10800, ! time interval between incoming real data, which will be the interval
40 between the lateral boundary condition file (in seconds)
41 input_from_file (max_dom) = T, ! whether nested run will have input files for domains other than 1
42 fine_input_stream (max_dom) = 0, ! field selection from nest input for its initialization
43 0: all fields are used! 2: only static and time-varying, masked land
44 surface fields are used. This requires the use of
45 io_form_auxinput2
46 history_interval (max_dom) = 60, ! history output file interval in minutes
47 frames_per_outfile (max_dom) = 1, ! number of output times per history output file,
48 used to split output into multiple files
49 into smaller pieces
50 restart = F, ! whether this run is a restart run
51 cycling = F, ! whether this run is a cycling run, if so, initializes look-up table for Thompson schemes only
52 restart_interval = 1440, ! restart output file interval in minutes
53 reset_simulation_start = F, ! whether to overwrite simulation_start_date with forecast start time
54 io_form_history = 2, ! 2 = netCDF
55 io_form_restart = 2, ! 2 = netCDF
56 io_form_input = 2, ! 2 = netCDF
57 io_form_boundary = 2, ! netCDF format
58 = 4, ! PHD5 format
59 = 5, ! GRIB1 format
60 = 10, ! GRIB2 format
61 = 11, ! pnetCDF format
62 ncd_nofill = .true., ! only a single write, not the write/read/write sequence
63 frames_per_emissfile = 12, ! number of times in each chemistry emission file.
64 io_style_emiss = 1, ! style to use for the chemistry emission files.
65 ! 0 = Do not read emissions from files.
66 ! 1 = Cycle between two 12 hour files (set frames_per_emissfile=12)
67 ! 2 = Dated files with length set by frames_per_emissfile
68 debug_level = 0, ! 50,100,200,300 values give increasing prints
69 diag_print = 0, ! print out time series of model diagnostics
70 0 = no print
71 1 = domain averaged 3-hourly hydrostatic surface pressure tendency
72 (Dpsfc/Dt), and dry-hydrostatic column pressure tendency (Dmu/Dt)
73 will appear in stdout file
74 2 = in addition to those above, domain averaged rainfall,
75 surface evaporation, and sensible and latent heat fluxes will be output
76 all_ic_times = .false., ! whether to write out wrfinput for all processing times
77 adjust_output_times = .false., ! adjust output times to the nearest hour
78 override_restart_timers = .false., ! whether to change the alarms from what is previously set
79 write_hist_at_0h_rst = .false., ! whether to output history file at the start of restart run
80 write_restart_at_0h = .false., ! whether to output restart file at the start of restart run
81 output_ready_flag = .true., ! asks the model to write-out an empty file with the name
82 'wrfoutReady_d<domain>_<date>.
83 Useful in production runs so that post-processing code can check on the
84 completeness of this file
85 force_use_old_data = .false., ! if set to .true., allow WRF Version 3 input data
86 =.false., stop when WRF model detects Version 3 input data
88 To choose between SI and WPS input to real for EM core:
89 auxinput1_inname = "met_em.d<domain>.<date>" ! Input to real from WPS (default since 3.0)
90 = "wrf_real_input_em.d<domain>.<date>" ! Input to real from SI
92 Other output options:
94 auxhist2_outname = "rainfall" ! file name for extra output! if not specified,
95 auxhist2_d<domain>_<date> will be used
96 also note that to write variables in output other
97 than the history file requires Registry.EM file change
98 auxhist2_interval (max_dom) = 10, ! interval in minutes
99 io_form_auxhist2 = 2, ! output in netCDF
100 frames_per_auxhist2 = 1000, ! number of output times in this file
102 For SST updating (used only with sst_update=1):
104 auxinput4_inname = "wrflowinp_d<domain>"
105 auxinput4_interval = 360 ! minutes generally matches time given by interval_seconds
106 io_form_auxinput4 = 2 ! IO format
108 nwp_diagnostics = 0 ! set to = 1 to add 7 history-interval max diagnostic fields
110 For additional regional climate surface fields
112 output_diagnostics = 0 ! set to = 1 to add 36 surface diagnostic arrays (max/min/mean/std)
113 auxhist3_outname = 'wrfxtrm_d<domain>_<date>' ! file name for added diagnostics
114 io_form_auxhist3 = 2 ! netcdf
115 auxhist3_interval = 1440 ! minutes between outputs (1440 gives daily max/min)
116 frames_per_auxhist3 = 1 ! output times per file
117 Note: do restart only at multiple of auxhist3_intervals
119 For observation nudging:
120 auxinput11_interval = 10 ! interval in minutes for observation data. It should be
121 set as or more frequently as obs_ionf (with unit of
122 coarse domain time step).
123 auxinput11_end_h = 6 ! end of observation time in hours.
125 Options for run-time IO:
127 iofields_filename (max_dom) = "my_iofields_list.txt",
128 (example: +:h:21:rainc, rainnc, rthcuten)
129 ignore_iofields_warning = .true., ! what to do when encountering an error in the user-specified files
130 .false., ! abort when encountering an error in iofields_filename file
132 nocolons = .false., ! whether to write file date with colons
133 .true., e.g. 2020-05-20_12_00_00
135 Additional settings when running WRFVAR:
137 write_input = t, ! write input-formatted data as output
138 inputout_interval (max_dom) = 180, ! interval in minutes when writing input-formatted data
139 input_outname = 'wrfinput_d<domain>_<date>' ! you may change the output file name
140 inputout_begin_y (max_dom) = 0
141 inputout_begin_mo = 0
142 inputout_begin_d (max_dom) = 0
143 inputout_begin_h (max_dom) = 3
144 inputout_begin_m (max_dom) = 0
145 inputout_begin_s (max_dom) = 0
146 inputout_end_y (max_dom) = 0
147 inputout_end_mo = 0
148 inputout_end_d (max_dom) = 0
149 inputout_end_h (max_dom) = 12
150 inputout_end_m (max_dom) = 0
151 inputout_end_s (max_dom) = 0 ! the above shows that the input-formatted data are output
152 starting from hour 3 to hour 12 in 180 min interval.
154 For automatic moving nests: requires special input data, and environment variable TERRAIN_AND_LANDUSE set at compile time
155 (This option will overwrite input_from_file for nest domains)
156 input_from_hires (max_dom) = .true.,
157 rsmas_data_path = "path-to-terrain-and-landuse-dataset"
159 &domains
160 time_step = 60, ! time step for integration in integer seconds
161 recommend 6*dx (in km) for typical real-data cases
162 time_step_fract_num = 0, ! numerator for fractional time step
163 time_step_fract_den = 1, ! denominator for fractional time step
164 Example, if you want to use 60.3 sec as your time step,
165 set time_step = 60, time_step_fract_num = 3, and
166 time_step_fract_den = 10
167 time_step_dfi = 60, ! time step for DFI, may be different from regular time_step
168 reasonable_time_step_ratio = 6., ! Any d01, real-data case with a time step ratio larger than this is stopped. Except for specific circumstances (e.g., using IEVA), this value should be no larger than 6 (default).
169 max_dom = 1, ! number of domains - set it to > 1 if it is a nested run
170 s_we (max_dom) = 1, ! start index in x (west-east) direction (leave as is)
171 e_we (max_dom) = 91, ! end index in x (west-east) direction (staggered dimension)
172 s_sn (max_dom) = 1, ! start index in y (south-north) direction (leave as is)
173 e_sn (max_dom) = 82, ! end index in y (south-north) direction (staggered dimension)
174 s_vert (max_dom) = 1, ! start index in z (vertical) direction (leave as is)
175 e_vert (max_dom) = 30, ! end index in z (vertical) direction (staggered dimension)
176 Note: this refers to full levels including surface and top
177 vertical dimensions need to be the same for all nests
178 Note: most variables are unstaggered (= staggered dim - 1)
179 dx (max_dom) = 10000, ! grid length in x direction; unit in meters.
180 dy (max_dom) = 10000, ! grid length in y direction; unit in meters.
181 ztop (max_dom) = 19000. ! used in mass model for idealized cases
182 grid_id (max_dom) = 1, ! domain identifier
183 parent_id (max_dom) = 0, ! id of the parent domain
184 i_parent_start (max_dom) = 0, ! starting LLC I-indices from the parent domain
185 j_parent_start (max_dom) = 0, ! starting LLC J-indices from the parent domain
186 parent_grid_ratio (max_dom) = 1, ! parent-to-nest domain grid size ratio: for real-data cases
187 the ratio has to be odd! for idealized cases,
188 the ratio can be even if feedback is set to 0.
189 parent_time_step_ratio (max_dom) = 1, ! parent-to-nest time step ratio! it can be different
190 from the parent_grid_ratio
191 feedback = 1, ! feedback from nest to its parent domain; 0 = no feedback
192 smooth_option = 2 ! smoothing option for parent domain, used only with feedback
193 option on. 0: no smoothing; 1: 1-2-1 smoothing; 2: smoothing-desmoothing (default)
195 Namelist variables specifically for the WPS input for real:
197 num_metgrid_soil_levels = 4 ! number of vertical soil levels or layers input
198 from WPS metgrid program
199 num_metgrid_levels = 27 ! number of vertical levels of 3d meteorological fields coming
200 from WPS metgrid program
201 interp_type = 2 ! vertical interpolation
202 1 = linear in pressure
203 2 = linear in log(pressure)
204 extrap_type = 2 ! vertical extrapolation of non-temperature fields
205 1 = extrapolate using the two lowest levels
206 2 = use lowest level as constant below ground
207 t_extrap_type = 2 ! vertical extrapolation for potential temperature
208 1 = isothermal
209 2 = -6.5 K/km lapse rate for temperature
210 3 = constant theta
211 use_levels_below_ground = .true. ! in vertical interpolation, use levels below input surface level
212 T = use input isobaric levels below input surface
213 F = extrapolate when WRF location is below input surface value
214 use_surface = .true. ! use the input surface level data in the vertical interp and extrap
215 T = use the input surface data
216 F = do not use the input surface data
217 lagrange_order = 2 ! vertical interpolation order
218 1 = linear
219 2 = quadratic
220 9 = cubic spline
221 zap_close_levels = 500 ! ignore isobaric level above surface if delta p (Pa) < zap_close_levels
222 lowest_lev_from_sfc = .false. ! place the surface value into the lowest eta location
223 ! T = use surface value as lowest eta (u,v,t,q)
224 ! F = use traditional interpolation
225 force_sfc_in_vinterp = 1 ! use the surface level as the lower boundary when interpolating
226 through this many eta levels
227 0 = perform traditional trapping interpolation
228 n = first n eta levels directly use surface level
229 maxw_horiz_pres_diff = 5000 ! Pressure threshold (Pa). For using the level of max winds, when the
230 pressure difference between neighboring values exceeds this maximum,
231 the variable is NOT inserted into the column for vertical interpolation.
232 trop_horiz_pres_diff = 5000 ! Pressure threshold (Pa). For using the tropopause level, when the
233 pressure difference between neighboring values exceeds this maximum,
234 the variable is NOT inserted into the column for vertical interpolation.
235 maxw_above_this_level = 30000 ! Minimum height (actually it is pressure in Pa) to allow using the
236 level of max wind information in real. With a value of 300 hPa, then
237 a max wind value at 500 hPa will be ignored.
238 use_maxw_level 0=do not use max wind speed level in vertical interpolation inside
239 of the real program, 1 = use level
240 use_trop_level as above, with tropopause level data
241 sfcp_to_sfcp = .false. ! Optional method to compute model's surface pressure when incoming
242 data only has surface pressure and terrain, but not SLP
243 smooth_cg_topo = .false. ! Smooth the outer rows and columns of domain 1's topography w.r.t.
244 the input data
245 use_tavg_for_tsk = .false. ! whether to use diurnally averaged surface temp as skin temp. The
246 diurnally averaged surface temp can be computed using WPS utility
247 avg_tsfc.exe. May use this option when SKINTEMP is not present.
248 aggregate_lu = .false. ! whether to aggregate the grass, shrubs, trees in dominant landuse;
249 default is false.
250 rh2qv_wrt_liquid = .true., ! whether to compute RH with respect to water (true) or ice (false)
251 rh2qv_method = 1, ! which method to use to computer mixing ratio from RH:
252 default is option 1, the old MM5 method; option 2 uses a WMO
253 recommended method (WMO-No. 49, corrigendum, August 2000) -
254 there is a difference between the two methods though small
255 interp_theta = .false. ! If set to .false., it will vertically interpolate temperature
256 instead of potential temperature, which may reduce bias when
257 compared with input data
258 hypsometric_opt = 2, ! = 1: default method
259 = 2: it uses an alternative way (less biased
260 when compared against input data) to compute height in program
261 real and pressure in model.
262 wif_input_opt = 0 ! = 1: option to process the Water Ice Friendly Aerosol input from metgrid for use with mp_physics=28
263 = 2: since V4.4, option to use black carbon aerosol category with mp_physics=28, as well as its radiative effect. Must include
264 file QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat during WPS
265 num_wif_levels = 30 ! number of levels in the Thompson Water Ice Friendly aerosols (mp_physic=28)
266 p_top_requested = 5000 ! p_top (Pa) to use in the model
267 vert_refine_fact = 1 ! vertical refinement factor for ndown, not used for concurrent vertical grid refinement
268 vert_refine_method (max_dom) = 0 ! vertical refinement method
269 0: no vertical refinement
270 1: integer vertical refinement (must set rebalance=1)
271 2: use specified or computed eta levels for vertical refinement (must set rebalance=1)
272 rebalance = 0 ! must be set to =1 if vertical nesting is used.
274 ts_buf_size = 200 ! size of time series buffer
275 max_ts_locs = 5, ! maximum number of time series locations
276 max_ts_level = 15, ! highest model level for time series output
277 tslist_unstagger_winds = .false., ! whether to interpolate 3D u and v to cell centers
279 Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number
280 of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request
281 only the number of levels (with e_vert), and the real program will compute values. There are
282 two methods selected with auto_levels_opt = 1 (old) or 2 (new). The old computation assumes
283 a known first several layers, then generates equi-height spaced levels up to the top of the model.
284 The new method has surface and upper stretching factors (dz_stretch_s and dz_stretch_u) to stretch levels according to
285 log p up to where it reaches the maximum thickness (max_dz) and starting from thickness dzbot.
286 The stretching transitions from dzstretch_s to dzstretch_u by the time the thickness reaches max_dz/2.
288 max_dz = 1000. ! maximum level thickness allowed (m)
289 auto_levels_opt = 1 ! old
290 = 2 ! new default (also set dzstretch_s, dzstretch_u, dzbot, max_dz)
291 dzbot = 50. ! thickness of lowest layer (m) for auto_levels_opt=2
292 dzstretch_s = 1.3 ! surface stretch factor for auto_levels_opt=2
293 dzstretch_u = 1.1 ! upper stretch factor for auto_levels_opt=2
295 eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946,
296 0.922, 0.894, 0.860, 0.817, 0.766,
297 0.707, 0.644, 0.576, 0.507, 0.444,
298 0.380, 0.324, 0.273, 0.228, 0.188,
299 0.152, 0.121, 0.093, 0.069, 0.048,
300 0.029, 0.014, 0.000,
301 eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
302 0.934, 0.909, 0.880, 0.845, 0.807,
303 0.765, 0.719, 0.672, 0.622, 0.571,
304 0.520, 0.468, 0.420, 0.376, 0.335,
305 0.298, 0.263, 0.231, 0.202, 0.175,
306 0.150, 0.127, 0.106, 0.088, 0.070,
307 0.055, 0.040, 0.026, 0.013, 0.000
309 = 0,2, ! this allows vertical nesting in the nest domain
310 Note that with vertical nesting one can only use RRTM and RRTMG radiation physics
312 An example to define vertical nested levels (in program real):
314 e_vert = 35, 45,
315 eta_levels(1:35) = 1., 0.993, 0.983, 0.97, 0.954, 0.934, 0.909, 0.88, 0.8406663, 0.8013327,
316 0.761999, 0.7226653, 0.6525755, 0.5877361, 0.5278192, 0.472514,
317 0.4215262, 0.3745775, 0.3314044, 0.2917579, 0.2554026, 0.2221162,
318 0.1916888, 0.1639222, 0.1386297, 0.1156351, 0.09525016, 0.07733481,
319 0.06158983, 0.04775231, 0.03559115, 0.02490328, 0.0155102, 0.007255059, 0.
320 eta_levels(36:81) = 1.0000, 0.9946, 0.9875, 0.9789, 0.9685, 0.9562, 0.9413, 0.9238, 0.9037, 0.8813, 0.8514,
321 0.8210, 0.7906, 0.7602, 0.7298, 0.6812, 0.6290, 0.5796, 0.5333, 0.4901, 0.4493, 0.4109,
322 0.3746, 0.3412, 0.3098, 0.2802, 0.2524, 0.2267, 0.2028, 0.1803, 0.1593, 0.1398, 0.1219,
323 0.1054, 0.0904, 0.0766, 0.0645, 0.0534, 0.0433, 0.0341, 0.0259, 0.0185, 0.0118, 0.0056, 0.
325 ideal_init_method method to compute alb in idealized cases in start_em
326 = 1, alb from phb (default); = 2, alb from t_init
328 Horizontal interpolation options, coarse grid to fine grid. The default is to use
329 the Smolarkiewicz "SINT" method. However, this is known to break with the
330 implementation inside of WRF for large refinement ratios (such as 15:1). For those
331 extreme (and quite rare occurrences), other schemes are available. For options
332 1, 3, 4, and 12, the FG lateral boundaries use the same horizontal scheme for the
333 lateral BC computations.
334 interp_method_type = 1 ! bi-linear interpolation
335 = 2 ! SINT, (default)
336 = 3 ! nearest neighbor - only to be used for
337 testing purposes
338 = 4 ! overlapping quadratic
339 =12 ! again for testing, uses SINT horizontal
340 interpolation, and same scheme for
341 computation of FG lateral boundaries
343 Variables specifically for the 3d ocean initialization with a single profile. Set
344 the ocean physics option to #2. Specify a number of levels. For each of those levels,
345 provide a depth (m) below the surface. At each depth provide a temperature (K) and
346 a salinity (ppt). The default is not to use the 3d ocean. Even when the 3d ocean is
347 activated, the user must specify a reasonable ocean. Currently, this is the only way
348 available to run the 3d ocean option.
350 &physics
351 sf_ocean_physics = activate ocean model (0=no, 1=1d mixed layer; 2=3D PWP, no bathymetry)
353 &domains
354 ocean_levels = 30,
355 ocean_z ; vertical profile of layer depths for ocean (in meters), e.g.:
356 = 5., 15., 25., 35., 45., 55.,
357 65., 75., 85., 95., 105., 115.,
358 125., 135., 145., 155., 165., 175.,
359 185., 195., 210., 230., 250., 270.,
360 290., 310., 330., 350., 370., 390.
361 ocean_t ; vertical profile of ocean temps, e.g.:
362 = 302.3493, 302.3493, 302.3493, 302.1055, 301.9763, 301.6818,
363 301.2220, 300.7531, 300.1200, 299.4778, 298.7443, 297.9194,
364 297.0883, 296.1443, 295.1941, 294.1979, 293.1558, 292.1136,
365 291.0714, 290.0293, 288.7377, 287.1967, 285.6557, 284.8503,
366 284.0450, 283.4316, 283.0102, 282.5888, 282.1674, 281.7461
367 ocean_s ; vertical profile of salinity, e.g.:
368 = 34.0127, 34.0127, 34.0127, 34.3217, 34.2624, 34.2632,
369 34.3240, 34.3824, 34.3980, 34.4113, 34.4220, 34.4303,
370 34.6173, 34.6409, 34.6535, 34.6550, 34.6565, 34.6527,
371 34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247,
372 34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247
374 Namelist variables for controlling the specified moving nest:
375 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
376 to the ARCHFLAGS. The maximum number of moves, max_moves, is set to 50
377 but can be modified in source code file frame/module_driver_constants.F.
378 num_moves = 0 ! total number of moves
379 move_id(max_moves) = 2,2,2,2, ! a list of nest domain id's, one per move
380 move_interval(max_moves) = 60,120,150,180, ! time in minutes since the start of this domain
381 move_cd_x(max_moves) = 1,1,0,-1, ! the number of parent domain grid cells to move in i direction
382 move_cd_y(max_moves) = 1,0,-1,1, ! the number of parent domain grid cells to move in j direction
383 positive is to move in increasing i and j direction, and
384 negative is to move in decreasing i and j direction.
385 0 means no move. The limitation now is to move only 1 grid cell
386 at each move.
388 Namelist variables for controlling the automatic moving nest:
389 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
390 and -DVORTEX_CENTER to the ARCHFLAGS. This option uses a mid-level vortex following algorithm to
391 determine the nest move. This option is experimental.
392 vortex_interval(max_dom) = 15 ! how often the new vortex position is computed
393 max_vortex_speed(max_dom) = 40 ! used to compute the search radius for the new vortex position
394 corral_dist(max_dom) = 8 ! how many coarse grid cells the moving nest is allowed to get
395 near the mother domain boundary
396 track_level = 50000 ! pressure value in Pa where the vortex is tracked
397 time_to_move(max_dom) = 0. ! time (in minutes) to start the moving nests
399 tile_sz_x = 0, ! number of points in tile x direction
400 tile_sz_y = 0, ! number of points in tile y direction
401 can be determined automatically
402 numtiles = 1, ! number of tiles per patch (alternative to above two items)
403 nproc_x = -1, ! number of processors in x for decomposition
404 nproc_y = -1, ! number of processors in y for decomposition
405 -1: code will do automatic decomposition
406 >1: for both: will be used for decomposition
408 Namelist variables for controlling the adaptive time step option:
409 use_adaptive_time_step = .false. ! T/F use adaptive time stepping.
410 step_to_output_time = .true. ! if adaptive time stepping, T/F modify the
411 time steps so that the exact history time is reached
412 target_cfl(max_dom) = 1.2,1.2 ! vertical and horizontal CFL <= to this value implies
413 no reason to reduce the time step, and to increase it
414 target_hcfl(max_dom) = .84,.84 ! horizontal CFL <= to this value implies
415 max_step_increase_pct(max_dom) = 5,51 ! percentage of previous time step to increase, if the
416 max(vert cfl, horiz cfl) <= target_cfl, then the time
417 will increase by max_step_increase_pct. Use something
418 large for nests (51% suggested)
419 starting_time_step(max_dom) = -1,-1 ! flag = -1 implies use 4*dx (defined in start_em),
420 starting_time_step = 100 means the starting time step
421 for the coarse grid is 100 s
422 max_time_step(max_dom) = -1,-1 ! flag = -1 implies max time step is 8*dx,
423 max_time_step = 100 means that the time step will not
424 exceed 100 s
425 min_time_step(max_dom) = -1,-1 ! flag = -1 implies max time step is 3*dx,
426 min_time_step = 100 means that the time step will not
427 be less than 100 s
428 adaptation_domain = 1 ! default, all fine grid domains adaptive dt driven by coarse-grid
429 2 = Fine grid domain #2 determines the fundamental adaptive dt.
431 &dfi_control
432 dfi_opt = 0 ! which DFI option to use (3 is recommended)
433 0 = no digital filter initialization
434 1 = digital filter launch (DFL)
435 2 = diabatic DFI (DDFI)
436 3 = twice DFI (TDFI)
437 dfi_nfilter = 7 ! digital filter type to use (7 is recommended)
438 0 = uniform
439 1 = Lanczos
440 2 = Hamming
441 3 = Blackman
442 4 = Kaiser
443 5 = Potter
444 6 = Dolph window
445 7 = Dolph
446 8 = recursive high-order
447 dfi_write_filtered_input = .true. ! whether to write wrfinput file with filtered
448 model state before beginning forecast
449 dfi_write_dfi_history = .false. ! whether to write wrfout files during filtering integration
450 dfi_cutoff_seconds = 3600 ! cutoff period, in seconds, for the filter
451 dfi_time_dim = 1000 ! maximum number of time steps for filtering period
452 this value can be larger than necessary
453 dfi_bckstop_year = 2004 ! four-digit year of stop time for backward DFI integration
454 dfi_bckstop_month = 03 ! two-digit month of stop time for backward DFI integration
455 dfi_bckstop_day = 14 ! two-digit day of stop time for backward DFI integration
456 dfi_bckstop_hour = 12 ! two-digit hour of stop time for backward DFI integration
457 dfi_bckstop_minute = 00 ! two-digit minute of stop time for backward DFI integration
458 dfi_bckstop_second = 00 ! two-digit second of stop time for backward DFI integration
459 dfi_fwdstop_year = 2004 ! four-digit year of stop time for forward DFI integration
460 dfi_fwdstop_month = 03 ! two-digit month of stop time for forward DFI integration
461 dfi_fwdstop_day = 13 ! two-digit month of stop time for forward DFI integration
462 dfi_fwdstop_hour = 12 ! two-digit month of stop time for forward DFI integration
463 dfi_fwdstop_minute = 00 ! two-digit month of stop time for forward DFI integration
464 dfi_fwdstop_second = 00 ! two-digit month of stop time for forward DFI integration
465 dfi_savehydmeteors = 0 ! option for radar DA: 0: set hydrometeors to 0 before DFI and
466 let them spin up in DFI; 1: keep them unchanged in DFI
468 &physics
470 Note: even the physics options can be different in different nest domains,
471 caution must be used as what options are sensible to use
473 chem_opt (max_dom) = 0, ! chemistry option - use WRF-Chem
474 mp_physics (max_dom) microphysics option
475 = 0, no microphysics
476 = 1, Kessler scheme
477 = 2, Lin et al. scheme
478 = 3, WSM 3-class simple ice scheme
479 = 4, WSM 5-class scheme
480 = 5, Ferrier (new Eta) microphysics, operational High-Resolution Window version
481 = 6, WSM 6-class graupel scheme
482 = 7, Goddard 4-ice scheme
483 = 8, Thompson scheme
484 = 9, Milbrandt-Yau 2-moment scheme
485 = 10, Morrison (2 moments)
486 = 11, CAM 5.1 microphysics
487 = 13, SBU_YLIN scheme
488 = 14, WDM 5-class scheme
489 = 16, WDM 6-class scheme
490 = 17, NSSL 2-moment 4-ice scheme (steady background CCN)
491 = 18, NSSL 2-moment 4-ice scheme with predicted CCN (better for idealized than real cases)
492 to set a global CCN value, use
493 nssl_cccn = 0.7e9 ; CCN for NSSL scheme (18).
494 Also sets same value to ccn_conc for mp_physics=18
495 = 19, NSSL 1-moment (7 class: qv,qc,qr,qi,qs,qg,qh; predicts graupel density)
496 = 21, NSSL 1-moment, (6-class), very similar to Gilmore et al. 2004
497 Can set intercepts and particle densities in physics namelist, e.g., nssl_cnor
498 For NSSL 1-moment schemes, intercept and particle densities can be set for snow,
499 graupel, hail, and rain. For the 1- and 2-moment schemes, the shape parameters
500 for graupel and hail can be set.
501 nssl_alphah = 0. ! shape parameter for graupel
502 nssl_alphahl = 2. ! shape parameter for hail
503 nssl_cnoh = 4.e5 ! graupel intercept
504 nssl_cnohl = 4.e4 ! hail intercept
505 nssl_cnor = 8.e5 ! rain intercept
506 nssl_cnos = 3.e6 ! snow intercept
507 nssl_rho_qh = 500. ! graupel density
508 nssl_rho_qhl = 900. ! hail density
509 nssl_rho_qs = 100. ! snow density
510 = 24, WSM 7-class scheme (separate hail and graupel categories)
511 = 26, WDM 7-class scheme (separate hail and graupel categories)
512 = 28, aerosol-aware Thompson scheme with water- and ice-friendly aerosol climatology
513 This option has two climatological aerosol input options:
514 use_aero_icbc = .F. : use constant values
515 use_aero_icbc = .T. : use climatological aerosol input from WPS
516 use_rap_aero_icbc = .false. ! Set to .true. to ingest real-time data containing aerosols (new in 4.4)
517 qna_update = 0 ! set to 1 to update time-varying sfc aerosol emission from climatology or
518 real-time data
519 with mp_physics = 28. Use with input file ‘wrfqnainp_d0*’
520 (must set auxinput17_interval and io_form_auxinput17; new in 4.4)
521 wif_fire_emit = .false. ! set to .true. to include biomass burning organic and black carbon
522 aerosols with mp_physics = 28 (new in 4.4)
523 wif_fire_inj = 1 ! (default) vertically distribute biomass burning emissions
524 in mp_physics = 28 (new in 4.4)
525 = 30, HUJI (Hebrew University of Jerusalem, Israel) spectral bin microphysics,
526 fast version
527 = 32, HUJI spectral bin microphysics, full version
528 = 38, Thompson scheme with 2-moment graupel/hail (new in 4.5)
529 = 40, Morrison (2 moments) with consideration of CESM-NCSU RCP4.5 climatological
530 aerosol
531 = 50, P3 1-ice category, 1-moment cloud water
532 = 51, P3 1-ice category plus double-moment cloud water
533 = 52, P3 2-ice categories plus double-moment cloud water
534 = 53, P3 1-ice category, 3-moment ice, plus double-moment cloud water
535 = 55, Jensen-ISHMAEL (Ice-Spheroids Habit Model with Aspect-ratio Evolution) scheme
536 predicting qc, qr, and three ice species with four predictive equations for each ice
537 = 56, NTU multi-moment scheme ! for ntu3m
538 ccnty = 1 ! for aerosol background type, marine, for ntu3m
539 2 ! continental clean (default), for ntu3m
540 3 ! continental average, for ntu3m
541 4 ! continental urban, for ntu3m
543 = 95, Ferrier (old Eta) microphysics
544 = 96, Madwrf
545 = 97, Goddard GCE scheme (also uses gsfcgce_hail, gsfcgce_2ice)
547 For non-zero mp_physics options, to keep Qv .GE. 0, and to set the other moisture
548 fields .LT. a critical value to zero
550 mp_zero_out = 0, ! no action taken, no adjustment to any moist field
551 = 1, ! except for Qv, all other moist arrays are set to zero
552 if they fall below a critical value
553 = 2, ! Qv is .GE. 0, all other moist arrays are set to zero
554 if they fall below a critical value
555 mp_zero_out_thresh = 1.e-8 ! critical value for moist array threshold, below which
556 moist arrays (except for Qv) are set to zero (kg/kg)
558 gsfcgce_hail = 0 ! for running gsfcgce microphysics with graupel
559 = 1 ! for running gsfcgce microphysics with hail
560 default value = 0
561 gsfcgce_2ice = 0 ! for running with snow, ice and graupel/hail
562 = 1 ! for running with only ice and snow
563 = 2 ! for running with only ice and graupel
564 (only used in very extreme situation)
565 default value = 0
566 gsfcgce_hail is ignored if gsfcgce_2ice is set to 1 or 2.
567 hail_opt = 0 ! hail switch for WSM6 and WDM6 : 0 - off, 1 - on
568 morr_rimed_ice = 1 ! hail switch for Morrison schemes (mp_physics=10 and 40: 0 - off,
569 1 - on (new in 4.0)
570 clean_atm_diag = 0 ! if set to =1, turns on clean sky diagnostics (for chem); default is 0=off
571 acc_phy_tend = 0 ! set to =1 to output 16 accumulated physics tendencies for potential temp,
572 water vaopr mixing ratio, and U/V wind components; default is 0=off (new in 4.4)
573 progn (max_dom) = 0 ! switch to use mix-activate scheme (Only for Morrison, WDM6, WDM5,
574 and NSSL_2MOMCCN/NSSL_2MOM
575 ccn_conc = 1.E8 ! CCN concentration, used by WDM schemes
577 no_mp_heating = 0 ! normal
578 = 1 ! turn off latent heating from a microphysics scheme
579 use_mp_re = 1 ! whether to use effective radii computed in mp schemes in RRTMG
580 0: do not use; 1: use effective radii
581 (The mp schemes that compute effective radii are 3,4,6,7,8,10,14,16,17-21,24,26,28,50-53,55)
583 force_read_thompson = .false. ! whether to read tables for mp_physics = 8,28
584 write_thompson_tables = .true. ! whether to read or compute tables for mp_phyiscs = 8,28
585 write_thompson_mp38table = .false. ! whether to read table (qr_acr_qg_mp38V1.dat) for mp_physics = 38
587 ra_lw_physics (max_dom) longwave radiation option
588 = 0, no longwave radiation
589 = 1, rrtm scheme
590 (Default values for GHG in V4.2: co2vmr=(280. + 90.*exp(0.02*(yr-2000)))*1.e-6
591 n2ovmr=319.e-9, ch4vmr=1774.e-9
592 Values used in previous versions: co2vmr=330.e-6, n2ovmr=0., ch4vmr=0.)
593 = 3, cam scheme
594 also requires levsiz, paerlev, cam_abs_dim1/2 (see below)
595 = 4, rrtmg scheme
596 (Default values for GHG in V4.2: co2vmr=(280. + 90.*exp(0.02*(yr-2000)))*1.e-6
597 n2ovmr=319.e-9, ch4vmr=1774.e-9,
598 cfc11=0.251e-9, cfc12=0.538e-9,
599 = 14, rrtmg-k scheme from KIAPS
600 = 24, fast rrtmg scheme for GPU and MIC (since 3.7)
601 (Default values for GHG in V4.2: co2vmr=(280. + 90.*exp(0.02*(yr-2000)))*1.e-6
602 n2ovmr=319.e-9, ch4vmr=1774.e-9,
603 cfc11=0.251e-9, cfc12=0.538e-9
604 = 5, Goddard longwave scheme
605 = 7, FLG (UCLA) scheme
606 = 31, Earth Held-Suarez forcing
607 = 99, GFDL (Eta) longwave (semi-supported)
608 also must use co2tf = 1 for ARW
610 ra_sw_physics (max_dom) shortwave radiation option
611 = 0, no shortwave radiation
612 = 1, Dudhia scheme
613 = 2, Goddard short wave
614 = 3, cam scheme
615 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
616 = 4, rrtmg scheme
617 = 14, rrtmg-k scheme from KIAPS
618 = 24, fast rrtmg scheme for GPU and MIC (since 3.7)
619 = 5, Goddard shortwave scheme
620 = 7, FLG (UCLA) scheme
621 = 99, GFDL (Eta) longwave (semi-supported)
622 also must use co2tf = 1 for ARW
624 radt (max_dom) = 30, ! minutes between radiation physics calls
625 recommend 1 min per km of dx (e.g. 10 for 10 km);
626 use the same value for all nests.
627 cldovrlp = 2, ! cloud overlapping option for RRTMG only. 1=random, 2=maximum-random (default),
628 3=maximum, 4=exponential, 5=exponential-random
629 idcor = 0, ! decorrelation length flag for cldovrlp=4 or 5
630 0 = constant decorrelation length, 2500 m;
631 1 = latitude-varying decorrelation length
632 ra_sw_eclipse = 0, ! eclipse effect on shortwave radiation. 0: off, 1: on.
633 Applies to RRTMG, Goddard, old Goddard and Dudhia schemes
634 ghg_input = 1, ! Option to read CAMtr_volume_mixing_ratio files of green house gas values,
635 as of v4.4, the default is SSP 2 with RCP 4.5 => SSP245
636 Used for CAM LW and SW, RRTM, RRTMG LW and SW, RRTMG_fast LW and SW
637 0 = do not read in the annual data
638 1 = read in time dependent data for CO2, N2O, CH4, CFC11, CFC12
640 co2tf CO2 transmission function flag only for GFDL radiation
641 = 0, ! read CO2 function data from pre-generated file
642 = 1, ! generate CO2 functions internally in the forecast
644 ra_call_offset radiation call offset
645 = 0 ! (=0, no offset; =-1, old offset)
646 swint_opt Interpolation of short-wave radiation based on the updated solar zenith angle
647 between SW call
648 = 0, ! no interpolation
649 = 1, ! use interpolation
650 = 2, ! Calls the Fast All-Sky Radiation Model for Solar applications (FARMS)
651 ! every model time step.
652 couple_farms = .false. ! True) uses FARMS SW radiation to drive the LSM, False) Uses the SW radiation from rad_sw_physics
653 cam_abs_freq_s = 21600 ! default CAM clearsky longwave absorption calculation frequency
654 (recommended minimum value to speed scheme up)
655 levsiz = 59 ! for CAM radiation input ozone levels, set automatically
656 paerlev = 29 ! for CAM radiation input aerosol levels, set automatically
657 cam_abs_dim1 = 4 ! for CAM absorption save array, set automatically
658 cam_abs_dim2 = value of e_vert for CAM 2nd absorption save array, set automatically
660 o3input = ozone input option for radiation (currently rrtmg only)
661 = 0, ! using profile inside the code
662 = 2, ! using CAM ozone data (ozone.formatted)
663 aer_opt = aerosol input option for radiation (currently rrtmg only)
664 = 0, ! none
665 = 1, ! using Tegen (1997) data,
666 = 2, ! using J. A. Ruiz-Arias method (see other aer_* options)
667 = 3, ! using G. Thompson's water/ice friendly climatological aerosol
668 alevsiz = 12 ! for Tegen aerosol input levels, set automatically
669 no_src_types = 6 ! for Tegen aerosols: organic and black carbon, sea salt, sulfalte, dust,
670 and stratospheric aerosol (volcanic ashes - currently 0), set automatically
672 The following aerosol options allow RRTMG and new Goddard radiation schemes to see it, but the aerosols are
673 constant during the model integration.
674 aer_aod550_opt (max_dom) = [1,2]
675 1 = input constant value for AOD at 550 nm from namelist.
676 In this case, the value is read from aer_aod550_val;
677 2 = input value from auxiliary input 15. It is a time-varying 2D grid in netcdf
678 wrf-compatible format. The default is aer_aod550_opt=1 and aer_aod550_val=0.12
679 aer_aod550_val (max_dom) = 0.12
680 aer_angexp_opt (max_dom) = [1,2,3] :
681 1 = input constant value for Angstrom exponent from namelist. In this case,
682 the value is read from aer_angexp_val;
683 2 = input value from auxiliary input 15, as in aer_aod550_opt;
684 3 = Angstrom exponent value estimated from the aerosol type defined in aer_type, and modulated
685 with the RH in WRF. Default operation is aer_angexp_opt = 1, and aer_angexp_val=1.3.
686 aer_angexp_val (max_dom) = 1.3
687 aer_ssa_opt (max_dom) = [1,2,3] similar to aer_angexp_opt.
688 aer_ssa_val (max_dom) = 0.85
689 aer_asy_opt (max_dom) = [1,2,3] similar to aer_angexp_opt.
690 aer_asy_val (max_dom) = 0.9
691 aer_type (max_dom) = [1,2,3] : 1 for rural, 2 is urban and 3 is maritime.
693 sf_sfclay_physics (max_dom) surface-layer option (old bl_sfclay_physics option)
694 = 0, no surface-layer
695 = 1, Revised MM5 Monin-Obukhov scheme (Jimenez, renamed in v3.6)
696 = 2, Monin-Obukhov (Janjic) scheme
697 = 3, NCEP Global Forecast System scheme
698 = 4, QNSE surface layer
699 = 5, MYNN surface layer
700 = 7, Pleim-Xiu surface layer
701 = 10, TEMF surface layer
702 = 91, Old MM5 scheme (previously option 1)
704 sf_surface_physics (max_dom) land-surface option (old bl_surface_physics option)
705 = 0, no surface temp prediction
706 = 1, thermal diffusion scheme
707 = 2, Unified Noah land-surface model
708 = 3, RUC land-surface model
709 = 4, Noah-MP land-surface model (see additional &noah_mp namelist)
710 = 5, Community Land Model version 4 (CLM4), adapted from CAM
711 = 6, Community Terrestrial Systems Model (CTSM). See doc/README.CTSM for info on running this option.
712 = 7, Pleim-Xiu LSM
713 = 8, Simplified Simple Biosphere Model (SSiB)
714 - can be used with Dudhia/RRTM, CAM or RRTMG radiation options
716 sf_urban_physics(max_dom) = 0, ! activate urban canopy model (in Noah and Noah-MP LSMs only)
717 = 0: no
718 = 1: Single-layer, UCM
719 = 2: Multi-layer, Building Environment Parameterization (BEP) scheme
720 (works only with YSU, MYJ and BouLac PBL)
721 = 3: Multi-layer, Building Environment Model (BEM) scheme
722 (works only with YSU, MYJ and BouLac PBL)
724 num_urban_ndm = 1! (= 2 if BEP or BEM active) maximum number of street dimensions (ndm in BEP or BEM header)
725 num_urban_ng = 1! (= 10 if BEP or BEM active) number of grid levels in the ground (ng_u in BEP or BEM header)
726 num_urban_nwr = 1! (= 10 if BEP or BEM active) number of grid levels in the walls or roof (nwr_u in BEP or BEM header)
727 num_urban_nz = 1! (= 18 if BEP or BEM active) maximum number of vertical levels in the urban grid (nz_um in BEP or BEM header)
728 num_urban_ngb = 1! (= 10 if BEM active) number of grid levels in the ground below building (ngb_u in BEM header)
729 num_urban_nf = 1! (= 10 if BEM active) number of grid levels in the floors (nf_u in BEM header)
730 num_urban_nbui = 1! (= 15 if BEM active) maximum number of types of buildings in an urban class (nbui_max in BEM header)
732 sf_surf_irr_scheme (max_dom) = 0: no irrigation scheme
733 = 1: channel scheme (in Noah only)
734 = 2: drip scheme
735 = 3: sprinkler scheme
737 irr_daily_amount (max_dom) = daily equivalent irrigation water amount in mm (float, e.g. 5.7)
738 irr_start_hour (max_dom) = UTC start hour for irrigation (int, e.g. 5)
739 irr_num_hours (max_dom) = number of consecutive hours for irrigation (int, e.g. 3)
740 irr_start_julianday (max_dom) = Julian start day of irrigation (int, e.g. 135)
741 irr_end_julianday (max_dom) = Julian end day of irrigation (int, e.g. 255)
742 irr_ph (max_dom) = 0: all the static IRRIGATION field is activated
743 = 1: pseudo random activation field: ensures repeatability across compilers
744 = 2: random activation field: results might depend on the compiler fortran RANDOM_SEED
745 irr_freq (max_dom) = irrigation frequency in days (int, e.g. 3)
747 bl_pbl_physics (max_dom) boundary-layer option
748 = 0, no boundary-layer
749 = 1, YSU scheme
750 = 2, Mellor-Yamada-Janjic TKE scheme
751 = 3, Hybrid EDMF GFS scheme
752 = 4, Eddy-diffusivity Mass Flux, Quasi-Normal Scale Elimination PBL
753 = 5, MYNN 2.5 level TKE scheme, works with
754 sf_sfclay_physics=1 or 2 as well as 5
755 (option 6 removed in 4.5, replaced by bl_mynn_closure options)
756 bl_mynn_closure = 2.5: level 2.5
757 2.6: level 2.6
758 3.0: level 3.0
759 = 7, ACM2 (Pleim) PBL
760 = 8, Bougeault and Lacarrere (BouLac) PBL
761 = 9, UW boundary layer scheme from CAM5 (CESM 1_0_1)
762 = 10, TEMF (Total Energy Mass Flux) scheme
763 sf_sfclay_physics=10
764 = 11, Shin-Hong 'scale-aware' PBL scheme
765 = 12, Grenier-Bretherton-McCaa scheme
766 = 16, TKE+TKE dissipation rate (epsilon) scheme
767 works with surface layer options 1,91,2,5.
768 = 17, TKE+TKE dissipation rate+TPE (temparature variance) scheme
769 = 99, MRF scheme
771 bldt (max_dom) = 0, ! minutes between boundary-layer physics calls
773 grav_settling (max_dom) gravitational settling of fog/cloud droplets (Now works for any PBL scheme)
774 = 0, No settling of cloud droplets
775 = 1, Settling from Dyunkerke 1991 (in atmos and at surface)
776 = 2, Fogdes (vegetation & wind speed dependent; Katata et al. 2008) at surface
777 and Dyunkerke in the atmos.
778 ysu_topdown_pblmix = 0, ! whether to turn on top-down, radiation-driven mixing (1=yes)
779 mfshconv (max_dom) = 1, ! whether to turn on new day-time EDMF QNSE (0=no)
780 topo_wind (max_dom) = 0, ! turn off,
781 = 1, turn on topographic surface wind correction from Jimenez
782 (YSU PBL only, and require extra input from geogrid)
783 = 2, turn on topographic surface wind correction from Mass (YSU PBL only)
785 tke_budget (max_dom) = 0, ! default off; = 1 adds MYNN tke budget terms to output (new in 4.5)
786 (replacing bl_mynn_tkebudget in prior versions)
787 bl_mynn_tkeadvect (max_dom) = .false., ! default off; = .true. do MYNN tke advection
788 icloud_bl option to couple the subgrid-scale clouds from the PBL scheme (MYNN only)
789 to radiation schemes
790 0: no coupling; 1: activate coupling to radiation (default)
791 bl_mynn_cloudmix (max_dom) = 0 ! default off; =1 activates mixing of qc and qi in MYNN
792 0: no mixing of qc & qi; 1: mixing activated (default).
793 Note qnc and qni are mixed when scalar_pblmix =1.
794 bl_mynn_mixlength option to change mixing length formulation in MYNN
795 0: original as in Nakanishi and Niino 2009,
796 1: RAP/HRRR (including BouLac in free atmosphere),
797 2: experimental (default; includes cloud-specific mixing length and a scale-aware mixing
798 length, following Ito et al. 2015, BLM). Option 2 has been well tested with
799 the edmf options.
800 bl_mynn_cloudpdf option to switch to different cloud PDFs to represent subgrid clouds
801 0: original (Sommeria and Deardorf 1977);
802 1: Kuwano et al 2010, similar to option 0, but uses resolved scale gradients
803 as opposed to higher order moments ;
804 2: from Chaboureau and Bechtold (2002, JAS, with mods, default)
805 bl_mynn_edmf (max_dom) = 1 ! activate mass-flux option in MYNN, 0: mass-flux option off
806 Related (hidden) options:
807 bl_mynn_edmf_mom (max_dom) = 1 ! - activates momentum transport in MYNN mass-flux scheme
808 (assuming bl_mynn_edmf > 0); 0=off (default)
809 0 no momentum transport; 1: momentum transport activated (default)
810 bl_mynn_edmf_tke (max_dom) = 0, ! default; = 1 activates TKE transport in MYNN mass-flux scheme
811 (assuming bl_mynn_edmf > 0)
812 0 no TKE transport (default);1: activate TKE transport
813 bl_mynn_output (max_dom) = 0, ! do not output extra arrays
814 1 allocate and output extra 3D arrays from MYNN PBL
815 bl_mynn_mixscalars (max_dom) = 0 ! off, 1: activate mixing of scalars
816 bl_mynn_mixqt (max_dom) = 0 ! mixing moisture species separately, 1: mixing total water
818 scalar_pblmix (max_dom) = 1 ! mix scalar fields consistent with PBL option (exch_h)
819 tracer_pblmix (max_dom) = 1 ! mix tracer fields consistent with PBL option (exch_h)
820 shinhong_tke_diag (max_dom) = 0 ! diagnostic TKE and mixing length from Shin-Hong PBL
821 opt_thcnd option to treat thermal conductivity in Noah LSM
822 = 1, original (default)
823 = 2, McCumber and Pielke for silt loam and sandy loam
824 sf_surface_mosaic option to mosaic landuse categories for Noah LSM
825 = 0 ! default; use dominant category only
826 = 1 ! use mosaic landuse categories
827 mosaic_cat = 3 ! number of mosaic landuse categories in a grid cell
828 mosaic_lu = 0, ! default ; set to = 1 to use mosaic landuse categories in RUC
829 mosaic_soil = 0, ! default ; set to = 1 to use mosaic soil categories in RUC
830 flag_sm_adj = 0, ! default ; set to = 1 to adjust soil moisture for RUC in real.exe
831 (e.g. in case RUC is initialized from GFS LSM data)
833 cu_physics (max_dom) cumulus option
834 = 0, no cumulus
835 = 1, Kain-Fritsch (new Eta) scheme
836 = 2, Betts-Miller-Janjic scheme
837 = 3, Grell-Freitas ensemble scheme
838 = 4, Scale-aware GFS Simplified Arakawa-Schubert (SAS) scheme
839 = 5, Grell 3D ensemble scheme
840 = 6, Modifed Tiedtke scheme
841 = 7, Zhang-McFarlane scheme from CAM5 (CESM 1_0_1)
842 = 10, Modified Kain-Fritsch scheme with trigger function based on PDFs
843 = 11, Multi-scale Kain-Fritsch scheme
844 = 14, KIM Simplified Arakawa-Schubert scheme (KSAS) across gray-zone resolutions
845 = 16, A newer Tiedtke scheme
846 = 94, 2015 GFS Simplified Arakawa-Schubert scheme (from HWRF, use with caution)
847 = 95, Previous GFS Simplified Arakawa-Schubert scheme (from HWRF, use with caution)
848 = 96, Previous NEW GFS simplified Arakawa-Schubert scheme from YSU
849 = 93, Grell-Devenyi ensemble scheme
850 = 99, previous Kain-Fritsch scheme
852 shcu_physics (max_dom) independent shallow cumulus option (not tied to deep convection)
853 = 0, no independent shallow cumulus
854 = 1, Grell 3D ensemble scheme (use with cu_physics=93 or 5) (PLACEHOLDER: SWITCH NOT YET IMPLEMENTED--use ishallow)
855 = 2, Park and Bretherton shallow cumulus from CAM5 (CESM 1_0_1)
856 = 3, GRIMS shallow cumulus from YSU group
857 = 4, NSAS shallow cululus scheme. Must be used with KIAPS SAS cumulus scheme (cu_physics = 14)
858 = 5, Deng cumulus scheme (including both shallow and deep convections) from PSU and WRF-Solar.
859 However the deep part isn't very active. Not recommended to use alone for deep convection
860 case. Could work well for grid sizes 3-9 km.
861 This scheme only works with MYJ PBL scheme, and should not be combined with icloud=3
862 This scheme can also work with MYNN PBL scheme, but one should turn off EDMF (bl_mynn_edmf=0)
864 ishallow = 0, ! = 1 turns on shallow convection, used with Grell 3D ensemble schemes (cu_physics = 3 or 5)
865 clos_choice = 0, ! closure choice (place holder only)
866 cu_diag (max_dom) = 0, ! additional t-averaged stuff for cu physics (cu_phys = 3, 5, 10, and 93 only)
867 kf_edrates (max_dom) = 0, ! Add entrainment/detrainment rates and convective timescale output variables for KF-based
868 cumulus schemes (cu_phys = 1, 11 and 99 only)
869 = 0, ! no output; = 1, additional output
870 convtrans_avglen_m = 30, ! averaging time for variables used by convective transport (call cu_phys options) and radiation routines (only cu_phys=3,5 and 93) (minutes)
871 cu_rad_feedback (max_dom) = .false. ! sub-grid cloud effect to the optical depth in radiation
872 currently it works only for GF, G3, GD and KF scheme
873 One also needs to set cu_diag = 1 for GF, G3, KF-CuP, and GD schemes
874 bmj_rad_feedback (max_dom) = .false. ! BMJ convective precipitation-derived sub-grid cloud effect to radiation
875 cudt (max_dom) = 0, ! minutes between cumulus physics calls
876 kfeta_trigger KF trigger option (cu_physics=1 only):
877 = 1, default option
878 = 2, moisture-advection based trigger (Ma and Tan [2009])
879 = 3, RH-dependent additional perturbation to option 1 (JMA)
880 cugd_avedx ; number of grid boxes over which subsidence is spread.
881 = 1, default, for large grid distances
882 = 3, for small grid distances (DX < 5 km)
883 nsas_dx_factor = 0, ! default option
884 = 1, NSAS grid-distance dependent option
885 ntiedtke_dx_opt = 0, default option
886 = 1, new Tiedtke grid-distance dependent option (scale-aware)
887 when grid size falls below 15 km
888 For KF-CuP scheme: recommended to use with cu_rad_feedback
889 shallowcu_forced_ra(max_dom) radiative impact of shallow Cu by a prescribed maximum cloud fraction
890 = .false., option off, default
891 = .true., radiative impact of shallow cu with a cloud fraction value of 0.36
892 numbins(max_dom) number of perturbations for potential temperature and mixing ratio in the CuP PDF,
893 should be an odd number (21 is a recommended value)
894 thBinSize(max_dom) bin size of potential temperature perturbation increment (0.01 K)
895 rBinSize(max_dom) bin size of mixing ratio perturbation increment (1.0e-4 kg/kg)
896 minDeepFreq(max_dom) minimum frequency required before deep convection is allowed (0.333)
897 minShallowFreq(max_dom) minimum frequency required before shallow convection is allowed (1.0e-2)
898 shcu_aerosols_opt(max_dom) whether aerosols in shcu: 0=none, 2=prognostic (run with WRF-Chem),
900 For MSKF scheme:
901 aercu_opt option to control aerosol interaction in MSKF and Morrison microphysics (mp_physics=40)
902 = 0: no aerosol interaction
903 = 1: aerosol interaction with MSKF only
904 = 2: aerosol interaction with both MSKF and Morrison
905 aercu_fct factor to multiply with aerosol amount
906 = 1. ! default value
907 no_src_types_cu = 1 ! number of aerosol species in global aerosol data: 10 for CESM input,
908 set automatically
909 alevsiz_cu = 1 ! number of levels in global aerosol data: 30 for CESM input,
910 set automatically
911 isfflx = 1, ! heat and moisture fluxes from the surface
912 (only works for sf_sfclay_physics = 1,5,7,11)
913 1 = with fluxes from the surface
914 0 = no flux from the surface
915 with bl_pbl_physics=0 this uses tke_drag_coefficient
916 and tke_heat_flux in vertical diffusion
917 2 = use drag from sf_sfclay_physics and heat flux from
918 tke_heat_flux with bl_pbl_physics=0
919 ideal_xland = 1, ! sets XLAND (1=land,2=water) for ideal cases with no input land-use
920 run-time switch for wrf.exe physics_init (default 1 as before)
921 ifsnow = 1, ! snow-cover effects
922 (only works for sf_surface_physics = 1)
923 1 = with snow-cover effect
924 0 = without snow-cover effect
925 icloud = 1, ! cloud effect to the optical depth in radiation
926 (only works for ra_sw_physics = 1,4 and ra_lw_physics = 1,4)
927 Since 3.6, this also controls the cloud fraction options
928 1 = with cloud effect, and use cloud fraction option 1
929 (Xu-Randall method)
930 0 = without cloud effect
931 2 = with cloud effect, and use cloud fraction option 2 (0/1 based
932 on threshold
933 3 = with cloud effect, and use cloud fraction option 3, based on
934 Sundqvist et al. (1989) (since 3.7)
935 insert_init_cloud = .false., ! Default
936 = .true., ! For use with cold start (zero cloud fields), this option will
937 use the icloud=3 scheme to produce initial cloud water and cloud
938 ice fields as well as ensure initial water vapor field matches
939 saturation, which retains clouds beyond the first few timesteps.
940 swrad_scat = 1. ! scattering tuning parameter (default 1. is 1.e-5 m2/kg)
941 (works for ra_sw_physics = 1 option only)
942 surface_input_source = 3, ! where landuse and soil category data come from:
943 1 = WPS/geogrid but with dominant categories recomputed
944 2 = GRIB data from another model (only possible
945 (VEGCAT/SOILCAT are in met_em files from WPS)
946 3 = use dominant land and soil categories from WPS/geogrid (default since 3.8)
948 num_soil_layers = 5, ! number of soil layers in land surface model
949 = 5: thermal diffusion scheme
950 = 4: Noah landsurface model
951 = 6 or 9: RUC landsurface model
952 = 10: CLM4 landsurface model
953 = 2: Pleim-Xu landsurface model
954 = 3: SSiB landsurface model
955 num_land_cat = 21, ! number of land categories in input data.
956 24 - for USGS (default); 20 for MODIS
957 28 - for USGS if including lake category
958 21 - for MODIS if including lake category (default since 3.8)
959 40 - for NCLD
960 use_wudapt_lcz = 0, ! Option to use WUDAPT LCZ urban landuse categories, along with standard
961 urban classes (31-33).
962 = 0: use the traditional 33 landuse classes (31-33 for urban)
963 = 1: use the WUDAPT LCZ urban landuse categories
964 * Note: If the number of urban category in the input files is
965 inconsistent with the namelist option, error messages will occur.
966 The method to create the LCZ data is described here: http://www.wudapt.org/
967 num_soil_cat = 16, ! number of soil categories in input data
969 pxlsm_smois_init(max_dom) = 1 ! PXLSM Soil moisture initialization option
970 0 - From analysis, 1 - From moisture availability
971 or SLMO in LANDUSE.TBL
972 pxlsm_modis_veg = 1, ! PX LSM LAI and VEGFRA
973 0 - Old PX method that uses PX landuse look-up table
974 1 - Use VEGFRA and LAI in wrflowinp_d0* file
975 Note: Values used are called VEGF_PX and LAI_PX in output.
977 maxiens = 1, ! Grell-Devenyi only
978 maxens = 3, ! G-D only
979 maxens2 = 3, ! G-D only
980 maxens3 = 16 ! G-D only
981 ensdim = 144 ! G-D only
982 These are recommended numbers. If you would like to use
983 any other number, consult the code, know what you are doing.
984 seaice_threshold = 100. ! tsk < seaice_threshold, if water point and 5-layer slab
985 scheme, set to land point and permanent ice; if water point
986 and Noah scheme, set to land point, permanent ice, set temps
987 from 2 m to surface, and set smois and sh2o. The default value has changed
988 from 271 to 100 K in v3.5.1 to avoid mixed-up use with fractional seaice input
989 Used by land model option 1,2,3,4 and 8
990 sst_update = 0 ! time-varying sea-surface temp (0=no, 1=yes). If selected real
991 puts SST, XICE, ALBEDO and VEGFRA in wrflowinp_d01 file, and wrf updates
992 these from it at same interval as boundary file. Also requires
993 namelists in &time_control: auxinput4_interval, auxinput4_end_h,
994 auxinput4_inname = "wrflowinp_d<domain>", and io_form_auxinput4
995 usemonalb = .false. ! use monthly albedo map instead of table value
996 (recommended for sst_update=1)
997 rdmaxalb = .true. ! use snow albedo from geogrid; false means using values from table
998 rdlai2d = .false. ! use LAI from input; false means using values from table
999 if sst_update=1, LAI will also be in wrflowinp file
1000 dust_emis = 0 ! Enable (0=no, 1=yes) surface dust emission scheme to enter mp_physics=28 QNIFA (ice-friendly aerosol variable)
1001 erosion_dim = 3 ! In conjunction with dust_emis=1, this value can only be set equal to 3 (erodibility information)
1002 bucket_mm = -1. ! bucket reset value for water accumulations (value in mm, -1.=inactive)
1003 bucket_J = -1. ! bucket reset value for energy accumulations (value in J, -1.=inactive)
1004 tmn_update = 0 ! update deep soil temperature (1, yes; 0, no)
1005 lagday = 150 ! days over which tmn is computed using skin temperature
1006 sst_skin = 0 ! calculate skin SST
1007 slope_rad (max_dom) = 0 ! slope effects for solar radiation (1=on, 0=off)
1008 topo_shading (max_dom) = 0 ! neighboring-point shadow effects for solar radiation (1=on, 0=off)
1009 shadlen = 25000. ! max shadow length in meters for topo_shading=1
1010 sf_ocean_physics = 0 ! activate ocean model (0=no, 1=1d mixed layer; 2=3D PWP, no bathymetry)
1011 oml_hml0 = 50 ! oml model can be initialized with a constant depth everywhere (m)
1012 < 0, oml is initialized with real-time ocean mixed depth
1013 = 0, oml is initialized with climatological ocean mixed depth
1014 oml_gamma = 0.14 ! oml deep water lapse rate (K m-1)
1015 oml_relaxation_time = 0. ! Relaxation time (in second) of mixed layer ocean model back to original values
1016 (an example value is 259200 sec. (3 days))
1017 omdt = 1. ! 3D PWP time step (min). It can be set to be the same as WRF time step
1018 in corresponding nested grids, but omdt should be no less than 1.0 minute.
1019 ocean_levels = 30 ! number of vertical levels in 3DPWP. Note that the depth of each ocean
1020 model layers is specified in OM_DEPTH in wrfinput_d01
1021 traj_opt = 0 ! Forward trajectory calculation (Lee and Chen 2013)
1022 num_traj = 1000 ! number of trajectories to be released
1023 isftcflx = 0 ! alternative Ck, Cd formulation for tropical storm application
1024 sf_sfclay=1 and 11
1025 0=default
1026 1=Donelan Cd + const z0q
1027 2=Donelan Cd + Garratt
1028 sf_sfclay=5
1029 (default) =0: z0, zt, and zq from COARE3.0 (Fairall et al 2003)
1030 =1: z0 from Davis et al (2008), zt & zq from COARE3.0
1031 =2: z0 from Davis et al (2008), zt & zq from Garratt (1992)
1032 fractional_seaice = 0 ! treat sea-ice as fractional field (1) or ice/no-ice flag (0)
1033 works for sf_sfclay_physics=1,2,3,4,5,7,and 91
1034 If fractional_seaice = 1, also set seaice_threshold = 0.
1035 seaice_albedo_opt = 0 ! option to set albedo over sea ice
1036 ! 0 = seaice albedo is a constant value from namelist option seaice_albedo_default
1037 ! 1 = seaice albedo is f(Tair,Tskin,Snow) following Mills (2011) for Arctic Ocean
1038 ! 2 = seaice albedo read in from input variable ALBSI
1039 seaice_albedo_default = 0.65 ! default value of seaice albedo for seaice_albedo_opt=0
1040 seaice_snowdepth_opt = 0 ! method for treating snow depth on sea ice
1041 ! 0 = snow depth on sea ice is bounded by seaice_snowdepth_min and seaice_snowdepth_max
1042 ! 1 = snow depth on sea ice read in from input array SNOWSI (bounded by
1043 ; seaice_snowdepth_min and seaice_snodepth_max)
1044 seaice_snowdepth_max = 1.E10 ! maximum allowed accumulation of snow (m) on sea ice
1045 seaice_snowdepth_min = 0.001 ! minimum snow depth (m) on sea ice
1046 seaice_thickness_opt = 0 ! option for treating seaice thickness
1047 ! 0 = seaice thickness is uniform value taken from namelist variable seaice_thickness_default
1048 ! 1 = seaice_thickness is read in from input variable ICEDEPTH
1049 seaice_thickness_default = 3.0 ! default value of seaice thickness for seaice_thickness_opt=0
1050 tice2tsk_if2cold = .false. ! set Tice to Tsk to avoid unrealistically low sea ice temperatures
1051 iz0tlnd = 0 ! thermal roughness length for sfclay (0 = old, 1 = veg dependent Chen-Zhang Czil,
1052 2 = Zilitinkevitch (czil=0.1))
1053 for mynn sfc (0=Zilitinkevitch (def),1=Chen-Zhang,2=mod Yang,3=const zt)
1054 mp_tend_lim = 10., ! limit on temp tendency from mp latent heating from radar data assimilation
1055 prec_acc_dt (max_dom) = 0., ! number of minutes in precipitation bucket - will add three
1056 new 2d output fields: prec_acc_c, prec_acc_nc and snow_acc_nc
1057 topo_wind (max_dom) = 0, ! 1 = improve effect of topography for surface winds.
1058 ua_phys = .false. ! Option to activate UA Noah changes: a different snow-cover physics in
1059 Noah, aimed particularly toward improving treatment of snow as it relates
1060 to the vegetation canopy. Also uses new columns added in VEGPARM.TBL
1061 do_radar_ref = 0, ! 1 = allows radar reflectivity to be computed using mp-scheme-specific
1062 parameters. Currently works for mp_physics = 2,4,6,7,8,10,14,16,24,26,28
1063 hailcast_opt (max_dom) = 0, ! 1 = 1-D hail growth model which predicts 1st-5th rank-ordered hail diameters, mean hail
1064 diameter and standard deviation of hail diameter. (Adams-Selin and Ziegler, MWR Dec 2016.)
1065 haildt (max_dom) = 0., ! seconds between WRF-HAILCAST calls (s)
1067 Namelist variables for lake module:
1069 sf_lake_physics(max_dom) = 1, ! lake model on/off
1070 lakedepth_default(max_dom) = 50, ! default lake depth (If there is no lake_depth information in the input data, then lake depth
1071 is assumed to be 50m)
1072 lake_min_elev(max_dom) = 5, ! minimum elevation of lakes. May be used to determine whether a water point is a lake in the absence of lake
1073 category. If the landuse type includes 'lake' (i.e. Modis_lake and USGS_LAKE), this variable is of no effects.
1074 use_lakedepth (max_dom) = 1, ! option to use lake depth data. Lake depth data is available from 3.6 geogrid program. If one didn't process
1075 the lake depth data, but this switch is set to 1, the program will stop and tell one to go back to geogrid
1076 program.
1077 = 0, do not use lake depth data.
1078 shalwater_z0 = 0, ! Shallow water roughness scheme on (1) or off (0). Only compatible with sf_sfclay_physics = 1.
1079 Bathymetry dataset from GEBCO Compilation Group. Please acknowledge the following in presentations and
1080 publications: GEBCO Compilation Group (2021) GEBCO 2021 Grid (doi:10.5285/c6612cbe-50b3-0cff-e053-6c86abc09f8f)
1081 shalwater_depth = 0.0, ! Set constant depth [m] (must be positive) for shallow water roughness scheme if no bathymetry data availabile.
1082 The scheme is intended for depths between 10.0 and 100.0 m. Any depths outside of this range will be
1083 rounded to the nearest limit value.
1084 lightning_option (max_dom) Lightning parameterization option to allow flash rate prediction without chemistry
1085 = 0 ! off
1086 = 1 ! PR92 based on maximum w, redistributes flashes within dBZ > 20 (for convection resolved runs; must also use
1087 do_radar_ref = 1, and mp_physics = 2,4,6,7,8,10,14, or 16)
1088 = 2 ! PR92 based on 20 dBZ top, redistributes flashes within dBZ > 20 (for convection resolved runs; must also use
1089 do_radar_ref = 1, and mp_physics = 2,4,6,7,8,10,14, or 16)
1090 = 3 ! Predicting the potential for lightning activity (based on Yair et al, 2010, J. Geophys. Res., 115, D04205, doi:10.1029/2008JD010868)
1091 = 11 ! PR92 based on level of neutral buoyancy from convective parameterization (for scales
1092 where a CPS is used, intended for use at 10 < dx < 50 km; must also use cu_physics = 5 or 93)
1093 lightning_dt (max_dom) = 0. ! time interval (seconds) for calling lightning parameterization. Default uses model time step
1094 lightning_start_seconds (max_dom) = 0. ! Start time for calling lightning parameterization. Recommends at least 10 minutes for spin-up.
1095 flashrate_factor (max_dom) = 1.0 ! Factor to adjust the predicted number of flashes. Recommends 1.0 for lightning_option = 11
1096 between dx=10 and 50 km. Manual tuning recommended for all other options independently
1097 for each nest.
1098 cellcount_method (max_dom) Method for counting storm cells. Used by CRM options (lightning_options=1,2).
1099 = 0, ! model determines method used
1100 = 1, ! tile-wide, appropriate for large domains
1101 = 2, ! domain-wide, appropriate for sing-storm domains
1102 cldtop_adjustment (max_dom) = 0. ! Adjustment from LNB in km. Used by lightning_option=11. Default is 0, but recommends 2 km
1103 iccg_method (max_dom) IC:CG partitioning method (IC: intra-cloud; CG: cloud-to-ground)
1104 = 0 ! Default method depending on lightning option, currently all options use iccg_method=2 by default
1105 = 1 ! Constant everywhere, set with namelist options iccg_prescribed (num|den)#, default is 0./1. (all CG).
1106 = 2 ! Coarsely prescribed 1995-1999 NLDN/OTD climatology based on Boccippio et al. (2001)
1107 = 3 ! Parameterization by Price and Rind (1993) based on cold-cloud depth
1108 = 4 ! Gridded input via arrays iccg_in_(num|den) from wrfinput for monthly mapped ratios.
1109 Points with 0/0 values use ratio defined by iccg_prescribed_(num|den)
1110 iccg_prescribed_num (max_dom) = 0. ! Numerator of user-specified prescribed IC:CG
1111 iccg_prescribed_den (max_dom) = 1. ! Denominator of user-specified prescribed IC:CG
1112 ltng_temp_upper = -45. ! Temperature (C) of upper peak of LNOx vertical distribution for IC lightning (used by &chem lnox_opt=2).
1113 ltng_temp_lower = -15. ! Temperatures (C) of lower peak of LNOx vertical distribution for both IC and CG lightning (used by &chem lnox_opt=2).
1115 Options for MAD-WRF - see doc/README.madwrf for usage information
1116 madwrf_opt = 0 ! MAD-WRF model: 0) Off, 1) Advect/diffuse passive hydrometeors (qcloud, qice, qsnow), 2) Nudge active microphysics hydrometeors to passive ones
1117 madwrf_dt_relax = 60.0 ! Relaxation time for hydrometeor nudging [s]
1118 madwrf_dt_nudge = 60.0 ! Temporal period for hydrometeor nudging [min]
1119 madwrf_cldinit = 0 ! Enhance cloud initialization: 0: Off, 1: On
1122 Options for wind turbine drag parameterization:
1124 windfarm_opt (max_dom) = 0 ! 1 = Simulates the effects of wind turbines in the atmospheric evolution
1125 windfarm_ij = 0 ! whether to use lat-lon or i-j coordinate as wind turbine locations
1126 ! 0 = The coordinate of the turbines are defined in terms of lat-lon
1127 ! 1 = The coordinate of the turbines are defined in terms of grid points
1128 windfarm_tke_factor = 0.25 ! Correction factor applied to the TKE coefficient (deafault is 0.25, Archer et al. 2020)
1131 Stochastic parameterization schemes:
1133 &stoch
1134 ; Random perturbation field (rand_perturb=1)
1135 rand_perturb (max_dom) = 1 ! Generate array with random perturbations for user-determined use, 1: on
1136 gridpt_stddev_rand_pert (max_dom) = 0.03 ! Standard deviation of random perturbation field at each gridpoint.
1137 Determines amplitude of random perturbations
1138 lengthscale_rand_pert (max_dom) = 500000.0 ! Perturbation lengthscale (in m).
1139 timescale_rand_pert (max_dom) = 21600.0 ! Temporal decorrelation of random field (in s).
1140 stddev_cutoff_rand_pert (max_dom) = 3.0 ! Cutoff tails of perturbation pattern above this threshold standard deviation.
1141 rand_pert_vertstruc = 0 ! Vertical structure for random perturbation field: 0=constant; 1=random phase with tilt
1142 iseed_rand_pert Seed for random number stream for rand_perturb. Will be
1143 combined with seed nens signifying ensemble member number and initial
1144 start time to ensure different random number streams for forecasts
1145 starting from different initial times and for different ensemble members.
1147 ; Stochastically perturbed physics tendencies (SPPT) (sppt=1)
1148 sppt (max_dom) = 0 ! Stochastically perturbed physics tendencies (SPPT), 0: off, 1: on
1149 gridpt_stddev_sppt (max_dom) = 0.5 ! Standard deviation of random perturbation field at each gridpoint.
1150 Determines amplitude of random perturbations
1151 lengthscale_sppt (max_dom) = 150000.0 ! Perturbation lengthscale (in m).
1152 timescale_sppt (max_dom) = 21600.0 ! Temporal decorrelation of random field (in s).
1153 sppt_vertstruc = 0. ! vertical structure for sppt, 0: constant, 1: random phase.
1154 stddev_cutoff_sppt (max_dom) = 2.0 ! Cutoff tails of perturbation pattern above this threshold standard deviation.
1155 iseed_sppt Seed for random number stream for sppt. Will be
1156 combined with seed nens signifying ensemble member number and initial
1157 start time to ensure different random number streams for forecasts
1158 starting from different initial times and for different ensemble members.
1160 ; Stochastic kinetic-energy backscatter scheme (SKEBS)(skebs=1):
1162 skebs (max_dom) = 0 ! stochastic kinetic-energy backscatter scheme, 0: off, 1: on
1163 tot_backscat_psi (max_dom) = 1.0E-05 ! Controls amplitude of rotational wind perturbations
1164 tot_backscat_t (max_dom) = 1.0E-06 ! Controls amplitude of potential temperature perturbations
1165 ztau_psi = 10800.0 ! decorr. time of noise for psi perturb
1166 ztau_t = 10800.0 ! decorr. time of noise for theta perturb
1167 rexponent_psi = -1.83 ! spectral slope of forcing for psi
1168 rexponent_t = -1.83 ! spectral slope of forcing for theta
1169 zsigma2_eps = 0.0833 ! variance of noise for psi perturb
1170 zsigma2_eta = 0.0833 ! variance of noise for theta perturb
1171 kminforc = 1 ! min. forcing wavenumber in lon. for psi perturb
1172 lminforc = 1 ! min. forcing wavenumber in lat. for psi perturb
1173 kminforct = 1 ! min. forcing wavenumber in lon. for theta perturb
1174 lminforct = 1 ! min. forcing wavenumber in lat. for theta perturb
1175 kmaxforc = 1000000 ! max. forcing wavenumber in lon. for psi perturb
1176 lmaxforc = 1000000 ! max. forcing wavenumber in lat. for psi perturb
1177 kmaxforct = 1000000 ! max. forcing wavenumber in lon. for theta perturb
1178 lmaxforct = 1000000 ! max. forcing wavenumber in lat. for theta perturb
1179 skebs_vertstruc = 0 ! Vertical structure for random perturbation field: 0=constant; 1=random phase with tilt
1181 Stochastically perturbed parameter scheme (SPP) (spp=1)
1182 sppt (max_dom) = 0 ! Stochastically perturbed parameter (SPP) scheme for
1183 ; GF convection scheme, MYNN boundary layer scheme and RUC LSM. 0: off, 1: on
1184 spp_conv (max_dom) = 0 ! Perturb parameters of GF convection scheme only
1185 gridpt_stddev_spp_conv (max_dom) = 0.3 ! Standard deviation of random perturbation field at each gridpoint.
1186 Determines amplitude of random perturbations
1187 lengthscale_spp_conv (max_dom) = 150000.0 ! Perturbation lengthscale (in m).
1188 timescale_spp_conv (max_dom) = 21600.0 ! Temporal decorrelation of random field (in s).
1189 stddev_cutoff_spp_conv (max_dom) = 3.0 ! Cutoff tails of perturbation pattern above this threshold standard deviation.
1190 iseed_spp_conv Seed for random number stream for spp_conv. Will be
1191 combined with seed nens signifying ensemble member number and initial
1192 start time to ensure different random number streams for forecasts
1193 starting from different initial times and for different ensemble members.
1194 spp_pbl (max_dom) = 0 ! Perturb parameters of MYNN PBL scheme only
1195 gridpt_stddev_spp_pbl (max_dom) = 0.15 ! Standard deviation of random perturbation field at each gridpoint.
1196 Determines amplitude of random perturbations
1197 lengthscale_spp_pbl (max_dom) = 70000.0 ! Perturbation lengthscale (in m).
1198 timescale_spp_pbl (max_dom) = 21600.0 ! Temporal decorrelation of random field (in s).
1199 stddev_cutoff_spp_pbl (max_dom) = 2.0 ! Cutoff tails of perturbation pattern above this threshold standard deviation.
1200 iseed_spp_pbl Seed for random number stream for spp_pbl . Will be
1201 combined with seed nens signifying ensemble member number and initial
1202 start time to ensure different random number streams for forecasts
1203 starting from different initial times and for different ensemble members.
1204 spp_lsm (max_dom) = 0 ! Perturb parameters of RUC LSM
1205 gridpt_stddev_spp_lsm (max_dom) = 0.3 ! Standard deviation of random perturbation field at each gridpoint.
1206 ; Determines amplitude of random perturbations
1207 lengthscale_spp_lsm (max_dom) = 50000.0 ! Perturbation lengthscale (in m).
1208 timescale_spp_lsm (max_dom) = 86400.0 ! Temporal decorrelation of random field (in s).
1209 stddev_cutoff_spp_lsm (max_dom) = 3.0 ! Cutoff tails of perturbation pattern above this threshold standard deviation.
1210 iseed_spp_lsm Seed for random number stream for spp_lsm . Will be
1211 combined with seed nens signifying ensemble member number and initial
1212 start time to ensure different random number streams for forecasts
1213 starting from different initial times and for different ensemble members.
1215 Multiple perturbations for WRF-Solar EPS (multi_perturb = 1)
1216 multi_perturb = 0 ! 1) WRF-Solar EPS: turns on stochastic perturbations tailored for solar energy applications
1217 spdt = -1.0 ! Frequency to update the stochastic perturbations [minutes]. A negative value indicates for every time step
1218 pert_farms = .false. ! Activates perturbations to the FARMS parameterization
1219 pert_farms_albedo = 0.0 ! 1.0) Perturbs the albedo, 0.0) no perturbations. Similar for other entries below
1220 pert_farms_aod = 0.0
1221 pert_farms_angexp = 0.0
1222 pert_farms_aerasy = 0.0
1223 pert_farms_qv = 0.0
1224 pert_farms_qc = 0.0
1225 pert_farms_qs = 0.0
1226 pert_deng = .false. ! Activates perturbations to Deng's shcu parameterization
1227 pert_deng_qv = 0.0
1228 pert_deng_qc = 0.0
1229 pert_deng_t = 0.0
1230 pert_deng_w = 0.0
1231 pert_mynn = .false. ! Activates perturbations to the MYNN PBL parameterization
1232 pert_mynn_qv = 0.0
1233 pert_mynn_qc = 0.0
1234 pert_mynn_t = 0.0
1235 pert_mynn_qke = 0.0
1236 pert_noah = .false. ! Activates perturbations to the Noah LSM
1237 pert_noah_qv = 0.0
1238 pert_noah_t = 0.0
1239 pert_noah_smois = 0.0
1240 pert_noah_tslb = 0.0
1241 pert_thom = .false. ! Activates perturbations to thompson microphysics
1242 pert_thom_qv = 0.0
1243 pert_thom_qc = 0.0
1244 pert_thom_qi = 0.0
1245 pert_thom_qs = 0.0
1246 pert_thom_ni = 0.0
1247 pert_cld3 = .false. ! Activates perturbations to clouds generated with icloud = 3
1248 pert_cld3_qv = 0.0
1249 pert_cld3_t = 0.0
1251 num_pert_3d = 15 ! Number of entries in STOCHPERT.TBL plus one (No need to modify)
1253 ; Stochastic Perturbations to the boundary conditions?| (perturb_bdy)
1254 perturb_bdy = 0 ! No boundary perturbations
1255 1 Use SKEBS pattern for boundary perturbations
1256 2 Use other user-provided pattern for boundary perturbations
1258 ; Stochastic perturbations to the boundary tendencies in WRF-CHEM (perturb_chem_bdy)
1259 perturb_chem_bdy Options for perturbing lateral boundaries of chemical tracers:
1260 0 = off; 1 = on with RAND_PERTURB pattern
1262 ; Common to all stochastic schemes
1263 nens =1 ! Seed for random number stream. For ensemble forecasts this parameter needs to be
1264 different for each member. The seed is a function of initial start time
1265 to ensure different random number streams for forecasts starting from
1266 different initial times. Changing this seed changes the random number
1267 streams for all activated stochastic parameterization schemes.
1270 Options for use with the Noah-MP Land Surface Model (sf_surface_physics=4):
1272 &noah_mp
1273 dveg = 4, ! Noah-MP Dynamic Vegetation option:
1274 1 = Off (LAI from table; FVEG = shdfac)
1275 2 = On (LAI predicted; FVEG calculated)
1276 3 = Off (LAI from table; FVEG calculated)
1277 4 = Off (LAI from table; FVEG = maximum veg. fraction)
1278 5 = On (LAI predicted; FVEG = maximum veg. fraction)
1279 6 = On (use FVEG = SHDFAC from input)
1280 7 = Off (use input LAI; use FVEG = SHDFAC from input)
1281 8 = Off (use input LAI; calculate FVEG)
1282 9 = Off (use input LAI; use maximum vegetation fraction)
1283 opt_crs = 1, ! Noah-MP Stomatal Resistance option:
1284 1 = Ball-Berry
1285 2 = Jarvis
1286 opt_sfc = 1 ! Noah-MP surface layer drag coefficient calculation
1287 1 = Monin-Obukhov
1288 2 = original Noah (Chen97)
1289 opt_btr = 1, ! Noah-MP Soil Moisture Factor for Stomatal Resistance
1290 1 = Noah
1291 2 = CLM
1292 3 = SSiB
1293 opt_run = 3, ! Noah-MP Runoff and Groundwater option
1294 1 = TOPMODEL with groundwater
1295 2 = TOPMODEL with equilibrium water table
1296 3 = original surface and subsurface runoff (free drainage) - default
1297 4 = BATS surface and subsurface runoff (free drainage)
1298 5 = Miguez-Macho & Fan groundwater scheme (Miguez-Macho et al. 2007 JGR; Fan et al. 2007 JGR)
1299 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1300 6 = Variable Infiltration Capacity Model surface runoff scheme (Wood et al., 1992, JGR)
1301 7 = Xiananjiang Infiltration and surface runoff scheme ((Jayawardena and Zhou, 2000)
1302 8 = Dynamic VIC surface runoff scheme (Liang and Xie, 2001)
1303 opt_infdv = 0, ! Noah-MP infiltration option in dynamic VIC runoff scheme (only works for opt_run=8)
1304 1 = Philip scheme
1305 2 = Green-Ampt scheme
1306 3 = Smith-Parlange scheme
1307 opt_frz = 1, ! Noah-MP Supercooled Liquid Water option
1308 1 = No iteration
1309 2 = Koren's iteration
1310 opt_inf = 1, ! Noah-MP Soil Permeability option
1311 1 = Linear effects, more permeable
1312 2 = Non-linear effects, less permeable
1313 opt_rad = 3, ! Noah-MP Radiative Transfer option
1314 1 = Modified two-stream (known to cause problems when vegetation fraction is small)
1315 2 = Two-stream applied to grid-cell
1316 3 = Two-stream applied to vegetated fraction
1317 opt_alb = 2, ! Noah-MP Ground Surface Albedo option
1318 1 = BATS
1319 2 = CLASS
1320 opt_snf = 1, ! Noah-MP Precipitation Partitioning between snow and rain
1321 1 = Jordan (1991)
1322 2 = BATS: Snow when SFCTMP < TFRZ+2.2
1323 3 = Snow when SFCTMP < TFRZ
1324 4 = Use WRF precipitation partitioning
1325 5 = Use wetbulb temperature (Wang et al., 2019 GRL)
1326 opt_tbot = 2, ! Noah-MP Soil Temperature Lower Boundary Condition
1327 1 = Zero heat flux
1328 2 = TBOT at 8 m from input file
1329 opt_stc = 1, ! Noah-MP Snow/Soil temperature time scheme
1330 1 = semi-implicit
1331 2 = full-implicit
1332 3 = semi-implicit where Ts uses snow cover fraction
1333 opt_gla = 1, ! Noah-MP glacier treatment option
1334 1 = includes phase change
1335 2 = slab ice (Noah)
1336 opt_rsf = 1, ! Noah-MP surface evaporation resistance option
1337 1 = Sakaguchi and Zeng, 2009
1338 2 = Sellers (1992)
1339 3 = adjusted Sellers to decrease RSURF for wet soil
1340 4 = option 1 for non-snow; rsurf = rsurf_snow for snow (set in MPTABLE)
1341 opt_soil = 1, ! Noah-MP options for defining soil properties
1342 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1343 1 = use input dominant soil texture
1344 2 = use input soil texture that varies with depth
1345 3 = use soil composition (sand, clay, orgm) and pedotransfer functions (OPT_PEDO)
1346 4 = use input soil properties (BEXP_3D, SMCMAX_3D, etc.) (not valid in WRF)
1347 opt_pedo = 1, ! Noah-MP options for pedotransfer functions (used when OPT_SOIL = 3)
1348 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1349 1 = Saxton and Rawls (2006)
1350 opt_crop = 0, ! options for crop model
1351 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1352 0 = No crop model, will run default dynamic vegetation
1353 1 = Liu, et al. 2016
1354 2 = Gecros (Genotype-by-Environment interaction on CROp growth Simulator) Yin and van Laar, 2005
1356 opt_irr = 0, ! options for irrigation scheme
1357 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1358 0 = No irrigation
1359 1 = Irrigation ON
1360 2 = irrigation trigger based on crop season Planting and harvesting dates
1361 3 = irrigation trigger based on LAI threshold
1362 opt_irrm = 0, ! options for irrigation method (only if opt_irr > 0)
1363 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1364 0 = method based on geo_em fractions (all three methods are ON)
1365 1 = sprinkler method
1366 2 = micro/drip irrigation
1367 3 = surface flooding
1368 opt_tdrn = 0, ! Noah-MP tile drainage option (currently only tested and works with opt_run=3)
1369 geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
1370 0 = No tile drainage
1371 1 = Simple drainage
1372 2 = Hooghoudt's equation based tile drainage
1373 soiltstep = 0.0, ! Noah-MP soil process timestep (seconds) for solving soil water and temperature
1374 0 = default, the same as main NoahMP model timestep
1375 N * dt_noahmp, typically 15min or 30min
1376 noahmp_output = 1, ! Noah-MP output levels
1377 1 = standard output
1378 3 = standard output with additional water and energy budget term output
1379 noahmp_acc_dt = 0.0, ! Noah-MP bucket reset time interval (minutes) between outputs for accumulation (works when noahmp_output=3)
1380 /
1382 &fdda
1383 grid_fdda (max_dom) = 1 ! grid-nudging fdda on (=0 off) for each domain
1384 = 2 ! spectral nudging
1385 gfdda_inname = "wrffdda_d<domain>" ! defined name in real
1386 gfdda_interval_m (max_dom) = 360 ! time interval (in min) between analysis times (must use minutes)
1387 gfdda_end_h (max_dom) = 6 ! time (in hours) to stop nudging after start of forecast
1388 io_form_gfdda = 2 ! analysis data io format (2 = netCDF)
1389 fgdt (max_dom) = 0 ! calculation frequency (minutes) for grid-nudging (0=every step)
1390 if_no_pbl_nudging_uv (max_dom) = 0 ! 1= no nudging of u and v in the pbl, 0=nudging in the pbl
1391 if_no_pbl_nudging_t (max_dom) = 0 ! 1= no nudging of temp in the pbl, 0=nudging in the pbl
1392 if_no_pbl_nudging_q (max_dom) = 0 ! 1= no nudging of qvapor in the pbl, 0=nudging in the pbl
1393 if_zfac_uv (max_dom) = 0 ! 0= nudge u and v in all layers, 1= limit nudging to levels above k_zfac_uv
1394 k_zfac_uv (max_dom) = 10 ! 10=model level below which nudging is switched off for u and v
1395 if_zfac_t (max_dom) = 0 ! 0= nudge temp in all layers, 1= limit nudging to levels above k_zfac_t
1396 k_zfac_t (max_dom) = 10 ! 10=model level below which nudging is switched off for temp
1397 if_zfac_q (max_dom) = 0 ! 0= nudge qvapor in all layers, 1= limit nudging to levels above k_zfac_q
1398 k_zfac_q (max_dom) = 10 ! 10=model level below which nudging is switched off for qvapor
1399 guv (max_dom) = 0.0003 ! nudging coefficient for u and v (sec-1)
1400 gt (max_dom) = 0.0003 ! nudging coefficient for temp (sec-1)
1401 gq (max_dom) = 0.00001 ! nudging coefficient for qvapor (sec-1)
1402 if_ramping = 0 ! 0= nudging ends as a step function, 1= ramping nudging down at end of period
1403 dtramp_min = 0 ! time (min) for ramping function
1404 grid_sfdda (max_dom) = 0 ! surface fdda switch
1405 0: off;
1406 1: nudging selected surface fields;
1407 2: FASDAS (flux-adjusting surface data assimilation system)
1408 sgfdda_inname = "wrfsfdda_d<domain>" ! defined name for sfc nudgingi in input file (from program obsgrid)
1409 sgfdda_end_h (max_dom) = 6 ! time (in hours) to stop sfc nudging after start of forecast
1410 sgfdda_interval_m (max_dom) = 180 ! time interval (in min) between sfc analysis times (must use minutes)
1411 io_form_sgfdda = 2 ! sfc analysis data io format (2 = netCDF)
1412 guv_sfc (max_dom) = 0.0003 ! nudging coefficient for sfc u and v (sec-1)
1413 gt_sfc (max_dom) = 0.0003 ! nudging coefficient for sfc temp (sec-1)
1414 gq_sfc (max_dom) = 0.00001 ! nudging coefficient for sfc qvapor (sec-1)
1415 rinblw (max_dom) = 0. ! radius of influence used to determine the confidence (or weights) for
1416 the analysis, which is based on the distance between the grid point to the nearest
1417 obs. The analysis without nearby observation is used at a reduced weight.
1419 pxlsm_soil_nudge(max_dom) = 1 ! PXLSM Soil nudging option (requires wrfsfdda file)
1421 The following are for spectral nudging:
1422 fgdtzero (max_dom) = 0, ! 1= nudging tendencies are set to zero in between fdda calls
1423 if_no_pbl_nudging_uv (max_dom) = 0, ! 1= no nudging of uv in the pbl, 0= nudging in the pbl
1424 if_no_pbl_nudging_t (max_dom) = 0, ! 1= no nudging of t in the pbl, 0= nudging in the pbl
1425 if_no_pbl_nudging_ph (max_dom) = 0, ! 1= no nudging of ph in the pbl, 0= nudging in the pbl
1426 if_no_pbl_nudging_q (max_dom) = 0, ! 1= no nudging of q in the pbl, 0= nudging in the pbl
1427 if_zfac_uv (max_dom) = 0, ! 0= nudge uv in all layers, 1= limit nudging to levels above k_zfac_uv
1428 k_zfac_uv (max_dom) = 0, ! 0= model level below which nudging is switched off for uv
1429 dk_zfac_uv (max_dom) = 1, ! depth in k between k_zfac_X to dk_zfac_X where nudging increases
1430 ; linearly to full strength
1431 if_zfac_t (max_dom) = 0, ! 0= nudge t in all layers, 1= limit nudging to levels above k_zfac_t
1432 k_zfac_t (max_dom) = 0, ! 0= model level below which nudging is switched off for t
1433 dk_zfac_t (max_dom) = 1, ! depth in k between k_zfac_X to dk_zfac_X where nudging increases
1434 linearly to full strength
1435 if_zfac_ph (max_dom) = 0, ! 0= nudge ph in all layers, 1= limit nudging to levels above k_zfac_ph
1436 k_zfac_ph (max_dom) = 0, ! 0= model level below which nudging is switched off for ph
1437 dk_zfac_ph (max_dom) = 1, ! depth in k between k_zfac_X to dk_zfac_X where nudging increases
1438 linearly to full strength
1439 if_zfac_q (max_dom) = 0, ! 0= nudge q in all layers, 1= limit nudging to levels above k_zfac_q
1440 k_zfac_q (max_dom) = 0, ! 0= model level below which nudging is switched off for q
1441 dk_zfac_q (max_dom) = 1, ! depth in k between k_zfac_X to dk_zfac_X where nudging increases
1442 gph (max_dom) = 0.0003,
1443 ktrop = 0, ! layer nominally representing tropopause for limiting nudging to q and t
1444 ; setting ktrop = 0 allows nudging to extend to the top of the atmosphere
1445 xwavenum (max_dom) = 3, ! top wave number to nudge in x direction
1446 ywavenum (max_dom) = 3, ! top wave number to nudge in y direction
1448 The following are for observation nudging:
1449 obs_nudge_opt (max_dom) = 1 ! obs-nudging fdda on (=0 off) for each domain
1450 also need to set auxinput11_interval and auxinput11_end_h
1451 in time_control namelist
1452 max_obs = 0 ! max number of observations used on a domain during any
1453 given time window
1454 fdda_start (max_dom) = 0 ! obs nudging start time in minutes
1455 fdda_end (max_dom) = 0 ! obs nudging end time in minutes
1456 obs_nudge_wind (max_dom) = 1 ! whether to nudge wind: (=0 off)
1457 obs_coef_wind (max_dom) = 0, ! nudging coefficient for wind, unit: s-1
1458 obs_nudge_temp (max_dom) = 0, ! set to = 1 to turn to nudge temperature (default = 0; off)
1459 obs_coef_temp (max_dom) = 0, ! nudging coefficient for temperature, unit: s-1
1460 obs_nudge_mois (max_dom) = 1 ! whether to nudge water vapor mixing ratio: (=0 off)
1461 obs_coef_mois (max_dom) = 0, ! nudging coefficient for water vapor mixing ratio, unit: s-1
1462 obs_nudge_pstr (max_dom) = 0 ! whether to nudge surface pressure (not used)
1463 obs_coef_pstr (max_dom) = 0. ! nudging coefficient for surface pressure, unit: s-1 (not used)
1464 obs_rinxy (max_dom) = 0. ! horizonal radius of influence in km
1465 obs_rinsig = 0 ! vertical radius of influence in eta
1466 obs_twindo (max_dom) = 0.66667 ! half-period time window over which an observation
1467 will be used for nudging (hours)
1468 obs_npfi = 0, ! freq in coarse grid timesteps for diagnostic prints
1469 obs_ionf (max_dom) = 2 ! freq in coarse grid timesteps for obs input and err calc
1470 obs_idynin = 0 ! for dynamic initialization using a ramp-down function to gradually
1471 turn off the FDDA before the pure forecast (=1 on)
1472 obs_dtramp = 0 ! time period in minutes over which the nudging is ramped down
1473 from one to zero.
1474 obs_prt_freq (max_dom) = 10, ! Frequency in obs index for diagnostic printout
1475 obs_prt_max = 1000, ! Maximum allowed obs entries in diagnostic printout
1476 obs_ipf_in4dob = .true. ! print obs input diagnostics (=.false. off)
1477 obs_ipf_errob = .true. ! print obs error diagnostics (=.false. off)
1478 obs_ipf_nudob = .true. ! print obs nudge diagnostics (=.false. off)
1479 obs_ipf_init = .true. ! Enable obs init warning messages
1481 obs_no_pbl_nudge_uv (max_dom) = 0 ! 1=no wind-nudging within pbl
1482 obs_no_pbl_nudge_t (max_dom) = 0 ! 1=no temperature-nudging within pbl
1483 obs_no_pbl_nudge_q (max_dom) = 0 ! 1=no moisture-nudging within pbl
1484 obs_sfc_scheme_horiz = 0 ! horizontal spreading scheme for surf obs;
1485 ; 0=wrf scheme, 1=original mm5 scheme
1486 obs_sfc_scheme_vert = 0 ! vertical spreading scheme for surf obs
1487 0=regime vif scheme, 1=original simple scheme
1488 obs_max_sndng_gap = 20 ! Max pressure gap between soundings, in cb
1489 obs_nudgezfullr1_uv = 50 ! Vert infl full weight height for lowest model level (LML) obs, regime 1, winds
1490 obs_nudgezrampr1_uv = 50 ! Vert infl ramp-to-zero height for LML obs, regime 1, winds
1491 obs_nudgezfullr2_uv = 50 ! Vert infl full weight height for LML obs, regime 2, winds
1492 obs_nudgezrampr2_uv = 50 ! Vert infl ramp-to-zero height for LML obs, regime 2, winds
1493 obs_nudgezfullr4_uv = -5000 ! Vert infl full weight height for LML obs, regime 4, winds
1494 obs_nudgezrampr4_uv = 50 ! Vert infl ramp-to-zero height for LML obs, regime 4, winds
1495 obs_nudgezfullr1_t = 50 ! Vert infl full weight height for LML obs, regime 1, temperature
1496 obs_nudgezrampr1_t = 50 ! Vert infl ramp-to-zero height for LML obs, regime 1, temperature
1497 obs_nudgezfullr2_t = 50 ! Vert infl full weight height for LML obs, regime 2, temperature
1498 obs_nudgezrampr2_t = 50 ! Vert infl ramp-to-zero height for LML obs, regime 2, temperature
1499 obs_nudgezfullr4_t = -5000 ! Vert infl full weight height for LML obs, regime 4, temperature
1500 obs_nudgezrampr4_t = 50 ! Vert infl ramp-to-zero height for LML obs, regime 4, temperature
1501 obs_nudgezfullr1_q = 50 ! Vert infl full weight height for LML obs, regime 1, moisture
1502 obs_nudgezrampr1_q = 50 ! Vert infl ramp-to-zero height for LML obs, regime 1, moisture
1503 obs_nudgezfullr2_q = 50 ! Vert infl full weight height for LML obs, regime 2, moisture
1504 obs_nudgezrampr2_q = 50 ! Vert infl ramp-to-zero height for LML obs, regime 2, moisture
1505 obs_nudgezfullr4_q = -5000 ! Vert infl full weight height for LML obs, regime 4, moisture
1506 obs_nudgezrampr4_q = 50 ! Vert infl ramp-to-zero height for LML obs, regime 4, moisture
1507 obs_nudgezfullmin = 50 ! Min depth through which vertical infl fcn remains 1.0
1508 obs_nudgezrampmin = 50 ! Min depth (m) through which vert infl fcn decreases from 1 to 0
1509 obs_nudgezmax = 3000 ! Max depth (m) in which vert infl function is nonzero
1510 obs_sfcfact = 1.0 ! Scale factor applied to time window for surface obs
1511 obs_sfcfacr = 1.0 ! Scale factor applied to horiz radius of influence for surface obs
1512 obs_dpsmx = 7.5 ! Max pressure change (cb) allowed within horiz radius of influence
1514 obs_scl_neg_qv_innov = 0 ! 1 = prevent to nudge toward negative QV
1515 /
1517 &scm
1518 scm_force = 1, ! switch for single column forcing (=0 off)
1519 scm_force_dx = 4000. ! DX for SCM forcing (in meters)
1520 num_force_layers = 8 ! number of SCM input forcing layers
1521 scm_lu_index = 2 ! SCM landuse category (2 is dryland, cropland and pasture)
1522 scm_isltyp = 4 ! SCM soil category (4 is silt loam)
1523 scm_vegfra = 0.5 ! SCM vegetation fraction
1524 scm_canwat = 0.0 ! SCM canopy water
1525 scm_lat = 36.605 ! SCM latitude
1526 scm_lon = -97.485 ! SCM longitude
1527 scm_th_adv = .true. ! turn on theta advection in SCM
1528 scm_wind_adv = .true. ! turn on wind advection in SCM
1529 scm_qv_adv = .true. ! turn on moisture advection in SCM
1530 scm_ql_adv = .true. ! turn on cloud liquid water advection in SCM
1531 scm_vert_adv = .true. ! turn on vertical advection in SCM
1532 num_force_soil_layers = 5, ! Number of SCM soil forcing layer
1533 scm_soilT_force = .false. ! Turn on soil temp forcing in SCM
1534 scm_soilq_force = .false. ! Turn on soil moisture forcing in SCM
1535 scm_force_th_largescale = .false. ! Turn on large scale theta forcing in SCM
1536 scm_force_qv_largescale = .false. ! Turn on large scale qv forcing in SCM
1537 scm_force_ql_largescale = .false. ! Turn on large scale cloud water forcing in SCM
1538 scm_force_wind_largescale = .false. ! Turn on large scale wind forcing in SCM
1540 &dynamics
1541 hybrid_opt = 2, ! default; Klemp cubic form with etac
1542 0 = original WRF terrain-following coordinate (through V3)
1543 etac = 0.2 ! znw(k) < etac, eta surfaces are isobaric, 0.2 is a good default (Pa/Pa)
1544 rk_ord = 3, ! time-integration scheme option:
1545 2 = Runge-Kutta 2nd order
1546 3 = Runge-Kutta 3rd order
1547 zadvect_implicit = 0, ! toggle off on [0/1] the IEVA scheme. Default off.
1548 w_crit_cfl = 1.2 ! default vertical courant number where vertical velocity damping begins (see below).
1549 ! when zadvect_implicit is on, value can be ~ 2.0
1550 zdamp (max_dom) = 5000., ! damping depth (m) from model top
1551 w_damping = 0, ! vertical velocity damping flag (for operational use)
1552 0 = without damping
1553 1 = with damping
1554 diff_opt(max_dom) = 0, ! turbulence and mixing option:
1555 0 = no turbulence or explicit
1556 spatial numerical filters (km_opt IS IGNORED).
1557 1 = evaluates 2nd order
1558 diffusion term on coordinate surfaces.
1559 uses kvdif for vertical diff unless PBL option
1560 is used. may be used with km_opt = 1 and 4.
1561 (= 1, recommended for real-data cases)
1562 2 = evaluates mixing terms in
1563 physical space (stress form) (x,y,z).
1564 turbulence parameterization is chosen
1565 by specifying km_opt.
1566 km_opt(max_dom) = 1, ! eddy coefficient option
1567 1 = constant (use khdif kvdif)
1568 2 = 1.5 order TKE closure (3D)
1569 3 = Smagorinsky first order closure (3D)
1570 Note: option 2 and 3 are not recommended for DX > 2 km
1571 4 = horizontal Smagorinsky first order closure
1572 (recommended for real-data cases)
1573 5 = SMS-3DTKE scale-adaptive LES/PBL scheme. It must be
1574 used with diff_opt = 2. PBL schemes must be turned off
1575 (bl_pbl_physics=0). It can work with sf_sfclay_physics = 1, 5, 91.
1576 damp_opt = 0, ! upper level damping flag
1577 0 = without damping
1578 1 = with diffusive damping, maybe used for real-data cases
1579 (dampcoef nondimensional ~0.01-0.1)
1580 2 = with Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .003; idealized case only
1581 not for real-data cases)
1582 3 = with w-Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .2;
1583 for real-data cases)
1584 use_theta_m = 1 ! 1: use theta_m=theta(1+1.61Qv)
1585 0: use dry theta in dynamics
1586 use_q_diabatic = 0 ! whether to include QV and QC tendencies in advection
1587 0 = default, old behavior
1588 1 = include QV and QC tendencies - this helps to produce correct solution
1589 in an idealized 'moist benchmark' test case (Bryan, 2014).
1590 In real data testing, time step needs to be reduce to maintain stable solution
1591 c_s (max_dom) = 0.25 ! Smagorinsky coeff
1592 c_k (max_dom) = 0.15 ! TKE coeff
1593 diff_6th_opt (max_dom) = 0, ! 6th-order numerical diffusion
1594 0 = no 6th-order diffusion (default)
1595 1 = 6th-order numerical diffusion (not recommended)
1596 2 = 6th-order numerical diffusion but prohibit up-gradient diffusion
1597 diff_6th_factor (max_dom) = 0.12, ! 6th-order numerical diffusion non-dimensional rate (max value 1.0
1598 corresponds to complete removal of 2dx wave in one timestep)
1599 diff_6th_slopeopt (max_dom) = 0 ! if set to =1, turns on 6th-order numerical diffusion - terrain-slope tapering. default is 0=off
1600 diff_6th_thresh (max_dom) = 0.10 ! slope threshold (m/m) that turns off 6th order diff in steep terrain
1601 dampcoef (max_dom) = 0., ! damping coefficient (see above)
1602 base_temp = 290., ! real-data, em ONLY, base sea-level temp (K)
1603 base_pres = 10^5 ! real-data, em ONLY, base sea-level pres (Pa), DO NOT CHANGE
1604 base_lapse = 50., ! real-data, em ONLY, lapse rate (K), DO NOT CHANGE
1605 iso_temp = 200., ! real-data, em ONLY, reference temp in stratosphere, US Standard atmosphere 216.5 K
1606 base_pres_strat = 0. ! real-data, em ONLY, base state pressure (Pa) at bottom of the stratosphere,
1607 US Standard atmosphere 55 hPa
1608 base_lapse_strat = -11. ! real-data, em ONLY, base state lapse rate ( dT / d(lnP) ) in stratosphere,
1609 approx to US Standard atmosphere -12 K
1610 use_baseparam_fr_nml = .f., ! whether to use base state parameters from the namelist
1611 use_input_w = .f., ! whether to use vertical velocity from input file
1612 khdif (max_dom) = 0, ; horizontal diffusion constant (m^2/s). A typical value should be 0.1*DX in meters.
1613 kvdif (max_dom) = 0, ! vertical diffusion constant (m^2/s). A typical value should be 100.
1614 smdiv (max_dom) = 0.1, ! divergence damping (0.1 is typical)
1615 emdiv (max_dom) = 0.01, ! external-mode filter coef for mass coordinate model
1616 (0.01 is typical for real-data cases)
1617 epssm (max_dom) = .1, ! time off-centering for vertical sound waves
1618 non_hydrostatic (max_dom) = .true., ! whether running the model in hydrostatic or non-hydro mode
1619 pert_coriolis (max_dom) = .false., ! Coriolis only acts on wind perturbation (idealized)
1620 top_lid (max_dom) = .false., ! Zero vertical motion at top of domain
1621 mix_full_fields(max_dom) = .true., ! used with diff_opt = 2; value of ".true." is recommended, except for
1622 highly idealized numerical tests; damp_opt must not be 1 if ".true."
1623 is chosen. .false. means subtract 1-d base-state profile before mixing
1624 mix_isotropic(max_dom) = 0 ! 0=anisotropic vertical/horizontal diffusion coeffs, 1=isotropic
1625 mix_upper_bound(max_dom) = 0.1 ! non-dimensional upper limit for diffusion coeffs
1626 tke_drag_coefficient(max_dom) = 0., ! surface drag coefficient (Cd, dimensionless) for diff_opt=2 only
1627 tke_heat_flux(max_dom) = 0., ! surface thermal flux (H/(rho*cp), K m/s) for diff_opt=2 only
1628 h_mom_adv_order (max_dom) = 5, ! horizontal momentum advection order (5=5th, etc.)
1629 v_mom_adv_order (max_dom) = 3, ! vertical momentum advection order
1630 h_sca_adv_order (max_dom) = 5, ! horizontal scalar advection order
1631 v_sca_adv_order (max_dom) = 3, ! vertical scalar advection order
1633 momentum_adv_opt(max_dom) = 1, ! advection options for momentum variables:
1634 1=original, 3 = 5th-order WENO
1635 advection options for scalar variables: 0=simple, 1=positive definite,
1636 2=monotonic, 3=5th order WENO, 4=5th-order WENO with positive definite filter
1637 moist_adv_opt (max_dom) = 1 ! for moisture
1638 moist_adv_dfi_opt (max_dom) = 0 ! positive-definite RK3 transport switch. Default is 0=off
1639 scalar_adv_opt (max_dom) = 1 ! for scalars
1640 chem_adv_opt (max_dom) = 1 ! for chem variables
1641 tracer_adv_opt (max_dom) = 1 ! for tracer variables (WRF-Chem activated)
1642 tke_adv_opt (max_dom) = 1 ! for tke
1643 phi_adv_z (max_dom) = 1 ! vertical advection option for geopotential; 1: original (default) 2: avoid double staggering of omega
1644 moist_mix2_off (max_dom) = .false. ! if set to T, deactivate 2nd-order horizontal mixing for moisture. default is F.
1645 chem_mix2_off (max_dom) = .false. ! if set to T, deactivate 2nd-order horizontal mixing for chem species. default is F.
1646 tracer_mix2_off (max_dom) = .false. ! if set to T, deactivate 2nd-order horizontal mixing for tracers. default is F.
1647 scalar_mix2_off (max_dom) = .false. ! if set to T, deactivate 2nd-order horizontal mixing for scalars. default is F.
1648 tke_mix2_off (max_dom) = .false. ! if set to T, deactivate 2nd-order horizontal mixing for tke. default is F.
1649 moist_mix6_off (max_dom) = .false. ! if set to T, deactivate 6th-order horizontal mixing for moisture. default is F.
1650 chem_mix6_off (max_dom) = .false. ! if set to T, deactivate 6th-order horizontal mixing for chem species. default is F.
1651 tracer_mix6_off (max_dom) = .false. ! if set to T, deactivate 6th-order horizontal mixing for tracers. default is F.
1652 scalar_mix6_off (max_dom) = .false. ! if set to T, deactivate 6th-order horizontal mixing for scalars. default is F.
1653 tke_mix6_off (max_dom) = .false. ! if set to T, deactivate 6th-order horizontal mixing for tke. default is F.
1656 time_step_sound (max_dom) = 4 / ! number of sound steps per time-step (0=set automatically)
1657 (if using a time_step much larger than 6*dx (in km),
1658 proportionally increase number of sound steps - also
1659 best to use even numbers)
1660 do_avgflx_em (max_dom) = 0, ! whether to output time-averaged mass-coupled advective velocities
1661 0 = no (default)
1662 1 = yes
1663 do_avgflx_cugd (max_dom) = 0, ! whether to output time-averaged convective mass-fluxes from Grell-Devenyi ensemble scheme
1664 0 = no (default)
1665 1 = yes (only takes effect if do_avgflx_em=1 and cu_physics= 93
1666 do_coriolis (max_dom) = .true., ! whether to do Coriolis calculations (idealized) (inactive)
1667 do_curvature (max_dom) = .true., ! whether to do curvature calculations (idealized) (inactive)
1668 do_gradp (max_dom) = .true., ! whether to do horizontal pressure gradient calculations (idealized) (inactive)
1669 fft_filter_lat = 91. ! the latitude above which the polar filter is turned on (degrees) - 45 degrees is a
1670 reasonable latitude to start using polar filters
1671 coupled_filtering = .true. ! T/F mu coupled scalar arrays are run through the polar filters
1672 pos_def = .false. ! T/F remove negative values of scalar arrays by setting minimum value to zero
1673 swap_pole_with_next_j = .false. ! T/F replace the entire j=1 (jds-1) with the values from j=2 (jds-2)
1674 actual_distance_average = .false. ! T/F average the field at each i location in the j-loop with a number of grid points based on a map-factor ratio
1675 gwd_opt (max_dom) = 0 ! for running without gravity wave drag
1676 = 1 ; for running with gravity wave drag (Choi and Hong 2015)
1677 = 3 ; NOAA/GSL gravity wave drag and turbulent orographic form drag
1678 gwd_diags = 0 ; set to 1 to output diagnostics for gwd_opt = 3
1679 sfs_opt (max_dom) = 0 ! nonlinear backscatter and anisotropy (NBA) off
1680 = 1 ! NBA1 using diagnostic stress terms (km_opt=2,3 for scalars)
1681 = 2 ! NBA2 using tke-based stress terms (km_opt=2 needed)
1682 m_opt (max_dom) = 0 ! no added output
1683 = 1 ! adds output of Mij stress terms when NBA is not used
1684 tracer_opt(max_dom) = 0 ! whether to turn on tracers
1685 rad_nudge = 1 ! nudge T to initial values in idealized tropical cyclone case
1687 &bdy_control
1688 spec_bdy_width = 5, ! total number of rows for specified boundary value nudging
1689 spec_zone = 1, ! number of points in specified zone (spec b.c. option)
1690 relax_zone = 4, ! number of points in relaxation zone (spec b.c. option)
1691 specified (max_dom) = .false., ! specified boundary conditions (only can be used for domain 1)
1692 the above 4 are used for real-data runs
1693 spec_exp = 0. ! exponential multiplier for relaxation zone ramp for specified=.t.
1694 (0.=linear ramp default, e.g. 0.33=~3*dx exp decay factor)
1695 constant_bc = .false. ! constant boundary condition used with DFI
1697 periodic_x (max_dom) = .false., ! periodic boundary conditions in x direction
1698 symmetric_xs (max_dom) = .false., ! symmetric boundary conditions at x start (west)
1699 symmetric_xe (max_dom) = .false., ! symmetric boundary conditions at x end (east)
1700 open_xs (max_dom) = .false., ! open boundary conditions at x start (west)
1701 open_xe (max_dom) = .false., ! open boundary conditions at x end (east)
1702 periodic_y (max_dom) = .false., ! periodic boundary conditions in y direction
1703 symmetric_ys (max_dom) = .false., ! symmetric boundary conditions at y start (south)
1704 symmetric_ye (max_dom) = .false., ! symmetric boundary conditions at y end (north)
1705 open_ys (max_dom) = .false., ! open boundary conditions at y start (south)
1706 open_ye (max_dom) = .false., ! open boundary conditions at y end (north)
1707 nested (max_dom) = .false., ! nested boundary conditions (must be used for nests)
1708 polar (max_dom) = .false., ! polar boundary condition
1709 (v=0 at polarward-most v-point)
1710 have_bcs_moist (max_dom) = .false., ! model run after ndown only: do not use microphysics variables in bdy file
1711 = .true. , ! use microphysics variables in bdy file
1712 have_bcs_scalar (max_dom) = .false., ! model run after ndown only: do not use scalar variables in bdy file
1713 = .true. , ! use scalar variables in bdy file
1714 multi_bdy_files = .false., ! F=default; T=model will run with split and individually named LBC files
1715 requires usage of bdy_inname = "wrfbdy_d<domain>_<date>"
1716 the real program can generate all LBC files in a single run
1717 only a single time period is in each separate LBC file
1719 &ideal
1720 ideal_case = 0 ! ideal case number. Option selects CASEs within module_initialize_ideal.F
1721 realcase ideal_case=0
1722 hill2d_x ideal_case=1
1723 quarter_ss ideal_case=2
1724 convrad ideal_case=3
1725 squall2d_x ideal_case=4
1726 squall2d_y ideal_case=5
1727 grav2d_x ideal_case=6
1728 b_wave ideal_case=7
1729 seabreeze2d_x ideal_case=8
1730 les ideal_case=9
1732 &tc ; controls for tc_em.exe ONLY, no impact on real, ndown, or model
1734 insert_bogus_storm = .false. ! T/F for inserting a bogus tropical storm (TC)
1735 remove_storm = .false. ! T/F for only removing the original TC
1736 num_storm = 1 ! Number of bogus TC
1737 latc_loc = -999. ! center latitude of the bogus TC
1738 lonc_loc = -999. ! center longitude of the bogus TC
1739 vmax_meters_per_second(max_bogus) = -999. ! vmax of bogus storm in meters per second
1740 rmax = -999. ! maximum radius outward from storm center
1741 vmax_ratio(max_bogus) = -999. ! ratio for representative maximum winds, 0.75 for 45 km grid, and
1742 0.9 for 15 km grid.
1743 rankine_lid = -999. ! top pressure limit for the tc bogus scheme
1745 &namelist_quilt This namelist record controls asynchronized I/O for MPI applications.
1747 nio_tasks_per_group = 0, ! default value is 0: no quilting; > 0 quilting I/O
1748 nio_groups = 1, ! default 1. May be set to higher value for nesting IO
1749 or history and restart IO
1752 &grib2:
1753 background_proc_id = 255, ! Background generating process identifier, typically defined
1754 by the originating center to identify the background data that
1755 was used in creating the data. This is octet 13 of Section 4
1756 in the grib2 message
1757 forecast_proc_id = 255, ! Analysis or generating forecast process identifier, typically
1758 defined by the originating center to identify the forecast process
1759 that was used to generate the data. This is octet 14 of Section
1760 4 in the grib2 message
1761 production_status = 255, ! Production status of processed data in the grib2 message.
1762 See Code Table 1.3 of the grib2 manual. This is octet 20 of
1763 Section 1 in the grib2 record
1764 compression = 40, ! The compression method to encode the output grib2 message.
1765 Only 40 for jpeg2000 or 41 for PNG are supported
1768 By default the pressure and height level data goes into stream 23 and 22, respectively. Using
1769 the vertical interpolation options requires the user to define an io_form and interval for
1770 the requested stream. See examples.namelist.
1772 &diags:
1773 p_lev_diags = 1, ! Vertically interpolate diagnostics to p-levels
1774 0=NO, 1=YES
1775 num_press_levels = 0, ! Number of pressure levels to interpolate to, for example,
1776 could be 2
1777 press_levels = 0, ! Which pressure levels (Pa) to interpolate to, for example
1778 could be 85000, 70000
1779 use_tot_or_hyd_p = 2 ! Which half level pressure to use: 1=total (p+pb); 2=hydrostatic
1780 (p_hyd). The p_hyd option is the default and less noisy. Total
1781 pressure is consistent with what is done in various post-proc
1782 packages.
1783 z_lev_diags = 0, ! Vertically interpolate diagnostics to z-levels
1784 0=NO, 1=YES
1785 num_z_levels = 2, ! Number of height levels to interpolate to
1786 z_levels = 0, ! List of height values (m) to interpolate data to.
1787 Positive numbers are for height above mean sea level (i.e. a flight level)
1788 Negative numbers are for levels above ground
1789 /
1791 AFWA diagnostics:
1792 &afwa
1793 afwa_diag_opt (max_dom) = 0, ! AFWA Diagnostic option, 1: on
1794 afwa_ptype_opt (max_dom) = 0, ! Precip type option, 1: on
1795 afwa_vil_opt (max_dom) = 0, ! Vert Int Liquid option, 1: on
1796 afwa_radar_opt (max_dom) = 0, ! Radar option, 1: on
1797 afwa_severe_opt (max_dom) = 0, ! Severe Wx option, 1: on
1798 afwa_icing_opt (max_dom) = 0, ! Icing option, 1: on
1799 afwa_vis_opt (max_dom) = 0, ! Visibility option, 1: on
1800 afwa_cloud_opt (max_dom) = 0, ! Cloud option, 1: on
1801 afwa_therm_opt (max_dom) = 0, ! Thermal indices option, 1: on
1802 afwa_turb_opt (max_dom) = 0, ! Turbulence option, 1: on
1803 afwa_buoy_opt (max_dom) = 0, ! Buoyancy option, 1: on
1804 afwa_ptype_ccn_tmp = 264.15, ! CCN temperature for precipitation type calculation
1805 afwa_ptype_tot_melt = 50, ! Total melting energy for precipitation type calculation
1810 Add an extra set of 3d arrays for vertical interpolation to the
1811 processing for the real program. Typically, this extra data set
1812 is to include monthly aerosol data. The vertical coordinate of the
1813 aerosol data is able to be separate from the input meteorological
1814 data. To introduce new data sets, mods are required in the Registry
1815 and in module_initialize_real.F. There is a space-holder/practice
1816 set-up for "GCA". The actual data set for Thompson mp=28 (WIF) that
1817 utilizes QNWFA and QNIFA (water and ice friendly aerosols) has
1818 been tested.
1820 &domains
1821 num_wif_levels = 30
1822 wif_input_opt = 1
1823 /
1824 num_gca_levels = 13
1825 gca_input_opt = 1
1828 physics suite specification, which specifies physics options
1830 mp_physics
1831 cu_physics
1832 ra_lw_physics
1833 ra_sw_physics
1834 bl_pbl_physics
1835 sf_sfclay_physics
1836 sf_surface_physics
1838 with a new namelist physics_suite = 'X'. Two suites are available:
1840 physics_suite = 'CONUS'
1842 or
1844 physics_suite = 'tropical'
1846 where 'CONUS' is equivalent to
1848 mp_physics = 8,
1849 cu_physics = 6,
1850 ra_lw_physics = 4,
1851 ra_sw_physics = 4,
1852 bl_pbl_physics = 2,
1853 sf_sfclay_physics = 2,
1854 sf_surface_physics = 2,
1856 and 'tropical' is equivalent to
1858 mp_physics = 6,
1859 cu_physics = 16,
1860 ra_lw_physics = 4,
1861 ra_sw_physics = 4,
1862 bl_pbl_physics = 1,
1863 sf_sfclay_physics = 91,
1864 sf_surface_physics = 2,
1866 One can use physics_suite, and overwrite one or more options by explicitly including the physics
1867 namelists.
1869 To overwrite an option for a nest, one can have
1871 &physics
1872 physics_suite = 'tropical'
1873 cu_physics = -1, -1, 0,
1874 ...
1876 here '-1' means to use option specified by the suite, and '0' modifies the cumulus option from the
1877 suite option 16 to 0 (turning cumulus off).
1880 Hybrid Vertical Coordinate (HVC) vs Terrain Following (TF)
1882 1. To turn ON the HVC
1883 &dynamics
1884 hybrid_opt = 2,
1886 2. To turn OFF the HVC
1887 &dynamics
1888 hybrid_opt = 0,
1892 Traditional / full fields model output
1894 The WRF model output may be modified at run-time to also include a stream that contains
1895 more traditional fields: temperature, pressure, geopotential height, etc. The flag needs to
1896 be activated in the namelist (diag_nwp2 = 1). In the registry.trad_fields, the output stream
1897 is "h1", so that stream needs to be explicitly requested in the namelist.input file.
1899 &diags
1900 diag_nwp2 = 1
1901 /
1903 &time_control
1904 io_form_auxhist1 = 2
1905 auxhist1_interval = 180, 30, 30,
1906 frames_per_auxhist1 = 1, 1, 1,
1907 auxhist1_outname = "wrf_trad_fields_d<domain>_<date>"
1908 /
1912 Solar diagnostics output
1914 For solar diagnostics, 2-D fields of variables relevant to solar forecasting are output
1915 when option solar_diagnostics is activated in diags section of namelist, as shown below.
1916 Also, if tslist is present when solar_diagnostics is activated, then these same variables
1917 are output to the time series files for the location(s) specified in tslist (see README.tslist).
1918 All variables are calculated in phys/module_diag_solar.F and defined in registry.solar_fields
1920 &diags
1921 solar_diagnostics = 1
1922 /