- Timestamp:
- Dec 13, 2019 1:23:23 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/land_surface_model_mod.f90
r4321 r4338 25 25 ! ----------------- 26 26 ! $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 27 31 ! Initialization of relative surface fractions revised 28 32 ! … … 7230 7234 !-- wind speed, the sea surface becomes aerodynamically smooth and 7231 7235 !-- 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. 7233 7238 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 ) ) & 7235 7240 + ( alpha_ch * surf_lsm_h%us(m)**2 / g ) 7236 7241 7237 7242 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 ) 7239 7244 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 ) 7241 7246 7242 7247
Note: See TracChangeset
for help on using the changeset viewer.