Bugfixes from IMS C++ code additions

scripts/exglobal_snow_analysis.py

@@ -24,7 +24,7 @@
     snow_anl.initialize()
 
     # Process IMS snow cover (if applicable)
-    if snow_anl.task_config.cyc == 0:
+    if snow_anl.task_config.DO_IMS_SCF:
         snow_anl.execute('scf_to_ioda')
 
     # Execute JEDI snow analysis

scripts/exglobal_snowens_analysis.py

@@ -30,7 +30,7 @@
     # stage ensemble mean backgrounds
 
     # Process IMS snow cover (if applicable)
-    if snow_ens_anl.task_config.cyc == 0:
+    if snow_ens_anl.task_config.DO_IMS_SCF:
         snow_ens_anl.execute('scf_to_ioda')
 
     # Execute JEDI snow analysis

ush/python/pygfs/task/snow_analysis.py

@@ -83,6 +83,9 @@ def __init__(self, config: Dict[str, Any]):
         expected_keys = ['scf_to_ioda', 'snowanlvar']
         self.jedi_dict = Jedi.get_jedi_dict(self.task_config.JEDI_CONFIG_YAML, self.task_config, expected_keys)
 
+        # Boolean to decide if IMS snow cover processing is done
+        self.task_config.DO_IMS_SCF = False
+
     @logit(logger)
     def initialize(self) -> None:
         """Initialize a global snow analysis
@@ -147,6 +150,10 @@ def initialize(self) -> None:
             ims_scf_to_ioda_staging_dict = parse_j2yaml(self.task_config.STAGE_IMS_SCF2IODA_YAML, self.task_config)
             FileHandler(ims_scf_to_ioda_staging_dict).sync()
             self.jedi_dict['scf_to_ioda'].initialize(self.task_config)
+            # Check if file exists
+            ims_file = ims_scf_to_ioda_staging_dict['copy_opt'][0][1]
+            if os.path.exists(ims_file):
+                self.task_config.DO_IMS_SCF = True
 
         # initialize JEDI variational application
         logger.info(f"Initializing JEDI variational DA application")

ush/python/pygfs/task/snowens_analysis.py

@@ -84,6 +84,9 @@ def __init__(self, config: Dict[str, Any]):
         expected_keys = ['scf_to_ioda', 'snowanlvar', 'esnowanlensmean']
         self.jedi_dict = Jedi.get_jedi_dict(self.task_config.JEDI_CONFIG_YAML, self.task_config, expected_keys)
 
+        # Boolean to decide if IMS snow cover processing is done
+        self.task_config.DO_IMS_SCF = False
+
     @logit(logger)
     def initialize(self) -> None:
         """Initialize a global snow ensemble analysis
@@ -165,6 +168,10 @@ def initialize(self) -> None:
             ims_scf_to_ioda_staging_dict = parse_j2yaml(self.task_config.STAGE_IMS_SCF2IODA_YAML, self.task_config)
             FileHandler(ims_scf_to_ioda_staging_dict).sync()
             self.jedi_dict['scf_to_ioda'].initialize(self.task_config)
+            # Check if file exists
+            ims_file = ims_scf_to_ioda_staging_dict['copy_opt'][0][1]
+            if os.path.exists(ims_file):
+                self.task_config.DO_IMS_SCF = True
 
         # initialize JEDI variational application
         logger.info(f"Initializing JEDI variational DA application")