Changeset 3196 for palm/trunk/SOURCE


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

set maximum value for aerodynamic resistance over horiztonal surfaces

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r3161 r3196  
    2525! -----------------
    2626! $Id$
     27! Added maximum aerodynamic resistance of 300.
     28!
     29! 3161 2018-07-23 09:30:10Z maronga
    2730! Increased roughness of asphalt surfaces to account for turbulence production
    2831! by traffic and other obstacles
     
    19151918!
    19161919!--    Make sure that the resistance does not drop to zero for neutral
    1917 !--    stratification.
    1918        IF ( surf%r_a(m) < 1.0_wp )  surf%r_a(m) = 1.0_wp
     1920!--    stratification. Also, set a maximum resistance to avoid the breakdown of
     1921!--    MOST for locations with zero wind speed
     1922       IF ( surf%r_a(m) <   1.0_wp )  surf%r_a(m) =   1.0_wp
     1923       IF ( surf%r_a(m) > 300.0_wp )  surf%r_a(m) = 300.0_wp       
    19191924!
    19201925!--    Second step: calculate canopy resistance r_canopy
  • 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.