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

Timestamp:

Aug 13, 2018 12:26:14 PM (7 years ago)

Author:

maronga

Message:

set maximum value for aerodynamic resistance over horiztonal surfaces

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Added maximum aerodynamic resistance of 300 for horiztonal surfaces.
+! 
+! 3176 2018-07-26 17:12:48Z suehring
 ! Bugfix, update virtual potential surface temparture, else heat fluxes on 
 ! roofs might become unphysical
@@ -7373 +7376 @@
 !                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
                 
-!--        make sure that the resistance does not drop to zero
+!--        Make sure that the resistance does not drop to zero
            IF ( surf_usm_h%r_a(m)        < 1.0_wp )                            &
                surf_usm_h%r_a(m)        = 1.0_wp
@@ -7380 +7383 @@
            IF ( surf_usm_h%r_a_window(m) < 1.0_wp )                            &
                surf_usm_h%r_a_window(m) = 1.0_wp
+              
+!
+!--        Make sure that the resistacne does not exceed a maxmium value in case
+!--        of zero velocities
+           IF ( surf_usm_h%r_a(m)        > 300.0_wp )                          &
+               surf_usm_h%r_a(m)        = 300.0_wp
+           IF ( surf_usm_h%r_a_green(m)  > 300.0_wp )                          &
+               surf_usm_h%r_a_green(m) = 300.0_wp
+           IF ( surf_usm_h%r_a_window(m) > 300.0_wp )                          &
+               surf_usm_h%r_a_window(m) = 300.0_wp                
+                
                 
 !--        factor for shf_eb
@@ -7579 +7593 @@
 !
 !--          Limit aerodynamic resistance
-             IF ( surf_usm_v(l)%r_a(m) < 1.0_wp )                              &
-                surf_usm_v(l)%r_a(m)        = 1.0_wp          
+             IF ( surf_usm_v(l)%r_a(m) < 1.0_wp )  surf_usm_v(l)%r_a(m) = 1.0_wp   
+             
                            
              f_shf         = rho_cp / surf_usm_v(l)%r_a(m)