Changeset 2321 for palm

Timestamp:

Jul 24, 2017 3:57:07 PM (7 years ago)

Author:

schwenkel

Message:

Bugfix for calculate ol via lookup table

File:

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

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: Correct index in lookup table for Obukhov length 
+! 
+! 2299 2017-06-29 10:14:38Z suehring
 ! Adjusted for allow separate spinups of LSM and atmosphere code
 ! 
@@ -1346 +1349 @@
 !--          correct direction
              li = li_bnd
-             IF ( rib_tab(li) - surf%rib(m) > 1.0_wp )  THEN
-                DO  WHILE ( rib_tab(li-1) - surf%rib(m) > 1.0_wp  .AND.  li > 1 )
+             IF ( rib_tab(li) - surf%rib(m) > 0.0_wp )  THEN
+                DO  WHILE ( rib_tab(li-1) - surf%rib(m) > 0.0_wp  .AND.  li > 0 )
                    li = li-1
                 ENDDO