anchor.diff:          
           ! physical parameters used in the different parameterizations
           double precision :: tmin, tmax, albr, albl, alb_smax, alb_smin, hcrit, rcrit, &
  -                            ceff, tstic, amp, clh, csh, &
  +                            ceff, tstic, amp, clh, csh, shf_enh, lhf_enh, &
                               pdd_sigma, pdd_b, &
                               itm_cc, itm_t, &
                               tau_a, tau_f, w_crit, mcrit, &
                               mmfac, Pmaxfrac, &
  -                            Ts_wt(2), afac,tmid
  +                            Ts_wt(2), afac,tmid 
           double precision :: tsticsub
           integer :: ksub, n_ksub
   
  
positive.code:     double precision, parameter :: clv  = 2.5e6_dp   !< latent heat of condensation [J/kg]
      double precision, parameter :: cap  = 1000.0_dp  !< specific heat capacitiy of air [J/(kg K)]
      double precision, parameter :: rhow = 1000.0_dp  !< density of water [kg/m3]
      double precision, parameter :: hsmax= 5.0_dp    !< maximum snow height [m]
  
  
      ! Define all parameters needed for the surface module
  
negative.code:             tsurf_new = tsurf
              do ksub = 1, n_ksub
  
                  ! Update the 2m temperature
                  if (.not. bnd%t2m) then  
                      where (land_ice_ocean == 2)             ! ice
                          t2m = Ts_wt(2)*tsurf_new + (1.0_dp-Ts_wt(2))*tatm
                      elsewhere
                          t2m = Ts_wt(1)*tsurf_new + (1.0_dp-Ts_wt(1))*tatm
                      end where
                  end if
  
                  ! bulk formulation of sensible heat flux (W/m^2)
                  if (.not. bnd%shf) then
                      shf = 0.0_dp
                      call sensible_heat_flux(tsurf_new, t2m, usurf, rhoatm, csh, shf_enh, cap, shf)
                  end if
  
                  ! bulk formulation of latent heat flux (W/m^2), only accounts for sublimation/deposition,
  