[MITgcm-support] KPP parameter

Martin Losch Martin.Losch at awi.de
Thu Nov 19 02:50:51 EST 2009


Hi Nicolas,

the parameter BVDcon in table 6.5. does not exist in the code,  
probably a typo, it should be BVScon and the description is not  
accurate: It's only a lower limit of the N2 in the computation of  
mixing due to convection below the mixed layer. I have not been able  
to find the corresponding description in the original KPP paper (Large  
etal 1994, Rev Geophys. 32(4)), but it's part of the computations of  
vertical mixing coefficients from shear instability and (I guess)  
static instability (S/R Ri_iwmix), for shear instabilities it's  
described in eq(25-28), for buoyancy it's analogous. Here's the code  
(diffus(i,ki,2) is the local rho0*N2):
> c  evaluate f of Brunt-Vaisala squared for convection, store in fcon
>
>             Rig   = max ( diffus(i,ki,2) , BVSQcon )
>             ratio = min ( (BVSQcon - Rig) / BVSQcon, c1 )
>             fcon  = c1 - ratio * ratio
>             fcon  = fcon * fcon * fcon
and then diffusivity=fcon*factor. I don't understand, why it has to be  
done for N2/(dV/dz)^2 and for N2, but maybe someone else can explain  
that (again, it's not documented in the original paper).

In general you should not have to change these parameters at all. KPP  
usually works fine with the default parameters,

Martin


On Nov 18, 2009, at 10:22 PM, Nicolas Grisouard wrote:

> Hi everyone,
>
> I am moving from lab-scale simulations (10x1 m^2) to ocean-scale  
> ones (400x4 km^2), still using non-hydrostatic simulations and I am  
> starting my long journey in the world of sub-grid scale  
> parametrizations. I have a first question about the KPP scheme.  
> According to the manual (Table 6.5), a parameter to set is the Brunt- 
> Vaisala frequency squared. The default value is negative, does that  
> correspond to a convective situation or is it just a matter of  
> convention? I have a strongly inhomogeneous Brunt-Vaisala profile  
> (presence of a thermocline), any idea about what to set?
>
> This is all... for the moment. ;)
>
> Thanks,
> Nicolas
>
> -- 
> Nicolas Grisouard
> PhD Student - ERES research group
> Laboratoire des Ecoulements Geophysiques et Industriels
> BP53
> GRENOBLE CEDEX 9 FRANCE
> tel : +33 (0) 476 825 037 - fax : +33 (0) 476 827 022
> http://nicolas.grisouard.free.fr
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list