Note, this is not a namelist.input file. Find what interests you, and cut and paste them to your own namelist.input file. For more information on these namelist parameters, please see run/README.namelist or Chapter 5 of the User's Guide. ** More options for real in namelist record &domains: p_top_requested = 5000 interp_type = 2 extrap_type = 2 t_extrap_type = 2 lowest_lev_from_sfc = .false. use_levels_below_ground = .true. use_surface = .true. lagrange_order = 2 force_sfc_in_vinterp = 1 zap_close_levels = 500 sfcp_to_sfcp = .false. adjust_heights = .false. smooth_cg_topo = .false. eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946, 0.922, 0.894, 0.860, 0.817, 0.766, 0.707, 0.644, 0.576, 0.507, 0.444, 0.380, 0.324, 0.273, 0.228, 0.188, 0.152, 0.121, 0.093, 0.069, 0.048, 0.029, 0.014, 0.000, ** Using sst_update option (add these to namelist records &time_control and &physics respectively): &time_control auxinput4_inname = "wrflowinp_d" auxinput4_interval = 360, 360, 360, io_form_auxinput4 = 2 &physics sst_update = 1, ** Using Noah-MP option (Use sf_surface_physics option = 4, and add &noah_mp namelist record) &physics sf_surface_physics = 4 &noah_mp dveg = 4, opt_crs = 1, opt_btr = 1, opt_run = 1, opt_sfc = 1, opt_frz = 1, opt_inf = 1, opt_rad = 3, opt_alb = 2, opt_snf = 1, opt_tbot = 2, opt_stc = 1, / ** Using lake model &physics sf_lake_physics = 1, 1, 1, lakedepth_default = 50., 50., 50., lake_min_elev = 5., 5., 5., ** Using stochastic backscatter scheme (new namelist record since v3.6) &stoch stoch_force_opt = 1, 1, 1, stoch_vertstruc_opt = 0, 0, 0, tot_backscat_psi = 1.E-05, 1.E-05, 1.E-05, tot_backscat_t = 1.E-06, 1.E-06, 1.E-06, nens = 1, ztau_psi = 10800.0, ztau_t = 10800.0, rexponent_psi =-1.83, rexponent_t =-1.83, zsigma2_eps = 0.0833, zsigma2_eta = 0.0833, kminforc = 1, kminforc = 1, kminforct = 1, lminforct = 1, kmaxforc = 1000000, lmaxforc = 1000000, kmaxforct = 1000000, lmaxforct = 1000000, perturb_bdy = 0, ** Using DFI options (note this is a separate namelist record): &dfi_control dfi_opt = 3, dfi_nfilter = 7, dfi_cutoff_seconds = 3600, dfi_write_filtered_input = .true. dfi_write_dfi_history = .false. dfi_bckstop_year = 2000, dfi_bckstop_month = 01, dfi_bckstop_day = 24, dfi_bckstop_hour = 10, dfi_bckstop_minute = 00, dfi_bckstop_second = 00, dfi_fwdstop_year = 2000, dfi_fwdstop_month = 01, dfi_fwdstop_day = 24, dfi_fwdstop_hour = 13, dfi_fwdstop_minute = 00, dfi_fwdstop_second = 00, / &domains time_step_dfi = 60 ** Using gridded nudging option (note this is a separate namelist record) for upperair nudging ** Upper air gridded nudging requires an file generated by the real program. Activating grid_fdda for the real program is adequate. &fdda grid_fdda = 1, 1, 1, gfdda_inname = "wrffdda_d", gfdda_end_h = 24, 24, 24, gfdda_interval_m = 360, 360, 360, fgdt = 0, 0, 0, if_no_pbl_nudging_uv = 0, 0, 0, if_no_pbl_nudging_t = 1, 1, 1, if_no_pbl_nudging_q = 1, 1, 1, if_zfac_uv = 0, 0, 0, k_zfac_uv = 10, 10, 10, if_zfac_t = 0, 0, 0, k_zfac_t = 10, 10, 10, if_zfac_q = 0, 0, 0, k_zfac_q = 10, 10, 10, guv = 0.0003, 0.0003, 0.0003, gt = 0.0003, 0.0003, 0.0003, gq = 0.0003, 0.0003, 0.0003, if_ramping = 1, dtramp_min = 60.0, io_form_gfdda = 2, ** Using gridded surface nudging option (note this is a separate namelist record) ** Note that upper-air and surface gridded nudging may be used together or separately. Surface nudging requires an input file generated by the obsgrid program. &fdda grid_sfdda = 1, 1, 1, sgfdda_inname = "wrfsfdda_d", sgfdda_end_h = 24, 24, 24, sgfdda_interval_m = 360, 360, 360, io_form_sgfdda = 2, guv_sfc = 0.0003, 0.0003, 0.0003, gt_sfc = 0.0003, 0.0003, 0.0003, gq_sfc = 0.0003, 0.0003, 0.0003, rinblw = 250., ** Using observation nudging option (note &fdda is a separate namelist record): &time_control auxinput11_interval_s = 180 , 180 , 180 auxinput11_end_h = 6 , 6 , 6 &fdda obs_nudge_opt = 1,1,1, max_obs = 150000, fdda_start = 0., 0., 0., fdda_end = 720.,720.,720., obs_nudge_wind = 1,1,1, obs_coef_wind = 6.E-4,6.E-4,6.E-4, obs_nudge_temp = 1,1,1, obs_coef_temp = 6.E-4,6.E-4,6.E-4, obs_nudge_mois = 1,1,1, obs_coef_mois = 6.E-4,6.E-4,6.E-4, obs_rinxy = 240.,240.,180., obs_rinsig = 0.1, obs_twindo = 0.6666667,0.6666667,0.6666667, obs_npfi = 10, obs_ionf = 2, 2, 2, obs_idynin = 0, obs_dtramp = 40., obs_prt_freq = 10, 10, 10, obs_prt_max = 10 obs_ipf_errob = .true. obs_ipf_nudob = .true. obs_ipf_in4dob = .true obs_no_pbl_nudge_uv = 0 obs_no_pbl_nudge_t = 0 obs_no_pbl_nudge_q = 0 obs_sfc_scheme_horiz = 0 obs_sfc_scheme_vert = 0 obs_max_sndng_gap = 20 obs_nudgezfullr1_uv = 50 obs_nudgezrampr1_uv = 50 obs_nudgezfullr2_uv = 50 obs_nudgezrampr2_uv = 50 obs_nudgezfullr4_uv = -5000 obs_nudgezrampr4_uv = 50 obs_nudgezfullr1_t = 50 obs_nudgezrampr1_t = 50 obs_nudgezfullr2_t = 50 obs_nudgezrampr2_t = 50 obs_nudgezfullr4_t = -5000 obs_nudgezrampr4_t = 50 obs_nudgezfullr1_q = 50 obs_nudgezrampr1_q = 50 obs_nudgezfullr2_q = 50 obs_nudgezrampr2_q = 50 obs_nudgezfullr4_q = -5000 obs_nudgezrampr4_q = 50 obs_nudgezfullmin = 50 obs_nudgezrampmin = 50 obs_nudgezmax = 3000 obs_sfcfact = 1.0 obs_sfcfacr = 1.0 obs_dpsmx = 7.5 / ** Using spectral nudging option &fdda grid_fdda = 2, 2, 2, gfdda_inname = "wrffdda_d", gfdda_end_h = 24, 24, 24, gfdda_interval_m = 360, 360, 360, fgdt = 0, 0, 0, fgdtzero = 0, 0, 0, if_no_pbl_nudging_uv = 0, 0, 0, if_no_pbl_nudging_t = 0, 0, 0, if_no_pbl_nudging_ph = 0, 0, 0, if_zfac_uv = 0, 0, 0, k_zfac_uv = 10, 10, 10, if_zfac_t = 0, 0, 0, k_zfac_t = 10, 10, 10, if_zfac_ph = 0, 0, 0, k_zfac_ph = 10, 10, 10, dk_zfac_uv = 1, 1, 1, dk_zfac_t = 1, 1, 1, dk_zfac_ph = 1, 1, 1, guv = 0.0003, 0.0003, 0.0003, gt = 0.0003, 0.0003, 0.0003, gph = 0.0003, 0.0003, 0.0003, xwavenum = 3, ywavenum = 3, if_ramping = 1, dtramp_min = 60.0, io_form_gfdda = 2, ** Using adaptive time step option (add these in namelist record &domains): use_adaptive_time_step = .true., step_to_output_time = .true., target_cfl = 1.2, 1.2, 1.2, target_hcfl = .84, .84, .84, max_step_increase_pct = 5, 51, 51, starting_time_step = -1, -1, -1, max_time_step = 360, 120, 40, min_time_step = 90, 30, 10, adaptation_domain = 1, ** Using automatic vortex-following option (tropical storm tracking only; add these in namelist record &domains): vortex_interval = 15, 15, 15, max_vortex_speed = 40, 40, 40, corral_dist = 8, 15, 15, track_level = 50000, ** Miscellaneous physics options for namelist record &physics: Topographic shading (only effective when grid sizes are a few kilometers) slope_rad = 1, 1, 1, topo_shading = 1, 1, 1, shadlen = 25000, Setting threashold value for defining seaice if seaice is not in the input file: seaice_threshold = 271, Switching off latent heating from a microphysics scheme (must also set cu_physics = 0): no_mp_heating = 0, ** Using precipiatiion bucket in a time interval (minutes): prec_acc_dt = 60. ** Using bucket accumulations for multi-year simulations (guideline: mean monthly accumulation) bucket_mm = 100. bucket_J = 1.e9 ** Optional gravitational settling of fog/cloud droplets (MYNN PBL only) grav_settling = 1, ; default 0 ** Using implicit gravity-wave damping option (add these in namelist record &dynamics): damp_opt = 3, zdamp = 5000., 5000., 5000., dampcoef = 0.2, 0.2, 0.2 ** Using expanded boundary zone and exponential decay option (add or modify these in namelist record &bdy_control): spec_bdy_width = 10, spec_zone = 1, relax_zone = 9, specified = .true., .false.,.false., spec_exp = 0.33 For a tropical channel configuration, set the following: specified = .true., .false.,.false., periodic_x = .true., .false.,.false., ** using io quilting option to improve output efficiency for large domain runs (note that this is a separate namelist record): &namelist_quilt nio_tasks_per_group = 2, nio_groups = 1, ** for tc bogusing: &tc insert_bogus_storm = .true., remove_storm = .true., num_storm = 1, latc_loc = 15., lonc_loc = -90., vmax_meters_per_second = 30, rmax = 50000, vmax_ratio = 0.5, / ** for regional climate surface diagnostics such as max/min/mean/std of T2/Q2/wind/rainfall between selected output times (e.g. daily) in auxhist3 &time_control output_diagnostics = 1 auxhist3_outname = 'wrfxtrm_d_' io_form_auxhist3 = 2 auxhist3_interval = 1440 frames_per_auxhist3 = 1 / ** for pressure-level (and some surface) diagnostics, output is on stream 23, where the listed pressure levels are in Pa. For the height level interpolation, the unit is stream 22. &time_control io_form_auxhist23 = 2, auxhist23_interval = 30, 30, 30, frames_per_auxhist23 = 1, 1, 1, auxhist23_outname = "PLEVS_d_" / &diags p_lev_diags = 1 num_press_levels = 4 press_levels = 85000, 70000, 50000, 25000 use_tot_or_hyd_p = 2 / ** for height-level diagnostics, output is on stream 22 (negative values for z_levels means AGL, so -500 is 500 m AGL, and 10000 is 10 km). &time_control io_form_auxhist22 = 2, auxhist22_interval = 30, 30, 30, frames_per_auxhist22 = 1, 1, 1, auxhist22_outname = "ZLEVS_d_" / &diags z_lev_diags = 1 num_z_levels = 2 z_levels = -500, 10000 / ** Using different flux formulation for tropical storm simulations (best for grid spacing less than 2 km) simple 1-D ocean mixed layer, or University of Miami 3DPWP ocean model (add these in namelist record &physics): isftcflx = 1, sf_ocean_physics = 0 (off), 1 (OML), 2 (3D PWP) oml_hml0 = 50, oml_gamma = 0.14 omdt = 1 Add these in the &domains section. Note that the example here is for a warm ocean: ocean_z is the depth of each layer (m), from the surface downward; the associated temperature is ocean_t (K); and the associated salinity is ocean_s (ppt). ocean_levels = 30, ocean_z = 5., 15., 25., 35., 45., 55., 65., 75., 85., 95., 105., 115., 125., 135., 145., 155., 165., 175., 185., 195., 210., 230., 250., 270., 290., 310., 330., 350., 370., 390. ocean_t = 302.3493, 302.3493, 302.3493, 302.1055, 301.9763, 301.6818, 301.2220, 300.7531, 300.1200, 299.4778, 298.7443, 297.9194, 297.0883, 296.1443, 295.1941, 294.1979, 293.1558, 292.1136, 291.0714, 290.0293, 288.7377, 287.1967, 285.6557, 284.8503, 284.0450, 283.4316, 283.0102, 282.5888, 282.1674, 281.7461 ocean_s = 34.0127, 34.0127, 34.0127, 34.3217, 34.2624, 34.2632, 34.3240, 34.3824, 34.3980, 34.4113, 34.4220, 34.4303, 34.6173, 34.6409, 34.6535, 34.6550, 34.6565, 34.6527, 34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247, 34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247 References for 3D ocean circulation model: Lee, C.-Y., and S. S. Chen, 2013: Stable Boundary Layer and Its Impact on Tropical Cyclone Structure in a Coupled Atmosphere-Ocean Model, Mon. Wea. Rev. Price, J. F., T. B. Sanford, and G. Z. Forristal, 1994: Forced stage response to a moving hurricane. J. Phy. Oceanogr., 24, 233-260. ** Using U. Miami Forward Lagrangian trajectory calculation (add it in namelist record &physics): &domain num_traj = 25, &physics traj_opt = 1, ** Vertical nesting The WRF model now supports vertical nesting for a coincident (online) model simulation (during a single model run, differing numbers of vertical levels may be used per domain). This is activated with a switch to turn on the option (vert_refine_method). The namelist array eta_levels is manually filled in for each domain. Below is an example for two domains. NOTE: The user is restricted to either the RRTM or RRTMG radiation schemes. &domains max_dom = 2, e_vert = 35, 45, eta_levels(1:35) = 1., 0.993, 0.983, 0.97, 0.954, 0.934, 0.909, 0.88, 0.8406663, 0.8013327, 0.761999, 0.7226653, 0.6525755, 0.5877361, 0.5278192, 0.472514, 0.4215262, 0.3745775, 0.3314044, 0.2917579, 0.2554026, 0.2221162, 0.1916888, 0.1639222, 0.1386297, 0.1156351, 0.09525016, 0.07733481, 0.06158983, 0.04775231, 0.03559115, 0.02490328, 0.0155102, 0.007255059, 0. eta_levels(36:81) = 1.0000, 0.9946, 0.9875, 0.9789, 0.9685, 0.9562, 0.9413, 0.9238, 0.9037, 0.8813, 0.8514, 0.8210, 0.7906, 0.7602, 0.7298, 0.6812, 0.6290, 0.5796, 0.5333, 0.4901, 0.4493, 0.4109, 0.3746, 0.3412, 0.3098, 0.2802, 0.2524, 0.2267, 0.2028, 0.1803, 0.1593, 0.1398, 0.1219, 0.1054, 0.0904, 0.0766, 0.0645, 0.0534, 0.0433, 0.0341, 0.0259, 0.0185, 0.0118, 0.0056, 0. vert_refine_method = 0, 2, ** Tropopause data level of max winds data, for program real.exe only When information (mostly NCEP supplied through GFS or NAM) is available in the Grib2 file, AND extracted with ungrib, the metgrid program inserts flag information into the data stream that is input by real. The real program is able to use the available u, v, T, height fields (each may be on a tropopause .OR. the level of max winds) in the vertical interpolation. To "tune" the data, the user may select a level below which the max wind inforamtion is ignored. The default is 300 hPa. The user may also select the level that when exceeded the trop and maxw fields are ignored (due to the horizontal pressure difference detecting a user-defined discontinuity, and the metgrid horizontal interpolation across the discontinuity would be suspect). By default, this option is turned off. The user may choose to individually activate the vertical interpolation of either the level of max winds .OR. the tropopause level. Default (all pressure units Pa): &domains maxw_horiz_pres_diff = 5000 trop_horiz_pres_diff = 5000 maxw_above_this_level = 30000 use_maxw_level = 0 (0=do not use level, 1 = use level) use_trop_level = 0 (0=do not use level, 1 = use level) ** Using aerosol option aer_opt = 2: &physics aer_opt = 2, aer_type = 1, aer_aod550_opt = 1, aer_aod550_val = 0.12, aer_angexp_opt = 1, aer_angexp_val = 1.3, aer_ssa_opt = 1, aer_ssa_val = 0.85, aer_asy_opt = 1, aer_asy_val = 0.90, ** Using Jimenez wind-farm scheme In the &physics namelist record, set the MAX_DOM value: windfarm = 1, 1, 1 Also in the directory with namelist.input, is the file defining the specifics of the wind turbine type: wind-turbine-1.tbl The location of wind turbines are specified as lat lon, and the turbine type: windturbines.txt ** Using stochastic schemes &stoch skebs = 1, 1, 1, rand_perturb = 1, 1, 1, ** Using kfcupscheme cu_physics (cumulus scheme option). This has been tested with the CAM radiation scheme and is not recommended with other radiation schemes. Regarding cu_rad_feedback, users want the parameterized clouds to affect radiation. Turning it off is only useful as a sensitivity study to determine the importance of that effect. Regarding shallowcu_forced_ra option, setting it to true will override the cloud fraction calculations to a prescribed maximum cloud fraction (a value of 0.36) which can be changed by the user for sensitivity testing purposes. &physics cu_physics = 10, 10, 10, ra_lw_physics = 3, 3, 3, ra_sw_physics = 3, 3, 3, cu_rad_feedback =.true.,.true.,.true., shallowcu_forced_ra = .false.,.false.,.false., numBins = 21, 21, 21, thBinSize = 0.1, 0.1, 0.1, rBinSize = 0.0001,0.0001,0.0001, minDeepFreq = 0.333, 0.333, 0.333, minShallowFreq = 0.01, 0.01, 0.01, ** To write out an input file during the model simulation. &time_control write_input = .TRUE., inputout_interval_m = 60, 60, 60 input_outname = "wrfinput_out_d_" inputout_begin_h = 3, 6, 6 inputout_end_h = 9, 6, 6 ** To use climatological data in Thompson microphysics option 28: in addition to process additional data in metgrid, use these: &domains wif_input_opt = 1 num_wif_levels = 30 &physics mp_physics = 28, 28, 28, use_aero_icbc = .true. ** The hybrid vertical coordinate (HVC) requires three pieces: 1. The WRF and real codes must be built with HVC activated (./configure -hyb) 2. The eta location (etac) at which the eta levels higher up in the atmosphere become isobaric must be defined (suggested default in Registry). 3. The run-time option to select the HVC (hybrid_opt, default is OFF): &dynamics hybrid_opt = 2, etac = 0.2, ** New mechanism to specify physics: &physics physics_suite = 'CONUS' which is equivalent to mp_physics = 8, cu_physics = 6, ra_lw_physics = 4, ra_sw_physics = 4, bl_pbl_physics = 2, sf_sfclay_physics = 2, sf_surface_physics = 2, &physics physics_suite = 'TROPICAL' which is equivalent to mp_physics = 6, cu_physics = 16, ra_lw_physics = 4, ra_sw_physics = 4, bl_pbl_physics = 1, sf_sfclay_physics = 91, sf_surface_physics = 2, To overwrite the cu_physics option for a second nest, set &physics physics_suite = 'tropical' cu_physics = -1, -1, 0, ...