Add esnowanl to enkfgfs cycle (#3283)

This PR adds the capability for this task to run for the enkfgfs
cycle, incorporating the rotating background directory structure used by
other enkfgfs tasks.

Resolves #3276 

---------

Co-authored-by: Rahul Mahajan <aerorahul@users.noreply.github.com>

Author: TravisElless-NOAA

jobs/JGLOBAL_SNOWENS_ANALYSIS

@@ -13,6 +13,13 @@ export GDUMP
 CDUMP=${RUN/enkf}
 export CDUMP
 
+export NMEM_ENS_MAX=${NMEM_ENS}
+if [[ "${RUN}" == "enkfgfs" ]]; then
+   NMEM_ENS=${NMEM_ENS_GFS}
+   export mem_offset=$((NMEM_ENS_GFS_OFFSET * cyc/6))
+else
+   export mem_offset=0
+fi
 ##############################################
 # Begin JOB SPECIFIC work
 ##############################################

parm/gdas/snow_stage_ens_update.yaml.j2

@@ -22,13 +22,17 @@ copy:
 # copy ensemble background files
 ######################################
 {% for mem in range(1, NMEM_ENS + 1) %}
+    {% set gmem = mem+mem_offset %}
+    {% if gmem > NMEM_ENS_MAX %}
+        {% set gmem = gmem-NMEM_ENS_MAX %}
+    {% endif %}
     # define variables
     # Declare a dict of search and replace terms to run on each template
     {% set tmpl_dict = {'${ROTDIR}':ROTDIR,
-                        '${RUN}':RUN,
+                        '${RUN}':"enkfgdas",
                         '${YMD}':previous_cycle | to_YMD,
                         '${HH}':previous_cycle | strftime("%H"),
-                        '${MEMDIR}':"mem" + '%03d' % mem} %}
+                        '${MEMDIR}':"mem" + '%03d' % gmem} %}
 
     # copy coupler file
 - ["{{ COM_ATMOS_RESTART_TMPL | replace_tmpl(tmpl_dict) }}/{{ current_cycle | to_fv3time }}.coupler.res", "{{ DATA }}/bkg/mem{{ '%03d' % mem }}/{{ current_cycle | to_fv3time }}.coupler.res"]

workflow/applications/gfs_cycled.py

@@ -313,8 +313,8 @@ def get_task_names(self):
                     task_names[run] += ['eobs', 'eupd']
                     task_names[run].append('echgres') if 'gdas' in run else 0
                     task_names[run] += ['ediag'] if options['lobsdiag_forenkf'] else ['eomg']
-                    task_names[run].append('esnowanl') if options['do_jedisnowda'] and 'gdas' in run else 0
 
+                task_names[run].append('esnowanl') if options['do_jedisnowda'] else 0
                 task_names[run].append('efcs') if 'gdas' in run else 0
                 task_names[run].append('epos') if 'gdas' in run else 0
                 task_names[run] += ['stage_ic', 'ecen', 'esfc', 'earc', 'cleanup']

workflow/rocoto/gfs_tasks.py

@@ -619,7 +619,7 @@ def snowanl(self):
     def esnowanl(self):
 
         deps = []
-        dep_dict = {'type': 'metatask', 'name': f'{self.run}_epmn', 'offset': f"-{timedelta_to_HMS(self._base['interval_gdas'])}"}
+        dep_dict = {'type': 'metatask', 'name': 'enkfgdas_epmn', 'offset': f"-{timedelta_to_HMS(self._base['interval_gdas'])}"}
         deps.append(rocoto.add_dependency(dep_dict))
         dep_dict = {'type': 'task', 'name': f"{self.run.replace('enkf', '')}_prep"}
         deps.append(rocoto.add_dependency(dep_dict))