Implement GloSea Operational dignostics

### Are there any linked Issues or Pull Requests?

_No response_

### Brief description

| STASH Code | Section | Item | Name (STASHmaster) | Units | PPFC | CF Standard Name | CF Long Name | Description / Help |
|---|---|---|---|---|---|---|---|---|
| m01s00i385 | 0 | 385 | SNOW LAYER DENSITY ON TILES | kg m⁻³ | 0 | — | — | JULES snow scheme variable (multi-layer scheme only). |
| m01s01i231 | 1 | 231 | DIFFUSE SW FLUX ON LEVELS | W m⁻² | 0 | — | — | Diffuse (scattered) component of the downwards SW flux, on radiation (rho) levels. 3-D field available on radiation timesteps. Domain profile: level 1 = surface, levels 2–ATMOS_TOP = rho levels 2–ATMOS_TOP, level ATMOS_TOP+1 = top of atmosphere. |
| m01s03i509 | 3 | 509 | Seaice sublimation (cats) | kg m⁻² s⁻¹ | 1528 | — | — | Sea-ice sublimation flux on thickness categories (pseudo levels). Same diagnostic as prognostic STASH (1,0,182). Only available when sea-ice thickness categories are used fully in surface exchange. |
| m01s08i258 | 8 | 258 | ACCUM SURFACE RUNOFF RATE | kg m⁻² s⁻¹ | 1532 | `surface_runoff_flux` | — | Gridbox-mean accumulated surface runoff rate, accumulated over the global river-routing timestep (typically 1 day / 86400 s). Diagnostic on atmosphere grid; only available when the global river-routing scheme is used. Prognostic equivalent: STASH (0,155). |
| m01s26i001 | 26 | 1 | RIVER WATER STORAGE | kg | 1902 | — | river water storage | River water storage in a river-routing gridbox, from the global river-routing scheme, on the river-routing grid. Requires a time profile matching the river-routing timestep with an appropriate offset. |
| m01s26i002 | 26 | 2 | GRIDBOX OUTFLOW | kg s⁻¹ | 1904 | — | gridbox outflow | River water outflow from a river-routing gridbox, from the global river-routing scheme, on the river-routing grid. Requires a time profile matching the river-routing timestep with an appropriate offset. |
| m01s26i004 | 26 | 4 | RIVER OUTFLOW (ATMOS GRID) | kg m⁻² s⁻¹ | 1901 | `water_flux_into_sea_water_from_rivers` | — | River water outflow from coastal points, from the global river-routing scheme, remapped onto the atmosphere grid for coupling freshwater fluxes to the ocean. Requires a time profile matching the river-routing timestep with an appropriate offset. |
| m01s30i310 | 30 | 310 | RESIDUAL MN MERID. CIRC. VSTARBAR | m s⁻¹ | 1075 | `northward_transformed_eulerian_mean_air_velocity` | — | Residual mean meridional circulation (northward, v*). Used as the eddy part in the transformed Eulerian mean equations, defined such that dw*/dz + dv*/dy = 0. |
| m01s30i311 | 30 | 311 | RESIDUAL MN MERID. CIRC. WSTARBAR | m s⁻¹ | 1076 | `upward_transformed_eulerian_mean_air_velocity` | — | Residual mean meridional circulation (upward, w*). Used as the eddy part in the transformed Eulerian mean equations, defined such that dw*/dz + dv*/dy = 0. |
| m01s30i312 | 30 | 312 | ELIASSEN-PALM FLUX (MERID. COMPNT) | — | 1077 | `northward_eliassen_palm_flux_in_air` | — | Meridional (northward) component of the Eliassen-Palm flux — involves northward eddy fluxes of heat and angular momentum. |

### Further details of the issue.

_No response_