From 8b43f540cf0a3beef216350073db5be51f73d68c Mon Sep 17 00:00:00 2001 From: Stuart Daines Date: Sun, 22 Dec 2024 12:13:49 +0000 Subject: [PATCH] Bugfix for cartesian grids netcdf load/save with isotopes - function cartesian_to_internal(grid::CartesianLinearGrid, griddata::AbstractArray) strip out missing that will have been added by internal_to_cartesian --- src/Grids.jl | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/Grids.jl b/src/Grids.jl index fb5fe26..9c99b00 100644 --- a/src/Grids.jl +++ b/src/Grids.jl @@ -745,8 +745,9 @@ end """ function cartesian_to_internal(grid::CartesianLinearGrid, griddata::AbstractArray) size(grid.linear_index) == size(griddata) || error("grid and data size mismatch") - - return griddata[grid.cartesian_index] + + # recreate vector to change type back from eg Union{Missing, Float64} to Float64, where missing will have been added by internal_to_cartesian + return [griddata[i] for i in grid.cartesian_index] end cartesian_to_internal(grid::CartesianArrayGrid, griddata::AbstractArray) = griddata