anchor.diff:              ! not for evaporation/condensation (would require estimate of liquid water content)
               subl = 0.0_dp
               evap = 0.0_dp
  +            lhf  = 0.0_dp
               if (.not. bnd%lhf) then  
  -                lhf  = 0.0_dp
                   where (land_ice_ocean == 2.0_dp)
                       where (tsurf < t0)
                           esat_sur = ei_sat(tsurf)
  
positive.code:                         !end where
                          lhf = subl*cls
                      end where
                  else where
                  ! over land latent heat flux exchange coeff is equal to sensible heat coeff
                      where (tsurf < t0)
                          esat_sur = ei_sat(tsurf)
  
negative.code: 
              ! bulk formulation of sensible heat flux (W/m^2)
              if (.not. bnd%shf) then
                  shf = 0.0_dp
                  where (land_ice_ocean == 2.0_dp) 
                      shf = csh*cap*rhoatm*usurf*(tsurf-tatm)
                  end where
  