Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r4653 r4660  
    2727! -----------------
    2828! $Id$
     29! Fix of wrong limitation of calculation of ueff (avoid r_a=inf) (J. Resler)
     30!
     31! 4653 2020-08-27 08:54:43Z pavelkrc
     32! Remove separate direction indices for urban and land faces (part of RTM v4.0)
     33!
     34! 4630 2020-07-30 14:54:34Z suehring
    2935! - Bugfix in resistance calculation - avoid potential divisions by zero
    3036! - Minor formatting adjustment
     
    80528058!--        to calculation of roughness relative to concrete. Note, wind velocity is limited
    80538059!--        to avoid division by zero. The nominator can become <= 0.0 for values z0 < 3*10E-4.
    8054            ueff        = MAX ( SQRT( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +                   &
    8055                                      ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +                   &
    8056                                      ( ( w(k,j,i) + w(k-1,j,i) ) * 0.5_wp )**2 ),                  &
    8057                                1.0_wp / 4.2_wp                                                     &
    8058                                * ( 4.0_wp / ( surf_usm_v(l)%z0(m) * d_roughness_concrete )         &
    8059                                - 11.8_wp ),                                                        &
    8060                                0.1_wp                                                              &
     8060           ueff        = MAX ( SQRT( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +                  &
     8061                                     ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +                  &
     8062                                     ( ( w(k,j,i) + w(k-1,j,i) ) * 0.5_wp )**2 ),                 &
     8063                               ( ( 4.0_wp + 0.1_wp )                                              &
     8064                                 / ( surf_usm_v(l)%z0(m) * d_roughness_concrete )                 &
     8065                                - 11.8_wp ) / 4.2_wp                                              &
    80618066                              )
    8062 
    8063            surf_usm_v(l)%r_a(m) = rho_cp / ( surf_usm_v(l)%z0(m) * d_roughness_concrete            &
     8067           surf_usm_v(l)%r_a(m) = rho_cp / ( surf_usm_v(l)%z0(m) * d_roughness_concrete           &
    80648068                                  * ( 11.8_wp + 4.2_wp * ueff )  - 4.0_wp  )
    80658069!
Note: See TracChangeset for help on using the changeset viewer.