Merge pull request #23 from boriskaus/pa-precision

Reduce precision for paraview files

Author: aelligp

examples/MTK_GMG_3D_Lanin.jl

@@ -140,9 +140,9 @@ MatParam     = (SetMaterialParams(Name="Air", Phase=0,
 Grid, Arrays, Tracers, Dikes, time_props = MTK_GMG_3D.MTK_GeoParams_3D(MatParam, Num, Dike_params, CartData_input=Data_3D); # start the main code
 
 Data_set3D_out = Data_3D;
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temperature[C]",  Array(Arrays.Tnew ));   # in MPa
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temp",         Array(Arrays.Tnew));
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Phases",       Array(Arrays.Phases));
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "MeltFraction", Array(Arrays.ϕ));
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temperature[C]",  Float32.(Array(Arrays.Tnew )));   # in MPa
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temp",         Float32.(Array(Arrays.Tnew)));
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Phases",       Int32.(Array(Arrays.Phases)));
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "MeltFraction", Float32.(Array(Arrays.ϕ)));
 save_GMG(joinpath(Num.SimName,"Lanin3D_MTK_final"), Data_set3D_out)
 write_paraview(Data_3D, joinpath(Num.SimName,"Lanin3D_MTK_final"))

examples/MTK_GMG_3D_example.jl

@@ -119,8 +119,8 @@ MatParam     = (SetMaterialParams(Name="Air", Phase=0,
 # Call the main code with the specified material parameters
 Grid, Arrays, Tracers, Dikes, time_props = MTK_GMG_3D.MTK_GeoParams_3D(MatParam, Num, Dike_params, CartData_input=Data_3D); # start the main code
 Data_set3D_out = Data_3D;
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temperature[C]",  Array(Arrays.Tnew ));   # in MPa
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temp",         Array(Arrays.Tnew));
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Phases",       Array(Arrays.Phases));
-Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "MeltFraction", Array(Arrays.ϕ));
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temperature[C]",  Float32.(Array(Arrays.Tnew )));   # in MPa
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Temp",         Float32.(Array(Arrays.Tnew)));
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "Phases",       Int32.(Array(Arrays.Phases)));
+Data_set3D_out = MTK_GMG.add_data_CartData(Data_set3D_out, "MeltFraction", Float32.(Array(Arrays.ϕ)));
 save_GMG(joinpath(Num.SimName,"Unzen3D_MTK_final"), Data_set3D_out)

src/MTK_GMG.jl

@@ -294,9 +294,9 @@ function MTK_save_output(Grid::GridData, Arrays::NamedTuple, Tracers::StructArra
                 Data_set3D  =   CartData(X/1e3,Y/1e3,Z/1e3, (Z=Z,))
             end
             # add datasets
-            Data_set3D = add_data_CartData(Data_set3D, "Temp",         Array(Arrays.Tnew));
-            Data_set3D = add_data_CartData(Data_set3D, "Phases",       Array(Arrays.Phases));
-            Data_set3D = add_data_CartData(Data_set3D, "MeltFraction", Array(Arrays.ϕ));
+            Data_set3D = add_data_CartData(Data_set3D, "Temp",         Float32.(Array(Arrays.Tnew)));
+            Data_set3D = add_data_CartData(Data_set3D, "Phases",       Int32.(Array(Arrays.Phases)));
+            Data_set3D = add_data_CartData(Data_set3D, "MeltFraction", Float64.(Array(Arrays.ϕ)));
 
             # Save output to CartData
             Num.pvd  = write_paraview(Data_set3D, name, pvd=Num.pvd,time=Num.time/SecYear/1e3);