Modern diag manager breaks ocean component with lack of scalar_axis dimension

[GFDL-Eric]

The ocean component for CM4 (the default model configuration for gfdlvitals) relies on the util/extract_ocean_scalar.mom6 routine to grab data from the ocean_scalar_annual.nc file.

The line here:

gfdlvitals/gfdlvitals/util/extract_ocean_scalar.py

Line 28 in 8623f2d

if len(fdata.variables[varname].shape) == 2:

assumes that the output from ocean_scalar_annual.nc will have 2 dimensions (implicitly, time and scalar_axis). Now that scalar_axis is no longer produced for the modern diag manager, the util/extract_ocean_scalar.mom6 routine needs a way to handle the one dimensional output.