Skip to content

Saving a subset of the COSMO REA dataset is giving a dataset that is not openable again #54

@felixcremer

Description

@felixcremer

I saved a subset of the COSMO-REA dataset and tried to reopen it and I get the following error:

julia> c = Cube("test/data/cosmo_rea_subset.zarr/")
name = "hurs"
(ds[name]).attrs["_ARRAY_DIMENSIONS"] = Any["Time", "rlat", "rlon"]
name = "vertices_longitude"
(ds[name]).attrs["_ARRAY_DIMENSIONS"] = Any["rlat", "rlon", "vertices"]
name = "pr"
(ds[name]).attrs["_ARRAY_DIMENSIONS"] = Any["Time", "rlat", "rlon"]
name = "vertices_latitude"
(ds[name]).attrs["_ARRAY_DIMENSIONS"] = Any["rlat", "rlon", "vertices"]
name = "rotated_latitude_longitude"
(ds[name]).attrs["_ARRAY_DIMENSIONS"] = Dict{String, Any}()
ERROR: MethodError: no method matching reverse(::Dict{String, Any})
The function `reverse` exists, but no method is defined for this combination of argument types.

Closest candidates are:
  reverse(::DimensionalData.Dimensions.Lookups.ForwardOrdered)
   @ DimensionalData ~/Documents/GeoZarr.jl/dev/DimensionalData/src/Lookups/methods.jl:19
  reverse(::Test.GenericString)
   @ Test ~/.julia/juliaup/julia-1.12.6+0.x64.linux.gnu/share/julia/stdlib/v1.12/Test/src/Test.jl:2240
  reverse(::DimensionalData.Dimensions.Lookups.Unordered)
   @ DimensionalData ~/Documents/GeoZarr.jl/dev/DimensionalData/src/Lookups/methods.jl:20
  ...

Stacktrace:
  [1] get_var_dims(ds::ZarrExt.ZarrDataset, name::String)
    @ ZarrExt ~/Documents/GeoZarr.jl/dev/YAXArrayBase/ext/ZarrExt.jl:33
  [2] (::YAXArrays.Datasets.var"#60#61"{ZarrExt.ZarrDataset, Vector{String}, Set{Tuple{String, Int64, Int64}}})(k::String)
    @ YAXArrays.Datasets ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:263
  [3] foreach
    @ ./abstractarray.jl:3188 [inlined]
  [4] collectdims(g::ZarrExt.ZarrDataset; force_datetime::Bool)
    @ YAXArrays.Datasets ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:262
  [5] collectdims
    @ ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:259 [inlined]
  [6] (::YAXArrays.Datasets.var"#93#94"{Bool, Tuple{}})(g::ZarrExt.ZarrDataset)
    @ YAXArrays.Datasets ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:467
  [7] open_dataset_handle(f::YAXArrays.Datasets.var"#93#94"{Bool, Tuple{}}, ds::ZarrExt.ZarrDataset)
    @ YAXArrayBase ~/Documents/GeoZarr.jl/dev/YAXArrayBase/src/datasets/datasetinterface.jl:23
  [8] open_dataset(g::String; skip_keys::Tuple{}, driver::Symbol, force_datetime::Bool, path::Nothing)
    @ YAXArrays.Datasets ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:465
  [9] open_dataset
    @ ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:462 [inlined]
 [10] Cube(s::String)
    @ YAXArrays.Datasets ~/Documents/GeoZarr.jl/dev/YAXArrays/src/DatasetAPI/Datasets.jl:1003

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions