Make needed updates to run forecast from GEFS

parm/config/config.base.emc.dyn

@@ -340,6 +340,7 @@ if [ $DOHYBVAR = "YES" ]; then
     export FHMIN_ENKF=3
     export FHMAX_ENKF=9
     export FHMAX_ENKF_GFS=120
+    export FHOUT_ENKF_GFS=3
     if [ $l4densvar = ".true." ]; then
         export FHOUT=1
         export FHOUT_ENKF=1

parm/config/config.resources

@@ -737,7 +737,7 @@ elif [ ${step} = "esfc" ]; then
 elif [ ${step} = "efcs" ]; then
 
     if [ ${CASE} = "C768" ]; then
-      export wtime_efcs="01:20:00"
+      export wtime_efcs="06:00:00"
     else
       export wtime_efcs="00:40:00"
     fi

scripts/exgdas_enkf_fcst.sh

@@ -108,7 +108,9 @@ export LEVS=${LEVS_ENKF:-${LEVS:-64}}
 
 # nggps_diag_nml
 export FHOUT=${FHOUT_ENKF:-3}
-
+if [[ $CDUMP == "gfs" ]]; then
+    export FHOUT=${FHOUT_ENKF_GFS:-${FHOUT_ENKF:${FHOUT:-3}}}
+fi
 # model_configure
 export DELTIM=${DELTIM_ENKF:-${DELTIM:-225}}
 export FHMAX=${FHMAX_ENKF:-9}

ush/forecast_det.sh

@@ -41,7 +41,7 @@ FV3_GFS_det(){
   # determine if restart IC exists to continue from a previous forecast
   RERUN="NO"
   filecount=$(find "${RSTDIR_ATM:-/dev/null}" -type f | wc -l)
-  if [ ${CDUMP} = "gfs" -a ${rst_invt1} -gt 0 -a ${FHMAX} -gt ${rst_invt1} -a ${filecount} -gt 10 ]; then
+  if [[ ( ${CDUMP} = "gfs" || ( ${RUN} = "gefs" && ${CDATE_RST} = "" )) && ${rst_invt1} -gt 0 && ${FHMAX} -gt ${rst_invt1} && ${filecount} -gt 10 ]]; then
     reverse=$(echo "${restart_interval[@]} " | tac -s ' ')
     for xfh in ${reverse} ; do
       yfh=$((xfh-(IAU_OFFSET/2)))
@@ -76,6 +76,7 @@ FV3_GFS_det(){
 
 FV3_GEFS_det(){
   echo "SUB ${FUNCNAME[0]}: Defining variables for FV3GEFS"
+  FV3_GFS_det
 }
 
 WW3_det(){

ush/forecast_postdet.sh

@@ -14,6 +14,7 @@
 FV3_GEFS_postdet(){
   echo SUB ${FUNCNAME[0]}: Linking input data for FV3 $RUN
   # soft link commands insert here
+  FV3_GFS_postdet
 }
 
 DATM_postdet(){
@@ -307,9 +308,9 @@ EOF
   # inline post fix files
   if [ $WRITE_DOPOST = ".true." ]; then
     $NLN $PARM_POST/post_tag_gfs${LEVS}             $DATA/itag
-    $NLN $PARM_POST/postxconfig-NT-GFS-TWO.txt      $DATA/postxconfig-NT.txt
-    $NLN $PARM_POST/postxconfig-NT-GFS-F00-TWO.txt  $DATA/postxconfig-NT_FH00.txt
-    $NLN $PARM_POST/params_grib2_tbl_new            $DATA/params_grib2_tbl_new
+    $NLN ${FLTFILEGFS:-$PARM_POST/postxconfig-NT-GFS-TWO.txt}           $DATA/postxconfig-NT.txt
+    $NLN ${FLTFILEGFSF00:-$PARM_POST/postxconfig-NT-GFS-F00-TWO.txt}    $DATA/postxconfig-NT_FH00.txt
+    $NLN ${POSTGRB2TBL:-$PARM_POST/params_grib2_tbl_new}                $DATA/params_grib2_tbl_new
   fi
 
   #------------------------------------------------------------------
@@ -519,11 +520,19 @@ EOF
       logi=logf${FH3}
       pgbi=GFSPRS.GrbF${FH2}
       flxi=GFSFLX.GrbF${FH2}
-      atmo=$memdir/${CDUMP}.t${cyc}z.atmf${FH3}.$affix
-      sfco=$memdir/${CDUMP}.t${cyc}z.sfcf${FH3}.$affix
-      logo=$memdir/${CDUMP}.t${cyc}z.logf${FH3}.txt
-      pgbo=$memdir/${CDUMP}.t${cyc}z.master.grb2f${FH3}
-      flxo=$memdir/${CDUMP}.t${cyc}z.sfluxgrbf${FH3}.grib2
+      if [[ $RUN == "gefs" ]]; then
+        atmo=$memdir/history/${CDUMP}.t${cyc}z.atmf${FH3}.$affix
+        sfco=$memdir/history/${CDUMP}.t${cyc}z.sfcf${FH3}.$affix
+        logo=$memdir/history/${CDUMP}.t${cyc}z.logf${FH3}.txt
+        pgbo=$memdir/master/${CDUMP}.t${cyc}z.master.grb2f${FH3}
+        flxo=$memdir/master/${CDUMP}.t${cyc}z.sfluxgrbf${FH3}.grib2
+      else
+        atmo=$memdir/${CDUMP}.t${cyc}z.atmf${FH3}.$affix
+        sfco=$memdir/${CDUMP}.t${cyc}z.sfcf${FH3}.$affix
+        logo=$memdir/${CDUMP}.t${cyc}z.logf${FH3}.txt
+        pgbo=$memdir/${CDUMP}.t${cyc}z.master.grb2f${FH3}
+        flxo=$memdir/${CDUMP}.t${cyc}z.sfluxgrbf${FH3}.grib2
+      fi
       eval $NLN $atmo $atmi
       eval $NLN $sfco $sfci
       eval $NLN $logo $logi
@@ -534,11 +543,19 @@ EOF
     done
   else
     for n in $(seq 1 $ntiles); do
-      eval $NLN nggps2d.tile${n}.nc       $memdir/nggps2d.tile${n}.nc
-      eval $NLN nggps3d.tile${n}.nc       $memdir/nggps3d.tile${n}.nc
-      eval $NLN grid_spec.tile${n}.nc     $memdir/grid_spec.tile${n}.nc
-      eval $NLN atmos_static.tile${n}.nc  $memdir/atmos_static.tile${n}.nc
-      eval $NLN atmos_4xdaily.tile${n}.nc $memdir/atmos_4xdaily.tile${n}.nc
+      if [[ $RUN == "gefs" ]]; then
+        eval $NLN nggps2d.tile${n}.nc       $memdir/history/nggps2d.tile${n}.nc
+        eval $NLN nggps3d.tile${n}.nc       $memdir/history/nggps3d.tile${n}.nc
+        eval $NLN grid_spec.tile${n}.nc     $memdir/history/grid_spec.tile${n}.nc
+        eval $NLN atmos_static.tile${n}.nc  $memdir/history/atmos_static.tile${n}.nc
+        eval $NLN atmos_4xdaily.tile${n}.nc $memdir/history/atmos_4xdaily.tile${n}.nc
+      else
+        eval $NLN nggps2d.tile${n}.nc       $memdir/nggps2d.tile${n}.nc
+        eval $NLN nggps3d.tile${n}.nc       $memdir/nggps3d.tile${n}.nc
+        eval $NLN grid_spec.tile${n}.nc     $memdir/grid_spec.tile${n}.nc
+        eval $NLN atmos_static.tile${n}.nc  $memdir/atmos_static.tile${n}.nc
+        eval $NLN atmos_4xdaily.tile${n}.nc $memdir/atmos_4xdaily.tile${n}.nc
+      fi
     done
   fi
 }
@@ -556,6 +573,20 @@ FV3_GFS_nml(){
   echo SUB ${FUNCNAME[0]}: FV3 name lists and model configure file created
 }
 
+FV3_GEFS_nml(){
+  # namelist output for a certain component
+  echo SUB ${FUNCNAME[0]}: Creating name lists and model configure file for FV3
+  if [ $machine = 'sandbox' ]; then
+    cd $SCRIPTDIR
+    echo "MAIN: !!!Sandbox mode, writing to current directory!!!"
+  fi
+  # Call child scripts in current script directory
+  source $SCRIPTDIR/parsing_namelists_FV3.sh
+  FV3_namelists
+  echo SUB ${FUNCNAME[0]}: FV3 name lists and model configure file created
+}
+
+
 DATM_nml(){
   source $SCRIPTDIR/parsing_namelists_DATM.sh
   DATM_namelists
@@ -599,14 +630,22 @@ data_out_GFS() {
           $NCP $file $memdir/RESTART/$file
         done
       fi
-    elif [ $CDUMP = "gfs" ]; then
-      $NCP $DATA/input.nml $ROTDIR/${CDUMP}.${PDY}/${cyc}/atmos/
+    elif [ $CDUMP = "gfs" ] || [ ${RUN} = "gefs" ]; then
+      if [ ${RUN} = "gefs" ]; then
+        $NCP $DATA/input.nml ${ROTDIR}/${RUN}.${PDY}/$cyc/${mem}/atmos/history/
+      else
+        $NCP $DATA/input.nml $ROTDIR/${CDUMP}.${PDY}/${cyc}/atmos/
+      fi
     fi
   fi
 
   echo "SUB ${FUNCNAME[0]}: Output data for FV3 copied"
 }
 
+data_out_GEFS() {
+  data_out_GFS
+}
+
 WW3_postdet() {
   echo "SUB ${FUNCNAME[0]}: Linking input data for WW3"
   COMPONENTwave=${COMPONENTwave:-${RUN}wave}

ush/forecast_predet.sh

@@ -208,8 +208,8 @@ FV3_GFS_predet(){
   print_freq=${print_freq:-6}
 
   #-------------------------------------------------------
-  if [ $CDUMP = "gfs" ] && [ $rst_invt1 -gt 0 ] && [ $MEMBER -lt 0 ]; then
-    RSTDIR_ATM=${RSTDIR:-$ROTDIR}/${CDUMP}.${PDY}/${cyc}/atmos/RERUN_RESTART
+  if [[ $CDUMP = "gfs" || $RUN = "gefs" ]] && [ $rst_invt1 -gt 0 ] && [ $MEMBER -lt 0 ]; then
+    RSTDIR_ATM=${RSTDIR_ATM:-${RSTDIR:-$ROTDIR}/${CDUMP}.${PDY}/${cyc}/atmos/RERUN_RESTART}
     if [ ! -d $RSTDIR_ATM ]; then mkdir -p $RSTDIR_ATM ; fi
     $NLN $RSTDIR_ATM RESTART
     # The final restart written at the end doesn't include the valid date
@@ -232,7 +232,7 @@ FV3_GFS_predet(){
 
   #-------------------------------------------------------
   # member directory
-  if [ $MEMBER -lt 0 ]; then
+  if [[ $MEMBER -lt 0 || $RUN = "gefs" ]]; then
     prefix=$CDUMP
     rprefix=$rCDUMP
     memchar=""
@@ -241,13 +241,19 @@ FV3_GFS_predet(){
     rprefix=enkf$rCDUMP
     memchar=mem$(printf %03i $MEMBER)
   fi
-  memdir=$ROTDIR/${prefix}.$PDY/$cyc/atmos/$memchar
+  memdir=${memdir:-$ROTDIR/${prefix}.$PDY/$cyc/atmos/$memchar}
   if [ ! -d $memdir ]; then mkdir -p $memdir; fi
+  if [[ $RUN == "gefs" ]]; then
+    if [ ! -d ${memdir}/history ]; then mkdir -p ${memdir}/history; fi
+    if [ $WRITE_DOPOST = ".true." ]; then
+      if [ ! -d ${memdir}/master ]; then mkdir -p ${memdir}/master; fi
+    fi
+  fi
 
   GDATE=$($NDATE -$assim_freq $CDATE)
   gPDY=$(echo $GDATE | cut -c1-8)
   gcyc=$(echo $GDATE | cut -c9-10)
-  gmemdir=$ROTDIR/${rprefix}.$gPDY/$gcyc/atmos/$memchar
+  gmemdir=${gmemdir:-$ROTDIR/${rprefix}.$gPDY/$gcyc/atmos/$memchar}
 
   if [[ "$DOIAU" = "YES" ]]; then
     sCDATE=$($NDATE -3 $CDATE)
@@ -266,6 +272,14 @@ FV3_GFS_predet(){
   echo "SUB ${FUNCNAME[0]}: pre-determination variables set"
 }
 
+FV3_GEFS_predet(){
+  echo "SUB ${FUNCNAME[0]}: Defining variables for FV3GFS"
+
+  FV3_GFS_predet
+
+  echo "SUB ${FUNCNAME[0]}: pre-determination variables set"
+}
+
 WW3_predet(){
   echo "SUB ${FUNCNAME[0]}: Defining variables for WW3"
   if [ $CDUMP = "gdas" ]; then

workflow/rocoto/workflow_tasks.py

@@ -1259,9 +1259,10 @@ def epos(self):
         def _get_eposgroups(epos):
             fhmin = epos['FHMIN_ENKF']
             fhmax = epos['FHMAX_ENKF']
+            fhout = epos['FHOUT_ENKF']
             if self.cdump == "gfs":
                 fhmax = epos['FHMAX_ENKF_GFS']
-            fhout = epos['FHOUT_ENKF']
+                fhout = epos['FHOUT_ENKF_GFS']
             fhrs = range(fhmin, fhmax + fhout, fhout)
 
             neposgrp = epos['NEPOSGRP']