Convert to MKArray when loading from files

moment_kinetics/ext/file_io_netcdf.jl

@@ -261,7 +261,7 @@ function load_variable(file_or_group::NCDataset, name::String)
         if isa(var, Char)
             var = (var == Char(true))
         end
-        return var
+        return MKArray(var)
     catch
         println("An error occured while loading $name")
         rethrow()
@@ -273,7 +273,7 @@ function load_slice(file_or_group::NCDataset, name::String, slices_or_indices...
     # file or a group).
     try
         var = file_or_group[name].var[slices_or_indices...]
-        return var
+        return MKArray(var)
     catch
         println("An error occured while loading $name")
         rethrow()

moment_kinetics/src/load_data.jl

@@ -168,7 +168,7 @@ function load_variable(file_or_group::HDF5.H5DataStore, name::String)
     # This overload deals with cases where fid is an HDF5 `File` or `Group` (`H5DataStore`
     # is the abstract super-type for both
     try
-        return read(file_or_group[name])
+        return MKArray(read(file_or_group[name]))
     catch
         println("An error occured while loading $name")
         rethrow()
@@ -183,7 +183,7 @@ function load_slice(file_or_group::HDF5.H5DataStore, name::String, slices_or_ind
     # This overload deals with cases where fid is an HDF5 `File` or `Group` (`H5DataStore`
     # is the abstract super-type for both
     try
-        return file_or_group[name][slices_or_indices...]
+        return MKArray(file_or_group[name][slices_or_indices...])
     catch
         println("An error occured while loading $name")
         rethrow()
@@ -414,7 +414,7 @@ function load_time_data(fid; printout=false)
 
     group = get_group(first(fid), "dynamic_data")
     time = load_variable(group, "time")
-    restarts_nt = [length(time)]
+    restarts_nt = MKArray([length(time)])
     for f ∈ fid[2:end]
         group = get_group(f, "dynamic_data")
         # Skip first point as this is a duplicate of the last point of the previous