add error when using tuples as fields and writing ASAGI

Author: boriskaus

src/IO_ASAGI.jl

@@ -88,12 +88,19 @@ function fields_to_namedtuple(fields::NamedTuple, selected_fields)
     ndim    =   length(size(fields[1]))
 
     s2      =   NamedTuple{names}(zeros(nfield))
-    
+
+    # check that they are all DiskArrays
+    for ifield in 1:nfield
+        if !isa(getproperty(fields,names[ifield]), Array)
+            @show typeof(getproperty(fields,names[ifield]))
+            error("Field $(names[ifield]) is not an Array but instead a $(typeof(getproperty(fields,names[ifield]))); only Arrays are supported")
+        end
+    end
+
     material = Array{typeof(s2),ndim}(undef, size(fields[1]))
     for I in eachindex(material)
         data_local = []
         for ifield in 1:nfield
-            
            push!(data_local,getproperty(fields,names[ifield])[I])
 
         end
@@ -167,6 +174,7 @@ function read_ASAGI(fname_asagi::String)
     read_fields_data = ()
     for ifield=1:numfields
         data_1 = zeros(types[ifield],nx,ny,nz)
+
         for I in CartesianIndices(data)
             loc = data[I]
             data_1[I] = getproperty(loc, cnames[ifield])

test/test_ASAGI_IO.jl

@@ -27,7 +27,8 @@ fname_asagi = write_ASAGI("test", Data, (:Sxx,))
 Data_ASAGI2 = read_ASAGI(fname_asagi)
 @test sum(Data_ASAGI2.fields.Sxx - Data.fields.Sxx) == 0
 
-fname_asagi = write_ASAGI("test", Data_tuple)
+# test that it errors if we use a tuple with non-scalar fields
+@test_throws "Field Stress is not an Array but instead a NTuple{9, Array{Float64, 3}}; only Arrays are supported" write_ASAGI("test", Data_tuple)
 
 # Cleanup
 foreach(rm, filter(endswith(".nc"), readdir()))