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 |