 In case of an explicit treatment of this term, the stability requirement may necessitate an unacceptable abridgement of the time increment.To avoid this, vertical turbulent diffusion is treated using the second order Crank–Nicolson method.On principle, advective terms can be computed using any suitable advection scheme.