anchor.diff:          character (len=256) :: name, boundary(30), alb_scheme
           character (len=3)   :: method
           integer             :: nx, n_ksub
  -        double precision    :: ceff,albr,albl,alb_smax,alb_smin,hcrit,rcrit,amp,csh,tmin,tmax,tstic,clh
  +        double precision    :: ceff,albr,albl,alb_smax,alb_smin,hcrit,rcrit,amp,csh,&
  +                               tmin,tmax,tstic,clh, shf_enh, lhf_enh
           double precision    :: pdd_sigma, pdd_b
           double precision    :: itm_cc, itm_t
           double precision    :: tau_a, tau_f, w_crit, mcrit
  
positive.code:             tsurf_new = tsurf
              do ksub = 1, n_ksub
  
                  ! Update the 2m temperature too 
                  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
  
                  ! 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,cap,shf)
                  end if
  
                  ! bulk formulation of latent heat flux (W/m^2), only accounts for sublimation/deposition,
  
negative.code:         par%tmid       = tmid
  
          par%n_ksub     = n_ksub
          par%shf_enh    = shf_enh
          par%lhf_enh    = lhf_enh
  
          return
  
  