Reading/writing F-order data w/o compression is broken #1279

from pathlib import Path
from shutil import rmtree

import numpy as np
import zarr


def get_multichanneled_data(dtype: type) -> np.ndarray:
    data: np.ndarray = np.zeros((3, 250, 200, 10), dtype=dtype)
    for h in range(10):
        for i in range(250):
            for j in range(200):
                data[0, i, j, h] = i * 256
                data[1, i, j, h] = j * 256
                data[2, i, j, h] = 100 * 256
    return data


def main():
    ds_path = Path("testoutput") / "test_zarr"
    rmtree(ds_path, ignore_errors=True)

    data = get_multichanneled_data(np.uint16)
    ds = zarr.create(
        shape=data.shape,
        chunks=(3, 32, 32, 32),
        dtype=data.dtype,
        compressor=(None),
        store=zarr.storage.FSStore(url=str(ds_path), mode="a"),
        order="F",
    )

    ds[:, :, :, :] = data

    ds_reopened = zarr.open_array(
        store=zarr.storage.FSStore(url=str(ds_path), mode="r")
    )

    written_data = ds_reopened[:, :, :, :]

    assert data.shape == written_data.shape
    assert np.array_equal(data, written_data)


if __name__ == "__main__":
    main()