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

Timestamp:

Sep 1, 2020 2:49:39 PM (4 years ago)

Author:

pavelkrc

Message:

Fix of wrong limitation of calculation of ueff in urban_surface_mod (avoid r_a=inf)

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Fix of wrong limitation of calculation of ueff (avoid r_a=inf) (J. Resler)
+!
+! 4653 2020-08-27 08:54:43Z pavelkrc
+! Remove separate direction indices for urban and land faces (part of RTM v4.0)
+!
+! 4630 2020-07-30 14:54:34Z suehring
 ! - Bugfix in resistance calculation - avoid potential divisions by zero
 ! - Minor formatting adjustment
@@ -8052 +8058 @@
 !--        to calculation of roughness relative to concrete. Note, wind velocity is limited
 !--        to avoid division by zero. The nominator can become <= 0.0 for values z0 < 3*10E-4.
-           ueff        = MAX ( SQRT( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +                   &
-                                     ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +                   &
-                                     ( ( w(k,j,i) + w(k-1,j,i) ) * 0.5_wp )**2 ),                  &
-                               1.0_wp / 4.2_wp                                                     &
-                               * ( 4.0_wp / ( surf_usm_v(l)%z0(m) * d_roughness_concrete )         &
-                               - 11.8_wp ),                                                        &
-                               0.1_wp                                                              &
+           ueff        = MAX ( SQRT( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +                  &
+                                     ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +                  &
+                                     ( ( w(k,j,i) + w(k-1,j,i) ) * 0.5_wp )**2 ),                 &
+                               ( ( 4.0_wp + 0.1_wp )                                              &
+                                 / ( surf_usm_v(l)%z0(m) * d_roughness_concrete )                 &
+                                - 11.8_wp ) / 4.2_wp                                              &
                               )
-
-           surf_usm_v(l)%r_a(m) = rho_cp / ( surf_usm_v(l)%z0(m) * d_roughness_concrete            &
+           surf_usm_v(l)%r_a(m) = rho_cp / ( surf_usm_v(l)%z0(m) * d_roughness_concrete           &
                                   * ( 11.8_wp + 4.2_wp * ueff )  - 4.0_wp  )
 !