source: lmdz_wrf/trunk/WRFV3/run/README.namelist @ 915

Last change on this file since 915 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

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