we are initialize empty zarrs with shape:
(31, 96, 1000, 1000) and chunks : (31, 96, 50, 50) because we need all data from dim 1 and 2 in each query.
Now we write into region again:

dataset.to_zarr(store,region=region,mode="r+",compute=True,consolidated=Trueifzarr_format<3elseFalse,zarr_format=zarr_format,safe_chunks=safe_chunks,    )

During writing along the first dimension day by day, the whole chunk (aka dataset) needs to be loaded to write into region. This causes a lot of traffic.

We would like to go into s3 object store with the process but we need to decrease the traffic. Do you know how to achive this?

The previous solution used append mode, which causes chunk sizes (1, 96, 50, 50) and very long queries against the dataset.

Best regards
Daniel