POP's convective mixing scheme is only applied below the boundary layer. I have a branch of the repository (on my personal fork) where I've added a new option to the convective mixing initialization named lnoOBL, with a default value of .true. implying we do not apply in the boundary layer. I don't know if
(a) this is a feature we want in CVMix (alternatively, GCMs can zero out the convective mixing contribution in the OBL), or
(b) if this is the best implementation.
Some shortfallings -- currently, lnoOBL is only checked if we are using Brunt-Vaisala convective mixing with BVsqr_convect >= 0. Once I have a better idea on the desired application / defaults / variable names I'll clean up the code and apply it regardless of convective mixing scheme.
POP's convective mixing scheme is only applied below the boundary layer. I have a branch of the repository (on my personal fork) where I've added a new option to the convective mixing initialization named
lnoOBL
, with a default value of.true.
implying we do not apply in the boundary layer. I don't know if(a) this is a feature we want in CVMix (alternatively, GCMs can zero out the convective mixing contribution in the OBL), or (b) if this is the best implementation.
Some shortfallings -- currently,
lnoOBL
is only checked if we are using Brunt-Vaisala convective mixing withBVsqr_convect
>= 0. Once I have a better idea on the desired application / defaults / variable names I'll clean up the code and apply it regardless of convective mixing scheme.