Make the GSI and ENKF analysis horizontal and vertical localization s…

…cales as

HAFS application/workflow config options. Also, the following default values
were given as follows: gsi.s_ens_h=150, gsi.s_ens_v=-0.5, enkf.corrlength=500,
enkf.lnsigcutoff=1.3 based on the HAFS DA community discussion.
Note:
*These default localization scales are currently set based on the HWRF D03 DA
settings and the HAFS DA community discussions. Further tests will be needed to
determine and fine-tune these localization scale settings for HAFS DA system in
the future.
*Beside, we plan to also explore using the hybens_info file to specify the
vertical profile of horizontal and vertical localization length scales for HAFS
EnVar and EnKF analysis.

parm/analysis/gsi/enkf.nml.tmp

@@ -2,10 +2,10 @@
    datestring="_datestring_",datapath="./",
    analpertwtnh=0.85,analpertwtsh=0.85,analpertwttr=0.85,
    covinflatemax=1.e2,covinflatemin=1,pseudo_rh=.true.,iassim_order=0,
-   corrlengthnh=400,corrlengthsh=400,corrlengthtr=400,
-   lnsigcutoffnh=0.5,lnsigcutoffsh=0.5,lnsigcutofftr=0.5,
-   lnsigcutoffpsnh=0.5,lnsigcutoffpssh=0.5,lnsigcutoffpstr=0.5,
-   lnsigcutoffsatnh=0.5,lnsigcutoffsatsh=0.5,lnsigcutoffsattr=0.5,
+   corrlengthnh=_corrlength_,corrlengthsh=_corrlength_,corrlengthtr=_corrlength_,
+   lnsigcutoffnh=_lnsigcutoff_,lnsigcutoffsh=_lnsigcutoff_,lnsigcutofftr=_lnsigcutoff_,
+   lnsigcutoffpsnh=_lnsigcutoff_,lnsigcutoffpssh=_lnsigcutoff_,lnsigcutoffpstr=_lnsigcutoff_,
+   lnsigcutoffsatnh=_lnsigcutoff_,lnsigcutoffsatsh=_lnsigcutoff_,lnsigcutoffsattr=_lnsigcutoff_,
    obtimelnh=1.e30,obtimelsh=1.e30,obtimeltr=1.e30,
    saterrfact=1.0,numiter=1,
    sprd_tol=1.e30,paoverpb_thresh=0.98,

parm/analysis/gsi/gsiparm.anl.tmp

@@ -229,8 +229,8 @@ SUPEROB_RADAR::
   n_ens=_N_ENS_,
   uv_hyb_ens=T,
   beta_s0=_BETA_S0_,
-  s_ens_h=300,
-  s_ens_v=-0.5,
+  s_ens_h=_S_ENS_H_,
+  s_ens_v=_S_ENS_V_,
   readin_localization=F,
   generate_ens=F,
   regional_ensemble_option=_REGIONAL_ENSEMBLE_OPTION_,

parm/hafs.conf

@@ -179,8 +179,12 @@ jend_nest_ens={grid/jend_nest}
 [gsi]
 use_bufr_nr=no            ;; use non-restricted version of observational bufr files
 grid_ratio_fv3_regional=1 ;; ratio of analysis grid to fv3 model grid in fv3 grid units
+s_ens_h=150               ;; homogeneous isotropic horizontal ensemble localization scale (km)
+s_ens_v=-0.5              ;; vertical localization scale
 
 [enkf]
+corrlength=500            ;; length for horizontal localization in km
+lnsigcutoff=1.3           ;; length for vertical localization
 
 [merge]
 

parm/hafs_holdvars.conf

@@ -91,6 +91,10 @@ npz_ens={forecast_ens/npz_ens}
 
 use_bufr_nr={gsi/use_bufr_nr}
 grid_ratio_fv3_regional={gsi/grid_ratio_fv3_regional}
+s_ens_h={gsi/s_ens_h}
+s_ens_v={gsi/s_ens_v}
+corrlength={enkf/corrlength}
+lnsigcutoff={enkf/lnsigcutoff}
 
 ccpp_suite_regional={forecast/ccpp_suite_regional}
 ccpp_suite_glob={forecast/ccpp_suite_glob}

parm/hafs_holdvars.txt

@@ -82,6 +82,10 @@ export vcoord_file_target_grid={vcoord_file_target_grid}
 
 export use_bufr_nr={use_bufr_nr}
 export grid_ratio_fv3_regional={grid_ratio_fv3_regional}
+export s_ens_h={s_ens_h}
+export s_ens_v={s_ens_v}
+export corrlength={corrlength}
+export lnsigcutoff={lnsigcutoff}
 
 export ccpp_suite_regional={ccpp_suite_regional}
 export ccpp_suite_glob={ccpp_suite_glob}

scripts/exhafs_analysis.sh

@@ -10,6 +10,8 @@ export DONST=${DONST:-"NO"}
 export LEVS=${LEVS:-65}
 export use_bufr_nr=${use_bufr_nr:-no}
 export grid_ratio_fv3_regional=${grid_ratio_fv3_regional:-1}
+export s_ens_h=${s_ens_h:-150}
+export s_ens_v=${s_ens_v:--0.5}
 export out_prefix=${out_prefix:-$(echo "${STORM}${STORMID}.${YMDH}" | tr '[A-Z]' '[a-z]')}
 
 export RUN_GSI_VR=${RUN_GSI_VR:-NO}
