The spatial variability in the land surface properties is often not captured by the resolution of land surface models. To overcome this limitation, most models subdivide the grid cells into fractions with homogeneous characteristics, for which the land processes are calculated separately. In reality, the fractions interact via the lateral exchange of water and heat, and the present manuscript details an approach to include these fluxes in the land component of the ICON modeling framework.