Ignore:
Timestamp:
Aug 13, 2018 12:26:14 PM (6 years ago)
Author:
maronga
Message:

set maximum value for aerodynamic resistance over horiztonal surfaces

File:
1 edited

Legend:

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

    r3176 r3196  
    2828! -----------------
    2929! $Id$
     30! Added maximum aerodynamic resistance of 300 for horiztonal surfaces.
     31!
     32! 3176 2018-07-26 17:12:48Z suehring
    3033! Bugfix, update virtual potential surface temparture, else heat fluxes on
    3134! roofs might become unphysical
     
    73737376!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
    73747377               
    7375 !--        make sure that the resistance does not drop to zero
     7378!--        Make sure that the resistance does not drop to zero
    73767379           IF ( surf_usm_h%r_a(m)        < 1.0_wp )                            &
    73777380               surf_usm_h%r_a(m)        = 1.0_wp
     
    73807383           IF ( surf_usm_h%r_a_window(m) < 1.0_wp )                            &
    73817384               surf_usm_h%r_a_window(m) = 1.0_wp
     7385             
     7386!
     7387!--        Make sure that the resistacne does not exceed a maxmium value in case
     7388!--        of zero velocities
     7389           IF ( surf_usm_h%r_a(m)        > 300.0_wp )                          &
     7390               surf_usm_h%r_a(m)        = 300.0_wp
     7391           IF ( surf_usm_h%r_a_green(m)  > 300.0_wp )                          &
     7392               surf_usm_h%r_a_green(m) = 300.0_wp
     7393           IF ( surf_usm_h%r_a_window(m) > 300.0_wp )                          &
     7394               surf_usm_h%r_a_window(m) = 300.0_wp               
     7395               
    73827396               
    73837397!--        factor for shf_eb
     
    75797593!
    75807594!--          Limit aerodynamic resistance
    7581              IF ( surf_usm_v(l)%r_a(m) < 1.0_wp )                              &
    7582                 surf_usm_v(l)%r_a(m)        = 1.0_wp         
     7595             IF ( surf_usm_v(l)%r_a(m) < 1.0_wp )  surf_usm_v(l)%r_a(m) = 1.0_wp   
     7596             
    75837597                           
    75847598             f_shf         = rho_cp / surf_usm_v(l)%r_a(m)
Note: See TracChangeset for help on using the changeset viewer.