Plans for BLOM/iHAMOCC for NorESM3.0

[TomasTorsvik]

Tentative plans for default (and optional) settings for BLOM/iHAMOCC towards NorESM3.0

When shall MCT support phase out?
Due to changes in the timing of releases of NorESM2.3 and NorESM2.5, fixes in the hybrid coordinates would be good to enter a tag 1.10 of Blom that still supports MCT and could still be used together with NorESM2.3. See also: #439

General

• Maintain, further develop, and document tools for flexible creation of new setups including grids, initial conditions,
  and boundary conditions
• Make pre-defined sets of output size, e.g. blom_outvol = small/medium/large simplified output choices for BLOM #252
  Started for iHAMOCC Output switches for output volumes of iHAMOCC #478 Add a user namelist for reduced output settings #498

BLOM-specific

NorESM 2.3+ (BLOM tag v1.10)

• fix negative tracers in hybrid coordinates due to virtual freshwater fluxes

NorESM 3.0 and beyond

• only mks (default) (remove cgs dependencies) Remove CGS units support #515
• Add new advection scheme (Compatible Piecewise Parbolic Method, CPPM) Compatible Piecewise Parabolic Method (CPPM) added as option for layer thickness and tracer transport #531
• Test all model configurations with NUOPC coupler
• Apply existing higher order vertical reconstruction for the remapping of diagnostic fields on to reference depths (hybrid coord. only, reconstruction does not exist for isopycnic coord.)

iHAMOCC-specific

Default settings:

• Turn on extended N-cycle
• Remove separation between euphotic zone and deep ocean processes
• Turn on new dust scheme
• Turn on M4AGO?
• Turn on DOM?

NorESM 2.3+ (BLOM tag 1.10)

• Reconsider output for Hamocc kwrbioz off - enable remin and primary prod throughout full water column #367 & to make it default (+ remove the old split) - see Removing process separation between euphotic and deep ocean regions in iHAMOCC? #399
• Cold start of preformed (silicate) tracer - handle all preformed tracers equally - refer to Initialization of preformed tracers #356 for the provision of climatologies of preformed tracers. Clean-up of restart option for preformed silicate (currently checks, if prefsilicate is present, which could be removed)
• Write output for limiting nutrient Nutrient limitation output #539
• Move towards reading N-dep files - only NOy and NHx reading (abondon ndep) Abandon usage of ndep in input files, always read NHx and NOy #492
• Potentially re-introduce N-cycle coupling through MCT Fix N-cycle coupling #506
• Sediment enhancement, optional; ideally would enter tag 1.8 Introducing optional sediment quality-based remineralization #502
• iHAMOCC restart remapping tool - initial version available through Initial version for iHAMOCC restart remapping #476
• Input files for iHAMOCC (N-depsoition, dust & riverine inputs) can be run with MCT (Provide new N-deposition files and used remapping sh-script #486, Add iHAMOCC support for dustflux on tnx0.5 #490, Add river input file for tn0.5v1 grid #505)
• Input files for new dust deposition Option to use dust/iron deposition data from Myriokefalitakis et al. 2018 (doi:10.5194/bg-15-6659-2018) #575 - still requires tnx0.5+ support

NorESM3.0 (envisaged BLOM tag 1.11 or 2.0) - due in INES II by June 2025

• Write output for rivers and any interpolated input fields introduced by Mariana (or for NorESM 3.0+, new CDEPS stream capability to read in forcing files needed by HAMOCC #405) Add output for dust and soluble iron deposition fluxes #572
• Output for sediment rate tendencies for the default model (potentially also for the water column): aerobic remin., denitrification and "sulfate red." in the sediment (and watercol - in the latter, extNcycle diags could be partially used, when a bit adjusted) - for the sediment: Sediment output #528
• New DOC code (first) as optional DOC code #570 DOC code of Jerry #571
• Boundary conditions: CDEPS versus MCT-compatible, iHAMOCC-based, see
  • iHAMOCC support for tnx0.5v1 grid #376
  • Preparing N-deposition input files for CMIP6 scenarios and historical runs #397
  • Unify N-deposition input reading #398
  • new requirements for dust deposition and riverine input for new stream functionality #404
  • new CDEPS stream capability to read in forcing files needed by HAMOCC #405
  • always obtain ndep from datm or cam using cmeps/cdeps #410
  • ndep forcing file is incorrect for most compsets using DATM using MCT #411 (closed)
    Current status
  • N-deposition in using NUOPC does not use the CDEPS stream capability. It only uses NDEP obtained from either CAM or DATM.
  • The only the technical CDEPS stream functionality in master is for stream_sss and stream_sst. This will serve as templates for further implementation of CDEPS/inline-stream based input of boundary conditions (e.g. fedep) (see new CDEPS stream capability to read in forcing files needed by HAMOCC #405)
  • iHAMOCC can be still run with NUOPC and iHAMOCC-based boundary conditions (i.e. they are read iHAMOCC internally).
    How we proceed: we aim at continuing the work as laid out in new CDEPS stream capability to read in forcing files needed by HAMOCC #405 (useful note: Problem with DMS couple flux when restarting NorESM2.3 with 1/4 BLOM. NorESM#652 (comment))

NorESM3.0+

• Remove/modify compsets for the extended N cycle (in cime_config/config_component.xml), once MCT is depreciated
• Nitrogen cycling interactive coupling through NUOPC

General iHAMOCC ideas

• Potentially enable to include weathering rates
• Offline sediment spinup tool/routine - started a draft: Offline sediment spinup #560

After code freeze

• code clean up