Changeset 3138

Timestamp:

Jul 17, 2018 8:21:20 AM (6 years ago)

Author:

maronga

Message:

bugfix for sea surfaces with constant_roughness = .F.

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: limit roughness lengths in case of sea surface with constant_roughness
+! = .F.
+! 
+! 3136 2018-07-16 14:48:21Z suehring
 ! Limit also roughness length for heat and moisture where necessary;
 ! Limit surface resistance to positive values
@@ -6899 +6903 @@
 
        USE control_parameters,                                                 &
-           ONLY: g, kappa, molecular_viscosity
+           ONLY: g, kappa, message_string, molecular_viscosity
 
        IMPLICIT NONE
@@ -6952 +6956 @@
                                  surf_lsm_h%us(m)
 
+  
+             IF ( surf_lsm_h%z0(m) >= surf_lsm_h%z_mo(m) )  THEN
+  
+                surf_lsm_h%z0(m) = 0.9_wp * surf_lsm_h%z_mo(m)
+             
+                WRITE( message_string, * ) 'z0 exceeds surface-layer height' //&
+                            ' at horizontal sea surface and is ' //            &
+                            'decreased appropriately at grid point (i,j) = ',  &
+                            surf_lsm_h%i(m), surf_lsm_h%j(m)
+                CALL message( 'land_surface_model_mod', 'PA0508',              &
+                              0, 0, 0, 6, 0 )
+             ENDIF
+  
+             IF ( surf_lsm_h%z0h(m) >= surf_lsm_h%z_mo(m) )  THEN
+  
+                surf_lsm_h%z0h(m) = 0.9_wp * surf_lsm_h%z_mo(m)
+             
+                WRITE( message_string, * ) 'z0h exceeds surface-layer height'//&
+                            ' at horizontal sea surface and is ' //            &
+                            'decreased appropriately at grid point (i,j) = ',  &
+                            surf_lsm_h%i(m), surf_lsm_h%j(m)
+                CALL message( 'land_surface_model_mod', 'PA0508',              &
+                              0, 0, 0, 6, 0 )
+             ENDIF
+             
+             IF ( surf_lsm_h%z0(m) >= surf_lsm_h%z_mo(m) )  THEN
+  
+                surf_lsm_h%z0q(m) = 0.9_wp * surf_lsm_h%z_mo(m)
+             
+                WRITE( message_string, * ) 'z0q exceeds surface-layer height'//&
+                            ' at horizontal sea surface and is ' //            &
+                            'decreased appropriately at grid point (i,j) = ',  &
+                            surf_lsm_h%i(m), surf_lsm_h%j(m)
+                CALL message( 'land_surface_model_mod', 'PA0508',              &
+                              0, 0, 0, 6, 0 )
+             ENDIF
+  
+                                 
           ENDIF
        ENDDO