Actions needed for next V4 dataset release

[ifenty]

Issues to resolve before next release (V4r5)

NetCDF Processing

• Prepare grid files FIRST, use when generating other granules
• assign flux or velocity quantities their own 'bounds' variable that provides the lat/lon coordinates for the two grid cell "corners"
• define collection summaries separately, assign keys to collection title or PODAAC shortname or ID
• fold processing of 1D and 3D time-invariant cases into the main NetCDF generation code
• deliver collection summary text to podaac at the end for their CMR
• save each variable's min/max within the granules of each collection when first constructing netcdf and sweep those at the end to get global min/max for collection/variable

Metadata

• descriptions of DFrI_TH and _SLT should mention contributions from GGL
• momentum budget terms
• possibly new sea-ice terms
• new ice-shelf terms
• Define EmPmR in oceqnet metadata

Code Change

• add missing sea-ice diagnosic terms
• add momentum budget tendency terms
• mask sea-ice velocity where there is no sea ice
• add accurate mixed layer depths (use GGL not Kara)

New Datasets

• GM streamfunction 'rotated' to lat-lon
• https://journals.ametsoc.org/view/journals/phoc/28/6/1520-0485_1998_028_1205_moadti_2.0.co_2.xml
• barotropic streamfunctions
• https://xgcm.readthedocs.io/en/latest/xgcm-examples/02_mitgcm.html?highlight=barotropic#Barotropic-Transport-Streamfunction
• basic climatology fields (T, S, rho; U, V, W; SSH; OBP; sea-ice concentration, sea-ice and snow thicknesses; sea-ice velocity; others?)
• momentum budget closure terms
• ocean heat content 1D time series (full depth, 0-700m & 0-2000m; global and basin)
• OBP anomaly, (OBPANOM) remove reference H * rho, possibly at run time.

Offline adjustment

• Melt rate: time-mean and spin-up
• Global mean SSH and OBP

Checking budget closure

• Temperature budget
• Salinity budget
• Momentum
• sea-ice volume
• sea-ice energy

Documentation

• Synopsis
• User guide
• -- add diagram of different air/sea heat fluxes (EXFqnet, oceQnet, TFLUX, etc.)
• Budget

Figures

(cost plots generally include V4r4 [using V4r5 data and error], V4r5 iteration 0, and V4r5)

• Cost bar plot
• Global mean SSH and OBP (against observation-based estimate)
• Total melt rate time-series (against observation-based estimate)
• Cost vs time
  • LSC SSH
  • OBP
  • Argo
  • SST
  • SSS (Aquarius, SMOS, and SMAP)
  • Sea-ice concentration
• Cost vs space
  • MDT
  • LSC SSH
  • OBP
  • Argo
  • SST
  • SSS (Aquarius, SMOS, and SMAP)
  • Sea-ice concentration
  • Climatology TS
  • Melt rate (Only V4r5 iteration 0 and final V4r5 release)
• gcmfaces plots

Other

• River Runoff and Iceberg Runoff should be in forcing directories, not init directories.
• Geothermal flux should also be in forcing directories, not init directories.

[ifenty]

DFrI_TH:

*. The vertical flux 3rd component from GM-Redi is due to Kwz tensor component
and is easy to process implicitly (vertical flux function of vertical gradient)
like any vertical diffusivity contribution; this is the case when implicitDiffusion=.TRUE.

*. with useGMRedi= True, the 3rd component (from Kwz) is added to "DFrE_???" when using implicitDiffusion=.True. (implictDiffusion = True with V4r4)

DFrE_TH:

• For any tracer, when using pkg/gmredi, it adds to vertical fluxes
  two diffusive components that are proportional to horizontal tracer gradient
  (Kwx & Kwy components of the GM-Redi tensor) and these components are always
  evaluated explicitly (i.e., forward in time)

• only be non-zero if pkg/gmredi activated. "

• with useGMRedi=TRUE, the first 2 components (from Kwx & Kwy) are accounted for in diagnostics "DFrE_???";

• These terms (Kwx & Kwy) are due to the Redi tensor but are increased by some GM
  contributions when using the Skew-flux form (GM_AdvForm = .FALSE., i.e, the default).

• with useGMRedi= True, the 3rd component (from Kwz) is added to "DFrE_???" when using implicitDiffusion=.False. (implictDiffusion = True with V4r4)