Changeset 3136 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Jul 16, 2018 2:48:21 PM (6 years ago)

Author:

suehring

Message:

Limit roughness for heat and moisture in case it exceeds the model-assumed surface-layer height; Limit surface resistance to positive values.

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Limit also roughness length for heat and moisture where necessary
+! 
+! 3123 2018-07-12 16:21:53Z suehring
 ! Correct working precision for INTEGER number
 ! 
@@ -4602 +4605 @@
 !
 !--    Check for consistent initialization.
-!--    Check if roughness length exceed surface-layer height and decrease 
-!--    local roughness length where necessary. 
+!--    Check if roughness length for momentum, or heat, exceed surface-layer
+!--    height and decrease local roughness length where necessary. 
        DO  m = 1, surf_usm_h%ns
           IF ( surf_usm_h%z0(m) >= surf_usm_h%z_mo(m) )  THEN
@@ -4616 +4619 @@
                             0, 0, 0, 6, 0 )
           ENDIF
+          IF ( surf_usm_h%z0h(m) >= surf_usm_h%z_mo(m) )  THEN
+          
+             surf_usm_h%z0h(m) = 0.9_wp * surf_usm_h%z_mo(m)
+             surf_usm_h%z0q(m) = 0.9_wp * surf_usm_h%z_mo(m)
+             
+             WRITE( message_string, * ) 'z0h exceeds surface-layer height ' // &
+                            'at horizontal urban surface and is ' //           &
+                            'decreased appropriately at grid point (i,j) = ',  &
+                            surf_usm_h%i(m), surf_usm_h%j(m)
+             CALL message( 'urban_surface_model_mod', 'PA0507',                &
+                            0, 0, 0, 6, 0 )
+          ENDIF          
        ENDDO
        
@@ -4630 +4645 @@
                             surf_usm_v(l)%j(m)+surf_usm_v(l)%joff
                 CALL message( 'urban_surface_model_mod', 'PA0503',             &
+                            0, 0, 0, 6, 0 )
+             ENDIF
+             IF ( surf_usm_v(l)%z0h(m) >= surf_usm_v(l)%z_mo(m) )  THEN
+          
+                surf_usm_v(l)%z0h(m) = 0.9_wp * surf_usm_v(l)%z_mo(m)
+                surf_usm_v(l)%z0q(m) = 0.9_wp * surf_usm_v(l)%z_mo(m)
+             
+                WRITE( message_string, * ) 'z0h exceeds surface-layer height '//&
+                            'at vertical urban surface and is ' //             &
+                            'decreased appropriately at grid point (i,j) = ',  &
+                            surf_usm_v(l)%i(m)+surf_usm_v(l)%ioff,             &
+                            surf_usm_v(l)%j(m)+surf_usm_v(l)%joff
+                CALL message( 'urban_surface_model_mod', 'PA0507',             &
                             0, 0, 0, 6, 0 )
              ENDIF