| 1 | Description of namelist variables |
|---|
| 2 | --------------------------------- |
|---|
| 3 | |
|---|
| 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) |
|---|
| 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. |
|---|
| 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 |
|---|
| 51 | into smaller pieces |
|---|
| 52 | restart = F, ; whether this run is a restart run |
|---|
| 53 | restart_interval = 1440, ; restart output file interval in minutes |
|---|
| 54 | io_form_history = 2, ; 2 = netCDF |
|---|
| 55 | io_form_restart = 2, ; 2 = netCDF |
|---|
| 56 | io_form_input = 2, ; 2 = netCDF |
|---|
| 57 | io_form_boundary = 2, ; netCDF format |
|---|
| 58 | = 4, ; PHD5 format |
|---|
| 59 | = 5, ; GRIB1 format |
|---|
| 60 | frames_per_emissfile = 12, ; Number of times in each chemistry emission file. |
|---|
| 61 | io_style_emiss = 1, ; Style to use for the chemistry emission files. |
|---|
| 62 | ; 0 = Do not read emissions from files. |
|---|
| 63 | ; 1 = Cycle between two 12 hour files (set frames_per_emissfile=12) |
|---|
| 64 | ; 2 = Dated files with length set by frames_per_emissfile |
|---|
| 65 | debug_level = 0, ; 50,100,200,300 values give increasing prints |
|---|
| 66 | |
|---|
| 67 | To choose between SI and WPS input to real: |
|---|
| 68 | auxinput1_inname = "met_em.d<domain>.<date>" ; Input to real from WPS |
|---|
| 69 | = "wrf_real_input_em.d<domain>.<date>" ; Input to real from SI |
|---|
| 70 | |
|---|
| 71 | Other output options: |
|---|
| 72 | |
|---|
| 73 | auxhist2_outname = "rainfall" ; file name for extra output; if not specified, |
|---|
| 74 | auxhist2_d<domain>_<date> will be used |
|---|
| 75 | also note that to write variables in output other |
|---|
| 76 | than the history file requires Registry.EM file change |
|---|
| 77 | auxhist2_interval (max_dom) = 10, ; interval in minutes |
|---|
| 78 | io_form_auxhist2 = 2, ; output in netCDF |
|---|
| 79 | |
|---|
| 80 | Additional ones when running 3DVAR: |
|---|
| 81 | |
|---|
| 82 | write_input = t, ; write input-formatted data as output |
|---|
| 83 | inputout_interval = 180, ; interval in minutes when writing input-formatted data |
|---|
| 84 | input_outname = 'wrf_3dvar_input_d<domain>_<date>' ; you may change the output file name |
|---|
| 85 | inputout_begin_y = 0 |
|---|
| 86 | inputout_begin_mo = 0 |
|---|
| 87 | inputout_begin_d = 0 |
|---|
| 88 | inputout_begin_h = 3 |
|---|
| 89 | inputout_begin_m = 0 |
|---|
| 90 | inputout_begin_s = 0 |
|---|
| 91 | inputout_end_y = 0 |
|---|
| 92 | inputout_end_mo = 0 |
|---|
| 93 | inputout_end_d = 0 |
|---|
| 94 | inputout_end_h = 12 |
|---|
| 95 | inputout_end_m = 0 |
|---|
| 96 | inputout_end_s = 0 ; the above shows that the input-formatted data are output |
|---|
| 97 | starting from hour 3 to hour 12 in 180 min interval. |
|---|
| 98 | |
|---|
| 99 | &domains |
|---|
| 100 | time_step = 60, ; time step for integration in integer seconds |
|---|
| 101 | recommend 6*dx (in km) for typical real-data cases |
|---|
| 102 | time_step_fract_num = 0, ; numerator for fractional time step |
|---|
| 103 | time_step_fract_den = 1, ; denominator for fractional time step |
|---|
| 104 | Example, if you want to use 60.3 sec as your time step, |
|---|
| 105 | set time_step = 60, time_step_fract_num = 3, and |
|---|
| 106 | time_step_fract_den = 10 |
|---|
| 107 | max_dom = 1, ; number of domains - set it to > 1 if it is a nested run |
|---|
| 108 | s_we (max_dom) = 1, ; start index in x (west-east) direction (leave as is) |
|---|
| 109 | e_we (max_dom) = 91, ; end index in x (west-east) direction (staggered dimension) |
|---|
| 110 | s_sn (max_dom) = 1, ; start index in y (south-north) direction (leave as is) |
|---|
| 111 | e_sn (max_dom) = 82, ; end index in y (south-north) direction (staggered dimension) |
|---|
| 112 | s_vert (max_dom) = 1, ; start index in z (vertical) direction (leave as is) |
|---|
| 113 | e_vert (max_dom) = 28, ; end index in z (vertical) direction (staggered dimension) |
|---|
| 114 | Note: this refers to full levels including surface and top |
|---|
| 115 | vertical dimensions need to be the same for all nests |
|---|
| 116 | Note: most variables are unstaggered (= staggered dim - 1) |
|---|
| 117 | dx (max_dom) = 10000, ; grid length in x direction, unit in meters |
|---|
| 118 | dy (max_dom) = 10000, ; grid length in y direction, unit in meters |
|---|
| 119 | ztop (max_dom) = 19000. ; used in mass model for idealized cases |
|---|
| 120 | grid_id (max_dom) = 1, ; domain identifier |
|---|
| 121 | parent_id (max_dom) = 0, ; id of the parent domain |
|---|
| 122 | i_parent_start (max_dom) = 0, ; starting LLC I-indices from the parent domain |
|---|
| 123 | j_parent_start (max_dom) = 0, ; starting LLC J-indices from the parent domain |
|---|
| 124 | parent_grid_ratio (max_dom) = 1, ; parent-to-nest domain grid size ratio: for real-data cases |
|---|
| 125 | the ratio has to be odd; for idealized cases, |
|---|
| 126 | the ratio can be even if feedback is set to 0. |
|---|
| 127 | parent_time_step_ratio (max_dom) = 1, ; parent-to-nest time step ratio; it can be different |
|---|
| 128 | from the parent_grid_ratio |
|---|
| 129 | feedback = 1, ; feedback from nest to its parent domain; 0 = no feedback |
|---|
| 130 | smooth_option = 0 ; smoothing option for parent domain, used only with feedback |
|---|
| 131 | option on. 0: no smoothing; 1: 1-2-1 smoothing; 2: smoothing-desmoothing |
|---|
| 132 | |
|---|
| 133 | Namelist variables specifically for the WPS input for real: |
|---|
| 134 | |
|---|
| 135 | num_metgrid_levels = 27 ; number of vertical levels of 3d meteorological fields coming |
|---|
| 136 | ; from WPS metgrid program |
|---|
| 137 | interp_type = 1 ; vertical interpolation |
|---|
| 138 | ; 1 = linear in pressure |
|---|
| 139 | ; 2 = linear in log(pressure) |
|---|
| 140 | lagrange_order = 1 ; vertical interpolation order |
|---|
| 141 | ; 1 = linear |
|---|
| 142 | ; 2 = quadratic |
|---|
| 143 | zap_close_levels = 500 ; ignore isobaric level above surface if delta p (Pa) < zap_close_levels |
|---|
| 144 | lowest_lev_from_sfc = .false. ; place the surface value into the lowest eta location |
|---|
| 145 | ; T = use surface value as lowest eta (u,v,t,q) |
|---|
| 146 | ; F = use traditional interpolation |
|---|
| 147 | force_sfc_in_vinterp = 1 ; use the surface level as the lower boundary when interpolating |
|---|
| 148 | ; through this many eta levels |
|---|
| 149 | ; 0 = perform traditional trapping interpolation |
|---|
| 150 | ; n = first n eta levels directly use surface level |
|---|
| 151 | p_top_requested = 5000 ; p_top (Pa) to use in the model |
|---|
| 152 | |
|---|
| 153 | Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number |
|---|
| 154 | of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request |
|---|
| 155 | only the number of levels (with e_vert), and the real program will compute values. The computation assumes |
|---|
| 156 | a known first several layers, then generates equi-height spaced levels up to the top of the model. |
|---|
| 157 | |
|---|
| 158 | eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946, |
|---|
| 159 | 0.922, 0.894, 0.860, 0.817, 0.766, |
|---|
| 160 | 0.707, 0.644, 0.576, 0.507, 0.444, |
|---|
| 161 | 0.380, 0.324, 0.273, 0.228, 0.188, |
|---|
| 162 | 0.152, 0.121, 0.093, 0.069, 0.048, |
|---|
| 163 | 0.029, 0.014, 0.000, |
|---|
| 164 | eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954, |
|---|
| 165 | 0.934, 0.909, 0.880, 0.845, 0.807, |
|---|
| 166 | 0.765, 0.719, 0.672, 0.622, 0.571, |
|---|
| 167 | 0.520, 0.468, 0.420, 0.376, 0.335, |
|---|
| 168 | 0.298, 0.263, 0.231, 0.202, 0.175, |
|---|
| 169 | 0.150, 0.127, 0.106, 0.088, 0.070, |
|---|
| 170 | 0.055, 0.040, 0.026, 0.013, 0.000 |
|---|
| 171 | |
|---|
| 172 | Namelist variables for controling the specified moving nest: |
|---|
| 173 | Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS |
|---|
| 174 | to the ARCHFLAGS. The maximum number of moves, max_moves, is set to 50 |
|---|
| 175 | but can be modified in source code file frame/module_driver_constants.F. |
|---|
| 176 | num_moves = 4 ; total number of moves |
|---|
| 177 | move_id = 2,2,2,2, ; a list of nest domain id's, one per move |
|---|
| 178 | move_interval = 60,120,150,180, ; time in minutes since the start of this domain |
|---|
| 179 | move_cd_x = 1,1,0,-1,; the number of parent domain grid cells to move in i direction |
|---|
| 180 | move_cd_y = 1,0,-1,1,; the number of parent domain grid cells to move in j direction |
|---|
| 181 | positive is to move in increasing i and j direction, and |
|---|
| 182 | negative is to move in decreasing i and j direction. |
|---|
| 183 | 0 means no move. The limitation now is to move only 1 grid cell |
|---|
| 184 | at each move. |
|---|
| 185 | |
|---|
| 186 | Namelist variables for controling the automatic moving nest: |
|---|
| 187 | Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS |
|---|
| 188 | and -DVORTEX_CENTER to the ARCHFLAGS. This option uses an mid-level vortex following algorthm to |
|---|
| 189 | determine the nest move. This option is experimental. |
|---|
| 190 | vortex_interval = 15 ; how often the new vortex position is computed |
|---|
| 191 | max_vortex_speed = 40 ; used to compute the search radius for the new vortex position |
|---|
| 192 | corral_dist = 8 ; how many coarse grid cells the moving nest is allowed to get |
|---|
| 193 | near the mother domain boundary |
|---|
| 194 | |
|---|
| 195 | &physics |
|---|
| 196 | |
|---|
| 197 | Note: even the physics options can be different in different nest domains, |
|---|
| 198 | caution must be used as what options are sensible to use |
|---|
| 199 | |
|---|
| 200 | chem_opt = 0, ; chemistry option - not yet available |
|---|
| 201 | mp_physics (max_dom) microphysics option |
|---|
| 202 | = 0, no microphysics |
|---|
| 203 | = 1, Kessler scheme |
|---|
| 204 | = 2, Lin et al. scheme |
|---|
| 205 | = 3, WSM 3-class simple ice scheme |
|---|
| 206 | = 4, WSM 5-class scheme |
|---|
| 207 | = 5, Ferrier (new Eta) microphysics |
|---|
| 208 | = 6, WSM 6-class graupel scheme |
|---|
| 209 | = 8, Thompson et al. scheme |
|---|
| 210 | = 98, NCEP 3-class simple ice scheme (to be removed) |
|---|
| 211 | = 99, NCEP 5-class scheme (to be removed) |
|---|
| 212 | |
|---|
| 213 | For non-zero mp_physics options, to keep Qv .GE. 0, and to set the other moisture |
|---|
| 214 | fields .LT. a critcal value to zero |
|---|
| 215 | |
|---|
| 216 | mp_zero_out = 0, ; no action taken, no adjustment to any moist field |
|---|
| 217 | = 1, ; except for Qv, all other moist arrays are set to zero |
|---|
| 218 | ; if they fall below a critical value |
|---|
| 219 | = 2, ; Qv is .GE. 0, all other moist arrays are set to zero |
|---|
| 220 | ; if they fall below a critical value |
|---|
| 221 | mp_zero_out_thresh = 1.e-8 ; critical value for moist array threshold, below which |
|---|
| 222 | ; moist arrays (except for Qv) are set to zero (kg/kg) |
|---|
| 223 | |
|---|
| 224 | ra_lw_physics (max_dom) longwave radiation option |
|---|
| 225 | = 0, no longwave radiation |
|---|
| 226 | = 1, rrtm scheme |
|---|
| 227 | = 3, cam scheme |
|---|
| 228 | also must set levsiz, paerlev, cam_abs_dim1/2 (see below) |
|---|
| 229 | = 99, GFDL (Eta) longwave (semi-supported) |
|---|
| 230 | also must use co2tf = 1 for ARW |
|---|
| 231 | |
|---|
| 232 | ra_sw_physics (max_dom) shortwave radiation option |
|---|
| 233 | = 0, no shortwave radiation |
|---|
| 234 | = 1, Dudhia scheme |
|---|
| 235 | = 2, Goddard short wave |
|---|
| 236 | = 3, cam scheme |
|---|
| 237 | also must set levsiz, paerlev, cam_abs_dim1/2 (see below) |
|---|
| 238 | = 99, GFDL (Eta) longwave (semi-supported) |
|---|
| 239 | also must use co2tf = 1 for ARW |
|---|
| 240 | |
|---|
| 241 | radt (max_dom) = 30, ; minutes between radiation physics calls |
|---|
| 242 | recommend 1 min per km of dx (e.g. 10 for 10 km) |
|---|
| 243 | |
|---|
| 244 | nrads (max_dom) = FOR NMM: number of fundamental timesteps between |
|---|
| 245 | calls to shortwave radiation; the value |
|---|
| 246 | is set in Registry.NMM but is overridden |
|---|
| 247 | by namelist value; radt will be computed |
|---|
| 248 | from this. |
|---|
| 249 | |
|---|
| 250 | nradl (max_dom) = FOR NMM: number of fundamental timesteps between |
|---|
| 251 | calls to longwave radiation; the value |
|---|
| 252 | is set in Registry.NMM but is overridden |
|---|
| 253 | by namelist value. |
|---|
| 254 | |
|---|
| 255 | co2tf CO2 transmission function flag only for GFDL radiation |
|---|
| 256 | = 0, read CO2 function data from pre-generated file |
|---|
| 257 | = 1, generate CO2 functions internally in the forecast |
|---|
| 258 | |
|---|
| 259 | ra_call_offset radiation call offset |
|---|
| 260 | = 0 (no offset), =-1 (old offset) |
|---|
| 261 | |
|---|
| 262 | cam_abs_freq_s = 21600 CAM clearsky longwave absorption calculation frequency |
|---|
| 263 | (recommended minimum value to speed scheme up) |
|---|
| 264 | levsiz = 59 for CAM radiation input ozone levels |
|---|
| 265 | paerlev = 29 for CAM radiation input aerosol levels |
|---|
| 266 | cam_abs_dim1 = 4 for CAM absorption save array |
|---|
| 267 | cam_abs_dim2 = e_vert for CAM 2nd absorption save array |
|---|
| 268 | |
|---|
| 269 | sf_sfclay_physics (max_dom) surface-layer option (old bl_sfclay_physics option) |
|---|
| 270 | = 0, no surface-layer |
|---|
| 271 | = 1, Monin-Obukhov scheme |
|---|
| 272 | = 2, Monin-Obukhov (Janjic) scheme |
|---|
| 273 | = 3, NCEP Global Forecast System scheme |
|---|
| 274 | |
|---|
| 275 | sf_surface_physics (max_dom) land-surface option (old bl_surface_physics option) |
|---|
| 276 | = 0, no surface temp prediction |
|---|
| 277 | = 1, thermal diffusion scheme |
|---|
| 278 | = 2, Noah land-surface model |
|---|
| 279 | = 3, RUC land-surface model |
|---|
| 280 | |
|---|
| 281 | bl_pbl_physics (max_dom) boundary-layer option |
|---|
| 282 | = 0, no boundary-layer |
|---|
| 283 | = 1, YSU scheme |
|---|
| 284 | = 2, Mellor-Yamada-Janjic TKE scheme |
|---|
| 285 | = 3, NCEP Global Forecast System scheme |
|---|
| 286 | = 99, MRF scheme (to be removed) |
|---|
| 287 | |
|---|
| 288 | bldt (max_dom) = 0, ; minutes between boundary-layer physics calls |
|---|
| 289 | |
|---|
| 290 | nphs (max_dom) = FOR NMM: number of fundamental timesteps between |
|---|
| 291 | calls to turbulence and microphysics; |
|---|
| 292 | the value is set in Registry.NMM but is |
|---|
| 293 | overridden by namelist value; bldt will |
|---|
| 294 | be computed from this. |
|---|
| 295 | |
|---|
| 296 | cu_physics (max_dom) cumulus option |
|---|
| 297 | = 0, no cumulus |
|---|
| 298 | = 1, Kain-Fritsch (new Eta) scheme |
|---|
| 299 | = 2, Betts-Miller-Janjic scheme |
|---|
| 300 | = 3, Grell-Devenyi ensemble scheme |
|---|
| 301 | = 4, Simplified Arakawa-Schubert scheme |
|---|
| 302 | = 99, previous Kain-Fritsch scheme |
|---|
| 303 | |
|---|
| 304 | cudt = 0, ; minutes between cumulus physics calls |
|---|
| 305 | |
|---|
| 306 | ncnvc (max_dom) = FOR NMM: number of fundamental timesteps between |
|---|
| 307 | calls to convection; the value is set in Registry.NMM |
|---|
| 308 | but is overridden by namelist value; cudt will be |
|---|
| 309 | computed from this. |
|---|
| 310 | |
|---|
| 311 | tprec (max_dom) = FOR NMM: number of hours in precipitation bucket |
|---|
| 312 | theat (max_dom) = FOR NMM: number of hours in latent heating bucket |
|---|
| 313 | tclod (max_dom) = FOR NMM: number of hours in cloud fraction average |
|---|
| 314 | trdsw (max_dom) = FOR NMM: number of hours in short wave buckets |
|---|
| 315 | trdlw (max_dom) = FOR NMM: number of hours in long wave buckets |
|---|
| 316 | tsrfc (max_dom) = FOR NMM: number of hours in surface flux buckets |
|---|
| 317 | pcpflg (max_dom) = FOR NMM: logical switch for precipitation assimilation |
|---|
| 318 | |
|---|
| 319 | isfflx = 1, ; heat and moisture fluxes from the surface |
|---|
| 320 | (only works for sf_sfclay_physics = 1) |
|---|
| 321 | 1 = with fluxes from the surface |
|---|
| 322 | 0 = no flux from the surface |
|---|
| 323 | ifsnow = 0, ; snow-cover effects |
|---|
| 324 | (only works for sf_surface_physics = 1) |
|---|
| 325 | 1 = with snow-cover effect |
|---|
| 326 | 0 = without snow-cover effect |
|---|
| 327 | icloud = 1, ; cloud effect to the optical depth in radiation |
|---|
| 328 | (only works for ra_sw_physics = 1 and ra_lw_physics = 1) |
|---|
| 329 | 1 = with cloud effect |
|---|
| 330 | 0 = without cloud effect |
|---|
| 331 | swrad_scat = 1. ; scattering tuning parameter (default 1. is 1.e-5 m2/kg) |
|---|
| 332 | surface_input_source = 1, ; where landuse and soil category data come from: |
|---|
| 333 | 1 = SI/gridgen |
|---|
| 334 | 2 = GRIB data from another model (only possible |
|---|
| 335 | (VEGCAT/SOILCAT are in wrf_real_input_em files from SI) |
|---|
| 336 | num_soil_layers = 5, ; number of soil layers in land surface model |
|---|
| 337 | = 5: thermal diffusion scheme |
|---|
| 338 | = 4: Noah landsurface model |
|---|
| 339 | = 6: RUC landsurface model |
|---|
| 340 | ucmcall = 0, ; activate urban canopy model (in Noah LSM only) (0=no, 1=yes) |
|---|
| 341 | |
|---|
| 342 | maxiens = 1, ; Grell-Devenyi only |
|---|
| 343 | maxens = 3, ; G-D only |
|---|
| 344 | maxens2 = 3, ; G-D only |
|---|
| 345 | maxens3 = 16 ; G-D only |
|---|
| 346 | ensdim = 144 ; G-D only |
|---|
| 347 | These are recommended numbers. If you would like to use |
|---|
| 348 | any other number, consult the code, know what you are doing. |
|---|
| 349 | seaice_threshold = 271 ; tsk < seaice_threshold, if water point and 5-layer slab |
|---|
| 350 | ; scheme, set to land point and permanent ice; if water point |
|---|
| 351 | ; and Noah scheme, set to land point, permanent ice, set temps |
|---|
| 352 | ; from 3 m to surface, and set smois and sh2o |
|---|
| 353 | sst_update = 0 ; time-varying sea-surface temp (0=no, 1=yes). If selected real |
|---|
| 354 | ; puts SST and VEGFRA in wrflowinp_d01 file, and wrf updates these from it |
|---|
| 355 | ; at same interval as boundary file. To read this, the time-control |
|---|
| 356 | ; namelist must include auxinput5_interval, auxinput5_end_h, and |
|---|
| 357 | ; auxinput5_inname = "wrflowinp_d<domain>" |
|---|
| 358 | |
|---|
| 359 | &fdda |
|---|
| 360 | grid_fdda (max_dom) = 1 ; grid-nudging fdda on (=0 off) for each domain |
|---|
| 361 | gfdda_inname = "wrffdda_d<domain>" ; defined name in real |
|---|
| 362 | gfdda_interval_m (max_dom) = 360 ; time interval (min) between analysis times |
|---|
| 363 | gfdda_end_h (max_dom) = 6 ; time (h) to stop nudging after start of forecast |
|---|
| 364 | io_form_gfdda = 2 ; analysis data io format (2 = netCDF) |
|---|
| 365 | fgdt (max_dom) = 0 ; calculation frequency (minutes) for grid-nudging (0=every step) |
|---|
| 366 | if_no_pbl_nudging_uv (max_dom) = 0 ; 0= no nudging of u and v in the pbl, 1=nudging in the pbl |
|---|
| 367 | if_no_pbl_nudging_t (max_dom) = 0 ; 0= no nudging of temp in the pbl, 1=nudging in the pbl |
|---|
| 368 | if_no_pbl_nudging_q (max_dom) = 0 ; 0= no nudging of qvapor in the pbl, 1=nudging in the pbl |
|---|
| 369 | if_zfac_uv (max_dom) = 0 ; 0= nudge u and v all layers, 1= limit nudging to levels above k_zfac_uv |
|---|
| 370 | k_zfac_uv (max_dom) = 10 ; 10=model level below which nudging is switched off for u and v |
|---|
| 371 | if_zfac_t (max_dom) = 0 ; 0= nudge temp all layers, 1= limit nudging to levels above k_zfac_t |
|---|
| 372 | k_zfac_t (max_dom) = 10 ; 10=model level below which nudging is switched off for temp |
|---|
| 373 | if_zfac_q (max_dom) = 0 ; 0= nudge qvapor all layers, 1= limit nudging to levels above k_zfac_q |
|---|
| 374 | k_zfac_q (max_dom) = 10 ; 10=model level below which nudging is switched off for qvapor |
|---|
| 375 | guv (max_dom) = 0.0003 ; nudging coefficient for u and v (sec-1) |
|---|
| 376 | gt (max_dom) = 0.0003 ; nudging coefficient for temp (sec-1) |
|---|
| 377 | gq (max_dom) = 0.0003 ; nudging coefficient for qvapor (sec-1) |
|---|
| 378 | if_ramping = 0 ; 0= nudging ends as a step function, 1= ramping nudging down at end of period |
|---|
| 379 | dtramp_min = 60.0 ; time (min) for ramping function, 60.0=ramping starts at last analysis time, |
|---|
| 380 | -60.0=ramping ends at last analysis time |
|---|
| 381 | |
|---|
| 382 | The following are for observation nudging: |
|---|
| 383 | obs_nudge_opt (max_dom) = 1 ; obs-nudging fdda on (=0 off) for each domain |
|---|
| 384 | also need to set auxinput11_interval and auxinput11_end_h |
|---|
| 385 | in time_control namelist |
|---|
| 386 | max_obs = 150000 ; max number of observations used on a domain during any |
|---|
| 387 | given time window |
|---|
| 388 | fdda_start = 0 ; obs nudging start time in minutes |
|---|
| 389 | fdda_end = 180 ; obs nudging end time in minutes |
|---|
| 390 | obs_nudge_wind (max_dom) = 1 ; whether to nudge wind: (=0 off) |
|---|
| 391 | obs_coef_wind = 6.E-4, ; nudging coefficient for wind, unit: s-1 |
|---|
| 392 | obs_nudge_temp = 1 ; whether to nudge temperature: (=0 off) |
|---|
| 393 | obs_coef_temp = 6.E-4, ; nudging coefficient for temperature, unit: s-1 |
|---|
| 394 | obs_nudge_mois = 1 ; whether to nudge water vapor mixing ratio: (=0 off) |
|---|
| 395 | obs_coef_mois = 6.E-4, ; nudging coefficient for water vapor mixing ratio, unit: s-1 |
|---|
| 396 | obs_nudge_pstr = 0 ; whether to nudge surface pressure (not used) |
|---|
| 397 | obs_coef_pstr = 0. ; nudging coefficient for surface pressure, unit: s-1 (not used) |
|---|
| 398 | obs_rinxy = 200., ; horizonal radius of influence in km |
|---|
| 399 | obs_rinsig = 0.1, ; vertical radius of influence in eta |
|---|
| 400 | obs_twindo = 40, ; half-period time window over which an observation |
|---|
| 401 | will be used for nudging |
|---|
| 402 | obs_npfi = 10, ; freq in coarse grid timesteps for diag prints |
|---|
| 403 | obs_ionf = 2 ; freq in coarse grid timesteps for obs input and err calc |
|---|
| 404 | obs_idynin = 0 ; for dynamic initialization using a ramp-down function to gradually |
|---|
| 405 | turn off the FDDA before the pure forecast (=1 on) |
|---|
| 406 | obs_dtramp = 40 ; time period in minutes over which the nudging is ramped down |
|---|
| 407 | from one to zero. |
|---|
| 408 | obs_ipf_in4dob = .true. ; print obs input diagnostics (=.false. off) |
|---|
| 409 | obs_ipf_errob = .true. ; print obs error diagnostics (=.false. off) |
|---|
| 410 | obs_ipf_nudob = .true. ; print obs nudge diagnostics (=.false. off) |
|---|
| 411 | / |
|---|
| 412 | |
|---|
| 413 | |
|---|
| 414 | &dynamics |
|---|
| 415 | dyn_opt = 2, ; dynamical core option: advanced research WRF core (Eulerian mass) |
|---|
| 416 | rk_ord = 3, ; time-integration scheme option: |
|---|
| 417 | 2 = Runge-Kutta 2nd order |
|---|
| 418 | 3 = Runge-Kutta 3rd order |
|---|
| 419 | diff_opt = 0, ; turbulence and mixing option: |
|---|
| 420 | 0 = no turbulence or explicit |
|---|
| 421 | spatial numerical filters (km_opt IS IGNORED). |
|---|
| 422 | 1 = evaluates 2nd order |
|---|
| 423 | diffusion term on coordinate surfaces. |
|---|
| 424 | uses kvdif for vertical diff unless PBL option |
|---|
| 425 | is used. may be used with km_opt = 1 and 4. |
|---|
| 426 | (= 1, recommended for real-data case when grid distance < 10 km) |
|---|
| 427 | 2 = evaluates mixing terms in |
|---|
| 428 | physical space (stress form) (x,y,z). |
|---|
| 429 | turbulence parameterization is chosen |
|---|
| 430 | by specifying km_opt. |
|---|
| 431 | km_opt = 1, ; eddy coefficient option |
|---|
| 432 | 1 = constant (use khdif kvdif) |
|---|
| 433 | 2 = 1.5 order TKE closure (3D) |
|---|
| 434 | 3 = Smagorinsky first order closure (3D) |
|---|
| 435 | Note: option 2 and 3 are not recommended for DX > 2 km |
|---|
| 436 | 4 = horizontal Smagorinsky first order closure |
|---|
| 437 | (recommended for real-data case when grid distance < 10 km) |
|---|
| 438 | damp_opt = 0, ; upper level damping flag |
|---|
| 439 | 0 = without damping |
|---|
| 440 | 1 = with diffusive damping, maybe used for real-data cases |
|---|
| 441 | (dampcoef nondimensional ~0.01-0.1) |
|---|
| 442 | 2 = with Rayleigh damping (dampcoef inverse time scale [1/s] e.g. .003; |
|---|
| 443 | not for real-data cases) |
|---|
| 444 | diff_6th_opt = 0, ; 6th-order numerical diffusion |
|---|
| 445 | 0 = no 6th-order diffusion (default) |
|---|
| 446 | 1 = 6th-order numerical diffusion |
|---|
| 447 | 2 = 6th-order numerical diffusion but prohibit up-gradient diffusion |
|---|
| 448 | diff_6th_factor = 0.12, ; 6th-order numerical diffusion non-dimensional rate (max value 1.0 |
|---|
| 449 | corresponds to complete removal of 2dx wave in one timestep) |
|---|
| 450 | dampcoef (max_dom) = 0., ; damping coefficient (see above) |
|---|
| 451 | zdamp (max_dom) = 5000., ; damping depth (m) from model top |
|---|
| 452 | w_damping = 0, ; vertical velocity damping flag (for operational use) |
|---|
| 453 | 0 = without damping |
|---|
| 454 | 1 = with damping |
|---|
| 455 | base_temp = 290., ; real-data, em ONLY, base sea-level temp (K) |
|---|
| 456 | base_pres = 10^5 ; real-data, em ONLY, base sea-level pres (Pa), DO NOT CHANGE |
|---|
| 457 | base_lapse = 50., ; real-data, em ONLY, lapse rate (K), DO NOT CHANGE |
|---|
| 458 | khdif (max_dom) = 0, ; horizontal diffusion constant (m^2/s) |
|---|
| 459 | kvdif (max_dom) = 0, ; vertical diffusion constant (m^2/s) |
|---|
| 460 | smdiv (max_dom) = 0.1, ; divergence damping (0.1 is typical) |
|---|
| 461 | emdiv (max_dom) = 0.01, ; external-mode filter coef for mass coordinate model |
|---|
| 462 | (0.01 is typical for real-data cases) |
|---|
| 463 | epssm (max_dom) = .1, ; time off-centering for vertical sound waves |
|---|
| 464 | non_hydrostatic (max_dom) = .true., ; whether running the model in hydrostatic or non-hydro mode |
|---|
| 465 | pert_coriolis (max_dom) = .false., ; Coriolis only acts on wind perturbation (idealized) |
|---|
| 466 | mix_full_fields(max_dom) = .true., ; used with diff_opt = 2; value of ".true." is recommended, except for |
|---|
| 467 | highly idealized numerical tests; damp_opt must not be 1 if ".true." |
|---|
| 468 | is chosen. .false. means subtract 1-d base-state profile before mixing |
|---|
| 469 | tke_drag_coefficient(max_dom) = 0., ; surface drag coefficient (Cd, dimensionless) for diff_opt=2 only |
|---|
| 470 | tke_heat_flux(max_dom) = 0., ; surface thermal flux (H/(rho*cp), K m/s) for diff_opt=2 only |
|---|
| 471 | h_mom_adv_order (max_dom) = 5, ; horizontal momentum advection order (5=5th, etc.) |
|---|
| 472 | v_mom_adv_order (max_dom) = 3, ; vertical momentum advection order |
|---|
| 473 | h_sca_adv_order (max_dom) = 5, ; horizontal scalar advection order |
|---|
| 474 | v_sca_adv_order (max_dom) = 3, ; vertical scalar advection order |
|---|
| 475 | pd_moist = F ; positive definite advection of moisture |
|---|
| 476 | pd_scalar = F ; positive definite advection of scalars |
|---|
| 477 | pd_chem = F ; positive definite advection of chem variables |
|---|
| 478 | pd_tke = F ; positive definite advection of tke |
|---|
| 479 | time_step_sound (max_dom) = 4 / ; number of sound steps per time-step (0=set automatically) |
|---|
| 480 | (if using a time_step much larger than 6*dx (in km), |
|---|
| 481 | proportionally increase number of sound steps - also |
|---|
| 482 | best to use even numbers) |
|---|
| 483 | |
|---|
| 484 | |
|---|
| 485 | &bdy_control |
|---|
| 486 | spec_bdy_width = 5, ; total number of rows for specified boundary value nudging |
|---|
| 487 | spec_zone = 1, ; number of points in specified zone (spec b.c. option) |
|---|
| 488 | relax_zone = 4, ; number of points in relaxation zone (spec b.c. option) |
|---|
| 489 | specified (max_dom) = .false., ; specified boundary conditions (only for domain 1) |
|---|
| 490 | the above 4 are used for real-data runs |
|---|
| 491 | |
|---|
| 492 | periodic_x (max_dom) = .false., ; periodic boundary conditions in x direction |
|---|
| 493 | symmetric_xs (max_dom) = .false., ; symmetric boundary conditions at x start (west) |
|---|
| 494 | symmetric_xe (max_dom) = .false., ; symmetric boundary conditions at x end (east) |
|---|
| 495 | open_xs (max_dom) = .false., ; open boundary conditions at x start (west) |
|---|
| 496 | open_xe (max_dom) = .false., ; open boundary conditions at x end (east) |
|---|
| 497 | periodic_y (max_dom) = .false., ; periodic boundary conditions in y direction |
|---|
| 498 | symmetric_ys (max_dom) = .false., ; symmetric boundary conditions at y start (south) |
|---|
| 499 | symmetric_ye (max_dom) = .false., ; symmetric boundary conditions at y end (north) |
|---|
| 500 | open_ys (max_dom) = .false., ; open boundary conditions at y start (south) |
|---|
| 501 | open_ye (max_dom) = .false., ; open boundary conditions at y end (north) |
|---|
| 502 | nested (max_dom) = .false., ; nested boundary conditions (inactive) |
|---|
| 503 | |
|---|
| 504 | |
|---|
| 505 | &namelist_quilt This namelist record controls asynchronized I/O for MPI applications. |
|---|
| 506 | |
|---|
| 507 | nio_tasks_per_group = 0, default value is 0: no quilting; > 0 quilting I/O |
|---|
| 508 | nio_groups = 1, default 1, don't change |
|---|
| 509 | |
|---|
| 510 | |
|---|
| 511 | miscelleneous in &domains: |
|---|
| 512 | tile_sz_x = 0, ; number of points in tile x direction |
|---|
| 513 | tile_sz_y = 0, ; number of points in tile y direction |
|---|
| 514 | can be determined automatically |
|---|
| 515 | numtiles = 1, ; number of tiles per patch (alternative to above two items) |
|---|
| 516 | nproc_x = -1, ; number of processors in x for decomposition |
|---|
| 517 | nproc_y = -1, ; number of processors in y for decomposition |
|---|
| 518 | -1: code will do automatic decomposition |
|---|
| 519 | >1: for both: will be used for decomposition |
|---|
