def METHOD_NAME():
    land_sea_mask = xr.DataArray(
        np.array([0.0, 1.0, 2.0]), dims=["x"], attrs={"units": None}
    )
    reference_sfc_temp = xr.DataArray(
        np.array([1.0, 1.0, 1.0]), dims=["x"], attrs={"units": "degK"}
    )
    model_sfc_temp = xr.DataArray(
        np.array([-1.0, -1.0, -1.0]), dims=["x"], attrs={"units": "degK"}
    )
    assert np.allclose(
        _sst_from_reference(reference_sfc_temp, model_sfc_temp, land_sea_mask),
        [1.0, -1.0, -1.0],
    )