diff --git a/scripts/exregional_nonvarcldanl.sh b/scripts/exregional_nonvarcldanl.sh
index cde59a0b1..4db1a1f9f 100755
--- a/scripts/exregional_nonvarcldanl.sh
+++ b/scripts/exregional_nonvarcldanl.sh
@@ -316,6 +316,9 @@ cat << EOF > gsiparm.anl
l_saturate_bkCloud=.true.,
l_qnr_from_qr=${l_qnr_from_qr},
n0_rain=100000000.0
+ i_T_Q_adjust=${i_T_Q_adjust},
+ l_rtma3d=${l_rtma3d},
+ i_precip_vertical_check=${i_precip_vertical_check},
/
EOF
diff --git a/ush/config.sh.RTMA_CONUS_3km b/ush/config.sh.RTMA_CONUS_3km
index bb095514b..dea4ba40a 100644
--- a/ush/config.sh.RTMA_CONUS_3km
+++ b/ush/config.sh.RTMA_CONUS_3km
@@ -3,6 +3,7 @@ source fix_rrfs_locations.sh #determine current platform
MACHINE=${PLATFORM}
ACCOUNT="nrtrr"
RESERVATION="rrfsdet"
+RESERVATION_POST="rrfsdet"
EXPT_BASEDIR="/home/rtrr/RTMA"
EXPT_SUBDIR="RTMA_CONUS"
@@ -14,6 +15,7 @@ if [[ -n $RESERVATION ]] ; then
PARTITION_GRAPHICS=kjet
PARTITION_ANALYSIS=kjet
PARTITION_WGRIB2=kjet
+ PARTITION_POST=kjet
QUEUE_DEFAULT="rt"
QUEUE_FCST="rt"
@@ -21,6 +23,9 @@ if [[ -n $RESERVATION ]] ; then
QUEUE_WGRIB2="rt"
QUEUE_GRAPHICS="rt"
QUEUE_HPSS="rt"
+ if [[ -n ${RESERVATION_POST} ]] ; then
+ QUEUE_POST="rt"
+ fi
fi
OBSPATH_NSSLMOSIAC=/public/data/radar/nssl/mrms/conus
@@ -85,9 +90,11 @@ readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_
beta1_inv=0.5 #beata_inv is 1-ensemble_wgt
ens_h=20 #110
ens_v=1 #3
-regional_ensemble_option=1 #1 for GDAS
+regional_ensemble_option=5 #1 for GDAS, 5 for RRFSDAS
grid_ratio_fv3=2.0 #fv3 resolution 3km, so analysis=3*2=6km
grid_ratio_ens=3 #analysis 6km, so ensemble=6*3=18km. GDAS ensemble is 20km
+##MEMO1: regional_ensemble_option=5 will be changed to 1 if there are not enough RRFS ensembles at the runtime.
+##MEMO2: grid_ratio_ens will be hardwired to 1 for RRFSDAS, so grid_ratio_ens here only affects GDAS.
i_en_perts_io=1 #0 or 1: original file 3: pre-processed ensembles
# &RAPIDREFRESH_CLDSURF
@@ -95,6 +102,9 @@ l_PBL_pseudo_SurfobsT=.false.
l_PBL_pseudo_SurfobsQ=.false
i_use_2mQ4B=0
i_use_2mT4B=0
+i_T_Q_adjust=1
+l_rtma3d=.true.
+i_precip_vertical_check=3
#-------------------------------------------------------------------------------------
#
BERROR_FN="berror.rtma" #under $FIX_GSI
@@ -116,8 +126,8 @@ EXTRN_MDL_LBCS_SEARCH_OFFSET_HRS="0"
BOUNDARY_LEN_HRS="1"
BOUNDARY_LONG_LEN_HRS="1"
-DATE_FIRST_CYCL="20220406"
-DATE_LAST_CYCL="20220626"
+DATE_FIRST_CYCL="20220801"
+DATE_LAST_CYCL="20221231"
CYCL_HRS=( "00" "12" )
CYCL_HRS_SPINSTART=("03" "15")
CYCL_HRS_PRODSTART=("99") # RTMA always warm-starts
@@ -144,6 +154,9 @@ for i in {0..23}; do FCST_LEN_HRS_CYCLES[$i]=0.0003; done
DA_CYCLE_INTERV="1" #need to be 1 to have hourly DA
RESTART_INTERVAL="1"
+USE_RRFSE_ENS="TRUE"
+CYCL_HRS_HYB_FV3LAM_ENS=("19" "20" "21" "22" "23" "00")
+
NCORES_RUN_ANAL=240
WTIME_RUN_ANAL="01:00:00"
HYBENSMEM_NMIN=66
@@ -221,11 +234,16 @@ OBSPATH="/lfs4/BMC/wrfruc/gge/rtma/wcoss_obs"
STMP="/lfs4/BMC/nrtrr/NCO_dirs/stmp" # Path to directory STMP that mostly contains input files.
PTMP="/lfs4/BMC/nrtrr/NCO_dirs/ptmp" # Path to directory STMP that mostly contains input files.
NWGES="/lfs4/BMC/nrtrr/NCO_dirs/nwges" # Path to directory NWGES that save boundary, cold initial, restart files
+if [[ ${regional_ensemble_option} == "5" ]]; then
+ RRFSE_NWGES="/lfs4/BMC/wrfruc/RRFSE/NCO_dirs/nwges" # Path to RRFSE directory NWGES that mostly contains ensemble restart files for GSI hybrid.
+ NUM_ENS_MEMBERS=30 # FV3LAM ensemble size for GSI hybrid analysis
+fi
# change settings if it is wcoss_dell_p3
if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
ACCOUNT="RTMA-T2O"
RESERVATION=""
+ RESERVATION_POST=""
EXPT_BASEDIR="/gpfs/dell2/emc/modeling/save/${USER}/RTMA-CONUS"
SERVICE_ACCOUNT=""
PARTITION_DEFAULT=""
@@ -238,6 +256,10 @@ if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
OBSPATH_NSSLMOSIAC="/gpfs/dell1/nco/ops/dcom/prod/ldmdata/obs/upperair/mrms/conus/MergedReflectivityQC"
AIRCRAFT_REJECT="/gpfs/dell2/emc/modeling/noscrub/Edward.Colon/FixData/obsuselist/amdar_reject_lists"
SFCOBS_USELIST="/gpfs/dell2/emc/modeling/noscrub/Edward.Colon/FixData/obsuselist/mesonet_uselists"
+
+ regional_ensemble_option=1
+ USE_RRFSE_ENS="FALSE"
+ CYCL_HRS_HYB_FV3LAM_ENS=("99")
WTIME_RUN_ANAL="00:30:00"
NNODES_RUN_ANAL="28"
@@ -262,7 +284,4 @@ if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
PTMP="/gpfs/dell2/ptmp/${USER}/ptmp"
NWGES="/gpfs/dell2/ptmp/${USER}/nwges"
fi
-# Wait longer for rrfsdas
-#if [[ "$MACHINE" == "jet" ]] && [[ $regional_ensemble_option -eq 5 ]] ; then
-# START_TIME_ANAL="5:00:00"
-#fi
+
diff --git a/ush/config.sh.RTMA_NA_3km b/ush/config.sh.RTMA_NA_3km
index 16617caa4..05e7c85c5 100644
--- a/ush/config.sh.RTMA_NA_3km
+++ b/ush/config.sh.RTMA_NA_3km
@@ -14,6 +14,7 @@ if [[ -n $RESERVATION ]] ; then
PARTITION_GRAPHICS=kjet
PARTITION_ANALYSIS=kjet
PARTITION_WGRIB2=kjet
+ PARTITION_POST=kjet
QUEUE_DEFAULT="rt"
QUEUE_FCST="rt"
@@ -88,9 +89,11 @@ readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_
beta1_inv=0.5 #beata_inv is 1-ensemble_wgt
ens_h=20 #110
ens_v=1 #3
-regional_ensemble_option=1 #1 for GDAS
+regional_ensemble_option=1 #1 for GDAS, 5 for RRFSDAS
grid_ratio_fv3=2.0 #fv3 resolution 3km, so analysis=3*2=6km
grid_ratio_ens=3 #analysis 6km, so ensemble=6*3=18km. GDAS ensemble is 20km
+##MEMO1: regional_ensemble_option=5 will be changed to 1 if there are not enough RRFS ensembles at the runtime.
+##MEMO2: grid_ratio_ens will be hardwired to 1 for RRFSDAS, so grid_ratio_ens here only affects GDAS.
i_en_perts_io=1 #0 or 1: original file 3: pre-processed ensembles
# &RAPIDREFRESH_CLDSURF
@@ -98,6 +101,9 @@ l_PBL_pseudo_SurfobsT=.false.
l_PBL_pseudo_SurfobsQ=.false
i_use_2mQ4B=0
i_use_2mT4B=0
+i_T_Q_adjust=1
+l_rtma3d=.true.
+i_precip_vertical_check=3
#-------------------------------------------------------------------------------------
#
BERROR_FN="berror.rtma" #under $FIX_GSI
diff --git a/ush/config_defaults.sh b/ush/config_defaults.sh
index e5476bcf2..0430cd84b 100644
--- a/ush/config_defaults.sh
+++ b/ush/config_defaults.sh
@@ -633,6 +633,9 @@ l_PBL_pseudo_SurfobsT=.false.
l_PBL_pseudo_SurfobsQ=.false.
i_use_2mQ4B=0
i_use_2mT4B=0
+i_T_Q_adjust=1
+l_rtma3d=.false.
+i_precip_vertical_check=0
#-----------------------------------------------------------------------
# HYBENSMEM_NMIN:
# Minimum number of ensemble members required a hybrid GSI analysis
diff --git a/ush/templates/FV3LAM_wflow.xml b/ush/templates/FV3LAM_wflow.xml
index a4359d81f..843b6cbda 100644
--- a/ush/templates/FV3LAM_wflow.xml
+++ b/ush/templates/FV3LAM_wflow.xml
@@ -1130,7 +1130,7 @@ MODULES_RUN_TASK_FP script.
{%- endif %}
- {%- else %}
+ {%- elif do_envar_radar_ref %}
@@ -1147,7 +1147,10 @@ MODULES_RUN_TASK_FP script.
+ {%- else %}
+
{%- endif %}
+
@@ -1234,6 +1237,7 @@ MODULES_RUN_TASK_FP script.
************************************************************************
************************************************************************
-->
+{% if not is_rtma %}
{{ restart_hrs_prod }}
@@ -1272,7 +1276,7 @@ MODULES_RUN_TASK_FP script.
-
+{%- endif %}