Commit 3db99ed
authored
Marine DA Hotfix for release 1.0 incr qc (#2075)
# Description
Copy/pasted from
[JCSDA-internal/soca/pull/1223](JCSDA-internal/soca#1223):
This PR adds increment quality control (QC) utilities to SOCA to enforce
physical consistency of ocean analysis increments prior to output. The
new functionality applies configurable bounds on temperature and
salinity increments, checks and corrects water column stability to
prevent density inversions, and constrains steric height contributions
to limit unrealistic SSH increments. Density calculations use the UNESCO
1983 equation of state, and the QC can be applied iteratively with
optional smoothing.
The QC is integrated into the ensemble analysis postprocessing workflow
via a new qcIncrement() call in AnalysisPostproc, ensuring that all
analysis increments satisfy physical constraints before being written
out.
This PR adds new functionality; changes in analysis increments are
expected when QC thresholds are active, particularly in regions or
cycles with large or unbalanced increments.
# Companion PRs
- NOAA-EMC/jcb-gdas#232
- NOAA-EMC/soca#3
# Issues
Resolves #
# Automated CI tests to run in Global Workflow
<!-- 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_hybatmsnowDA <!-- JEDI snow cycled DA !-->
- [ ] C96_gcafs_cycled <!-- JEDI aerosol cycled DA !-->
- [ ] C48mx500_3DVarAOWCDA <!-- JEDI low-res marine 3DVar cycled DA !-->
- [ ] C48mx500_hybAOWCDA <!-- JEDI marine hybrid envar cycled DA !-->
- [ ] C96C48_ufsgsi_hybatmDA <!-- JEDI atm Var with GSI EnKF cycled DA
!-->
- [ ] C48_ufsenkf_atmDA <!-- JEDI atm EnKF cycled DA !-->
- [ ] C96C48_hybatmDA <!-- GSI atm cycled DA !-->1 parent b747e5c commit 3db99ed
2 files changed
Lines changed: 2 additions & 2 deletions
- src/mains/AnalysisPostproc.h+119-46
- src/soca/Utils/CMakeLists.txt+3
- src/soca/Utils/incrqc/CMakeLists.txt+9
- src/soca/Utils/incrqc/README.md+205
- src/soca/Utils/incrqc/include/soca_diagb_utils.h+284
- src/soca/Utils/incrqc/include/soca_incr_qc.h+151
- src/soca/Utils/incrqc/include/soca_incr_qc_utils.h+479
- src/soca/Utils/incrqc/include/soca_physlight_utils.h+210
- src/soca/Utils/incrqc/include/soca_utils.h+210
0 commit comments