Ignore:
Timestamp:
Jun 1, 2007 3:25:22 PM (17 years ago)
Author:
raasch
Message:

preliminary uncomplete changes for ocean version

File:
1 edited

Legend:

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

    r77 r94  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Adjustment of mixing length calculation for the ocean version.
     7! This is also a bugfix, because the height above the topography is now
     8! used instead of the height above level k=0.
    79!
    810! Former revisions:
     
    114116             ENDIF
    115117             IF ( adjust_mixing_length  .AND.  prandtl_layer )  THEN
    116                 l(k)  = MIN( l(k),  kappa * zu(k) / phi_m )
    117                 ll(k) = MIN( ll(k), kappa * zu(k) / phi_m )
     118                l(k)  = MIN( l(k),  kappa * &
     119                                    ( zu(k) - zw(nzb_s_inner(j,i)) ) / phi_m )
     120                ll(k) = MIN( ll(k), kappa * &
     121                                    ( zu(k) - zw(nzb_s_inner(j,i)) ) / phi_m )
    118122             ENDIF
    119123
Note: See TracChangeset for help on using the changeset viewer.