anchor.diff:              now%evap = 0.0_dp
               now%lhf  = 0.0_dp
               call latent_heat_flux(now%tsurf, now%wind, now%qq, now%sp, now%rhoa, now%mask, &
  -                                  par%clh, par%lhf_enh, eps, cls, clv, now%lhf, now%subl, now%evap)
  +                                  par%clh, eps, cls, clv, now%lhf, now%subl, now%evap)
           end if
   
           !> 3. surface_physics::longwave_upward \n
  
positive.code:         double precision    :: tmax      !< maximum temperature for which albedo decline becomes effective ("slater") [K]
          double precision    :: tstic     !< time step [s]
          double precision    :: tsticsub  !< sub-time step [s]
          double precision    :: shf_enh   !< sensible heat flux enhancement factor for positive u*(tsurf-t2m) [no unit]
          double precision    :: lhf_enh   !< latent heat flux enhancement factor over land (mask==1) [no unit]
          double precision    :: tau_a     !< dry albedo decline for "isba" albedo scheme [1/day]
          double precision    :: tau_f     !< wet albedo decline for "isba" albedo scheme [1/day]
          double precision    :: w_crit    !< critical liquid water content for "isba" albedo scheme [kg/m2]
  
negative.code:         par%tmid       = tmid
  
          par%n_ksub     = n_ksub
  
          ! initialize sub-daily time step tsticsub
          par%tsticsub = par%tstic / dble(par%n_ksub)
  