anchor.diff:          shf = coeff*cap*rhoa*wind*(ts-ta)
       end subroutine
   
  -    elemental subroutine latent_heat_flux(ts, wind, shum, sp, rhoatm, lmask, clh, enh, eps, cls, clv, lhf, subl, evap)
  +    elemental subroutine latent_heat_flux(ts, wind, shum, sp, rhoatm, mask, clh, enh, eps, cls, clv, lhf, subl, evap)
           double precision, intent(in)  :: ts, shum, sp, rhoatm, wind
           double precision, intent(in)  :: clh, eps, cls, clv, enh
  -        integer, intent(in)  :: lmask
  +        integer, intent(in)  :: mask
           double precision, intent(out) :: lhf, subl, evap
           double precision :: esat_sur, shum_sat, coeff
   
           subl = 0.0_dp
           evap = 0.0_dp
           lhf  = 0.0_dp
  -        if (lmask == 1) then
  +        if (mask == 1) then
               coeff = enh*clh
           else
               coeff = clh
  
positive.code:                     bnd%melt    = .TRUE. 
                  case("refr")
                      bnd%refr    = .TRUE. 
                  case("massbal")
                      bnd%massbal = .TRUE. 
                  case("acc")
                      bnd%acc     = .TRUE. 
                  case("lhf")
  
negative.code:         par%tmin       = tmin           ! minimum albedo-affecting temperature
          par%tmax       = tmax           ! maximum albedo-affecting temperature
          par%tstic      = tstic
          par%alb_scheme = alb_scheme
          par%tau_a      = tau_a
          par%tau_f      = tau_f
          par%w_crit     = w_crit
  