From 2fd43d134f9d8479eb37c3800f6bcb715ceb7c3f Mon Sep 17 00:00:00 2001 From: Kate Friedman Date: Mon, 8 May 2023 16:53:21 -0400 Subject: [PATCH] Redo v16.3 GSI script updates to scripts/exglobal_atmos_analysis.sh (#1535) Updates to scripts/exglobal_atmos_analysis.sh are being put in again after a CRTM issue was resolved. These updates are needed to resolve an issue with the GSI assimilating too low of a number of observations. Also update gsi_ver to 20230112 in fix.ver Refs #1494 #1550 --- parm/config/config.anal | 2 +- parm/config/config.base.emc.dyn | 4 ++-- parm/config/config.base.nco.static | 4 ++-- scripts/exglobal_atmos_analysis.sh | 12 +++++------- versions/fix.ver | 2 +- 5 files changed, 11 insertions(+), 13 deletions(-) diff --git a/parm/config/config.anal b/parm/config/config.anal index f68ec28efe..e3a17f9c6a 100644 --- a/parm/config/config.anal +++ b/parm/config/config.anal @@ -29,7 +29,7 @@ fi # Set parameters specific to L127 if [[ ${LEVS} = "128" ]]; then export GRIDOPTS="nlayers(63)=1,nlayers(64)=1," - export SETUP="gpstop=55,nsig_ext=56,${SETUP:-}" + export SETUP="gpstop=55,nsig_ext=45,${SETUP:-}" fi # Set namelist option for LETKF diff --git a/parm/config/config.base.emc.dyn b/parm/config/config.base.emc.dyn index e5e29d1440..95a071dc9b 100644 --- a/parm/config/config.base.emc.dyn +++ b/parm/config/config.base.emc.dyn @@ -354,13 +354,13 @@ export MAKE_NSSTBUFR="@MAKE_NSSTBUFR@" export MAKE_ACFTBUFR="@MAKE_ACFTBUFR@" # Analysis increments to zero in CALCINCEXEC -export INCREMENTS_TO_ZERO="'liq_wat_inc','icmr_inc'" +export INCREMENTS_TO_ZERO="'liq_wat_inc','icmr_inc','rwmr_inc','snmr_inc','grle_inc'" # Write analysis files for early cycle EnKF export DO_CALC_INCREMENT_ENKF_GFS="YES" # Stratospheric increments to zero -export INCVARS_ZERO_STRAT="'sphum_inc','liq_wat_inc','icmr_inc'" +export INCVARS_ZERO_STRAT="'sphum_inc','liq_wat_inc','icmr_inc','rwmr_inc','snmr_inc','grle_inc'" export INCVARS_EFOLD="5" # Swith to generate netcdf or binary diagnostic files. If not specified, diff --git a/parm/config/config.base.nco.static b/parm/config/config.base.nco.static index e3702852f4..a2620e5ab1 100644 --- a/parm/config/config.base.nco.static +++ b/parm/config/config.base.nco.static @@ -228,10 +228,10 @@ if [ $DONST = "YES" ]; then export FNTSFA=" "; fi export nst_anl=.true. # Analysis increments to zero in CALCINCEXEC -export INCREMENTS_TO_ZERO="'liq_wat_inc','icmr_inc'" +export INCREMENTS_TO_ZERO="'liq_wat_inc','icmr_inc','rwmr_inc','snmr_inc','grle_inc'" # Stratospheric increments to zero -export INCVARS_ZERO_STRAT="'sphum_inc','liq_wat_inc','icmr_inc'" +export INCVARS_ZERO_STRAT="'sphum_inc','liq_wat_inc','icmr_inc','rwmr_inc','snmr_inc','grle_inc'" export INCVARS_EFOLD="5" # Swith to generate netcdf or binary diagnostic files. If not specified, diff --git a/scripts/exglobal_atmos_analysis.sh b/scripts/exglobal_atmos_analysis.sh index ccedbdabda..6ab24523c8 100755 --- a/scripts/exglobal_atmos_analysis.sh +++ b/scripts/exglobal_atmos_analysis.sh @@ -55,8 +55,8 @@ lupp=${lupp:-".true."} cnvw_option=${cnvw_option:-".false."} # Observation usage options -cao_check=${cao_check:-".false."} -ta2tb=${ta2tb:-".false."} +cao_check=${cao_check:-".true."} +ta2tb=${ta2tb:-".true."} # Diagnostic files options lobsdiag_forenkf=${lobsdiag_forenkf:-".false."} @@ -423,9 +423,7 @@ ${NLN} ${RTMFIX}/NPOESS.VISsnow.EmisCoeff.bin ./crtm_coeffs/NPOESS.VISsnow.Emis ${NLN} ${RTMFIX}/NPOESS.VISwater.EmisCoeff.bin ./crtm_coeffs/NPOESS.VISwater.EmisCoeff.bin ${NLN} ${RTMFIX}/FASTEM6.MWwater.EmisCoeff.bin ./crtm_coeffs/FASTEM6.MWwater.EmisCoeff.bin ${NLN} ${RTMFIX}/AerosolCoeff.bin ./crtm_coeffs/AerosolCoeff.bin -${NLN} ${RTMFIX}/CloudCoeff.bin ./crtm_coeffs/CloudCoeff.bin -#$NLN $RTMFIX/CloudCoeff.GFDLFV3.-109z-1.bin ./crtm_coeffs/CloudCoeff.bin - +${NLN} ${RTMFIX}/CloudCoeff.GFDLFV3.-109z-1.bin ./crtm_coeffs/CloudCoeff.bin ############################################################## # Observational data @@ -702,7 +700,7 @@ cat > gsiparm.anl << EOF iguess=-1, tzr_qc=${TZR_QC}, oneobtest=.false.,retrieval=.false.,l_foto=.false., - use_pbl=.false.,use_compress=.true.,nsig_ext=12,gpstop=50.,commgpstop=45.,commgpserrinf=1.0, + use_pbl=.false.,use_compress=.true.,nsig_ext=45,gpstop=50.,commgpstop=45.,commgpserrinf=1.0, use_gfs_nemsio=.false.,use_gfs_ncio=.true.,sfcnst_comb=.true., use_readin_anl_sfcmask=${USE_READIN_ANL_SFCMASK}, lrun_subdirs=${lrun_subdirs}, @@ -755,7 +753,7 @@ cat > gsiparm.anl << EOF ${OBSQC} / &OBS_INPUT - dmesh(1)=145.0,dmesh(2)=150.0,dmesh(3)=100.0,dmesh(4)=25.0,time_window_max=3.0, + dmesh(1)=145.0,dmesh(2)=150.0,dmesh(3)=100.0,dmesh(4)=50.0,time_window_max=3.0, ${OBSINPUT} / OBS_INPUT:: diff --git a/versions/fix.ver b/versions/fix.ver index 0046ed5eb8..02b0953a94 100644 --- a/versions/fix.ver +++ b/versions/fix.ver @@ -12,7 +12,7 @@ export gdas_fv3jedi_ver=20220805 export gdas_gsibec_ver=20221031 export gldas_ver=20220920 export glwu_ver=20220805 -export gsi_ver=20221128 +export gsi_ver=20230112 export lut_ver=20220805 export mom6_ver=20220805 export orog_ver=20220805