| 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. 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 | |
|---|
| 78 | To 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 | |
|---|
| 82 | To 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 | |
|---|
| 86 | Other 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 | |
|---|
| 96 | For 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 | |
|---|
| 102 | Options 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 | |
|---|
| 109 | Additional 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 | |
|---|
| 163 | Namelist 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 | |
|---|
| 217 | Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number |
|---|
| 218 | of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request |
|---|
| 219 | only the number of levels (with e_vert), and the real program will compute values. The computation assumes |
|---|
| 220 | a 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 | |
|---|
| 236 | Namelist 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 | |
|---|
| 250 | Namelist 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 | |
|---|
| 270 | Namelist 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 | |
|---|
| 602 | Options 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 | |
|---|
| 617 | Options 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 | |
|---|
| 667 | The 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 | |
|---|
| 680 | The 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 |
|---|
![(please configure the [header_logo] section in trac.ini)](/LMDZ_WRF/chrome/site/your_project_logo.png){kind=logo}
