From 417ed5e686ae73af7be65350fda20037c5934242 Mon Sep 17 00:00:00 2001
From: Anna Shlyaeva <annash@ucar.edu>
Date: Wed, 22 Jan 2025 14:21:42 -0700
Subject: [PATCH] Plot postprocessed ice concentration increment (#1446)

# Description

Add plots to offline vrfy to output postprocessed (aka what get added to
CICE restart) sea ice increment.

# Companion PRs

https://github.com/NOAA-EMC/jcb-gdas/pull/64
https://github.com/NOAA-EMC/global-workflow/pull/3235

# Issues

Resolves #1402

# Automated CI tests to run in Global Workflow

Since it's offline vrfy, no tests would cover it currently.

<!-- Which Global Workflow CI tests are required to adequately test this
PR? -->
- [ ] atm_jjob <!-- JEDI atm single cycle DA !-->
- [ ] C96C48_ufs_hybatmDA <!-- JEDI atm cycled DA !-->
- [ ] C96C48_hybatmaerosnowDA  <!-- JEDI aero/snow cycled DA !-->
- [ ] C48mx500_3DVarAOWCDA <!-- JEDI low-res marine 3DVar cycled DA !-->
- [ ] C48mx500_hybAOWCDA <!-- JEDI marine hybrid envar cycled DA !-->
- [ ] C96C48_hybatmDA <!-- GSI atm cycled DA !-->

---------

Co-authored-by: Anna Shlyaeva <anna.v.shlyaeva@noaa.gov>
---
 parm/jcb-gdas                                         |  2 +-
 .../exgdas_global_marine_analysis_vrfy_manual.py      | 11 ++++++++++-
 2 files changed, 11 insertions(+), 2 deletions(-)

diff --git a/parm/jcb-gdas b/parm/jcb-gdas
index 10a744aa2..9e069cc58 160000
--- a/parm/jcb-gdas
+++ b/parm/jcb-gdas
@@ -1 +1 @@
-Subproject commit 10a744aa2a8103e523aba54bebc44a965c2ec466
+Subproject commit 9e069cc585e8e6b09643c350cd0cc9c7d6909119
diff --git a/utils/soca/fig_gallery/exgdas_global_marine_analysis_vrfy_manual.py b/utils/soca/fig_gallery/exgdas_global_marine_analysis_vrfy_manual.py
index 7c8efd0a6..ccb37a55c 100644
--- a/utils/soca/fig_gallery/exgdas_global_marine_analysis_vrfy_manual.py
+++ b/utils/soca/fig_gallery/exgdas_global_marine_analysis_vrfy_manual.py
@@ -146,7 +146,16 @@
                                                'hs_h': [-0.1, 0.1]},
                               colormap='seismic',
                               projs=['North', 'South'],
-                              comout=os.path.join(comout, 'vrfy', 'incr'))]   # sea ice increment
+                              comout=os.path.join(comout, 'vrfy', 'incr')),   # sea ice increment
+                   plotConfig(grid_file=grid_file,
+                              data_file=os.path.join(comout, f'{RUN}.t'+cyc+'z.ice.incr.postproc.nc'),
+                              lats=np.arange(-60, 60, 10),
+                              variables_horiz={'aice_h': [-0.2, 0.2],
+                                               'hi_h': [-0.5, 0.5],
+                                               'hs_h': [-0.1, 0.1]},
+                              colormap='seismic',
+                              projs=['North', 'South'],
+                              comout=os.path.join(comout, 'vrfy', 'incr.postproc'))]   # sea ice increment after postprocessing
     configs.extend(config_incr)