From 73d8f20969311a70f9ec1f57dd61bd2982acba8d Mon Sep 17 00:00:00 2001 From: "catherine.thomas" Date: Wed, 13 Nov 2024 14:58:23 -0500 Subject: [PATCH 1/3] Enable tapering of atm ens perts at the model top This commit turns on the tapering of atmospheric ensemble perturbations at the model top in the EnKF to reduce the upper layer ensemble spread. This commit also includes a minor update of changing the IC tag date for the staged HR4 ICs. Resolves #3096 --- parm/config/gfs/config.eupd | 2 +- parm/config/gfs/config.stage_ic | 2 +- scripts/exgdas_enkf_update.sh | 3 ++- 3 files changed, 4 insertions(+), 3 deletions(-) diff --git a/parm/config/gfs/config.eupd b/parm/config/gfs/config.eupd index 2ff48240aed..7b0c8994c2a 100644 --- a/parm/config/gfs/config.eupd +++ b/parm/config/gfs/config.eupd @@ -28,7 +28,7 @@ export lnsigcutoff=2.75 # ignored if modelspace_vloc=.true. export lobsdiag_forenkf=".true." # use jacobian. must be .true. if modelspace_vloc=".true." # need to specify .true. setting since config.anal sets to .false. - +export taperanalperts=".true." export NAM_ENKF="smoothparm=35," echo "END: config.eupd" diff --git a/parm/config/gfs/config.stage_ic b/parm/config/gfs/config.stage_ic index 7aa0c25f321..67b38fd9702 100644 --- a/parm/config/gfs/config.stage_ic +++ b/parm/config/gfs/config.stage_ic @@ -15,7 +15,7 @@ export STAGE_IC_YAML_TMPL="${PARMgfs}/stage/master_gfs.yaml.j2" # Set ICSDIR (if not defined) if [[ -z "${ICSDIR}" ]] ; then - ic_ver="20240610" + ic_ver="20240829" if (( NMEM_ENS > 0 )) ; then ensic="${CASE_ENS}" diff --git a/scripts/exgdas_enkf_update.sh b/scripts/exgdas_enkf_update.sh index e924274d395..3880872dff1 100755 --- a/scripts/exgdas_enkf_update.sh +++ b/scripts/exgdas_enkf_update.sh @@ -77,6 +77,7 @@ write_spread_diag=${write_spread_diag:-".false."} cnvw_option=${cnvw_option:-".false."} netcdf_diag=${netcdf_diag:-".true."} modelspace_vloc=${modelspace_vloc:-".false."} # if true, 'vlocal_eig.dat' is needed +taperanalperts=${taperanalperts:-".false."} IAUFHRS_ENKF=${IAUFHRS_ENKF:-6} NMEM_ENS_MAX=${NMEM_ENS:-80} if [ "${RUN}" = "enkfgfs" ]; then @@ -287,7 +288,7 @@ cat > enkf.nml << EOFnml univaroz=.false.,adp_anglebc=.true.,angord=4,use_edges=.false.,emiss_bc=.true., letkf_flag=${letkf_flag},nobsl_max=${nobsl_max},denkf=${denkf},getkf=${getkf}., nhr_anal=${IAUFHRS_ENKF},nhr_state=${IAUFHRS_ENKF}, - lobsdiag_forenkf=$lobsdiag_forenkf, + lobsdiag_forenkf=$lobsdiag_forenkf,taperanalperts=$taperanalperts, write_spread_diag=$write_spread_diag, modelspace_vloc=$modelspace_vloc, use_correlated_oberrs=${use_correlated_oberrs}, From a245732c8d2b82e82efdb85f637ba0a231382af4 Mon Sep 17 00:00:00 2001 From: Rahul Mahajan Date: Wed, 13 Nov 2024 16:56:12 -0500 Subject: [PATCH 2/3] Update scripts/exgdas_enkf_update.sh --- scripts/exgdas_enkf_update.sh | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/scripts/exgdas_enkf_update.sh b/scripts/exgdas_enkf_update.sh index 3880872dff1..752cb07a6b6 100755 --- a/scripts/exgdas_enkf_update.sh +++ b/scripts/exgdas_enkf_update.sh @@ -288,7 +288,7 @@ cat > enkf.nml << EOFnml univaroz=.false.,adp_anglebc=.true.,angord=4,use_edges=.false.,emiss_bc=.true., letkf_flag=${letkf_flag},nobsl_max=${nobsl_max},denkf=${denkf},getkf=${getkf}., nhr_anal=${IAUFHRS_ENKF},nhr_state=${IAUFHRS_ENKF}, - lobsdiag_forenkf=$lobsdiag_forenkf,taperanalperts=$taperanalperts, + lobsdiag_forenkf=${lobsdiag_forenkf},taperanalperts=${taperanalperts}, write_spread_diag=$write_spread_diag, modelspace_vloc=$modelspace_vloc, use_correlated_oberrs=${use_correlated_oberrs}, From d47da43d940cb5af38f2f9003974151f44eb39ab Mon Sep 17 00:00:00 2001 From: CatherineThomas-NOAA <59020064+CatherineThomas-NOAA@users.noreply.github.com> Date: Wed, 13 Nov 2024 18:08:48 -0500 Subject: [PATCH 3/3] Revert config.stage_ic --- parm/config/gfs/config.stage_ic | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/parm/config/gfs/config.stage_ic b/parm/config/gfs/config.stage_ic index 67b38fd9702..7aa0c25f321 100644 --- a/parm/config/gfs/config.stage_ic +++ b/parm/config/gfs/config.stage_ic @@ -15,7 +15,7 @@ export STAGE_IC_YAML_TMPL="${PARMgfs}/stage/master_gfs.yaml.j2" # Set ICSDIR (if not defined) if [[ -z "${ICSDIR}" ]] ; then - ic_ver="20240829" + ic_ver="20240610" if (( NMEM_ENS > 0 )) ; then ensic="${CASE_ENS}"