@@ -664,7 +664,7 @@ def initialize_zarr_store(root: Any, big_dataset: "Dataset") -> None:
664664 # Create or append to "dates" dataset.
665665 if "dates" not in root :
666666 full_length = len (big_dataset .dates )
667- root .create_dataset ("dates" , data = np .array ([], dtype = "datetime64[s]" ), chunks = (full_length ,))
667+ root .create_dataset ("dates" , data = np .array ([], dtype = "datetime64[s]" ), chunks = (full_length ,), shape = ( 0 ,) )
668668
669669 if "data" not in root :
670670 dims = (1 , len (big_dataset .variables ), ensembles , big_dataset .shape [- 1 ])
@@ -681,12 +681,15 @@ def initialize_zarr_store(root: Any, big_dataset: "Dataset") -> None:
681681 k ,
682682 data = v ,
683683 compressor = None ,
684+ shape = v .shape ,
684685 )
685686
686687 # Create spatial coordinate datasets if missing.
687688 if "latitudes" not in root or "longitudes" not in root :
688- root .create_dataset ("latitudes" , data = big_dataset .latitudes , compressor = None )
689- root .create_dataset ("longitudes" , data = big_dataset .longitudes , compressor = None )
689+ root .create_dataset ("latitudes" , data = big_dataset .latitudes , compressor = None , shape = big_dataset .latitudes .shape )
690+ root .create_dataset (
691+ "longitudes" , data = big_dataset .longitudes , compressor = None , shape = big_dataset .longitudes .shape
692+ )
690693 for k , v in big_dataset .metadata ().items ():
691694 if k not in root .attrs :
692695 root .attrs [k ] = v
0 commit comments