K-Profile Parametrization

[navidcy]

We might need to implement the K-Profile boundary layer turbulence closure so that we can use it to compare the CATKE performance.

Ideally, we'd like to calibrate the K-Profile parametrisation using OceanTurbulenceParameterEstimation.jl alongside with CATKE closure and then compare how the two perform under various scenarios.

@glwagner, would it be easy to implement KPP, similarly as done in OceanTurb.jl?

cc @rafferrari

[glwagner]

It is possible, but not easy.

I suggest that we implement a more modern turbulence closure if we would like a second one for comparisons with CATKE. Implementing any "production" turbulence closure will be a significant effort, and we would like it to be worth both the time to implement, and the cost of code maintenance into the future. This will also set us up for comparisons with other ocean models like MOM6 or NEMO.

I suggest:

• Harcourt 2015
• ePBL with surface waves (Reichl and Li 2019)
• The OSMOSIS scheme, the details of which I don't think are yet published, but is implemented in CVmix.

[glwagner]

I'm closing this issue because I'm judging that it's not of current, timely relevance to Oceananigans development. If you would like to make it a higher priority or if you think the issue was closed in error please feel free to re-open.