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