.. _gsi_namelist: C. GSI Namelist: Name, Default Value, Explanation ================================================== The following are lists and explanations of the GSI namelist variables. You can also find them in the source code **gsimod.F90**. ==================== ============= ========================================================= Variable name Default value Description ==================== ============= ========================================================= **SETUP** **General control namelist** gencode 80 source generation code factqmin 1 weighting factor for negative moisture constraint factqmax 1 weighting factor for supersaturated moisture constraint clip_supersaturation .false. flag to remove supersaturation during each outer loop factv 1 weighting factor for negative visibility constraint factl factp factg factw10m facthowv deltim 1200 model timestep dtphys 3600 physics timestep biascor -1 background error bias correction coefficient bcoption 1 0=ibc (no bias correction to bkg); \ 1= sbc(original implementation) diurnalbc 0 1= diurnal bias; 0= persistent bias ==================== ============= ========================================================= ================= ============= ================================================================================= Variable name Default value Description ================= ============= ================================================================================= **SETUP** **General control namelist** niter(0:50) 0,... Maximum number of inner loop iterations for each outer loop niter_no_qc(0:50) 1000000 Inner loop iteration at which to turn on variational quality control miter 1 number of outer loops qoption 1 option for moisture analysis variable; 1:q/qsatg 2:normalized RH cwoption pseudo_q2 .false. breed between q1/q2 options, that is, (q1/sig(q)) nhr_assimilation 6 assimilation time interval (currently 6 hours for global, 3 hours for regional ) min_offset 3 time of analysis in assimilation window iout_iter 220 output file number for iteration information npredp 6 number of predictors for precipitation bias correction retrieval .false. logical to turn off or on the SST physical retrieval nst_gsi 0 indicator to control the Tr Analysis mode: \ 0 = no nst info ingsi at all; \ 1 = input nst info, but used for monitoring only \ 2 = input nst info, and used in CRTM simulation, but no Tr analysis \ 3 = input nst info, and used in CRTM simulation and Tr analysis is on nst_tzr 0 indicator to control the Tzr_QC mode: \ 0 = no Tz retrieval; \ 1 = Do Tz retrieval and applied to QC nstinfo 0 number of nst variables fac_dtl 0 index to apply diurnal thermocline layer or not: 0 = no; 1 = yes fac_tsl 0 index to apply thermal skin layer or not: 0 = no; 1 = yes. nst_tzr tzr_bufrsave .false. logical to turn off or on the bufr Tz retrieval file true=on diag_rad .true. logical to turn off or on the diagnostic radiance file (true=on) diag_pcp .true. logical to turn off or on the diagnostic precipitation file (true=on) diag_conv .true. logical to turn off or on the diagnostic conventional file (true=on) diag_ozone .true. logical to turn off or on the diagnostic ozone file (true=on) diag_aero .false. logical to turn off or on the diagnostic aerosol file (true=on) diag_co .false. logical to turn off or on the diagnostic carbon monoxide file (true=on) iguess 1 flag for guess solution (currently not working) \ -1 do not use guess file \ 0 write only guess file \ 1 read and write guess file \ 2 read only guess file write_diag .false., ... logical to write out diagnostic files for outer iteration reduce_diag .false. namelist logical to produce reduced radiance diagnostic files ================= ============= ================================================================================= ===================== ============= ========================================================================================================== Variable name Default value Description ===================== ============= ========================================================================================================== **SETUP** **General control namelist** oneobtest .false. one observation test flag true=on sfcmodel .false. if true, then use boundary layer forward model for surface temperature data. dtbduv_on .true. logical for switching on (.true.) sensitivity of uv winds to microwave brightness temperatures. if true, use d(microwave brightness temperature)/d(uv wind) in inner loop ifact10 0 flag for recomputing 10m wind factor \ = 1 compute using GFS surface physics \ = 2 compute using MM5 surface physics \ = 0 or any other value - DO NOT recompute - use value from guess file l_foto .false. option for First-Order Time extrapolation to observation offtime_data .false. if true, then allow use of obs files with ref time different from analysis time. default value = .false., in which case analysis fails if observation file reference time is different from analysis time. npred_conv_max 0 maximum number of conventional observation bias correction coefficients id_bias_ps 0 prepbufr id to have conv_bias added for testing id_bias_t 0 prepbufr id to have conv_bias added for testing id_bias_spd 120 prepbufr id to have conv_bias added for testing conv_bias_ps 0 magnitude of ps bias(mb) conv_bias_t 0 magnitude of t bias(deg K) conv_bias_spd 0 magnitude of spd bias(m/sec) id_bias_pm2_5 conv_bias_pm2_5 id_bias_pm10 conv_bias_pm10 stndev_conv_ps 1.0 stndev_conv_t 1.0 stndev_conv_spd 1.0 use_pbl .false. Logical flag to include PBL effects in tendency model. use_compress .false. option to turn on the use of compressibility factors in geopotential heights nsig_ext 13 number of layers above the model top which are necessary to compute the bending angle for gpsro gpstop 30.0 maximum height for gpsro data assimilation. Reject anything above this height. (km) perturb_obs .false. logical flag to perturb observation (true=on) perturb_fact 1 magnitude factor for observation perturbation oberror_tune .false. logical to tune (=true) oberror preserve_restart_date .false. if true, then do not update regional restart file date. crtm_coeffs_path ./ path of directory w/ CRTM coeffs files berror_stats berror_stats filename if other than "berror_stats" newpc4pred .false. option for additional preconditioning for pred coeff adp_anglebc .false. option to perform variational angle bias correction angord 0 order of polynomial for variational angle bias correction ===================== ============= ========================================================================================================== ============== ============= ===================================================================================================== Variable name Default value Description ============== ============= ===================================================================================================== **SETUP** **General control namelist** passive_bc .false. option to turn on bias correction for passive (monitored) channels use_edges .true. option to exclude radiance data on scan edges biaspredvar 0.1 set background error variance for radiance bias coeffs lobsdiagsave .false. write out additional observation diagnostics l4dvar .false. turn 4D-Var on/off (default=off=3D-Var) lbicg .false. use B-precond w/ bi-conjugate gradient for minimization lsqrtb .false. Use sqrt(B) preconditioning lcongrad .false. Use conjugate gradient/Lanczos minimizer lbfgsmin .false. Use L-BFGS minimizer ltlint .false. Use TL inner loop (ie TL intall) nhr_obsbin -1 length of observation bins nhr_subwin -1 length of weak constraint 4d-Var sub-window intervals nwrvecs -1 Number of precond vectors (Lanczos) or pairs of vectors (QN) being saved iorthomax 0 max number of vectors used for orthogonalization of various CG options ladtest .false. Run adjoint test ladtest_obs .false. if true, doing the adjoint check for the observation operators lgrtest .false. Run gradient test lobskeep .false. keep obs from first outer loop for subsequent OL lsensrecompute .false. does adjoint by recomputing forward solution jsiga -1 calculate approximate analysis errors from lanczos for jiter=jsiga ltcost .false. calculate true cost when using Lanczos (this is very expensive) lobsensfc .false. compute forecast sensitivity to observations lobsensjb .false. compute Jb sensitivity to observations lobsensincr .false. compute increment sensitivity to observations lobsensadj .false. use adjoint of approx. Hessian to compute obs sensitivity lobsensmin .false. use minimisation to compute obs sensitivity iobsconv 0 compute convergence test in observation space \ =1 at final point, =2 at every iteration idmodel .false. uses identity model when running 4D-Var (test purposes) iwrtinc .false. when .t., writes out increments instead of analysis jiterstart 1 first outloop iteration number jiterend 1 last outloop iteration number lobserver .false. when .t., calculate departure vectors only lanczosave .false. save lanczos vectors for forecast sensitivity computation llancdone .false. use to tell adjoint that Lanczos vecs have been pre-computed lferrscale .false. Something related to forecast error print_diag_pcg .false. logical turn on of printing of GMAO diagnostics in pcgsoi.f90 tsensible .false. option to use sensible temperature as the analysis variable. Works only for twodvar_regional=.true. lgschmidt .false. option for re-biorthogonalization of the gradx and grady set from pcgsoi when twodvar_regional=.true. ============== ============= ===================================================================================================== ==================== ============= ======================================================================================================================================= Variable name Default value Description ==================== ============= ======================================================================================================================================= **SETUP** **General control namelist** lread_obs_save .false. option to write out collective obs selection info lread_obs_skip .false. option to read in collective obs selection info use_gfs_ozone .false. option to read in gfs ozone and interpolate to regional model domain check_gfs_ozone_date .false. option to date check gfs ozone before interpolating to regional model domain regional_ozone .false. option to turn on ozone in regional analysis lwrite_predterms .false. option to write out actual predictor terms instead of predicted bias to the radiance diagnostic files lwrite_peakwt .false. option to writ out the approximate pressure of the peak of the weighting function for satellite data to the radiance diagnostic files use_gfs_nemsio .false. option to use nemsio to read global model NEMS/GFS first guess liauon .false. treat 4dvar CV as tendency perturbation (default=false) use_prepb_satwnd .false. allow using satwnd’s from prepbufr (historical) file l4densvar .false. logical to turn on ensemble 4dvar ens4d_nstarthr 3 start hour for ensemble perturbations (generally should match min_offset) use_gfs_stratosphere When true, a guess gfs valid at the same time as the nems-nmmb guess is used to replace the upper levels with gfs values. The purpose of this is to allow direct use of gdas derived sat radiance bias correction coefs. pblend0 152 The nems-nmmb vertical coordinate is smoothly merged with gfs above this level. Below this level, is original nems-nmmb. pblend1 79.0 The nems-nmmb vertical coordinate is smoothly merged with gfs below this level. Above this level,is gfs. step_start 1.e-4 initial stepsize in minimization diag_precon .false. if true do preconditioning lrun_subdirs .false. logical to toggle use of subdirectires at runtime for pe specific files emiss_bc .false. option to turn on emissivity bias predictor upd_pred 1 bias update indicator for radiance bias correction; 1.0=bias correction coefficients evolve use_reflectivity .false. option of using reflectivity lnested_loops .false. allow for nested resolution outer/inner loops lwrite4danl .false. logical to write out 4d analysis states if 4dvar or 4denvar mode lsingleradob .false. logical for single radiance observation assimilation. Uses existing bufr file and rejects all radiances that don't fall within a tight threshold around oblat/oblon (SINGLEOB_TEST) ssmis_method 1 choose method for SSMIS noise reduction 0=no smoothing 1=default ssmis_precond 0.01 weighting factor for SSMIS preconditioning (if not using newpc4pred) R_option .false. Option to use variable correlation length for lcbas based on data density - follows Hayden and Purser (1995) (twodvar_regional only) thin4d ==================== ============= ======================================================================================================================================= =================== ============= ==================================================================================================== Variable name Default value Description =================== ============= ==================================================================================================== **GRIDOPTS** **Grid setup variables, including regional specific variables** jcap 62 spectral resolution of the analysis jcap_b 62 spectral resolution of background (model guess field) nsig 42 number of sigma levels nlat 96 number of latitudes nlon 384 number of longitudes hybrid logical hybrid data file flag true=hybrid nlat_regional 0 Number of y grid point in whole regional domain nlon_regional 0 Number of x grid point in whole regional domain diagnostic_reg .false. logical for regional debugging update_regsfc .false. logical to write out updated surface fields to the regional analysis file (default = false) netcdf .false. if true, then wrf files are in netcdf format, otherwise wrf files are in binary format. regional .false. logical for regional GSI run wrf_nmm_regional .false. logical for input from WRF NMM nems_nmmb_regional .false. logical for input from NEMS NMMB wrf_mass_regional .false. logical for input from WRF MASS-CORE (ARW) twodvar_regional .false. logical for regional 2d-var analysis filled_grid .false. logical to fill in points on WRF-NMM E-grid half_grid .false. logical to use every other row of WRF-NMM E-Grid nvege_type 24 number of types of vegetation; old=24, IGBP=20 nlayers(100) 1 number of sub-layers to break indicated model layer into prior to calling radiative transfer model cmaq_regional .false. Background input is from CMAQ model nmmb_reference_grid H =’H’, then analysis grid covers H grid domain \ = ’V’, then analysis grid covers V grid domain grid_ratio_nmmb sqrt(2) ratio of analysis grid to nmmb model grid in nmmb model grid units. grid_ratio_wrfmass 1.0 ratio of analysis grid to wrf mass grid in wrf grid units jcap_gfs spectral truncation used to transform high wavenumber spectral coefficients to a coarser resolution grid,when use_gfs_ozone = .true. or use_gfs_stratosphere = .true. jcap_cut use_sp_eqspac .false. if .true., then ensemble grid is equal spaced, staggered 1/2 grid unit off poles. if .false., then gaussian grid assumed for ensemble (global only) =================== ============= ==================================================================================================== ============= ============= ============================================================================================================ Variable name Default value Description ============= ============= ============================================================================================================ **BKGERR** **Background error related variables** vs 1/1.5 scale factor for vertical correlation lengths for background error nhscrf 3 number of horizontal scales for recursive filter hzscl(3) 1, 1, 1 scale factor for horizontal smoothing, n=1,number of scales (3 for now) \ specifies factor by which to reduce horizontal scales (i.e. 2 would then apply 1/2 of the horizontal scale) hswgt(3) 1/3, 1/3, 1/3 empirical weights to apply to each horizontal scale norh 2 order of interpolation in smoothing ndeg 4 degree of smoothing in recursive filters noq 3 1/4 of accuracy in compact finite differencing bw 0 factor in background error calculation norsp 0 order of interpolation for smooth polar cascade routine default is norsp=0, in which case norh is used with original polar cascade interpolation (global only). fstat .false. logical to separate f from balance projection pert_berr .false. logical to turn on random inflation/deflation of background error tuning parameters pert_berr_fct 0 factor for increasing/decreasing berror parameters, this is multiplied by random number bkgv_flowdep .false. flag to turn on flow dependence to background error variances bkgv_rewgtfct 0 factor used to perform flow dependent reweighting of error variances bkgv_write .false. flag to turn on=.true. /off=.false. generation of binary file with reweighted variances fpsproj .true. controls full nsig projection to surface pressure fut2ps controls the projection from unbalance T to surface pressure adjustozvar adjusts ozone variances in the stratosphere based on guess field cwcoveqqcov sets cw Bcov to be the same as B-cov(q) (presently glb default) ============= ============= ============================================================================================================ ===================== ============= ========================================================================================================= Variable name Default value Description ===================== ============= ========================================================================================================= **ANBKGERR** **Anisotropic background error related variables** anisotropic .false. if true, then use anisotropic background error covariance ancovmdl 0 covariance model settings - 0: pt-based, 1: ensemble based triad4 .true. for 2d variables, if true, use blended triad algorithm ifilt_ord 4 filter order for anisotropic filters npass 1 2×npass = number of factors in background error normal 200 number of random vectors to use for filter normalization ( if < 0 then slightly slower, but results independent of number of processors) binom .true. if true, weight correlation lengths of factors using binomial distribution, with shortest scales on outside, longest scales on inside. This can help to produce smoother correlations in the presence of strong anisotropy ngauss 3 number of Gaussians to add together in each factor rgauss 0 multipliers on reference aspect tensor for each Gaussian factor anhswgt 1.0 empirical weights to apply to each gaussian an_vs 1 scale factor for background error vertical scales (temporary carry over from isotropic inhomogeneous option) grid_ratio 2.0 ratio of coarse to fine grid in fine grid units grid_ratio_p 0 ratio of coarse to fine grid in fine grid units for polar patches nord_f2a 4 order of interpolation for transfer operators between filter grid and analysis grid an_flen_u 1 coupling parameter for connecting horizontal wind to background error an_flen_t 1 coupling parameter for connecting grad(potential temperature) to background error an_flen_z 1 coupling parameter for connecting grad(terrain) to background error rtma_subdomain_option .false. if true, then call alternative code which calls recursive filter directly from subdomain mode, bypassing transition to/from horizontal slabs. This is mainly to improve efficiency for 2d rtma analysis. at the moment, this only works for twodvar_regional=.true. rtma_subdomain_option will be forced to false when twodvar_regional=.false. lreadnorm .false. if true, then read normalization from fixed files nsmooth 0 number of 1-2-1 smoothing passes before and after background error application nsmooth_shapiro 0 number of 2nd moment preserving (shapiro) smoothing passes before and after background error application. \ NOTE: default for nsmooth and nsmooth_shapiro is 0. \ if both are > 0, then nsmooth will be forced to zero. afact0 0.0 anistropy effect parameter, the range must be in 0.0-1.0. covmap .false. if true, covariance map would be drawn ===================== ============= ========================================================================================================= ===================== ============= ========================================================================================================== Variable name Default value Description ===================== ============= ========================================================================================================== **JCOPTS** **Constraint term in cost function (Jc)** ljcdfi .false. if .false., uses original formulation based on wind, temp, and ps tends \ when .t. uses digital filter initialization of increments (4dvar) alphajc 10.0 parameter for digital filter switch_on_derivatives .false., … if true, then compute horizontal derivatives of all state variables (to be used eventually for time derivatives, dynamic constraints and observation forward models that need horizontal derivatives) tendsflag .false. if true, compute time tendencies ljcpdry .false. when .t. uses dry pressure constraint on increment bamp_jcpdry 0.0 parameter for pdry_jc eps_eer -1.0 Errico-Ehrendofer parameter for q-term in energy norm ljc4tlevs .false. when true and in 4D mode, apply any weak constraints over all time levels instead of just at a single time ===================== ============= ========================================================================================================== ============== ============= ============================================================================================================= Variable name Default value Description ============== ============= ============================================================================================================= **STRONGOPTS** **Strong dynamic constraint** reg_tlnmc_type 1 =1 for 1st version of regional strong constraint \ =2 for 2nd version of regional strong constraint tlnmc_option 0 integer flag for strong constraint (various capabilities for hybrid): \ =0: no TLNMC \ =1: TLNMC for 3DVAR mode \ =2: TLNMC on total increment for single time level only (for 3D EnVar) or if 4D EnVar mode, TLNMC applied \ to increment in center of window \ =3: TLNMC on total increment over all time levels (if in 4D EnVar mode) \ =4: TLNMC on static contribution to increment ONLY for any EnVar mode nstrong 0 if > 0, then number of iterations of implicit normal mode initialization to apply for each inner loop iteration period_max 1000000.0 cutoff period for gravity waves included in implicit normal mode initialization (units = hours) period_width 1.0 defines width of transition zone from included to excluded gravity waves nvmodes_keep 0 number of vertical modes to use in implicit normal mode initialization baldiag_full .false. flag to toggle balance diagnostics for the full fields baldiag_inc .false. flag to toggle balance diagnostics for the analysis increment ============== ============= ============================================================================================================= ================= ============= ============================================================================================== Variable name Default value Description ================= ============= ============================================================================================== **OBSQC** **Observation quality control variables** \ Parameters used for gross error checks are set in file convinfo (ermin, ermax, \ ratio) Parameters below used for nonlinear (variational) quality control dfact 0 factor for duplicate observation at same location for conventional data dfact1 3.0 time factor for duplicate observation at same location for conventional data erradar_inflate 1 radar error inflation factor tdrerr_inflate .false. logical for tdr obs error inflation tdrgross_fact 1 factor applied to tdr gross error oberrflg .false. logical for reading in new observation error table (if set to true) vadfile ’none’ character(10) variable holding name of VAD wind bufr file noiqc .false. logical flag to bypass OI QC (if set to true) c_varqc 1 constant number to control variance qc turning on speed blacklst .false. logical for reading in raob blacklist (if set to true) use_poq7 .false. Logical to toggle accept (.true.) or reject (.false.) SBUV/2 ozone observations flagged with profile ozone quality mark hilbert_curve .false. option for hilbert-curve based cross-validation. works only with twodvar_regional=.true. tcp_refps 1000.0 reference pressure for tcps oberr calculation (mb) tcp_width 50.0 parameter for tcps oberr inflation (width, mb) tcp_ermin 0.75 parameter for tcps oberr inflation (minimum oberr, mb) tcp_ermax 5.0 parameter for tcps oberr inflation (maximum oberr, mb) qc_noirjaco3 .false. controls whether to use O3 Jac from IR instruments qc_noirjaco3_pole .false. controls wheter to use O3 Jac from IR instruments near poles qc_satwnds .true. allow bypass sat-winds qc normally removing lots of mid-tropo obs njqc vqc aircraft_t_bc_pof .false. logical for aircraft temperature bias correction, pof is used for predictor aircraft_t_bc .false. logical for aircraft temperature bias correction aircraft_t_bc_ext .false. logical for reading aircraft temperature bias correction from external file buddycheck_t .false. When true, run buddy check algorithm on temperature observations buddydiag_save .false. When true, output files containing buddy check QC info for all obs run through the buddy check biaspredt 1 berror var for temperature bias correction coefficients upd_aircraft .true. indicator if update bias at 06Z & 18Z cleanup_tail .false. logical to remove tail number no longer used ================= ============= ============================================================================================== ================= ============= ============================================================================================================ Variable name Default value Description ================= ============= ============================================================================================================ **OBS_INPUT** **Controls input data** dfile ``' '`` input observation file name dtype ``' '`` observation type dplat ``' '`` satellite (platform) id (for satellite data) dsis ``' '`` sensor/instrument/satellite flag from satinfo files dthin ``' '`` satellite group dval ``' '`` relative value of each profile within group relative weight for observation = dval/sum(dval) within grid box dmesh(max(dthin)) thinning mesh for each group \ mesh size (km) for radiance thinning grid (used in satthin) dsfcalc ``' '`` specifies method to determine surface fields within a FOV. when equal to one, integrate model fields over FOV. when not one, bilinearly interpolate model fields to FOV center. time_window_max 3 upper limit on time window for all input data ext_sonde .false. logical for extended forward model on sonde data l_foreaft_thin .false. separate TDR fore/aft scan for thinning ================= ============= ============================================================================================================ ================= ============= ================================================================================================== Variable name Default value Description ================= ============= ================================================================================================== **SINGLEOB_TEST** **Single observation test case setup** maginnov 1 magnitude of innovation for one observation magoberr 1 magnitude of observational error oneob_type ``' '`` observation type (t, u, v, etc.) oblat 0 observation latitude oblon 0 observation longitude obpres 1000.0 observation pressure (hPa) obdattim 2000010100 observation date (YYYYMMDDHH) obhourset 0 observation delta time from analysis time pctswitch .false. if .true. innovation & oberr are relative (%) of background value (level ozone only) obchan 0 if > 0, selects the channel number. If <= zero, it will use all channels that pass qc in setuprad. ================= ============= ================================================================================================== ================= ============= ===================================================================================== Variable name Default value Description ================= ============= ===================================================================================== **SUPEROB_RADAR** **Level 2 bufr file to radar wind superobs** del_azimuth 5.0 azimuth range for superob box (default 5 degrees) del_elev 0.25 elevation angle range for superob box (default .05 degrees) del_range 5000.0 radial range for superob box (default 5 km) del_time 0.5 1/2 time range for superob box (default .5 hours) elev_angle_max 5.0 max elevation angle (default of 5 deg) minnum 50 minimum number of samples needed to make a superob range_max 100000.0 max radial range in meters to use in constructing superobs (default 100km) l2superob_only .false. if true, then process level 2 data creating superobs, then quit. (added for easier retrospective testing, since level 2 bufr files are very large and hard to work with) ================= ============= ===================================================================================== ==================== ================ =============================================================== Variable name Default value Description ==================== ================ =============================================================== **LAG_DATA\|** **Lagrangian data assimilation related variables** lag_accur 1.0e-6 Accuracy used to decide whether or not a balloon is on the grid infile_lag inistate_lag.dat File containing the initial position of the balloon lag_stepduration 900.0 Duration of one time step for the propagation model lag_nmax_bal 1000 Maximum number of balloons at starting time lag_vorcore_stderr_a 2.0e3 Observation error for vorcore balloon lag_vorcore_stderr_b 0.0 error = b + a*timestep(in hours) ==================== ================ =============================================================== ======================== ============= ================================================================================================= Variable name Default value Description ======================== ============= ================================================================================================= **HYBRID_ENSEMBLE** **Parameters for use with hybrid ensemble option** l_hyb_ens .false. if true, then turn on hybrid ensemble option uv_hyb_ens .false. if true, then ensemble perturbation wind variables are u,v, otherwise, ensemble \ perturbation wind variables are stream, pot. Functions. q_hyb_ens .false. if true, then use specific humidity ensemble perturbations, otherwise, use relative humidity aniso_a_en .false. if true, then use anisotropic localization of hybrid ensemble control variable a_en. generate_ens .true. if true, then generate internal ensemble based on existing background error n_ens 0 number of ensemble members. nlon_ens 0 number of longitudes on ensemble grid (may be different from analysis grid nlon) nlat_ens 0 number of latitudes on ensemble grid (may be different from analysis grid nlat) jcap_ens 0 for global spectral model, spectral truncation pseudo_hybens .false. if true, turn on pseudo ensemble hybrid for HWRF merge_two_grid_ensperts .false. if true, merge ensemble perturbations from two forecast domains to analysis domain (one way \ to deal with hybrid DA for HWRF moving nest) regional_ensemble_option 0 integer, used to select type of ensemble to read in for regional application. Currently takes \ values from 1 to 4 \ =1: use GEFS internally interpolated to ensemble grid. \ =2: ensembles are WRF NMM format \ =3: ensembles are ARW netcdf format. \ =4: ensembles are NEMS NMMB format. full_ensemble .false. if true, first ensemble perturbation on first guess istead of on ens mean betaflg .false. if true, use vertical weighting on beta1_inv and beta2_inv, for regional coef_bw 0.9 fraction of weight given to the vertical boundaries when betaflg is true pwgtflg .false. if true, use vertical integration function on ensemble contribution of Psfc jcap_ens_test 0 for global spectral model, test spectral truncation (to test dual resolution) beta1_inv 1 1/beta1, the default weight given to static background error covariance if (.not. readin_beta) \ 0 <= beta1_inv <= 1, tuned for optimal performance \ =1, then ensemble information turned off \ =0, then static background turned off the weights are applied per vertical level such that : \ betas_inv(:) = beta1_inv , vertically varying weights given to static ; \ betae_inv(:) = 1 - beta1_inv , vertically varying weights given ensemble derived covariance. \ If (readin_beta) then betas_inv and betae_inv are read from a file and beta1_inv is not used. ======================== ============= ================================================================================================= ===================== ============= =========================================================================================================== Variable name Default value Description ===================== ============= =========================================================================================================== **HYBRID_ENSEMBLE** **Parameters for use with hybrid ensemble option** s_ens_h 2828 homogeneous isotropic horizontal ensemble localization scale (km) s_ens_v 30 vertical localization scale (grid units for now) s_ens_h, s_ens_v, and beta1_inv are tunable parameters. use_gfs_ens .true. controls use of global ensemble: .t. use GFS (default); .f. uses user-defined ens readin_localization .false. flag to read (.true.)external localization information file readin_beta .false. flag to read (.true.) the vertically varying beta parameters betas_inv and betae_inv from a file. eqspace_ensgrid .false. if .true., then ensemble grid is equal spaced, staggered 1/2 grid unit off ploes. \ if .false., then gaussian grid assumed for ensemble (global only) use_localization_grid .false. if true, then use extra lower res gaussian grid for horizontal localization (global runs only–allows \ possiblity for non-gaussian ensemble grid) grid_ratio_ens 1 for regional runs, ratio of ensemble grid resolution to analysis grid resolution \ default value = 1 (dual resolution off) oz_univ_static .false. if true, decouple ozone from other variables and defaults to static B (ozone only) write_ens_sprd .false. writing global ensemble spread in byte addressable format for plotting with grads enspreproc .false. flag to read(.true.) pre-processed ensemble data already i_en_perts_io 0 flag to read in ensemble perturbations in ensemble grid. \ This is to speed up RAP/HRRR hybrid runs because the same ensemble perturbations are used in 6 cycles \ =0: No ensemble perturbations IO (default) \ =2: skip get_gefs_for_regional and read in ensemble perturbations from saved files. l_ens_in_diff_time .false. if use ensembles that are available at different time from analysis time. \ =false: only ensembles available at analysis time can be used for hybrid. (default) \ =true: ensembles available time can be different from analysis time in hybrid analysis ensemble_path ===================== ============= =========================================================================================================== ================================= ============= ============================================================================== Variable name Default value Description ================================= ============= ============================================================================== **rapidrefresh_cldsurf** **Options for cloud analysis and surface enhancement for RR application** dfi_radar_latent_heat_time_period 30.0 DFI forward integration window in minutes metar_impact_radius 10.0 metar cloud observation impact radius in grid number metar_impact_radius_lowCloud 4.0 impact radius for METAR cloud observation that indicate low cloud base l_gsd_terrain_match_surfTobs .false. if .true., GSD terrain match for surface temperature observation l_sfcobserror_ramp_t .false. namelist logical for adjusting surface temperature observation error l_sfcobserror_ramp_q .false. namelist logical for adjusting surface moisture observation error l_PBL_pseudo_SurfobsT .false. if .true. produce pseudo-obs in PBL layer based on surface obs T l_PBL_pseudo_SurfobsQ .false. if .true. produce pseudo-obs in PBL layer based on surface obs Q l_PBL_pseudo_SurfobsUV .false. if .true. produce pseudo-obs in PBL layer based on surface obs UV pblH_ration 0.75 percent of the PBL height within which to add pseudo-obs pps_press_incr 30hPa pressure increase for each additional pseudo-obs on top of previous level l_gsd_limit_ocean_q .false. if .true. do GSD limitation of Q over ocean l_pw_hgt_adjust .false. if .true. do GSD PW adjustment for model vs. obs station height l_limit_pw_innov .false. if .true. do GSD limitation of PW obs max_innov_pct 0.1 sets limit of PW ob to a percent of the background value (0-1) l_cleanSnow_WarmTs .false. if .true. do GSD limitation of using retrieved snow over warn area \ (Ts > r_cleanSnow_WarmTs_threshold) l_conserve_thetaV .false. if .true. conserve thetaV during moisture adjustment in cloud analysis r_cleanSnow_WarmTs_threshold 8.0 threshold for using retrieved snow over warn area i_conserve_thetaV_iternum 3 iteration number for conserving thetaV during moisture adjustment l_gsd_soilTQ_nudge .false. if .true. do GSD GOES cloud building l_cld_bld .false. if .true. do GSD soil T and Q nudging based on the lowest t analysis increment cld_bld_hgt 1200m sets limit below which GOES cloud building occurs build_cloud_frac_p 0.95 sets the threshold for building clouds from satellite clear_cloud_frac_p 0.1 sets the threshold for clearing clouds from satellite nesdis_npts_rad 1 NESDIS cloud product impact radiu (grid points) iclean_hydro_withRef 1 if =1, then clean hydrometeors if the grid point has no echo and maxref=0 iclean_hydro_withRef_allcol 0 if =1, then clean whole column hydrometeors if the observed max ref =0 and satellite cloud shows clean ================================= ============= ============================================================================== ======================== ============= ========================================================================== Variable name Default value Description ======================== ============= ========================================================================== **rapidrefresh_cldsurf** **Options for cloud analysis and surface enhancement for RR application** l_use_2mq4b 0 background used for calculate surface moisture observation innovation \ =0 Use Q from the 1st model level. (default) \ =1 use 2m Q as part of background i_use_2mt4b 0 background used for calculate surface temperature observation innovation \ =0 Use T from the 1st model level. (default) \ =1 use 2m T as part of background i_gsdcldanal_type 0 options for how GSD cloud analysis should be conducted \ =0. no cloud analysis (default) \ =1. cloud analysis after var analysis \ =5. skip cloud analysis and NETCDF file update i_gsdsfc_uselist 0 options for how to use surface observation use or rejection list \ =0 . EMC method (default) \ =1 . GSD method i_lightpcp 0 options for how to deal with light precipitation \ =0 . don’t add light precipitation (default) \ =1 . add light precipitation in warm section i_sfct_gross 0 if use extended threshold for surface T gross check \ =0 use threshold from convinfo (default) \ =1 for cold surface, threshold for gross check is enlarged to bring more \ large negative innovation into analysis. ======================== ============= ========================================================================== =============== ================== ============================================================================= Variable name Default value Description =============== ================== ============================================================================= **CHEM** **Chemistry data assimilation** berror_chem .false. if berror file is supplied for chemistry oneobtest_chem .false. single observation test for chemistry maginnov_chem 30.0 if oneobtest_chem=T magnitude of innovation for chemistry magoberr_chem 2.0 if oneobtest_chem=T magnitude of observation error for chemistry oneob_type_chem pm2_5 if oneobtest_chem=T type of chemical observation oblat_chem 45.0 if oneobtest_chem=T latitude of the observation oblon_chem 270.0 if oneobtest_chem=T longitude of the observation obpres_chem 1000.0 if oneobtest_chem=T pressure of the observation diag_incr .false. if user wishes to output to a binary file increment elev_tolerance 500.0 for surface chemical observation sometimes elevation (elev_obs) of the \ measurement is available (sometimes not). tunable_error 0.5 tuning parameter to specify representativeness error for in-situ observations in_fname cmaq_input.bin name of background file for cmaq out_fname cmaq_output.bin name analysis file for cmaq incr_fname chem_increment.bin if diag_incr=T name of the binary dump for pm2_5 laeroana_gocart .false. when true, do chem analysis with wrfchem and modis l_aoderr_table aod_qa_limit luse_deepblue aero_ratios wrf_pm2_5 =============== ================== =============================================================================