add $HYBENSMEM_NMIN which sets the minimum number of ensembles requir…

…ed for a hybrid GSI analysis (NOAA-EMC#194)

 add a variable HYBENSMEM_NMIN
 GSI will do hybrid analysis if $HYBENSMEM_NMIN or more ensembles are available
 This allow GSI to do hybrid analysis even we miss a few members of 80 GDAS ensembles.

scripts/exregional_run_analysis.sh

@@ -281,10 +281,13 @@ ifhyb=.false.
 memname='atmf009'
 nummem=$(more filelist03 | wc -l)
 nummem=$((nummem - 3 ))
-if [[ ${nummem} -eq 80 ]]; then
+if [[ ${nummem} -ge ${HYBENSMEM_NMIN} ]]; then
   print_info_msg "$VERBOSE" "Do hybrid with ${memname}"
   ifhyb=.true.
   print_info_msg "$VERBOSE" " Cycle ${YYYYMMDDHH}: GSI hybrid uses ${memname} with n_ens=${nummem}" 
+else
+  print_info_msg "$VERBOSE" " Cycle ${YYYYMMDDHH}: GSI does pure 3DVAR."
+  print_info_msg "$VERBOSE" " Hybrid needs at least ${HYBENSMEM_NMIN} ${memname} ensembles, only ${nummem} available"
 fi
 
 #

ush/config.sh.3DRTMA_dev1

@@ -209,6 +209,7 @@ DA_CYCLE_INTERV="1" #need to be 1 to have hourly DA
 RESTART_INTERVAL="1"
 
 NCORES_RUN_ANAL=240
+HYBENSMEM_NMIN=66
 HALO_BLEND=10
 #LAYOUT_X="17"
 #LAYOUT_Y="10"

ush/config.sh.RRFS_NA_3km

@@ -114,6 +114,7 @@ RESTART_INTERVAL="1 2"
 
 NCORES_RUN_ANAL=240
 WTIME_RUN_ANAL="01:00:00"
+HYBENSMEM_NMIN=66
 HALO_BLEND=10
 
 NNODES_MAKE_ICS="10"

ush/config.sh.RRFS_dev1

@@ -164,6 +164,7 @@ DA_CYCLE_INTERV="1"
 RESTART_INTERVAL="1 2 3 4 5 6"
 
 NCORES_RUN_ANAL=240
+HYBENSMEM_NMIN=66
 HALO_BLEND=10
 
 if [[ -n $RESERVATION ]] ; then

ush/config.sh.RTMA_NA_3km

@@ -166,6 +166,7 @@ RESTART_INTERVAL="1"
 
 NCORES_RUN_ANAL=240
 WTIME_RUN_ANAL="01:00:00"
+HYBENSMEM_NMIN=66
 HALO_BLEND=10
 
 NNODES_MAKE_ICS="10"

ush/config_defaults.sh

@@ -581,7 +581,10 @@ l_PBL_pseudo_SurfobsQ=.false.
 i_use_2mQ4B=0
 i_use_2mT4B=0
 #-----------------------------------------------------------------------
+# HYBENSMEM_NMIN:
+#    Minimum number of ensemble members required a hybrid GSI analysis 
 #
+HYBENSMEM_NMIN=80
 ANAVINFO_FN="anavinfo.fv3lam_hrrr"
 CONVINFO_FN="convinfo.rrfs"
 BERROR_FN="rap_berror_stats_global_RAP_tune" #under $FIX_GSI