diff --git a/sorc/gdas.cd b/sorc/gdas.cd
index 15113adfdb..417ed5e686 160000
--- a/sorc/gdas.cd
+++ b/sorc/gdas.cd
@@ -1 +1 @@
-Subproject commit 15113adfdbf2500ec2d5099fc9b62b21fbfcc9b8
+Subproject commit 417ed5e686ae73af7be65350fda20037c5934242
diff --git a/ush/python/pygfs/task/marine_analysis.py b/ush/python/pygfs/task/marine_analysis.py
index dd48c9f80f..d4ee08ba32 100644
--- a/ush/python/pygfs/task/marine_analysis.py
+++ b/ush/python/pygfs/task/marine_analysis.py
@@ -281,6 +281,8 @@ def _prep_checkpoint(self: Task) -> None:
         soca2cice_param = AttrDict({
             "ocn_ana": f"./Data/ocn.3dvarfgat_pseudo.an.{self.task_config.MARINE_WINDOW_MIDDLE_ISO}.nc",
             "ice_ana": f"./Data/ice.3dvarfgat_pseudo.an.{self.task_config.MARINE_WINDOW_MIDDLE_ISO}.nc",
+            "ocn_inc": f"./Data/ocn.3dvarfgat_pseudo.incr.{self.task_config.MARINE_WINDOW_MIDDLE_ISO}.nc",
+            "ice_inc": f"./Data/ice.3dvarfgat_pseudo.incr.{self.task_config.MARINE_WINDOW_MIDDLE_ISO}.nc",
             "ice_rst": ice_rst_ana,
             "fcst_begin": fcst_begin
         })
@@ -288,11 +290,10 @@ def _prep_checkpoint(self: Task) -> None:
 
         # render the SOCA to CICE YAML file for the Arctic and Antarctic
         logger.info("render the SOCA to CICE YAML file for the Arctic and Antarctic")
-        varchgyamls = ['soca_2cice_global.yaml']
-        for varchgyaml in varchgyamls:
-            soca2cice_config = parse_j2yaml(path=os.path.join(self.task_config.MARINE_JCB_GDAS_ALGO, f'{varchgyaml}.j2'),
-                                            data=soca2cice_param)
-            soca2cice_config.save(os.path.join(self.task_config.DATA, varchgyaml))
+        varchgyaml = 'soca_2cice_global.yaml'
+        soca2cice_config = parse_j2yaml(path=os.path.join(self.task_config.MARINE_JCB_GDAS_ALGO, f'{varchgyaml}.j2'),
+                                        data=soca2cice_param)
+        soca2cice_config.save(os.path.join(self.task_config.DATA, varchgyaml))
 
     @logit(logger)
     def variational(self: Task) -> None:
@@ -388,6 +389,10 @@ def list_all_files(dir_in, dir_out, wc='*', fh_list=[]):
             post_file_list.append([os.path.join(anl_dir, 'Data', f'{domain}.3dvarfgat_pseudo.an.{mdate}.nc'),
                                    os.path.join(com_ocean_analysis, f'{RUN}.t{cyc}z.{domain}ana.nc')])
 
+        # Copy soca2cice ice increment
+        post_file_list.append([os.path.join(anl_dir, 'Data', f'ice.soca2cice.incr.{bdate}.nc'),
+                              os.path.join(com_ocean_analysis, f'{RUN}.t{cyc}z.ice.incr.postproc.nc')])
+
         # Copy of the ssh diagnostics
         if nmem_ens > 2:
             for string in ['ssh_steric_stddev', 'ssh_unbal_stddev', 'ssh_total_stddev', 'steric_explained_variance']: