| blob | 9dc6b8e87d066833e9bf5c3f081705ec8d02ceff |
1 Note, this is not a namelist.input file. Find what interests you, and cut and paste them
2 to your own namelist.input file. For more information on these namelist parameters,
3 please see run/README.namelist or Chapter 5 of the User's Guide.
5 ** More options for real in namelist record &domains:
7 p_top_requested = 5000
8 interp_type = 2
9 extrap_type = 2
10 t_extrap_type = 2
11 lowest_lev_from_sfc = .false.
12 use_levels_below_ground = .true.
13 use_surface = .true.
14 lagrange_order = 2
15 force_sfc_in_vinterp = 1
16 zap_close_levels = 500
17 sfcp_to_sfcp = .false.
18 adjust_heights = .false.
19 smooth_cg_topo = .false.
20 eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946,
21 0.922, 0.894, 0.860, 0.817, 0.766,
22 0.707, 0.644, 0.576, 0.507, 0.444,
23 0.380, 0.324, 0.273, 0.228, 0.188,
24 0.152, 0.121, 0.093, 0.069, 0.048,
25 0.029, 0.014, 0.000,
27 ** Using sst_update option (add these to namelist records &time_control and
28 &physics respectively):
30 &time_control
31 auxinput4_inname = "wrflowinp_d<domain>"
32 auxinput4_interval = 360, 360, 360,
33 io_form_auxinput4 = 2
35 &physics
36 sst_update = 1,
38 ** Using qna_update option (add these to namelist records &time_control and
39 &physics respectively):
41 &time_control
42 auxinput17_inname = "wrfqnainp_d<domain>"
43 auxinput17_interval = 360, 360, 360,
44 io_form_auxinput17 = 2
46 &physics
47 qna_update = 1,
49 ** Using Noah-MP option (Use sf_surface_physics option = 4, and add
50 &noah_mp namelist record)
52 &physics
53 sf_surface_physics = 4
55 &noah_mp
56 dveg = 4,
57 opt_crs = 1,
58 opt_btr = 1,
59 opt_run = 1,
60 opt_sfc = 1,
61 opt_frz = 1,
62 opt_inf = 1,
63 opt_rad = 3,
64 opt_alb = 2,
65 opt_snf = 1,
66 opt_tbot = 2,
67 opt_stc = 1,
68 /
70 ** Using UCM urban model
71 &physics
72 sf_urban_physics = 1, 1, 1,
74 ** Using BEP urban model
75 &physics
76 sf_urban_physics = 2, 2, 2,
78 ** Using BEM urban model
79 &physics
80 sf_urban_physics = 3, 3, 3,
82 ** Using lake model
83 &physics
84 sf_lake_physics = 1, 1, 1,
85 lakedepth_default = 50., 50., 50.,
86 lake_min_elev = 5., 5., 5.,
88 ** Using shallow water roughness (constant depth)
89 &physics
90 shalwater_z0 = 1, 1, 1,
91 shalwater_depth = 40.0,
93 ** Using shallow water roughness (with bathymetry data)
94 &physics
95 shalwater_z0 = 1, 1, 1,
96 shalwater_depth = 0.0,
98 ** Using stochastic backscatter scheme (new namelist record since v3.6)
100 &stoch
101 stoch_force_opt = 1, 1, 1,
102 stoch_vertstruc_opt = 0, 0, 0,
103 tot_backscat_psi = 1.E-05, 1.E-05, 1.E-05,
104 tot_backscat_t = 1.E-06, 1.E-06, 1.E-06,
105 nens = 1,
106 ztau_psi = 10800.0,
107 ztau_t = 10800.0,
108 rexponent_psi =-1.83,
109 rexponent_t =-1.83,
110 zsigma2_eps = 0.0833,
111 zsigma2_eta = 0.0833,
112 kminforc = 1,
113 kminforc = 1,
114 kminforct = 1,
115 lminforct = 1,
116 kmaxforc = 1000000,
117 lmaxforc = 1000000,
118 kmaxforct = 1000000,
119 lmaxforct = 1000000,
120 perturb_bdy = 0,
123 ** Using DFI options (note this is a separate namelist record):
125 &dfi_control
126 dfi_opt = 3,
127 dfi_nfilter = 7,
128 dfi_cutoff_seconds = 3600,
129 dfi_write_filtered_input = .true.
130 dfi_write_dfi_history = .false.
131 dfi_bckstop_year = 2000,
132 dfi_bckstop_month = 01,
133 dfi_bckstop_day = 24,
134 dfi_bckstop_hour = 10,
135 dfi_bckstop_minute = 00,
136 dfi_bckstop_second = 00,
137 dfi_fwdstop_year = 2000,
138 dfi_fwdstop_month = 01,
139 dfi_fwdstop_day = 24,
140 dfi_fwdstop_hour = 13,
141 dfi_fwdstop_minute = 00,
142 dfi_fwdstop_second = 00,
143 /
145 &domains
146 time_step_dfi = 60
148 ** Using gridded nudging option (note this is a separate namelist record) for
149 upperair nudging
151 ** Upper air gridded nudging requires an file generated by the real
152 program. Activating grid_fdda for the real program is adequate.
154 &fdda
155 grid_fdda = 1, 1, 1,
156 gfdda_inname = "wrffdda_d<domain>",
157 gfdda_end_h = 24, 24, 24,
158 gfdda_interval_m = 360, 360, 360,
159 fgdt = 0, 0, 0,
160 if_no_pbl_nudging_uv = 0, 0, 0,
161 if_no_pbl_nudging_t = 1, 1, 1,
162 if_no_pbl_nudging_q = 1, 1, 1,
163 if_zfac_uv = 0, 0, 0,
164 k_zfac_uv = 10, 10, 10,
165 if_zfac_t = 0, 0, 0,
166 k_zfac_t = 10, 10, 10,
167 if_zfac_q = 0, 0, 0,
168 k_zfac_q = 10, 10, 10,
169 guv = 0.0003, 0.0003, 0.0003,
170 gt = 0.0003, 0.0003, 0.0003,
171 gq = 0.0003, 0.0003, 0.0003,
172 if_ramping = 1,
173 dtramp_min = 60.0,
174 io_form_gfdda = 2,
176 ** Using gridded surface nudging option (note this is a separate namelist
177 record)
179 ** Note that upper-air and surface gridded nudging may be used together
180 or separately. Surface nudging requires an input file generated by the
181 obsgrid program.
183 &fdda
184 grid_sfdda = 1, 1, 1,
185 sgfdda_inname = "wrfsfdda_d<domain>",
186 sgfdda_end_h = 24, 24, 24,
187 sgfdda_interval_m = 360, 360, 360,
188 io_form_sgfdda = 2,
189 guv_sfc = 0.0003, 0.0003, 0.0003,
190 gt_sfc = 0.0003, 0.0003, 0.0003,
191 gq_sfc = 0.0003, 0.0003, 0.0003,
192 rinblw = 250., 250., 250.,
194 ** Using observation nudging option (note &fdda is a separate namelist record):
196 &time_control
197 auxinput11_interval_s = 180 , 180 , 180
198 auxinput11_end_h = 6 , 6 , 6
200 &fdda
201 obs_nudge_opt = 1,1,1,
202 max_obs = 150000,
203 fdda_start = 0., 0., 0.,
204 fdda_end = 720.,720.,720.,
205 obs_nudge_wind = 1,1,1,
206 obs_coef_wind = 6.E-4,6.E-4,6.E-4,
207 obs_nudge_temp = 1,1,1,
208 obs_coef_temp = 6.E-4,6.E-4,6.E-4,
209 obs_nudge_mois = 1,1,1,
210 obs_coef_mois = 6.E-4,6.E-4,6.E-4,
211 obs_rinxy = 240.,240.,180.,
212 obs_rinsig = 0.1,
213 obs_twindo = 0.6666667,0.6666667,0.6666667,
214 obs_npfi = 10,
215 obs_ionf = 2, 2, 2,
216 obs_idynin = 0,
217 obs_dtramp = 40.,
218 obs_prt_freq = 10, 10, 10,
219 obs_prt_max = 10
220 obs_ipf_errob = .true.
221 obs_ipf_nudob = .true.
222 obs_ipf_in4dob = .true
223 obs_no_pbl_nudge_uv = 0
224 obs_no_pbl_nudge_t = 0
225 obs_no_pbl_nudge_q = 0
226 obs_sfc_scheme_horiz = 0
227 obs_sfc_scheme_vert = 0
228 obs_max_sndng_gap = 20
229 obs_nudgezfullr1_uv = 50
230 obs_nudgezrampr1_uv = 50
231 obs_nudgezfullr2_uv = 50
232 obs_nudgezrampr2_uv = 50
233 obs_nudgezfullr4_uv = -5000
234 obs_nudgezrampr4_uv = 50
235 obs_nudgezfullr1_t = 50
236 obs_nudgezrampr1_t = 50
237 obs_nudgezfullr2_t = 50
238 obs_nudgezrampr2_t = 50
239 obs_nudgezfullr4_t = -5000
240 obs_nudgezrampr4_t = 50
241 obs_nudgezfullr1_q = 50
242 obs_nudgezrampr1_q = 50
243 obs_nudgezfullr2_q = 50
244 obs_nudgezrampr2_q = 50
245 obs_nudgezfullr4_q = -5000
246 obs_nudgezrampr4_q = 50
247 obs_nudgezfullmin = 50
248 obs_nudgezrampmin = 50
249 obs_nudgezmax = 3000
250 obs_sfcfact = 1.0
251 obs_sfcfacr = 1.0
252 obs_dpsmx = 7.5
253 /
255 ** Using spectral nudging option
257 &fdda
258 grid_fdda = 2, 2, 2,
259 gfdda_inname = "wrffdda_d<domain>",
260 gfdda_end_h = 24, 24, 24,
261 gfdda_interval_m = 360, 360, 360,
262 fgdt = 0, 0, 0,
263 fgdtzero = 0, 0, 0,
264 if_no_pbl_nudging_uv = 0, 0, 0,
265 if_no_pbl_nudging_t = 0, 0, 0,
266 if_no_pbl_nudging_ph = 0, 0, 0,
267 if_zfac_uv = 0, 0, 0,
268 k_zfac_uv = 10, 10, 10,
269 if_zfac_t = 0, 0, 0,
270 k_zfac_t = 10, 10, 10,
271 if_zfac_ph = 0, 0, 0,
272 k_zfac_ph = 10, 10, 10,
273 dk_zfac_uv = 1, 1, 1,
274 dk_zfac_t = 1, 1, 1,
275 dk_zfac_ph = 1, 1, 1,
276 guv = 0.0003, 0.0003, 0.0003,
277 gt = 0.0003, 0.0003, 0.0003,
278 gph = 0.0003, 0.0003, 0.0003,
279 xwavenum = 3,
280 ywavenum = 3,
281 if_ramping = 1,
282 dtramp_min = 60.0,
283 io_form_gfdda = 2,
285 ** Using adaptive time step option (add these in namelist record &domains):
287 use_adaptive_time_step = .true.,
288 step_to_output_time = .true.,
289 target_cfl = 1.2, 1.2, 1.2,
290 target_hcfl = .84, .84, .84,
291 max_step_increase_pct = 5, 51, 51,
292 starting_time_step = -1, -1, -1,
293 max_time_step = 360, 120, 40,
294 min_time_step = 90, 30, 10,
295 adaptation_domain = 1,
297 ** Using automatic vortex-following option (tropical storm tracking only;
298 add these in namelist record &domains):
300 vortex_interval = 15, 15, 15,
301 max_vortex_speed = 40, 40, 40,
302 corral_dist = 8, 15, 15,
303 track_level = 50000,
305 ** Miscellaneous physics options for namelist record &physics:
307 Topographic shading (only effective when grid sizes are a few kilometers)
309 slope_rad = 1, 1, 1,
310 topo_shading = 1, 1, 1,
311 shadlen = 25000,
313 Setting threshold value for defining seaice if seaice is not in the input file:
315 seaice_threshold = 271,
317 Switching off latent heating from a microphysics scheme (must also set cu_physics = 0):
319 no_mp_heating = 0,
321 ** Using precipiatiion bucket in a time interval (minutes):
323 prec_acc_dt = 60.
325 ** Using bucket accumulations for multi-year simulations (guideline: mean monthly accumulation)
327 bucket_mm = 100.
328 bucket_J = 1.e9
330 ** Optional gravitational settling of fog/cloud droplets (MYNN PBL only)
332 grav_settling = 1, ; default 0
334 ** Using implicit gravity-wave damping option (add these in namelist record &dynamics):
336 damp_opt = 3,
337 zdamp = 5000., 5000., 5000.,
338 dampcoef = 0.2, 0.2, 0.2
340 ** Using expanded boundary zone and exponential decay option (add or modify these in
341 namelist record &bdy_control).
342 spec_zone is ALWAYS = 1
343 relax_zone is ALWAYS = spec_bdy_width - spec_zone
345 &bdy_control
346 spec_bdy_width = 10,
347 specified = .true.,
348 spec_exp = 0.33
350 For a tropical channel configuration, set the following:
352 specified = .true.,
353 periodic_x = .true.,
355 ** Using split lateral boundary files
356 The run-time flag multi_bdy_files must be set to TRUE (default is false), and the
357 lateral boundary files must have a date associated. When using the split LBC option,
358 there is ALWAYS and ONLY a single time LBC time in each file.
360 multi_bdy_files = .true.
362 Also requires
363 &time_control
364 bdy_inname = "wrfbdy_d<domain>_<date>"
366 ** using io quilting option to improve output efficiency for large domain runs
367 (note that this is a separate namelist record):
369 &namelist_quilt
370 nio_tasks_per_group = 2,
371 nio_groups = 1,
373 ** for tc bogusing:
374 Helpful hints
375 1) The TC bogus *MUST* only be run with a single processor: serial
376 (with or without a nest build option), or as a DM build with np=1.
377 2) Run the TC program for only the initial model time. Remember to copy
378 your metgrid file to a safe location and then let the TC scheme
379 generate the auxiliary files.
380 3) The TC program does not handle soil moisture or soil temperature
381 correctly. After running the TC program, use an NCL script to put the
382 modified fields back into the original metgrid file (remember that you
383 are keeping a pristine copy elsewhere). The modified fields to copy
384 are: "RH", "TT", "UU", "VV", "GHT", "PRES", "PMSL", "PSFC".
385 4) The TC program runs quickly, only a few seconds. It is a good idea to
386 process multiple runs, and vary the TC initialization. It would be
387 reasonable to select the maximum wind speed (vmax_meters_per_second)
388 from 30 to 60 by 5, select the radius of maximum winds ((m), rmax) =
389 50000 to 200000 by 50000, and select a tuning parameter for the
390 vortex (vmax_ratio) = 0.5 to 0.9 by 0.1. Look at the resulting
391 bogus storms (sea level pressure, surface wind speeds, and
392 overall storm size), and choose the best.
393 5) Remember to also consider starting the TC program at a later or
394 earlier time during your parameter space exploration.
395 6) The bogus storm will organize better if the initialization of the
396 tropical storm is entirely over water, and if the storm is able to
397 develop over water for several days.
398 7) After the metgrid file is finalized, and the real program is run,
399 use DFI in WRF. Removing and then introducing an entire bogus typhoon
400 is the definition of unbalanced. Running DFI for 2-3 hours back and forth
401 will be a benefit. For a successful implementation of a TC bogus storm in
402 WRF, the storm should not oscillate in size, should not be radiating
403 massive amounts of gravity waves, and should not rapidly weaken.
405 &tc
406 insert_bogus_storm = .true.,
407 remove_storm = .false.,
408 num_storm = 1,
409 latc_loc = 15.,
410 lonc_loc = -90.,
411 vmax_meters_per_second = 30,
412 rmax = 50000,
413 vmax_ratio = 0.5,
414 /
416 ** for regional climate surface diagnostics such as max/min/mean/std of T2/Q2/wind/rainfall
417 between selected output times (e.g. daily) in auxhist3
419 &time_control
420 output_diagnostics = 1
421 auxhist3_outname = 'wrfxtrm_d<domain>_<date>'
422 io_form_auxhist3 = 2
423 auxhist3_interval = 1440
424 frames_per_auxhist3 = 1
425 /
428 ** for pressure-level (and some surface) diagnostics, output is on stream 23,
429 where the listed pressure levels are in Pa. For the height level interpolation,
430 the unit is stream 22.
432 &time_control
433 io_form_auxhist23 = 2,
434 auxhist23_interval = 30, 30, 30,
435 frames_per_auxhist23 = 1, 1, 1,
436 auxhist23_outname = "PLEVS_d<domain>_<date>"
437 /
439 &diags
440 p_lev_diags = 1
441 num_press_levels = 4
442 press_levels = 85000, 70000, 50000, 25000
443 use_tot_or_hyd_p = 2
444 /
447 ** for height-level diagnostics, output is on stream 22 (negative values
448 for z_levels means AGL, so -500 is 500 m AGL, and 10000 is 10 km).
450 &time_control
451 io_form_auxhist22 = 2,
452 auxhist22_interval = 30, 30, 30,
453 frames_per_auxhist22 = 1, 1, 1,
454 auxhist22_outname = "ZLEVS_d<domain>_<date>"
455 /
457 &diags
458 z_lev_diags = 1
459 num_z_levels = 2
460 z_levels = -500, 10000
461 /
464 ** For solar diagnostics, 2-D fields of variables relevant to solar forecasting are output
465 when option solar_diagnostics is activated in diags section of namelist, as shown below.
466 Also, if tslist is present when solar_diagnostics is activated, then these same variables
467 are output to the time series files for the location(s) specified in tslist.
468 All variables are calculated in phys/module_diag_solar.F and defined in registry.solar_fields
470 &diags
471 solar_diagnostics = 1
472 /
475 ** Using different flux formulation for tropical storm simulations
476 (best for grid spacing less than 2 km)
477 simple 1-D ocean mixed layer, or University of Miami 3DPWP ocean model
478 (add these in namelist record &physics):
480 isftcflx = 1,
481 sf_ocean_physics = 0 (off), 1 (OML), 2 (3D PWP)
482 oml_hml0 = 50,
483 oml_gamma = 0.14
484 omdt = 1
486 Add these in the &domains section. Note that the example here is for
487 a warm ocean: ocean_z is the depth of each layer (m), from the surface
488 downward; the associated temperature is ocean_t (K); and the associated
489 salinity is ocean_s (ppt).
491 ocean_levels = 30,
492 ocean_z = 5., 15., 25., 35., 45., 55.,
493 65., 75., 85., 95., 105., 115.,
494 125., 135., 145., 155., 165., 175.,
495 185., 195., 210., 230., 250., 270.,
496 290., 310., 330., 350., 370., 390.
497 ocean_t = 302.3493, 302.3493, 302.3493, 302.1055, 301.9763, 301.6818,
498 301.2220, 300.7531, 300.1200, 299.4778, 298.7443, 297.9194,
499 297.0883, 296.1443, 295.1941, 294.1979, 293.1558, 292.1136,
500 291.0714, 290.0293, 288.7377, 287.1967, 285.6557, 284.8503,
501 284.0450, 283.4316, 283.0102, 282.5888, 282.1674, 281.7461
502 ocean_s = 34.0127, 34.0127, 34.0127, 34.3217, 34.2624, 34.2632,
503 34.3240, 34.3824, 34.3980, 34.4113, 34.4220, 34.4303,
504 34.6173, 34.6409, 34.6535, 34.6550, 34.6565, 34.6527,
505 34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247,
506 34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247
509 References for 3D ocean circulation model:
510 Lee, C.-Y., and S. S. Chen, 2013: Stable Boundary Layer and Its Impact on Tropical Cyclone Structure in a Coupled Atmosphere-Ocean Model, Mon. Wea. Rev.
511 Price, J. F., T. B. Sanford, and G. Z. Forristal, 1994: Forced stage response to a moving hurricane. J. Phy. Oceanogr., 24, 233-260.
515 ** Using the ACOM Forward Lagrangian trajectory calculation:
517 &domains
518 num_traj = 25,
520 &physics
521 traj_opt = 1,
522 dm_has_traj = .true., ..true., .true.
525 For domain #1, the file must "wrfinput_traj_d01" exist in the working directory. Similarly for domain 2, 3, etc. Each domain
526 has a separate file for a namelist.
528 &traj_default
529 traj_def%start_time = '2000-01-24_12:00:00',
530 traj_def%stop_time = '2000-01-25_12:00:00',
531 traj_def%dyn_name(1:6) = 'p', 'T', 'z', 'u', 'v', 'w',
532 traj_def%hyd_name(1) = 'QVAPOR',
533 /
535 &traj_spec
536 traj_type%start_time = '2000-01-24_12:00:00', '2000-01-24_12:00:00',
537 '2000-01-24_12:00:00', '2000-01-24_12:00:00',
538 '2000-01-24_12:00:00', '2000-01-24_12:00:00',
539 '2000-01-24_12:00:00', '2000-01-24_12:00:00',
540 '2000-01-24_12:00:00', '2000-01-24_12:00:00',
541 '2000-01-24_12:00:00',
542 traj_type%stop_time = '2000-01-25_12:00:00', '2000-01-25_12:00:00',
543 '2000-01-25_12:00:00', '2000-01-25_12:00:00',
544 '2000-01-25_12:00:00', '2000-01-25_12:00:00',
545 '2000-01-25_12:00:00', '2000-01-25_12:00:00',
546 '2000-01-25_12:00:00', '2000-01-25_12:00:00',
547 '2000-01-25_12:00:00',
548 traj_type%lev = 60., 60., 60., 60., 60., 60., 60., 60., 60., 60., 60.,
549 traj_type%lon = -79.88470, -79.74551, -79.60422, -79.46072,
550 -79.31503, -79.16708, -79.01682, -78.86417,
551 -78.70911, -78.55151, -78.39142,
552 traj_type%lat = 29.18063, 29.70515, 30.23069, 30.75718,
553 31.28461, 31.81292, 32.34208, 32.87204,
554 33.40276, 33.93421, 34.46631,
555 /
559 ** Vertical nesting
561 The WRF model now supports vertical nesting for a coincident (online)
562 model simulation (during a single model run, differing numbers of
563 vertical levels may be used per domain). This is activated with a switch
564 to turn on the option (vert_refine_method). The namelist array eta_levels
565 is manually filled in for each domain. Below is an example for two
566 domains.
568 NOTE: The user is restricted to using the RRTM LW and Dudhia SW radiation schemes.
569 NOTE: The user is restricted from using the hybrid vertical coordinate.
571 &domains
572 max_dom = 2,
573 e_vert = 35, 45,
574 eta_levels(1:35) = 1., 0.993, 0.983, 0.97, 0.954, 0.934, 0.909, 0.88, 0.8406663, 0.8013327,
575 0.761999, 0.7226653, 0.6525755, 0.5877361, 0.5278192, 0.472514,
576 0.4215262, 0.3745775, 0.3314044, 0.2917579, 0.2554026, 0.2221162,
577 0.1916888, 0.1639222, 0.1386297, 0.1156351, 0.09525016, 0.07733481,
578 0.06158983, 0.04775231, 0.03559115, 0.02490328, 0.0155102, 0.007255059, 0.
579 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,
580 0.8210, 0.7906, 0.7602, 0.7298, 0.6812, 0.6290, 0.5796, 0.5333, 0.4901, 0.4493, 0.4109,
581 0.3746, 0.3412, 0.3098, 0.2802, 0.2524, 0.2267, 0.2028, 0.1803, 0.1593, 0.1398, 0.1219,
582 0.1054, 0.0904, 0.0766, 0.0645, 0.0534, 0.0433, 0.0341, 0.0259, 0.0185, 0.0118, 0.0056, 0.
583 vert_refine_method = 0, 2,
587 ** Tropopause data level of max winds data, for program real.exe only
589 When information (mostly NCEP supplied through GFS or NAM) is available
590 in the Grib2 file, AND extracted with ungrib, the metgrid program inserts
591 flag information into the data stream that is input by real. The real
592 program is able to use the available u, v, T, height fields (each may be on
593 a tropopause .OR. the level of max winds) in the vertical interpolation.
594 To "tune" the data, the user may select a level below which the max wind
595 inforamtion is ignored. The default is 300 hPa. The user may also
596 select the level that when exceeded the trop and maxw fields are ignored
597 (due to the horizontal pressure difference detecting a user-defined
598 discontinuity, and the metgrid horizontal interpolation across the
599 discontinuity would be suspect). By default, this option is turned off.
600 The user may choose to individually activate the vertical interpolation of
601 either the level of max winds .OR. the tropopause level.
603 Default (all pressure units Pa):
604 &domains
605 maxw_horiz_pres_diff = 5000
606 trop_horiz_pres_diff = 5000
607 maxw_above_this_level = 30000
608 use_maxw_level = 0 (0=do not use level, 1 = use level)
609 use_trop_level = 0 (0=do not use level, 1 = use level)
614 ** Using aerosol option aer_opt = 2:
616 &physics
617 aer_opt = 2,
618 aer_type = 1,
619 aer_aod550_opt = 1,
620 aer_aod550_val = 0.12,
621 aer_angexp_opt = 1,
622 aer_angexp_val = 1.3,
623 aer_ssa_opt = 1,
624 aer_ssa_val = 0.85,
625 aer_asy_opt = 1,
626 aer_asy_val = 0.90,
629 ** Using Jimenez wind-farm scheme
631 In the &physics namelist record, set the MAX_DOM value:
632 windfarm = 1, 1, 1
634 Also in the directory with namelist.input, is the file defining the
635 specifics of the wind turbine type: wind-turbine-1.tbl
637 The location of wind turbines are specified as lat lon, and
638 the turbine type: windturbines.txt
640 ** Using stochastic schemes
641 &stoch
642 skebs = 1, 1, 1,
643 rand_perturb = 1, 1, 1,
646 ** Using kfcupscheme cu_physics (cumulus scheme option). This has been tested with
647 the CAM radiation scheme and is not recommended with other radiation schemes.
648 Regarding cu_rad_feedback, users want the parameterized clouds to affect radiation.
649 Turning it off is only useful as a sensitivity study to determine the importance of
650 that effect. Regarding shallowcu_forced_ra option, setting it to true will override
651 the cloud fraction calculations to a prescribed maximum cloud fraction (a value of 0.36)
652 which can be changed by the user for sensitivity testing purposes.
654 &physics
655 cu_physics = 10, 10, 10,
656 ra_lw_physics = 3, 3, 3,
657 ra_sw_physics = 3, 3, 3,
658 cu_rad_feedback =.true.,.true.,.true.,
659 shallowcu_forced_ra = .false.,.false.,.false.,
660 numBins = 21, 21, 21,
661 thBinSize = 0.1, 0.1, 0.1,
662 rBinSize = 0.0001,0.0001,0.0001,
663 minDeepFreq = 0.333, 0.333, 0.333,
664 minShallowFreq = 0.01, 0.01, 0.01,
667 ** To write out an input file during the model simulation.
669 &time_control
670 write_input = .TRUE.,
671 inputout_interval_m = 60, 60, 60
672 input_outname = "wrfinput_out_d<domain>_<date>"
673 inputout_begin_h = 3, 6, 6
674 inputout_end_h = 9, 6, 6
676 ** To use climatological data in Thompson microphysics option 28: in addition to
677 process additional data in metgrid, use these:
679 &domains
680 wif_input_opt = 1
681 num_wif_levels = 30
683 &physics
684 mp_physics = 28, 28, 28,
685 use_aero_icbc = .true.
688 ** The hybrid vertical coordinate (HVC) requires three pieces:
689 1. The WRF and real codes must be built with HVC activated (./configure -hyb)
690 2. The eta location (etac) at which the eta levels higher up in the atmosphere
691 become isobaric must be defined (suggested default in Registry).
692 3. The run-time option to select the HVC (hybrid_opt, default is OFF):
693 &dynamics
694 hybrid_opt = 2,
695 etac = 0.2,
697 ** New mechanism to specify physics:
699 &physics
700 physics_suite = 'CONUS'
702 which is equivalent to
704 mp_physics = 8,
705 cu_physics = 6,
706 ra_lw_physics = 4,
707 ra_sw_physics = 4,
708 bl_pbl_physics = 2,
709 sf_sfclay_physics = 2,
710 sf_surface_physics = 2,
712 &physics
713 physics_suite = 'TROPICAL'
715 which is equivalent to
717 mp_physics = 6,
718 cu_physics = 16,
719 ra_lw_physics = 4,
720 ra_sw_physics = 4,
721 bl_pbl_physics = 1,
722 sf_sfclay_physics = 91,
723 sf_surface_physics = 2,
726 To overwrite the cu_physics option for a second nest, set
728 &physics
729 physics_suite = 'tropical'
730 cu_physics = -1, -1, 0,
731 ...
733 ** New way to automatically choose levels (default auto_levels_opt=2)
734 max_dz = 1000.
735 auto_levels_opt = 2
736 dzbot = 50.
737 dzstretch_s = 1.3
738 dzstretch_u = 1.1
740 ** Using the scale-adaptive SMS-3DTKE subgrid turbulent mixing scheme:
741 &physics
742 bl_pbl_physics = 0,
743 sf_sfclay_physics = 1(or 5, 91),
744 /
745 &dynamics
746 diff_opt = 2,
747 km_opt = 5,
748 /
751 ** Using lightning option
752 &physics
753 lightning_option = 2,
754 iccg_method = 1,
755 iccg_prescribed_num = 2,
756 iccg_prescribed_den = 1,
757 cldtop_adjustment = 0,
758 lightning_dt = 5, ; default: time_step
759 lightning_start_seconds = 600,
760 flashrate_factor = 1.0,
761 ltng_temp_upper = -45., ; used by lnox_opt=2
762 ltng_temp_lower = -15., ; used by lnox_opt=2
763 /
765 &chem
766 lnox_opt = 2,
767 n_ic = 200,
768 n_cg = 200,
769 lnox_passive = .false.,
770 /
773 ** Using the irrigation parameterizations:
774 &physics
775 sf_surf_irr_scheme = 3,
776 irr_daily_amount = 5.7,
777 irr_start_hour = 2,
778 irr_num_hours = 3,
779 irr_start_julianday = 135,
780 irr_end_julianday = 350,
781 irr_ph = 0,
782 irr_freq = 1,
783 /