Changeset 4338 for palm


Ignore:
Timestamp:
Dec 13, 2019 1:23:23 PM (5 years ago)
Author:
suehring
Message:

To avoid divisions by zero, add security factor in calculation of roughness length over water surfaces.

File:
1 edited

Legend:

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

    r4321 r4338  
    2525! -----------------
    2626! $Id$
     27! To avoid divisions by zero, add security factor in calculation of roughness
     28! length over water surfaces.
     29!
     30! 4321 2019-12-04 10:26:38Z pavelkrc
    2731! Initialization of relative surface fractions revised
    2832!
     
    72307234!--          wind speed, the sea surface becomes aerodynamically smooth and
    72317235!--          the roughness scales with the viscosity. At high wind speed, the
    7232 !--          Charnock relation is used.
     7236!--          Charnock relation is used. Add a security factor of 1E-8 to avoid
     7237!--          divisions by zero.
    72337238             surf_lsm_h%z0(m)  = ( 0.11_wp * molecular_viscosity /             &
    7234                                  surf_lsm_h%us(m) )                            &
     7239                                 ( surf_lsm_h%us(m) + 1E-8_wp ) )              &
    72357240                               + ( alpha_ch * surf_lsm_h%us(m)**2 / g )
    72367241
    72377242             surf_lsm_h%z0h(m) = 0.40_wp * molecular_viscosity /               &
    7238                                  surf_lsm_h%us(m)
     7243                                 ( surf_lsm_h%us(m) + 1E-8_wp )
    72397244             surf_lsm_h%z0q(m) = 0.62_wp * molecular_viscosity /               &
    7240                                  surf_lsm_h%us(m)
     7245                                 ( surf_lsm_h%us(m) + 1E-8_wp )
    72417246
    72427247 
Note: See TracChangeset for help on using the changeset viewer.