Large-eddy simulations were employed over half-ice and half-water surfaces, with varying surface temperatures, wind speeds, directions, as to test if the atmospheric interaction with the heterogeneous surface can be predicted via a heterogeneity Richardson number. This dataset was used to determine that surface heat fluxes over ice, water, and the aggregate surface seem to be captured reasonably well by the wind direction and the heterogeneity Richardson number, but the mean wind and turbulent kinetic energy (TKE) profiles were not, suggesting that not only the difference in stability between the two surface, but also the individual stabilities over each surface influence the dynamics. Information on this LES model can be found in these two papers: https://doi.org/10.1175/JAS3826.1 and https://doi.org/10.1063/1.1839152 Corresponding paper: Fogarty, J., Bou-Zeid, E. The Atmospheric Boundary Layer Above the Marginal Ice Zone: Scaling, Surface Fluxes, and Secondary Circulations. Boundary-Layer Meteorol (2023). https://doi.org/10.1007/s10546-023-00825-x