Add code changes for the implementation of the low-resolution B

[guillaumevernieres]

As the title says.

Dependencies

• Configuration for the low res marine B jcb-gdas#66
• Add capability to run diff resolutions for marine anl and background global-workflow#3238
• Update of the static files in glopara: Update to the marine JEDI fix files global-workflow#3239

Work done in this PR

Interp to low-resolution:

• variance partitioning
• preparation of the perturbations for the envar
• hybrid weights
• ...

Testing

Visual check showing the variance explained by the steric height (30 members)

left is native resolution, right is on the 1/2 deg grid.

"Science" note: The explained variance above is almost the opposite to what I would have expected, and what we saw in the previous offline ensemble, or even what we see in tendencies. Check if there's a bug in the explained variance calculation.

Automated CI tests to run in Global Workflow

• atm_jjob
• C96C48_ufs_hybatmDA
• C96C48_hybatmaerosnowDA
• C48mx500_3DVarAOWCDA
• C48mx500_hybAOWCDA
• C96C48_hybatmDA

[emcbot]

Automated GW-GDASApp Testing Results:
Machine: hera

Start: Mon Jan 13 21:38:37 UTC 2025 on hfe10
---------------------------------------------------
Build:                                 *SUCCESS*
Build: Completed at Mon Jan 13 22:32:01 UTC 2025
---------------------------------------------------
Tests: ctest -j12 -R gdasapp -E atm_jjob|C96C48_ufs_hybatmDA|C96C48_hybatmaerosnowDA|C96C48_hybatmDA
Tests:                                 *SUCCESS*
Tests: Completed at Mon Jan 13 22:55:24 UTC 2025
Tests: 100% tests passed, 0 tests failed out of 53

[shlyaeva]

Did not test, but I reviewed all the code changes, they look good to me.

[emcbot]

Automated GW-GDASApp Testing Results:
Machine: hera

Start: Wed Jan 15 21:37:00 UTC 2025 on hfe09
---------------------------------------------------
Build:                                 *SUCCESS*
Build: Completed at Wed Jan 15 22:26:38 UTC 2025
---------------------------------------------------
Tests: ctest -j12 -R gdasapp -E atm_jjob|C96C48_ufs_hybatmDA|C96C48_hybatmaerosnowDA|C96C48_hybatmDA
Tests:                                  *Failed*
Tests: Failed at Wed Jan 15 22:47:43 UTC 2025
Tests: 92% tests passed, 4 tests failed out of 53
	2042 - test_gdasapp_C48mx500_3DVarAOWCDA_gdas_marineanlchkpt_202103250000 (Failed)
	2043 - test_gdasapp_C48mx500_3DVarAOWCDA_gdas_marineanlfinal_202103250000 (Failed)
	2054 - test_gdasapp_C48mx500_hybAOWCDA_gdas_marineanlchkpt_202103250000 (Failed)
	2055 - test_gdasapp_C48mx500_hybAOWCDA_gdas_marineanlfinal_202103250000 (Failed)
Tests: see output at /scratch1/NCEPDEV/da/role.jedipara/CI/GDASApp/workflow/PR/1441/global-workflow/sorc/gdas.cd/build/log.ctest

[emcbot]

Automated GW-GDASApp Testing Results:
Machine: hera

Start: Mon Feb  3 23:46:06 UTC 2025 on hfe07
---------------------------------------------------
Build:                                 *SUCCESS*
Build: Completed at Tue Feb  4 00:34:26 UTC 2025
---------------------------------------------------
Tests: ctest -j12 -R gdasapp -E atm_jjob|C96C48_ufs_hybatmDA|C96C48_hybatmaerosnowDA|C96C48_hybatmDA
Tests:                                 *SUCCESS*
Tests: Completed at Tue Feb  4 01:39:01 UTC 2025
Tests: 100% tests passed, 0 tests failed out of 54

[emcbot]

Automated GW-GDASApp Testing Results:
Machine: hercules

