Prognostic continuous isopycnal does not rearrange entire model state

[ashao]

The 'convective adjustment' performed prior to building the continuous isopycnal coordinate only rearranges layer thicknesses, temperatures, and salinities. If a water column is statically unstable, the regridding thus unphysically breaks the relationships between tracers, velocities, and h, T, and S. To fix this, the convective adjustment thus needs to be made aware of all prognostic quantities associated with a layer.

While this is straightforward to fix for quantities at thickness points, I'm not sure how this would be done for quantities at U/V points. For example, consider two columns: J and J+1 where column J needs to be flipped upside down to be statically stable. How should the meridional-velocities V(J) be rearranged?

In general, this is probably better handled by parameterizing the mixing of both momentum and tracers by convection separately from, but before regridding is done.

[Hallberg-NOAA]

@ashao , it is not clear to me what parts of the code you are discussing. All of the changes to the water column due to mixing (convective or otherwise) in diabatic_driver are provided to the tracers and are properly accounted for with the velocities as well. There is no re-sorting or flipping of layers in diabatic_driver. If you are referring to some other part of the code, please be more explicit about what routines or files are being discussed with this issue. Thank you.