Fix rounding issue in output depths

src/Setup/output.jl

@@ -25,8 +25,8 @@ function CryoGridOutput(sol::TSol, ts=sol.t) where {TSol <: ODESolution}
     log = get_log(sol, ts)
     ts_datetime = Dates.epochms2datetime.(round.(ts*1000.0))
     # Helper functions for mapping variables to appropriate DimArrays by grid/shape.
-    withdims(var::Var{name,T,OnGrid{Edges}}, arr, i) where {name,T} = DimArray(arr, (Z(setup.meta[i].grids[varname(var)]u"m"),Ti(ts_datetime)))
-    withdims(var::Var{name,T,OnGrid{Cells}}, arr, i) where {name,T} = DimArray(arr, (Z(setup.meta[i].grids[varname(var)]u"m"),Ti(ts_datetime)))
+    withdims(var::Var{name,T,OnGrid{Edges}}, arr, i) where {name,T} = DimArray(arr, (Z(round.(setup.meta[i].grids[varname(var)], digits=5)u"m"),Ti(ts_datetime)))
+    withdims(var::Var{name,T,OnGrid{Cells}}, arr, i) where {name,T} = DimArray(arr, (Z(round.(setup.meta[i].grids[varname(var)], digits=5)u"m"),Ti(ts_datetime)))
     withdims(var::Var, arr, i) = DimArray(nestedview(arr), (Ti(ts_datetime),))
     layerstates = NamedTuple()
     for (i,node) in enumerate(setup.strat.nodes)