Start: Mon Feb  3 12:19:57 CST 2025 on hercules-login-1.hpc.msstate.edu
---------------------------------------------------
Build:                                 *SUCCESS*
Build: Completed at Mon Feb  3 13:00:49 CST 2025
---------------------------------------------------
Tests: ctest -j12 -R gdasapp -E atm_jjob|C96C48_ufs_hybatmDA|C96C48_hybatmaerosnowDA|C96C48_hybatmDA
Tests:                                  *Failed*
Tests: Failed at Mon Feb  3 19:55:29 CST 2025
Tests: 98% tests passed, 1 tests failed out of 54
	2039 - test_gdasapp_C48mx500_3DVarAOWCDA_gdas_waveinit_202103241800 (Timeout)
Tests: see output at /work2/noaa/da/role-da/CI/hercules/GDASApp/workflow/PR/1441/global-workflow/sorc/gdas.cd/build/log.ctest

[guillaumevernieres]

Weird, it looks like it used the wrong commit of the g-w branch. I merged develop before resubmitting the ci but that;s not what was cloned here: /work2/noaa/da/role-da/CI/hercules/GDASApp/workflow/PR/1441/global-workflow

[emcbot]

Automated GW-GDASApp Testing Results:
Machine: orion

Start: Tue Feb  4 08:04:59 AM CST 2025 on orion-login-1.hpc.msstate.edu
---------------------------------------------------
Build:                                 *SUCCESS*
Build: Completed at Tue Feb  4 09:11:58 AM CST 2025
---------------------------------------------------
Tests: ctest -j12 -R gdasapp -E atm_jjob|C96C48_ufs_hybatmDA|C96C48_hybatmaerosnowDA|C96C48_hybatmDA
Tests:                                 *SUCCESS*
Tests: Completed at Tue Feb  4 09:39:22 AM CST 2025
Tests: 100% tests passed, 0 tests failed out of 54

[emcbot]

Automated GW-GDASApp Testing Results:
Machine: hercules

Start: Tue Feb  4 08:05:32 CST 2025 on hercules-login-1.hpc.msstate.edu
---------------------------------------------------
Build:                                 *SUCCESS*
Build: Completed at Tue Feb  4 08:46:00 CST 2025
---------------------------------------------------
Tests: ctest -j12 -R gdasapp -E atm_jjob|C96C48_ufs_hybatmDA|C96C48_hybatmaerosnowDA|C96C48_hybatmDA
Tests:                                 *SUCCESS*
Tests: Completed at Tue Feb  4 10:09:04 CST 2025
Tests: 100% tests passed, 0 tests failed out of 54

[RussTreadon-NOAA]

GDASApp ctests including C48mx500_3DVarAOWCDA and C48mx500_hybAOWCDA pass on Hera, Hercules, and Orion.

Approve.

[AndrewEichmann-NOAA]

looks ok

[shlyaeva]

Here are the changes that I needed to do for the recenter task to succeed with a zero exit status:

diff --git a/parm/soca/berror/soca_ensrecenter.yaml b/parm/soca/berror/soca_ensrecenter.yaml
index b700601..b7d8e20 100644
--- a/parm/soca/berror/soca_ensrecenter.yaml
+++ b/parm/soca/berror/soca_ensrecenter.yaml
@@ -60,5 +60,5 @@ output increment:
   date: '{{ MARINE_WINDOW_END_ISO }}'
   exp: trash
   type: incr
-  output file: 'ocn.recenter.incr.%mem%.nc'
+  output file: 'recenter.incr.%mem%.nc'
   pattern: '%mem%'

and

diff --git a/utils/soca/gdas_ens_handler.h b/utils/soca/gdas_ens_handler.h
index 741055c..dd5fe9c 100644
--- a/utils/soca/gdas_ens_handler.h
+++ b/utils/soca/gdas_ens_handler.h
@@ -169,7 +171,30 @@ namespace gdasapp {
           postProcIncr.setToZero(incr);

           // Save the increments used to initialize the ensemble forecast
-          result = postProcIncr.save(mom6_incr, i+1);
+          result = postProcIncr.save(mom6_incr, i+1, {"ocn"});

It fixes two issues: missing ice file (that we don't use in this usecase anyway) and wrong ocean file name.
I didn't test beyond this task.