anchor.diff:                          !end where
                           lhf = subl*cls
                       end where
  -                else where
  +                else where (land_ice_ocean == 1.0_dp)
                   ! over land latent heat flux exchange coeff is equal to sensible heat coeff
                       where (tsurf < t0)
                           esat_sur = ei_sat(tsurf)
  
positive.code: 
              ! bulk formulation of sensible heat flux (W/m^2)
              if (.not. bnd%shf) then
                  where (land_ice_ocean == 2.0_dp) 
                      shf = csh*cap*rhoatm*usurf*(tsurf-tatm)
                  end where
  
negative.code:             ! 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  
                  where (land_ice_ocean == 2.0_dp)
                      where (tsurf < t0)
                          esat_sur = ei_sat(tsurf)
  