Skip to content

Commit

Permalink
Update namelist descriptions in registry.var and var/README.namelist …
Browse files Browse the repository at this point in the history 
…(#228)

TYPE: text only

KEYWORDS: WRFDA, README.namelist

SOURCE: internal

DESCRIPTION OF CHANGES:

1. in registry.var, add some descriptions and change a few namelist orders
    in the new radar_da record.
2. update var/README.namelist

LIST OF MODIFIED FILES:
M Registry/registry.var
M var/README.namelist

TESTS CONDUCTED:
none
  • Loading branch information
@jamiebresch
jamiebresch committed Mar 31, 2017
1 parent db53d72 commit ce975f6
Show file tree
Hide file tree
Showing 2 changed files with 65 additions and 20 deletions.
20 changes: 10 additions & 10 deletions Registry/registry.var
View file
Original file line number Diff line number Diff line change
Expand Up @@ -256,8 +256,8 @@ rconfig logical read_lanczos namelist,wrfvar6 1 .false. - "rea
rconfig logical write_lanczos namelist,wrfvar6 1 .false. - "write_lanczos" "" ""
rconfig logical orthonorm_gradient namelist,wrfvar6 1 .false. - "orthonorm_gradient" "" ""
rconfig integer cv_options namelist,wrfvar7 1 5 - "cv_options" "" ""
rconfig integer cloud_cv_options namelist,wrfvar7 1 0 - "cloud_cv_options" "" ""
rconfig logical use_cv_w namelist,wrfvar7 1 .false. - "use_cv_w" "if activate w control variable" ""
rconfig integer cloud_cv_options namelist,wrfvar7 1 0 - "cloud_cv_options" "0: off, 1: qt, 3: specified qc,qr,qi,qs,qg BE" ""
rconfig logical use_cv_w namelist,wrfvar7 1 .false. - "use_cv_w" "if activate w control variable when cloud_cv_options=3" ""
rconfig real as1 namelist,wrfvar7 3*max_outer_iterations -1.0 - "as1" "" ""
rconfig real as2 namelist,wrfvar7 3*max_outer_iterations -1.0 - "as2" "" ""
rconfig real as3 namelist,wrfvar7 3*max_outer_iterations -1.0 - "as3" "" ""
Expand Down Expand Up @@ -428,8 +428,8 @@ rconfig real pseudo_z namelist,wrfvar15 1 1.0 - "ps
rconfig real pseudo_val namelist,wrfvar15 1 1.0 - "pseudo_val" "" ""
rconfig real pseudo_err namelist,wrfvar15 1 1.0 - "pseudo_err" "" ""
rconfig real pseudo_elv namelist,wrfvar15 1 -999.99 - "pseudo_elv" "pseudo ob elevation, used by pseudo tpw/ztd" "m"
rconfig integer ep_para_read namelist,wrfvar16 1 0 - "ep_para_read" "" ""
rconfig integer rden_bin namelist,wrfvar16 1 1 - "rden_bin" "" ""
rconfig integer ep_para_read namelist,wrfvar16 1 0 - "ep_para_read" "how the ensemble perturbations are read" "0: serial read, 1: parallel read"
rconfig integer rden_bin namelist,wrfvar16 1 1 - "rden_bin" "number of bins/batches to read ensemble perturbations in parallel" ""
rconfig integer alphacv_method namelist,wrfvar16 1 2 - "alphacv_method" "" ""
rconfig integer ensdim_alpha namelist,wrfvar16 1 0 - "ensdim_alpha" "" ""
rconfig integer alpha_truncation namelist,wrfvar16 1 0 - "alpha_truncation" "" ""
Expand All @@ -440,7 +440,7 @@ rconfig logical alpha_vertloc namelist,wrfvar16 1 .false. - "al
rconfig logical alpha_hydrometeors namelist,wrfvar16 1 .false. - "alpha_hydrometeors" "" ""
rconfig logical hybrid_dual_res namelist,wrfvar16 1 .false. - "hybrid_dual_res" "" ""
rconfig integer dual_res_upscale_opt namelist,wrfvar16 1 3 - "dual_res_upscale_opt" "" ""
rconfig logical use_4denvar namelist,wrfvar16 1 .false. - "4D-En-Var" "" ""
rconfig logical use_4denvar namelist,wrfvar16 1 .false. - "4D-En-Var" "switch for activating 4D-Ensemble-Var" ""
rconfig character analysis_type namelist,wrfvar17 1 "3D-VAR" - "analysis_type" "" ""
rconfig integer sensitivity_option namelist,wrfvar17 1 -1 - "sensitivity_option" "" ""
rconfig logical adj_sens namelist,wrfvar17 1 .false. - "adj_sens" "" ""
Expand All @@ -450,15 +450,15 @@ rconfig character documentation_url namelist,wrfvar20 1 "http://www.mm
rconfig character time_window_min namelist,wrfvar21 1 "2002-08-02_21:00:00.0000" - "time_window_min" "" ""
rconfig character time_window_max namelist,wrfvar22 1 "2002-08-03_03:00:00.0000" - "time_window_max" "" ""
rconfig integer radar_non_precip_opt namelist,radar_da 1 0 - "radar_non_precip_opt" "" "0: off, 1: KNU scheme"
rconfig integer cloudbase_calc_opt namelist,radar_da 1 1 - "cloudbase_calc_opt" "" "1: KNU scheme, 2: NCAR scheme"
rconfig real radar_non_precip_rf namelist,radar_da 1 -999.99 - "radar_non_precip_rf" "rf value used to indicate non-precip ob" "dBZ"
rconfig real radar_non_precip_rh_w namelist,radar_da 1 95 - "radar_non_precip_rh_w" "RH wrt water for non_precip rqv" "%"
rconfig real radar_non_precip_rh_i namelist,radar_da 1 85 - "radar_non_precip_rh_i" "RH wrt ice for non_precip rqv" "%"
rconfig real radar_non_precip_rh_w namelist,radar_da 1 95.0 - "radar_non_precip_rh_w" "RH wrt water for non_precip rqv" "%"
rconfig real radar_non_precip_rh_i namelist,radar_da 1 85.0 - "radar_non_precip_rh_i" "RH wrt ice for non_precip rqv" "%"
rconfig integer cloudbase_calc_opt namelist,radar_da 1 1 - "cloudbase_calc_opt" "" "1: KNU scheme, 2: NCAR scheme"
rconfig real radar_saturated_rf namelist,radar_da 1 25.0 - "radar_saturated_rf" "rf value used to indicate saturated rqv" "dBZ"
rconfig real radar_rqv_thresh1 namelist,radar_da 1 40.0 - "radar_rqv_thresh1" "rf value used to scale down rqv" "dBZ"
rconfig real radar_rqv_thresh2 namelist,radar_da 1 50.0 - "radar_rqv_thresh2" "rf value used to scale down rqv" "dBZ"
rconfig real radar_rqv_rh1 namelist,radar_da 1 85 - "radar_rqv_rh1" "RH for (radar_saturated_rf < rf < radar_rqv_thresh1)" "%"
rconfig real radar_rqv_rh2 namelist,radar_da 1 95 - "radar_rqv_rh2" "RH for (radar_rqv_thresh1 < rf < radar_rqv_thresh2)" "%"
rconfig real radar_rqv_rh1 namelist,radar_da 1 85.0 - "radar_rqv_rh1" "RH for (radar_saturated_rf < rf < radar_rqv_thresh1)" "%"
rconfig real radar_rqv_rh2 namelist,radar_da 1 95.0 - "radar_rqv_rh2" "RH for (radar_rqv_thresh1 < rf < radar_rqv_thresh2)" "%"
rconfig real radar_rqv_h_lbound namelist,radar_da 1 -999.0 - "radar_rqv_h_lbound" "height lower bound for rqv" "m"
rconfig real radar_rqv_h_ubound namelist,radar_da 1 -999.0 - "radar_rqv_h_ubound" "height upper bound for rqv" "m"
rconfig logical jcdfi_use namelist,perturbation 1 .false. - "jcdfi_use" "JcDFI on/off" ""
Expand Down
65 changes: 55 additions & 10 deletions var/README.namelist
View file
Original file line number Diff line number Diff line change
Expand Up @@ -86,9 +86,8 @@ Description of WRFDA namelist variables
use_radarobs = .false. ; .true.: Assimilate radar data; see User's Guide for more information
use_radar_rv = .false. ; Assimilate radar velocity observations
use_radar_rf = .false. ; Assimilate radar reflectivity using original reflectivity operator (total mixing ratio)
use_radar_rqv = .false. ; Assimilate radar reflectivity using estimated humidity from radar reflectivity
use_radar_rhv = .false. ; Assimilate radar reflectivity using rainwater and ice mixing ratios
use_3dvar_phy = .true. ; Partition hydrometeors via the moist explicit scheme (warm rain process)
use_radar_rhv = .false. ; Assimilate retrieved hydrometeors (qr, qs, qg) from radar reflectivity
use_radar_rqv = .false. ; Assimilate estimated humidity (qv) from radar reflectivity
use_rainobs = .false. ; .true.: Assimilate precipitation data; see User's Guide for more information
thin_rainobs = .true. ; perform thinning on preciptation data
use_airsretobs = .true.
Expand Down Expand Up @@ -162,7 +161,7 @@ Description of WRFDA namelist variables
&wrfvar7
cv_options = 5 ; 3: NCEP Background Error model
; 5: NCAR Background Error model (default)
; 6: Multivariate background error statistics
; 6: Moisture-multivariate background error statistics
; 7: New NCAR Background Error model (CV7)
; For cv_options= 5, 6, or 7, be.dat must be generated by GEN_BE. See User's Guide for details.
;
Expand Down Expand Up @@ -202,6 +201,20 @@ Description of WRFDA namelist variables
len_scaling4 (max_outer_iterations) = 1.0 ; tuning factor of scale-length for control variable 4.
len_scaling5 (max_outer_iterations) = 1.0 ; tuning factor of scale-length for control variable 5.
je_factor = 1.0 ; ensemble covariance weighting factor
cloud_cv_options = 0 ; 0: no cloud control variables
; 1: Total water (qt) is the moist/cloud control variable.
; Cloud water and rainwater are combined with water vapor as total water.
; A warm-rain scheme is used to partition the total water into water vapor,
; cloud water and rainwater during minimization.
; 2: Moist control variable is pseudo relative humidity.
; The additional 5 cloud control variables are cloud water, rainwater, cloud ice,
; snow and graupel. Both horizontal and vertical correlations are considered.
; This option requires the background error statistics information of the
; aforementioned 5 cloud variables to be included in the be.dat file.
; 3: Moist control variable is pseudo relative humidity. The additional 5 cloud
; control variables are cloud water, rainwater, cloud ice, snow and graupel.
; Only horizontal correlation is considered. Length scales and variances are
; hard-coded in var/da/da_setup_structures/da_setup_be_regional.inc.
/
&wrfvar8 ; not used
/
Expand Down Expand Up @@ -250,24 +263,37 @@ Description of WRFDA namelist variables
sfc_assi_options = 1 ; 1: surface observations will be assimilated based on the
; lowest model level first guess. Observations are not used
; when the height difference of the elevation of the observing
; site and the lowest model level height is larger than 100m.
; site and the lowest model level height is larger than the
; height (in meters) specified in namelist max_stheight_diff.
; 2: surface observations will be assimilated based on surface
; similarity theory in PBL. Innovations are computed based
; on 10-m wind, 2-m temperature and 2-m moisture.
max_stheight_diff = 100.0 ; threshold in meters. See sfc_assi_options=1.
obs_err_inflate = .false. ; .true.: inflate obs error by exp(|Zdiff|/stn_ht_diff_scale) for
; SYNOP t, p, q with sfc_assi_options=1 only.
; Zdiff is the height difference between model and station that is
; calculated in WRFDA.
stn_ht_diff_scale = 200.0 ; height in meters used for calculating obs error inflation factors
; when obs_err_inflate = .true.
psfc_from_slp = .false. ; .false.: use reported Psfc.
; This is the recommended (default) option.
; .true.: when sfc_assi_options=1, re-calculates Psfc from SLP
; when the observation elevation is below the lowest
; model level height. This is the behavior before V3.8.
calculate_cg_cost_fn = .false. ; .true.: Compute and write out cost function and gradient of
; each iteration into files called cost_fn and grad_fn
; for diagnostic purpose.
; .false.: Only the initial and final cost functions are computed
calculate_cg_cost_fn = .false. ; .true.: Compute and write out cost function of each iteration
; into a file cost_fn for diagnostic purpose.
; .false.: Only the initial and final cost functions are computed
; and output in cost_fn.
; The cost function values in rsl.out.0000 are estimates.
; The actual cost function of each iteration is not needed for
; minimization. For production runs, set calculate_cg_cost_fn=
; .false. to reduce run time
; .false. to reduce run time.
write_detail_grad_fn = .false. ; .true.: calculate and write out detailed gradient components of
; each iteration into a file grad_fn for diagnostic purpose.
; This will increase runtime, especially for 4DVAR.
; The summed gradient info is available in rsl.out.0000.
; Keeping write_detail_grad_fn=false (default) is the best practice.

seed_array1 = 1 ; First integer for seeding the random function when put_rand_seed=.true.
seed_array2 = 1 ; Second integer for seeding the random function when put_rand_seed=.true.
; It is not necessary to change both seeds to get different perturbations
Expand Down Expand Up @@ -440,6 +466,25 @@ Description of WRFDA namelist variables
; Also used for ob_format=2 to check if the obs
; are within the specified time window.
/
&radar_da
; tuning parameters for use_radar_rqv=.true. only
; except that radar_non_precip_opt applies to both use_radar_rhv=.true. and use_radar_rqv=.true.
radar_non_precip_opt = 0 ; 0: do not consider non-precipitating radar echo
; 1: KNU scheme for assimilating non-precipitating radar echo
radar_non_precip_rf = -999.99 ; reflectivity (dBZ) value used to indicate non-precip ob
radar_non_precip_rh_w = 95.0 ; RH (%) wrt water for non_precip rqv
radar_non_precip_rh_i = 85.0 ; RH (%) wrt ice for non_precip rqv
cloudbase_calc_opt = 1 ; how model LCL is calculated
; 1: KNU scheme
; 2: NCAR scheme
radar_saturated_rf = 25.0 ; reflectivity (dBZ) value used to indicate saturated rqv
radar_rqv_thresh1 = 40.0 ; reflectivity (dBZ) threshold for scaling down rqv
radar_rqv_thresh2 = 50.0 ; reflectivity (dBZ) threshold for scaling down rqv
radar_rqv_rh1 = 95.0 ; RH (%) for (radar_saturated_rf < rf < radar_rqv_thresh1)
radar_rqv_rh2 = 85.0 ; RH (%) for (radar_rqv_thresh1 < rf < radar_rqv_thresh2)
radar_rqv_h_lbound = -999.0 ; height (m) lower bound for rqv
radar_rqv_h_ubound = -999.0 ; height (m) upper bound for rqv
/
&perturbation ; settings related to the 4D-Var penalty term option, which controls the
; high-frequency gravity waves using a digital filter.
jcdfi_use = .true. ; .true.: include JcDF term in the cost function
Expand Down

0 comments on commit ce975f6

Please sign in to comment.