1 Description of namelist variables
2 ---------------------------------
4 For WRF-NMM users, please see Chapter 5 of the WRF-NMM User's Guide for
5 information on NMM specific settings (http://www.dtcenter.org/wrf-nmm/users)
8 Note: variables followed by (max_dom) indicate that this variable needs to
9 be defined for the nests when max_dom > 1.
12 run_days = 1, ; run time in days
13 run_hours = 0, ; run time in hours
14 Note: if it is more than 1 day, one may use both run_days and run_hours
15 or just run_hours. e.g. if the total run length is 36 hrs, you may
16 set run_days = 1, and run_hours = 12, or run_days = 0, and run_hours = 36
17 run_minutes = 0, ; run time in minutes
18 run_seconds = 0, ; run time in seconds
19 start_year (max_dom) = 2001, ; four digit year of starting time
20 start_month (max_dom) = 06, ; two digit month of starting time
21 start_day (max_dom) = 11, ; two digit day of starting time
22 start_hour (max_dom) = 12, ; two digit hour of starting time
23 start_minute (max_dom) = 00, ; two digit minute of starting time
24 start_second (max_dom) = 00, ; two digit second of starting time
25 Note: the start time is used to name the first wrfout file.
26 It also controls the start time for nest domains, and the time to restart
27 tstart (max_dom) = 00, ; FOR NMM: starting hour of the forecast
28 end_year (max_dom) = 2001, ; four digit year of ending time
29 end_month (max_dom) = 06, ; two digit month of ending time
30 end_day (max_dom) = 12, ; two digit day of ending time
31 end_hour (max_dom) = 12, ; two digit hour of ending time
32 end_minute (max_dom) = 00, ; two digit minute of ending time
33 end_second (max_dom) = 00, ; two digit second of ending time
34 It also controls when the nest domain integrations end
35 All start and end times are used by real.exe.
37 Note that one may use either run_days/run_hours etc. or
38 end_year/month/day/hour etc. to control the length of
39 model integration. But run_days/run_hours
40 takes precedence over the end times.
41 Program real.exe uses start and end times only.
43 interval_seconds = 10800, ; time interval between incoming real data, which will be the interval
44 between the lateral boundary condition file
45 input_from_file (max_dom) = T, ; whether nested run will have input files for domains other than 1
46 fine_input_stream (max_dom) = 0, ; field selection from nest input for its initialization
47 0: all fields are used; 2: only static and time-varying, masked land
48 surface fields are used.
49 history_interval (max_dom) = 60, ; history output file interval in minutes
50 frames_per_outfile (max_dom) = 1, ; output times per history output file, used to split output files
52 restart = F, ; whether this run is a restart run
53 restart_interval = 1440, ; restart output file interval in minutes
54 reset_simulation_start = F, ; whether to overwrite simulation_start_date with forecast start time
55 io_form_history = 2, ; 2 = netCDF
56 io_form_restart = 2, ; 2 = netCDF
57 io_form_input = 2, ; 2 = netCDF
58 io_form_boundary = 2, ; netCDF format
61 frames_per_emissfile = 12, ; Number of times in each chemistry emission file.
62 io_style_emiss = 1, ; Style to use for the chemistry emission files.
63 ; 0 = Do not read emissions from files.
64 ; 1 = Cycle between two 12 hour files (set frames_per_emissfile=12)
65 ; 2 = Dated files with length set by frames_per_emissfile
66 debug_level = 0, ; 50,100,200,300 values give increasing prints
68 To choose between SI and WPS input to real for EM core:
69 auxinput1_inname = "met_em.d<domain>.<date>" ; Input to real from WPS
70 = "wrf_real_input_em.d<domain>.<date>" ; Input to real from SI
72 To choose between SI and WPS input to real for NMM core:
73 auxinput1_inname = "met_nm.d<domain>.<date>" ; Input to real from WPS
74 = "wrf_real_input_nm.d<domain>.<date>" ; Input to real from SI
78 auxhist2_outname = "rainfall" ; file name for extra output; if not specified,
79 auxhist2_d<domain>_<date> will be used
80 also note that to write variables in output other
81 than the history file requires Registry.EM file change
82 auxhist2_interval (max_dom) = 10, ; interval in minutes
83 io_form_auxhist2 = 2, ; output in netCDF
85 For SST updating (used only with sst_update=1):
87 auxinput4_inname = "wrflowinp_d<domain>"
88 auxinput4_interval = 360 ; minutes generally matches time given by interval_seconds
90 Additional settings when running WRFVAR:
92 write_input = t, ; write input-formatted data as output
93 inputout_interval = 180, ; interval in minutes when writing input-formatted data
94 input_outname = 'wrfinput_d<domain>_<date>' ; you may change the output file name
106 inputout_end_s = 0 ; the above shows that the input-formatted data are output
107 starting from hour 3 to hour 12 in 180 min interval.
110 time_step = 60, ; time step for integration in integer seconds
111 recommend 6*dx (in km) for typical real-data cases
112 time_step_fract_num = 0, ; numerator for fractional time step
113 time_step_fract_den = 1, ; denominator for fractional time step
114 Example, if you want to use 60.3 sec as your time step,
115 set time_step = 60, time_step_fract_num = 3, and
116 time_step_fract_den = 10
117 max_dom = 1, ; number of domains - set it to > 1 if it is a nested run
118 s_we (max_dom) = 1, ; start index in x (west-east) direction (leave as is)
119 e_we (max_dom) = 91, ; end index in x (west-east) direction (staggered dimension)
120 s_sn (max_dom) = 1, ; start index in y (south-north) direction (leave as is)
121 e_sn (max_dom) = 82, ; end index in y (south-north) direction (staggered dimension)
122 s_vert (max_dom) = 1, ; start index in z (vertical) direction (leave as is)
123 e_vert (max_dom) = 28, ; end index in z (vertical) direction (staggered dimension)
124 Note: this refers to full levels including surface and top
125 vertical dimensions need to be the same for all nests
126 Note: most variables are unstaggered (= staggered dim - 1)
127 dx (max_dom) = 10000, ; grid length in x direction, unit in meters
128 dy (max_dom) = 10000, ; grid length in y direction, unit in meters
129 ztop (max_dom) = 19000. ; used in mass model for idealized cases
130 grid_id (max_dom) = 1, ; domain identifier
131 parent_id (max_dom) = 0, ; id of the parent domain
132 i_parent_start (max_dom) = 0, ; starting LLC I-indices from the parent domain
133 j_parent_start (max_dom) = 0, ; starting LLC J-indices from the parent domain
134 parent_grid_ratio (max_dom) = 1, ; parent-to-nest domain grid size ratio: for real-data cases
135 the ratio has to be odd; for idealized cases,
136 the ratio can be even if feedback is set to 0.
137 parent_time_step_ratio (max_dom) = 1, ; parent-to-nest time step ratio; it can be different
138 from the parent_grid_ratio
139 feedback = 1, ; feedback from nest to its parent domain; 0 = no feedback
140 smooth_option = 0 ; smoothing option for parent domain, used only with feedback
141 option on. 0: no smoothing; 1: 1-2-1 smoothing; 2: smoothing-desmoothing
143 Namelist variables specifically for the WPS input for real:
145 num_soil_layers_in = 4 ; number of vertical soil levels or layers input
146 ; from WPS metgrid program
147 num_metgrid_levels = 27 ; number of vertical levels of 3d meteorological fields coming
148 ; from WPS metgrid program
149 interp_type = 2 ; vertical interpolation
150 ; 1 = linear in pressure
151 ; 2 = linear in log(pressure)
152 extrap_type = 2 ; vertical extrapolation of non-temperature fields
153 ; 1 = extrapolate using the two lowest levels
154 ; 2 = use lowest level as constant below ground
155 t_extrap_type = 2 ; vertical extrapolation for potential temperature
157 ; 2 = -6.5 K/km lapse rate for temperature
159 use_levels_below_ground = .true. ; in vertical interpolation, use levels below input surface level
160 ; T = use input isobaric levels below input surface
161 ; F = extrapolate when WRF location is below input surface value
162 use_surface = .true. ; use the input surface level data in the vertical interp and extrap
163 ; T = use the input surface data
164 ; F = do not use the input surface data
165 lagrange_order = 1 ; vertical interpolation order
168 zap_close_levels = 500 ; ignore isobaric level above surface if delta p (Pa) < zap_close_levels
169 lowest_lev_from_sfc = .false. ; place the surface value into the lowest eta location
170 ; T = use surface value as lowest eta (u,v,t,q)
171 ; F = use traditional interpolation
172 force_sfc_in_vinterp = 1 ; use the surface level as the lower boundary when interpolating
173 ; through this many eta levels
174 ; 0 = perform traditional trapping interpolation
175 ; n = first n eta levels directly use surface level
176 sfcp_to_sfcp = .false. ; Optional method to compute model's surface pressure when incoming
177 ; data only has surface pressure and terrain, but not SLP
178 smooth_cg_topo = .false. ; Smooth the outer rows and columns of domain 1's topography w.r.t.
180 p_top_requested = 5000 ; p_top (Pa) to use in the model
181 ptsgm = 42000. ; FOR NMM: defines the pressure interface dividing
182 ; the terrain following portion of the hybrid vertical
183 ; coordinate (p > ptsgm) and the purely
184 ; isobaric portion of the vertical coordinate (p < ptsgm)
186 Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number
187 of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request
188 only the number of levels (with e_vert), and the real program will compute values. The computation assumes
189 a known first several layers, then generates equi-height spaced levels up to the top of the model.
191 eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946,
192 0.922, 0.894, 0.860, 0.817, 0.766,
193 0.707, 0.644, 0.576, 0.507, 0.444,
194 0.380, 0.324, 0.273, 0.228, 0.188,
195 0.152, 0.121, 0.093, 0.069, 0.048,
197 eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
198 0.934, 0.909, 0.880, 0.845, 0.807,
199 0.765, 0.719, 0.672, 0.622, 0.571,
200 0.520, 0.468, 0.420, 0.376, 0.335,
201 0.298, 0.263, 0.231, 0.202, 0.175,
202 0.150, 0.127, 0.106, 0.088, 0.070,
203 0.055, 0.040, 0.026, 0.013, 0.000
205 Namelist variables for controling the specified moving nest:
206 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
207 to the ARCHFLAGS. The maximum number of moves, max_moves, is set to 50
208 but can be modified in source code file frame/module_driver_constants.F.
209 num_moves = 4 ; total number of moves
210 move_id(max_moves) = 2,2,2,2, ; a list of nest domain id's, one per move
211 move_interval(max_moves) = 60,120,150,180, ; time in minutes since the start of this domain
212 move_cd_x(max_moves) = 1,1,0,-1,; the number of parent domain grid cells to move in i direction
213 move_cd_y(max_moves) = 1,0,-1,1,; the number of parent domain grid cells to move in j direction
214 positive is to move in increasing i and j direction, and
215 negative is to move in decreasing i and j direction.
216 0 means no move. The limitation now is to move only 1 grid cell
219 Namelist variables for controling the automatic moving nest:
220 Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
221 and -DVORTEX_CENTER to the ARCHFLAGS. This option uses an mid-level vortex following algorthm to
222 determine the nest move. This option is experimental.
223 vortex_interval(max_dom) = 15 ; how often the new vortex position is computed
224 max_vortex_speed(max_dom) = 40 ; used to compute the search radius for the new vortex position
225 corral_dist(max_dom) = 8 ; how many coarse grid cells the moving nest is allowed to get
226 near the mother domain boundary
227 track_level = 50000 ; pressure value in Pa where the vortex is tracked
229 tile_sz_x = 0, ; number of points in tile x direction
230 tile_sz_y = 0, ; number of points in tile y direction
231 can be determined automatically
232 numtiles = 1, ; number of tiles per patch (alternative to above two items)
233 nproc_x = -1, ; number of processors in x for decomposition
234 nproc_y = -1, ; number of processors in y for decomposition
235 -1: code will do automatic decomposition
236 >1: for both: will be used for decomposition
238 Namelist variables for controlling the adaptive time step option:
239 These options are only valid for the ARW core.
240 use_adaptive_time_step = .false. ; T/F use adaptive time stepping, ARW only
241 step_to_output_time = .true. ; if adaptive time stepping, T/F modify the
242 time steps so that the exact history time is reached
243 target_cfl(max_dom) = 1.2,1.2 ; vertical and horizontal CFL <= to this value implies
244 no reason to reduce the time step, and to increase it
245 max_step_increase_pct(max_dom) = 5,51 ; percentage of previous time step to increase, if the
246 max(vert cfl, horiz cfl) <= target_cfl, then the time
247 will increase by max_step_increase_pct. Use something
248 large for nests (51% suggested)
249 starting_time_step(max_dom) = -1,-1 ; flag = -1 implies use 6 * dx (defined in start_em),
250 starting_time_step = 100 means the starting time step
251 for the coarse grid is 100 s
252 max_time_step(max_dom) = -1,-1 ; flag = -1 implies max time step is 3 * starting_time_step,
253 max_time_step = 100 means that the time step will not
255 min_time_step(max_dom) = -1,-1 ; flag = -1 implies max time step is 0.5 * starting_time_step,
256 min_time_step = 100 means that the time step will not
260 dfi_opt = 0 ; which DFI option to use (3 is recommended)
261 ; 0 = no digital filter initialization
262 ; 1 = digital filter launch (DFL)
263 ; 2 = diabatic DFI (DDFI)
264 ; 3 = twice DFI (TDFI)
265 dfi_nfilter = 7 ; digital filter type to use (7 is recommended)
274 ; 8 = recursive high-order
275 dfi_write_filtered_input = .true. ; whether to write wrfinput file with filtered
276 ; model state before beginning forecast
277 dfi_write_dfi_history = .false. ; whether to write wrfout files during filtering integration
278 dfi_cutoff_seconds = 3600 ; cutoff period, in seconds, for the filter
279 dfi_time_dim = 1000 ; maximum number of time steps for filtering period
280 ; this value can be larger than necessary
281 dfi_bckstop_year = 2004 ; four-digit year of stop time for backward DFI integration
282 dfi_bckstop_month = 03 ; two-digit month of stop time for backward DFI integration
283 dfi_bckstop_day = 14 ; two-digit day of stop time for backward DFI integration
284 dfi_bckstop_hour = 12 ; two-digit hour of stop time for backward DFI integration
285 dfi_bckstop_minute = 00 ; two-digit minute of stop time for backward DFI integration
286 dfi_bckstop_second = 00 ; two-digit second of stop time for backward DFI integration
287 dfi_fwdstop_year = 2004 ; four-digit year of stop time for forward DFI integration
288 dfi_fwdstop_month = 03 ; two-digit month of stop time for forward DFI integration
289 dfi_fwdstop_day = 13 ; two-digit month of stop time for forward DFI integration
290 dfi_fwdstop_hour = 12 ; two-digit month of stop time for forward DFI integration
291 dfi_fwdstop_minute = 00 ; two-digit month of stop time for forward DFI integration
292 dfi_fwdstop_second = 00 ; two-digit month of stop time for forward DFI integration
296 Note: even the physics options can be different in different nest domains,
297 caution must be used as what options are sensible to use
299 chem_opt = 0, ; chemistry option - not yet available
300 mp_physics (max_dom) microphysics option
303 = 2, Lin et al. scheme
304 = 3, WSM 3-class simple ice scheme
305 = 4, WSM 5-class scheme
306 = 5, Ferrier (new Eta) microphysics
307 = 6, WSM 6-class graupel scheme
308 = 7, Goddard GCE scheme (also uses gsfcgce_hail, gsfcgce_2ice)
310 = 10, Morrison (2 moments)
312 For non-zero mp_physics options, to keep Qv .GE. 0, and to set the other moisture
313 fields .LT. a critcal value to zero
315 mp_zero_out = 0, ; no action taken, no adjustment to any moist field
316 = 1, ; except for Qv, all other moist arrays are set to zero
317 ; if they fall below a critical value
318 = 2, ; Qv is .GE. 0, all other moist arrays are set to zero
319 ; if they fall below a critical value
320 mp_zero_out_thresh = 1.e-8 ; critical value for moist array threshold, below which
321 ; moist arrays (except for Qv) are set to zero (kg/kg)
323 gwd_opt = 0 ; for running without gravity wave drag
324 = 1 ; for running the WRF-ARW with its gravity wave drag
325 = 2 ; for running the WRF-NMM with its gravity wave drag
327 gsfcgce_hail = 0 ; for running gsfcgce microphysics with graupel
328 = 1 ; for running gsfcgce microphysics with hail
330 gsfcgce_2ice = 0 ; for running with snow, ice and graupel/hail
331 = 1 ; for running with only ice and snow
332 = 2 ; for running with only ice and graupel
333 (only used in very extreme situation)
335 gsfcgce_hail is ignored if gsfcgce_2ice is set to 1 or 2.
337 no_mp_heating = 0 ; normal
338 = 1 ; turn off latent heating from a microphysics scheme
340 ra_lw_physics (max_dom) longwave radiation option
341 = 0, no longwave radiation
344 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
346 = 31, Earth Held-Suarez forcing
347 = 99, GFDL (Eta) longwave (semi-supported)
348 also must use co2tf = 1 for ARW
350 ra_sw_physics (max_dom) shortwave radiation option
351 = 0, no shortwave radiation
353 = 2, Goddard short wave
355 also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
357 = 99, GFDL (Eta) longwave (semi-supported)
358 also must use co2tf = 1 for ARW
360 radt (max_dom) = 30, ; minutes between radiation physics calls
361 recommend 1 min per km of dx (e.g. 10 for 10 km)
363 nrads (max_dom) = FOR NMM: number of fundamental timesteps between
364 calls to shortwave radiation; the value
365 is set in Registry.NMM but is overridden
366 by namelist value; radt will be computed
369 nradl (max_dom) = FOR NMM: number of fundamental timesteps between
370 calls to longwave radiation; the value
371 is set in Registry.NMM but is overridden
374 co2tf CO2 transmission function flag only for GFDL radiation
375 = 0, read CO2 function data from pre-generated file
376 = 1, generate CO2 functions internally in the forecast
378 ra_call_offset radiation call offset
379 = 0 (no offset), =-1 (old offset)
381 cam_abs_freq_s = 21600 CAM clearsky longwave absorption calculation frequency
382 (recommended minimum value to speed scheme up)
383 levsiz = 59 for CAM radiation input ozone levels
384 paerlev = 29 for CAM radiation input aerosol levels
385 cam_abs_dim1 = 4 for CAM absorption save array
386 cam_abs_dim2 = value of e_vert for CAM 2nd absorption save array
388 sf_sfclay_physics (max_dom) surface-layer option (old bl_sfclay_physics option)
389 = 0, no surface-layer
390 = 1, Monin-Obukhov scheme
391 = 2, Monin-Obukhov (Janjic) scheme
392 = 3, NCEP Global Forecast System scheme (NMM only)
393 = 7, Pleim-Xiu surface layer (ARW only)
395 sf_surface_physics (max_dom) land-surface option (old bl_surface_physics option)
396 = 0, no surface temp prediction
397 = 1, thermal diffusion scheme
398 = 2, Unified Noah land-surface model
399 = 3, RUC land-surface model
400 = 7, Pleim-Xiu LSM (ARW)
402 bl_pbl_physics (max_dom) boundary-layer option
403 = 0, no boundary-layer
405 = 2, Mellor-Yamada-Janjic TKE scheme
406 = 3, NCEP Global Forecast System scheme (NMM only)
407 = 7, ACM2 (Pleim) PBL (ARW)
408 = 99, MRF scheme (to be removed)
410 bldt (max_dom) = 0, ; minutes between boundary-layer physics calls
412 nphs (max_dom) = FOR NMM: number of fundamental timesteps between
413 calls to turbulence and microphysics;
414 the value is set in Registry.NMM but is
415 overridden by namelist value; bldt will
416 be computed from this.
418 cu_physics (max_dom) cumulus option
420 = 1, Kain-Fritsch (new Eta) scheme
421 = 2, Betts-Miller-Janjic scheme
422 = 3, Grell-Devenyi ensemble scheme
423 = 4, Simplified Arakawa-Schubert scheme (NMM only)
424 = 5, Grell 3D ensemble scheme
425 = 99, previous Kain-Fritsch scheme
427 cudt = 0, ; minutes between cumulus physics calls
429 ncnvc (max_dom) = FOR NMM: number of fundamental timesteps between
430 calls to convection; the value is set in Registry.NMM
431 but is overridden by namelist value; cudt will be
434 tprec (max_dom) = FOR NMM: number of hours in precipitation bucket
435 theat (max_dom) = FOR NMM: number of hours in latent heating bucket
436 tclod (max_dom) = FOR NMM: number of hours in cloud fraction average
437 trdsw (max_dom) = FOR NMM: number of hours in short wave buckets
438 trdlw (max_dom) = FOR NMM: number of hours in long wave buckets
439 tsrfc (max_dom) = FOR NMM: number of hours in surface flux buckets
440 pcpflg (max_dom) = FOR NMM: logical switch for precipitation assimilation
442 isfflx = 1, ; heat and moisture fluxes from the surface
443 (only works for sf_sfclay_physics = 1)
444 1 = with fluxes from the surface
445 0 = no flux from the surface
446 with bl_pbl_physics=0 this uses tke_drag_coefficient
447 and tke_heat_flux in vertical diffusion
448 2 = use drag from sf_sfclay_physics and heat flux from
449 tke_heat_flux with bl_pbl_physics=0
450 ifsnow = 0, ; snow-cover effects
451 (only works for sf_surface_physics = 1)
452 1 = with snow-cover effect
453 0 = without snow-cover effect
454 icloud = 1, ; cloud effect to the optical depth in radiation
455 (only works for ra_sw_physics = 1 and ra_lw_physics = 1)
456 1 = with cloud effect
457 0 = without cloud effect
458 swrad_scat = 1. ; scattering tuning parameter (default 1. is 1.e-5 m2/kg)
459 surface_input_source = 1, ; where landuse and soil category data come from:
461 2 = GRIB data from another model (only possible
462 (VEGCAT/SOILCAT are in met_em files from WPS)
463 num_soil_layers = 5, ; number of soil layers in land surface model
464 = 5: thermal diffusion scheme
465 = 4: Noah landsurface model
466 = 6: RUC landsurface model
467 = 2: Pleim-Xu landsurface model
468 sf_urban_physics(max_dom) = 0, ; activate urban canopy model (in Noah LSM only)
470 = 1: Single-layer, Noah UCM (Hiroyuki Kusaka)
471 = 2: Multi-layer, BEP scheme (Alberto Martilli)
472 num_land_cat = 24, ; number of land categories in input data
473 num_soil_cat = 16, ; number of soil categories in input data
475 pxlsm_smois_init(max_dom) = 1 ; PXLSM Soil moisture initialization option
476 0 - From analysis, 1 - From MAVAIL
478 maxiens = 1, ; Grell-Devenyi only
479 maxens = 3, ; G-D only
480 maxens2 = 3, ; G-D only
481 maxens3 = 16 ; G-D only
482 ensdim = 144 ; G-D only
483 These are recommended numbers. If you would like to use
484 any other number, consult the code, know what you are doing.
485 seaice_threshold = 271 ; tsk < seaice_threshold, if water point and 5-layer slab
486 ; scheme, set to land point and permanent ice; if water point
487 ; and Noah scheme, set to land point, permanent ice, set temps
488 ; from 3 m to surface, and set smois and sh2o
489 sst_update = 0 ; time-varying sea-surface temp (0=no, 1=yes). If selected real
490 ; puts SST, XICE, ALBEDO and VEGFRA in wrflowinp_d01 file, and wrf updates
491 ; these from it at same interval as boundary file. To read this, the time-control
492 ; namelist must include auxinput4_interval, auxinput4_end_h, and
493 ; auxinput4_inname = "wrflowinp_d<domain>"
494 usemonalb = .true. ; use monthly albedo map instead of LANDUSE.TBL value
495 ; (must be used for NMM and recommended for sst_update=1)
496 slope_rad = 0 ; slope effects for ra_sw_physics=1 (1=on, 0=off)
497 topo_shading = 0 ; neighboring-point shadow effects for ra_sw_physics=1 (1=on, 0=off)
498 shadlen = 25000. ; max shadow length in meters for topo_shading=1
499 omlcall = 0 ; activate simple ocean mixed layer model (0=no, 1=yes)
500 oml_hml0 = 50 ; oml model can be initialized with a constant depth everywhere (m)
501 oml_gamma = 0.14 ; oml deep water lapse rate (K m-1)
502 isftcflx = 0 ; alternative Ck, Cd formulation for tropical storm application (0=default, 1=new)
503 fractional_seaice = 0 ; fractional_seaice (1=on, 0=off)
506 grid_fdda (max_dom) = 1 ; grid-nudging fdda on (=0 off) for each domain
507 = 2: spectral nudging
508 gfdda_inname = "wrffdda_d<domain>" ; defined name in real
509 gfdda_interval_m (max_dom) = 360 ; time interval (min) between analysis times
510 gfdda_end_h (max_dom) = 6 ; time (h) to stop nudging after start of forecast
511 io_form_gfdda = 2 ; analysis data io format (2 = netCDF)
512 fgdt (max_dom) = 0 ; calculation frequency (minutes) for grid-nudging (0=every step)
513 if_no_pbl_nudging_uv (max_dom) = 0 ; 1= no nudging of u and v in the pbl, 0=nudging in the pbl
514 if_no_pbl_nudging_t (max_dom) = 0 ; 1= no nudging of temp in the pbl, 0=nudging in the pbl
515 if_no_pbl_nudging_q (max_dom) = 0 ; 1= no nudging of qvapor in the pbl, 0=nudging in the pbl
516 if_zfac_uv (max_dom) = 0 ; 0= nudge u and v all layers, 1= limit nudging to levels above k_zfac_uv
517 k_zfac_uv (max_dom) = 10 ; 10=model level below which nudging is switched off for u and v
518 if_zfac_t (max_dom) = 0 ; 0= nudge temp all layers, 1= limit nudging to levels above k_zfac_t
519 k_zfac_t (max_dom) = 10 ; 10=model level below which nudging is switched off for temp
520 if_zfac_q (max_dom) = 0 ; 0= nudge qvapor all layers, 1= limit nudging to levels above k_zfac_q
521 k_zfac_q (max_dom) = 10 ; 10=model level below which nudging is switched off for qvapor
522 guv (max_dom) = 0.0003 ; nudging coefficient for u and v (sec-1)
523 gt (max_dom) = 0.0003 ; nudging coefficient for temp (sec-1)
524 gq (max_dom) = 0.0003 ; nudging coefficient for qvapor (sec-1)
525 if_ramping = 0 ; 0= nudging ends as a step function, 1= ramping nudging down at end of period
526 dtramp_min = 60.0 ; time (min) for ramping function, 60.0=ramping starts at last analysis time,
527 -60.0=ramping ends at last analysis time
529 The following are for spectral nudging:
530 fgdtzero (max_dom) = 0,
531 if_no_pbl_nudging_ph = 0,
532 if_zfac_ph (max_dom) = 0,
533 k_zfac_ph (max_dom) = 10,
534 dk_zfac_uv (max_dom) = 1, ; depth in k between k_zfac_X to dk_zfac_X where nudging increases
535 linearly to full strength
536 dk_zfac_t (max_dom) = 1,
537 dk_zfac_ph (max_dom) = 1,
538 gph (max_dom) = 0.0003,
539 xwavenum (max_dom) = 3, ; top wave number to nudge in x direction
540 ywavenum (max_dom) = 3, ; top wave number to nudge in y direction
542 The following are for observation nudging:
543 obs_nudge_opt (max_dom) = 1 ; obs-nudging fdda on (=0 off) for each domain
544 also need to set auxinput11_interval and auxinput11_end_h
545 in time_control namelist
546 max_obs = 150000 ; max number of observations used on a domain during any
548 fdda_start = 0 ; obs nudging start time in minutes
549 fdda_end = 180 ; obs nudging end time in minutes
550 obs_nudge_wind (max_dom) = 1 ; whether to nudge wind: (=0 off)
551 obs_coef_wind = 6.E-4, ; nudging coefficient for wind, unit: s-1
552 obs_nudge_temp = 1 ; whether to nudge temperature: (=0 off)
553 obs_coef_temp = 6.E-4, ; nudging coefficient for temperature, unit: s-1
554 obs_nudge_mois = 1 ; whether to nudge water vapor mixing ratio: (=0 off)
555 obs_coef_mois = 6.E-4, ; nudging coefficient for water vapor mixing ratio, unit: s-1
556 obs_nudge_pstr = 0 ; whether to nudge surface pressure (not used)
557 obs_coef_pstr = 0. ; nudging coefficient for surface pressure, unit: s-1 (not used)
558 obs_rinxy = 200., ; horizonal radius of influence in km
559 obs_rinsig = 0.1, ; vertical radius of influence in eta
560 obs_twindo (max_dom) = 0.66667 ; half-period time window over which an observation
561 will be used for nudging (hours)
562 obs_npfi = 10, ; freq in coarse grid timesteps for diag prints
563 obs_ionf (max_dom) = 2 ; freq in coarse grid timesteps for obs input and err calc
564 obs_idynin = 0 ; for dynamic initialization using a ramp-down function to gradually
565 turn off the FDDA before the pure forecast (=1 on)
566 obs_dtramp = 40 ; time period in minutes over which the nudging is ramped down
568 obs_nobs_prt (max_dom) = 10, ; Number of current obs to print grid coord. info.
569 obs_ipf_in4dob = .true. ; print obs input diagnostics (=.false. off)
570 obs_ipf_errob = .true. ; print obs error diagnostics (=.false. off)
571 obs_ipf_nudob = .true. ; print obs nudge diagnostics (=.false. off)
572 obs_ipf_init = .true. ; Enable obs init warning messages
576 scm_force = 1, ; switch for single column forcing (=0 off)
577 scm_force_dx = 4000. ; DX for SCM forcing (in meters)
578 num_force_layers = 8 ; number of SCM forcing layers
579 scm_lu_index = 2 ; SCM landuse category (2 is dryland, cropland and pasture)
580 scm_isltyp = 4 ; SCM soil category (4 is silt loam)
581 scm_ivgtyp = 2 ; SCM landuse category
582 scm_vegfra = 0.5 ; SCM vegetation fraction
583 scm_canwat = 0.0 ; SCM canopy water
584 scm_lat = 37.600 ; SCM latitude
585 scm_lon = -96.700 ; SCM longitude
586 scm_th_adv = .true. ; turn on theta advection in SCM
587 scm_wind_adv = .true. ; turn on wind advection in SCM
588 scm_qv_adv = .true. ; turn on moisture advection in SCM
589 scm_vert_adv = .true. ; turn on vertical advection in SCM
592 rk_ord = 3, ; time-integration scheme option:
593 2 = Runge-Kutta 2nd order
594 3 = Runge-Kutta 3rd order
595 diff_opt = 0, ; turbulence and mixing option:
596 0 = no turbulence or explicit
597 spatial numerical filters (km_opt IS IGNORED).
598 1 = evaluates 2nd order
599 diffusion term on coordinate surfaces.
600 uses kvdif for vertical diff unless PBL option
601 is used. may be used with km_opt = 1 and 4.
602 (= 1, recommended for real-data cases)
603 2 = evaluates mixing terms in
604 physical space (stress form) (x,y,z).
605 turbulence parameterization is chosen
606 by specifying km_opt.
607 km_opt = 1, ; eddy coefficient option
608 1 = constant (use khdif kvdif)
609 2 = 1.5 order TKE closure (3D)
610 3 = Smagorinsky first order closure (3D)
611 Note: option 2 and 3 are not recommended for DX > 2 km
612 4 = horizontal Smagorinsky first order closure
613 (recommended for real-data cases)
614 damp_opt = 0, ; upper level damping flag
616 1 = with diffusive damping, maybe used for real-data cases
617 (dampcoef nondimensional ~0.01-0.1)
618 2 = with Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .003;
619 not for real-data cases)
620 3 = with w-Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .05;
622 diff_6th_opt = 0, ; 6th-order numerical diffusion
623 0 = no 6th-order diffusion (default)
624 1 = 6th-order numerical diffusion (not recommended)
625 2 = 6th-order numerical diffusion but prohibit up-gradient diffusion
626 diff_6th_factor = 0.12, ; 6th-order numerical diffusion non-dimensional rate (max value 1.0
627 corresponds to complete removal of 2dx wave in one timestep)
628 dampcoef (max_dom) = 0., ; damping coefficient (see above)
629 zdamp (max_dom) = 5000., ; damping depth (m) from model top
630 w_damping = 0, ; vertical velocity damping flag (for operational use)
633 base_temp = 290., ; real-data, em ONLY, base sea-level temp (K)
634 base_pres = 10^5 ; real-data, em ONLY, base sea-level pres (Pa), DO NOT CHANGE
635 base_lapse = 50., ; real-data, em ONLY, lapse rate (K), DO NOT CHANGE
636 khdif (max_dom) = 0, ; horizontal diffusion constant (m^2/s)
637 kvdif (max_dom) = 0, ; vertical diffusion constant (m^2/s)
638 smdiv (max_dom) = 0.1, ; divergence damping (0.1 is typical)
639 emdiv (max_dom) = 0.01, ; external-mode filter coef for mass coordinate model
640 (0.01 is typical for real-data cases)
641 epssm (max_dom) = .1, ; time off-centering for vertical sound waves
642 non_hydrostatic (max_dom) = .true., ; whether running the model in hydrostatic or non-hydro mode
643 pert_coriolis (max_dom) = .false., ; Coriolis only acts on wind perturbation (idealized)
644 top_lid (max_dom) = .false., ; Zero vertical motion at top of domain
645 mix_full_fields(max_dom) = .true., ; used with diff_opt = 2; value of ".true." is recommended, except for
646 highly idealized numerical tests; damp_opt must not be 1 if ".true."
647 is chosen. .false. means subtract 1-d base-state profile before mixing
648 mix_isotropic(max_dom) = 0 ; 0=anistropic vertical/horizontal diffusion coeffs, 1=isotropic
649 mix_upper_bound(max_dom) = 0.1 ; non-dimensional upper limit for diffusion coeffs
650 tke_drag_coefficient(max_dom) = 0., ; surface drag coefficient (Cd, dimensionless) for diff_opt=2 only
651 tke_heat_flux(max_dom) = 0., ; surface thermal flux (H/(rho*cp), K m/s) for diff_opt=2 only
652 h_mom_adv_order (max_dom) = 5, ; horizontal momentum advection order (5=5th, etc.)
653 v_mom_adv_order (max_dom) = 3, ; vertical momentum advection order
654 h_sca_adv_order (max_dom) = 5, ; horizontal scalar advection order
655 v_sca_adv_order (max_dom) = 3, ; vertical scalar advection order
657 ; advection options for scalar variables: 0=simple, 1=positive definite, 2=monotonic
658 moist_adv_opt (max_dom) = 1 ; for moisture
659 scalar_adv_opt (max_dom) = 1 ; for scalars
660 chem_adv_opt (max_dom) = 1 ; for chem variables
661 tke_adv_opt (max_dom) = 1 ; for tke
663 time_step_sound (max_dom) = 4 / ; number of sound steps per time-step (0=set automatically)
664 (if using a time_step much larger than 6*dx (in km),
665 proportionally increase number of sound steps - also
666 best to use even numbers)
667 do_coriolis (max_dom) = .true., ; whether to do Coriolis calculations (idealized) (inactive)
668 do_curvature (max_dom) = .true., ; whether to do curvature calculations (idealized) (inactive)
669 do_gradp (max_dom) = .true., ; whether to do horizontal pressure gradient calculations (idealized) (inactive)
670 fft_filter_lat = 45. ; the latitude above which the polar filter is turned on
673 spec_bdy_width = 5, ; total number of rows for specified boundary value nudging
674 spec_zone = 1, ; number of points in specified zone (spec b.c. option)
675 relax_zone = 4, ; number of points in relaxation zone (spec b.c. option)
676 specified (max_dom) = .false., ; specified boundary conditions (only can be used for domain 1)
677 the above 4 are used for real-data runs
678 spec_exp = 0. ; exponential multiplier for relaxation zone ramp for specified=.t.
679 (0.=linear ramp default, e.g. 0.33=~3*dx exp decay factor)
681 periodic_x (max_dom) = .false., ; periodic boundary conditions in x direction
682 symmetric_xs (max_dom) = .false., ; symmetric boundary conditions at x start (west)
683 symmetric_xe (max_dom) = .false., ; symmetric boundary conditions at x end (east)
684 open_xs (max_dom) = .false., ; open boundary conditions at x start (west)
685 open_xe (max_dom) = .false., ; open boundary conditions at x end (east)
686 periodic_y (max_dom) = .false., ; periodic boundary conditions in y direction
687 symmetric_ys (max_dom) = .false., ; symmetric boundary conditions at y start (south)
688 symmetric_ye (max_dom) = .false., ; symmetric boundary conditions at y end (north)
689 open_ys (max_dom) = .false., ; open boundary conditions at y start (south)
690 open_ye (max_dom) = .false., ; open boundary conditions at y end (north)
691 nested (max_dom) = .false., ; nested boundary conditions (must be used for nests)
692 polar = .false., ; polar boundary condition
693 (v=0 at polarward-most v-point)
696 &namelist_quilt This namelist record controls asynchronized I/O for MPI applications.
698 nio_tasks_per_group = 0, default value is 0: no quilting; > 0 quilting I/O
699 nio_groups = 1, default 1, don't change
703 background_proc_id = 255, ; Background generating process identifier, typically defined
704 by the originating center to identify the background data that
705 was used in creating the data. This is octet 13 of Section 4
707 forecast_proc_id = 255, ; Analysis or generating forecast process identifier, typically
708 defined by the originating center to identify the forecast process
709 that was used to generate the data. This is octet 14 of Section
710 4 in the grib2 message
711 production_status = 255, ; Production status of processed data in the grib2 message.
712 See Code Table 1.3 of the grib2 manual. This is octet 20 of
713 Section 1 in the grib2 record
714 compression = 40, ; The compression method to encode the output grib2 message.
715 Only 40 for jpeg2000 or 41 for PNG are supported