Changeset 3868 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Apr 8, 2019 11:52:36 AM (6 years ago)

Author:

suehring

Message:

More strict limitation of roughness length when it is in the order of the vertical grid spacing

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! More strict limitation of roughness length when it is in the order of the 
+! vertical grid spacing 
+! 
+! 3856 2019-04-03 11:06:59Z suehring
 ! Bugfix in lsm_init in case no surface-fractions are provided
 ! 
@@ -4879 +4883 @@
 !--    necessary. 
        DO  m = 1, surf_lsm_h%ns
-          IF ( surf_lsm_h%z0(m) >= surf_lsm_h%z_mo(m) )  THEN
+          IF ( surf_lsm_h%z0(m) > 0.5_wp * surf_lsm_h%z_mo(m) )  THEN
           
-             surf_lsm_h%z0(m) = 0.9_wp * surf_lsm_h%z_mo(m)
+             surf_lsm_h%z0(m) = 0.5_wp * surf_lsm_h%z_mo(m)
              
              WRITE( message_string, * ) 'z0 exceeds surface-layer height ' //  &
@@ -4890 +4894 @@
                             0, 0, 0, 6, 0 )
           ENDIF
-          IF ( surf_lsm_h%z0h(m) >= surf_lsm_h%z_mo(m) )  THEN
+          IF ( surf_lsm_h%z0h(m) > 0.5_wp * surf_lsm_h%z_mo(m) )  THEN
           
-             surf_lsm_h%z0h(m) = 0.9_wp * surf_lsm_h%z_mo(m)
-             surf_lsm_h%z0q(m) = 0.9_wp * surf_lsm_h%z_mo(m)
+             surf_lsm_h%z0h(m) = 0.5_wp * surf_lsm_h%z_mo(m)
+             surf_lsm_h%z0q(m) = 0.5_wp * surf_lsm_h%z_mo(m)
              
              WRITE( message_string, * ) 'z0h exceeds surface-layer height ' // &
@@ -4906 +4910 @@
        DO  l = 0, 3
           DO  m = 1, surf_lsm_v(l)%ns
-             IF ( surf_lsm_v(l)%z0(m) >= surf_lsm_v(l)%z_mo(m) )  THEN
+             IF ( surf_lsm_v(l)%z0(m) > 0.5_wp * surf_lsm_v(l)%z_mo(m) )  THEN
           
-                surf_lsm_v(l)%z0(m) = 0.9_wp * surf_lsm_v(l)%z_mo(m)
+                surf_lsm_v(l)%z0(m) = 0.5_wp * surf_lsm_v(l)%z_mo(m)
              
                 WRITE( message_string, * ) 'z0 exceeds surface-layer height '//&
@@ -4918 +4922 @@
                             0, 0, 0, 6, 0 )
              ENDIF
-             IF ( surf_lsm_v(l)%z0h(m) >= surf_lsm_v(l)%z_mo(m) )  THEN
+             IF ( surf_lsm_v(l)%z0h(m) > 0.5_wp * surf_lsm_v(l)%z_mo(m) )  THEN
           
-                surf_lsm_v(l)%z0h(m) = 0.9_wp * surf_lsm_v(l)%z_mo(m)
-                surf_lsm_v(l)%z0q(m) = 0.9_wp * surf_lsm_v(l)%z_mo(m)
+                surf_lsm_v(l)%z0h(m) = 0.5_wp * surf_lsm_v(l)%z_mo(m)
+                surf_lsm_v(l)%z0q(m) = 0.5_wp * surf_lsm_v(l)%z_mo(m)
              
                 WRITE( message_string, * ) 'z0h exceeds surface-layer height '//&