Changeset 3135

Timestamp:

Jul 16, 2018 2:43:04 PM (6 years ago)

Author:

gronemeier

Message:

Bugfix: add +1 to c_3 according to Sogachev et al. (2012)

File:

• palm/trunk/SOURCE/model_1d_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/model_1d_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: add +1 to c_3 according to Sogachev et al. (2012)
+! 
+! 3126 2018-07-13 15:13:54Z gronemeier
 ! Bugfix: define c_0 = 0.03^0.25 = 0.416179145
 ! 
@@ -532 +535 @@
                                               * l1d(k) / lambda
                    ENDIF
-                   c_3 = ( c_1 - c_2 ) * alpha_buoyancy
+                   c_3 = ( c_1 - c_2 ) * alpha_buoyancy + 1.0_wp
                    te_diss(k) = ( km1d(k) *                                    &
                                   ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2  &