Fourth minor update to v0.60 (v0.64)
Several small changes in cvmix_kpp, including some that change answers:
-
Updated bulk Richardson computation to use (d-dr) rather than d, where dr is the depth of the center of the surface layer. Since the surface layer depth is assumed to be surf_layer_ext_d, this means the bulk Richardson number is computed using (1-0.5_surf_layer_ext)*d instead of d in the numerator.
-
Change default behavior for computing the unresolved shear if you pass in N or N^2: instead of using N at the deeper interface as the cell center value as well, the default behavior is to average the two interfaces to get the cell center value. The old behavior is preserved by setting lavg_N_or_Nsqr = .false. in cvmix_init_kpp()
-
The way N^2 is averaged has changed; rather than N_cntr = sqrt(avg(N^2[above,below])), we now use N_cntr = avg(sqrt(N^2)[above,below])
-
Added a new parameter, minVtsqr, that can be set in cvmix_init_kpp to enforce a minimum value on the unresolved shear. Default is 1e-10 m^2/s^2
-
Added two new parameters, minOBLdepth and maxOBLdepth, that can be set in cvmix_init_kpp to enforce a minimum and / or maximum value on the OBL depth. maxOBLdepth = 0 => no maximum; by default both parameters are set to 0.