source: (kappa.f90) [ c_92_rp1 release ]
This module provides algorithmns for vertical mixing coefficients and declares the Type (T_vmix_coefs) which contains the parameters needed for those algorithms.
The main functions are VerticalDiffusivity and VerticalViscosity They are overloaded for 0 through 4 dimensional fields of Richardson number, but they can be extended to allow for more parameters, like N. I tried to make these elemental functions, but I did not like the idea of
The Richardson number that is available is currently simply based on the gradient of buoyancy, which is presumed to be done in terms of potential density. This needs to be improved.