@@ -514,6 +516,8 @@ sed -e "s/_MITER_/${MITER:-2}/g" \
     -e "s/_REDUCE_DIAG_/${REDUCE_DIAG:-.false.}/g" \
     -e "s/_L_HYB_ENS_/${L_HYB_ENS:-.false.}/g" \
     -e "s/_N_ENS_/${N_ENS:-80}/g" \
+    -e "s/_S_ENS_H_/${s_ens_h:-150}/g" \
+    -e "s/_S_ENS_V_/${s_ens_v:--0.5}/g" \
     -e "s/_BETA_S0_/${BETA_S0:-0.2}/g" \
     -e "s/_GRID_RATIO_ENS_/${GRID_RATIO_ENS:-1}/g" \
     -e "s/_REGIONAL_ENSEMBLE_OPTION_/${REGIONAL_ENSEMBLE_OPTION:-1}/g" \

scripts/exhafs_analysis_vr.sh

@@ -27,6 +27,8 @@ export COMINhafs=${COMgfs:-/gpfs/dell1/nco/ops/com/gfs/para}
 export DONST=${DONST:-"NO"}
 export use_bufr_nr=${use_bufr_nr:-no}
 export grid_ratio_fv3_regional=${grid_ratio_fv3_regional:-1}
+export s_ens_h=${s_ens_h:-150}
+export s_ens_v=${s_ens_v:--0.5}
 export out_prefix=${out_prefix:-$(echo "${STORM}${STORMID}.${YMDH}" | tr '[A-Z]' '[a-z]')}
 
 export RUN_GSI_VR=${RUN_GSI_VR:-NO}
@@ -331,6 +333,8 @@ sed -e "s/_MITER_/${MITER:-2}/g" \
     -e "s/_REDUCE_DIAG_/${REDUCE_DIAG:-.false.}/g" \
     -e "s/_L_HYB_ENS_/${L_HYB_ENS:-.false.}/g" \
     -e "s/_N_ENS_/${N_ENS:-80}/g" \
+    -e "s/_S_ENS_H_/${s_ens_h:-150}/g" \
+    -e "s/_S_ENS_V_/${s_ens_v:--0.5}/g" \
     -e "s/_BETA_S0_/${BETA_S0:-0.2}/g" \
     -e "s/_GRID_RATIO_ENS_/${GRID_RATIO_ENS:-1}/g" \
     -e "s/_REGIONAL_ENSEMBLE_OPTION_/${REGIONAL_ENSEMBLE_OPTION:-1}/g" \

scripts/exhafs_enkf.sh

@@ -50,6 +50,8 @@ fi
 export PARMgsi=${PARMgsi:-${PARMhafs}/analysis/gsi}
 export ldo_enscalc_option=${ldo_enscalc_option:-1}
 export nens=${ENS_SIZE:-40}
+export corrlength=${corrlength:-500}
+export lnsigcutoff=${lnsigcutoff:-1.3}
 
 # Diagnostic files options
 netcdf_diag=${netcdf_diag:-".true."}
@@ -233,6 +235,8 @@ ${NLN} ${COMgfs}/gdas.$PDYprior/${hhprior}/${atmos}gdas.t${hhprior}z.abias_pc
 ${NCP} ${PARMgsi}/enkf.nml.tmp ./
 
 sed -e "s/_datestring_/${CDATE}/g" \
+    -e "s/_corrlength_/${corrlength:-500}/g" \
+    -e "s/_lnsigcutoff_/${lnsigcutoff:-1.3}/g" \
     -e "s/_nlons_/$((${npx_ens:-$npx}-1))/g" \
     -e "s/_nlats_/$((${npy_ens:-$npy}-1))/g" \
     -e "s/_nlevs_/${npz_ens:-$npz}/g" \

scripts/exhafs_enkf_hx.sh

@@ -11,6 +11,8 @@ export COMINhafs=${COMgfs:-/gpfs/dell1/nco/ops/com/gfs/para}
 export DONST=${DONST:-"NO"}
 export use_bufr_nr=${use_bufr_nr:-no}
 export grid_ratio_fv3_regional=${grid_ratio_fv3_regional:-1}
+export s_ens_h=${s_ens_h:-150}
+export s_ens_v=${s_ens_v:--0.5}
 export out_prefix=${out_prefix:-$(echo "${STORM}${STORMID}.${YMDH}" | tr '[A-Z]' '[a-z]')}
 
 if [ ${ENSDA} = YES ]; then
@@ -483,6 +485,8 @@ sed -e "s/_MITER_/${MITER:-2}/g" \
     -e "s/_REDUCE_DIAG_/${REDUCE_DIAG:-.false.}/g" \
     -e "s/_L_HYB_ENS_/${L_HYB_ENS:-.false.}/g" \
     -e "s/_N_ENS_/${N_ENS:-80}/g" \
+    -e "s/_S_ENS_H_/${s_ens_h:-150}/g" \
+    -e "s/_S_ENS_V_/${s_ens_v:--0.5}/g" \
     -e "s/_BETA_S0_/${BETA_S0:-0.2}/g" \
     -e "s/_GRID_RATIO_ENS_/${GRID_RATIO_ENS:-1}/g" \
     -e "s/_REGIONAL_ENSEMBLE_OPTION_/${REGIONAL_ENSEMBLE_OPTION:-1}/g" \