Ignore:
Timestamp:
Jul 17, 2018 10:14:30 PM (3 years ago)
Author:
maronga
Message:

modified calculation of surface resistance in land surface scheme

File:
1 edited

Legend:

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

    r3142 r3143  
    2525! -----------------
    2626! $Id$
     27! Modified calculation of the surface resistance r_s
     28!
     29! 3142 2018-07-17 15:27:45Z suehring
    2730! Minor bugfix for last commit.
    2831!
     
    21512154
    21522155!
    2153 !--    Calculate the true surface resistance
     2156!--    Calculate the true surface resistance. ABS is used here to avoid negative
     2157!--    values that can occur for very small fluxes due to the artifical addition
     2158!--    of 1.0E-20.
    21542159       IF ( .NOT.  humidity )  THEN
    21552160          surf%r_s(m) = 1.0E10_wp
    21562161       ELSE
    2157           surf%r_s(m) = - rho_lv * ( surf%qv1(m) - q_s + dq_s_dt               &
    2158                           * surf_t_surface%var_1d(m) - dq_s_dt *               &
    2159                             surf_t_surface_p%var_1d(m) ) /                     &
    2160                             (surf%qsws(m) + 1.0E-20)  - surf%r_a(m)
     2162          surf%r_s(m) = ABS(rho_lv / (f_qsws + 1.0E-20_wp) - surf%r_a(m))
    21612163       ENDIF
    2162 !
    2163 !--    Limit resistance resistance
    2164        surf%r_s(m) = MAX( surf%r_s(m), 1.0_wp )
    21652164!
    21662165!--    Calculate change in liquid water reservoir due to dew fall or
Note: See TracChangeset for help on using the changeset viewer.