anchor.diff:          par%tmid       = tmid
   
           par%n_ksub     = n_ksub
  -        par%shf_enh    = shf_enh
  -        par%lhf_enh    = lhf_enh
   
           ! initialize sub-daily time step tsticsub
           par%tsticsub = par%tstic / dble(par%n_ksub)
  
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:         !! bulk formulation of sensible heat flux (W/m^2)
          if (.not. bnd%shf) then
              now%shf = 0.0_dp
              call sensible_heat_flux(now%tsurf, now%t2m, now%wind, now%rhoa, par%csh, cap, now%shf)
          end if
  
          !> 2. surface_physics::latent_heat_flux \n
  