gross checks

The intent here is to have tools which can do quick gross checking of model components. Ideally all components, but for now it is just the ice (CICE) and ocean (MOM6), which have output in NetCDF4. The scripts are fundamentally agnostic as to source, but the interface to reading the data is not as variables have different names between the two (TLAT vs. latitude, for instance).

The README should give sufficient information. There is also a sbatch-submittable script (from Hera) trial.csh which gives an example of calling and running the gross checking scripts.

Files are: [gross|wholesale]_[ice|ocean].py -- python3 scripts specialized to check either a single forecast (gross) or an entire suite of 168 forecasts (wholesale) from either the ice or ocean.

Whenever an extreme is exceeded, the i,j coordinates (python-convention -- 0 to N-1 rather than 1 to N) and lat-long of the point are given, as well as its value and the allowed extremum.

In addition to absolute limits, there are also checks for the low being sufficiently low and the high being sufficiently high. While sea ice thickness won't be > 1.5 meters everywhere, it should be this thick somewhere. While the ocean isn't colder than 1 C everywhere, some points should be this cold or colder.

[icesubset|ocean].[extremes|vhigh] -- dictionary with bounds for the standard icesubset or ocean surface parameters (3d to be added)

• .extremes files give bounds which should never be exceeded, either because they're impossible (negative precipitation rates, speeds > light, salinity > salt saturation, ...) or because they're outside what's been observed (e.g. ocean speeds > 2.5 m/s).
• .vhigh gives bounds which are very high -- so extreme that they're exceeded only on the order of 1 grid point in a million. These can be quite different. For example, the limit for ice u component of velocity is 1.75 m/s for .extreme, but only 0.75 m/s for .vhigh
• ... more as desired. The intention is that there be a progression of bounds checks with narrowing tolerances.
  • secondary is to add regionalized exceptions. e.g., the North Equatorial Currents can exceed 1.5 m/s, as may the major western boundary currents. But outside those modest areas, such speeds are rare or impossible.