anchor.diff:              qcold = 0.0_dp
           end where
   
  -        ! 1) ablation: melt (m/s); potential melt resulting from available melt energy
  +        !> 4. Ablation: melt (m/s); potential melt resulting from available melt energy
           if (.not. bnd%melt) then
               ! potential melt
               now%melt = qmelt/(rhow*clm)
  
positive.code: 
      end subroutine
  
      elemental subroutine sensible_heat_flux(ts,ta,wind,rhoa,csh,enh,cap,shf)
          double precision, intent(in) :: ts, ta, wind, rhoa
          double precision, intent(in) :: csh, cap, enh
          double precision, intent(out) :: shf
          double precision :: coeff
  
          ! enhancement over land
  
negative.code:         end if
      end subroutine diurnal_cycle
  
      !> calculate snow albedo according to 
      elemental subroutine albedo_isba(alb,sf,melt,tstic,tau,tau_a,tau_f,w_crn,mcrit,alb_smin,alb_smax)
  
          double precision, intent(in out) :: alb
  