Ignore:
Timestamp:
Jul 16, 2018 2:43:04 PM (3 years ago)
Author:
gronemeier
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/model_1d_mod.f90

    r3126 r3135  
    2525! -----------------
    2626! $Id$
     27! Bugfix: add +1 to c_3 according to Sogachev et al. (2012)
     28!
     29! 3126 2018-07-13 15:13:54Z gronemeier
    2730! Bugfix: define c_0 = 0.03^0.25 = 0.416179145
    2831!
     
    532535                                              * l1d(k) / lambda
    533536                   ENDIF
    534                    c_3 = ( c_1 - c_2 ) * alpha_buoyancy
     537                   c_3 = ( c_1 - c_2 ) * alpha_buoyancy + 1.0_wp
    535538                   te_diss(k) = ( km1d(k) *                                    &
    536539                                  ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2  &
Note: See TracChangeset for help on using the changeset viewer.