diff --git a/CMakeLists.txt b/CMakeLists.txt
index f2dcc98f82..8290227be6 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -426,6 +426,11 @@ set(CORE_SOURCE
         src/backend/PatchRecord.cpp
         src/backend/PatchRecordComponent.cpp
         src/backend/Writable.cpp
+        src/backend/scientific_defaults/ScientificDefaults.cpp
+        src/backend/scientific_defaults/ScientificDefaults_auxiliary.cpp
+        src/backend/scientific_defaults/ProcessParsedAttribute.cpp
+        src/backend/scientific_defaults/AttributeReader.cpp
+        src/backend/scientific_defaults/ConfigAttribute.cpp
         src/auxiliary/OneDimensionalBlockSlicer.cpp
         src/helper/list_series.cpp
         src/snapshots/ContainerImpls.cpp
diff --git a/include/openPMD/IO/AbstractIOHandler.hpp b/include/openPMD/IO/AbstractIOHandler.hpp
index 9b7735b5ba..b54b661c78 100644
--- a/include/openPMD/IO/AbstractIOHandler.hpp
+++ b/include/openPMD/IO/AbstractIOHandler.hpp
@@ -26,6 +26,7 @@
 #include "openPMD/IterationEncoding.hpp"
 #include "openPMD/config.hpp"
 #include "openPMD/version.hpp"
+#include <ostream>
 
 #if openPMD_HAVE_MPI
 #include <mpi.h>
@@ -81,6 +82,66 @@ enum class FlushLevel
     CreateOrOpenFiles
 };
 
+std::ostream &operator<<(std::ostream &, FlushLevel);
+
+namespace flush_level
+{
+    inline constexpr auto global_flushpoint(FlushLevel fl)
+    {
+        switch (fl)
+        {
+        case FlushLevel::UserFlush:
+            return true;
+        case FlushLevel::InternalFlush:
+        case FlushLevel::SkeletonOnly:
+        case FlushLevel::CreateOrOpenFiles:
+            return false;
+        }
+        return false; // unreachable
+    }
+    // same as global_flushpoint for now, but we will soon introduce
+    // immediate_flush
+    inline constexpr auto write_datasets(FlushLevel fl)
+    {
+        switch (fl)
+        {
+        case FlushLevel::UserFlush:
+            return true;
+        case FlushLevel::InternalFlush:
+        case FlushLevel::SkeletonOnly:
+        case FlushLevel::CreateOrOpenFiles:
+            return false;
+        }
+        return false; // unreachable
+    }
+    inline constexpr auto write_attributes(FlushLevel fl)
+    {
+        switch (fl)
+        {
+        case FlushLevel::UserFlush:
+        case FlushLevel::InternalFlush:
+            return true;
+        case FlushLevel::SkeletonOnly:
+        case FlushLevel::CreateOrOpenFiles:
+            return false;
+        }
+        return false; // unreachable
+    }
+    inline constexpr auto flush_hierarchy(FlushLevel fl)
+    {
+        switch (fl)
+        {
+        case FlushLevel::UserFlush:
+        case FlushLevel::InternalFlush:
+        case FlushLevel::SkeletonOnly:
+            return true;
+        case FlushLevel::CreateOrOpenFiles:
+            return false;
+        }
+        return false; // unreachable
+    }
+} // namespace flush_level
+
 enum class OpenpmdStandard
 {
     v_1_0_0,
@@ -121,6 +182,7 @@ namespace internal
      * To be used for reading
      */
     FlushParams const defaultFlushParams{};
+    FlushParams const publicFlush{FlushLevel::UserFlush};
 
     struct ParsedFlushParams;
 
diff --git a/include/openPMD/IO/AbstractIOHandlerImpl.hpp b/include/openPMD/IO/AbstractIOHandlerImpl.hpp
index d45ce1bdcc..fdb8af3599 100644
--- a/include/openPMD/IO/AbstractIOHandlerImpl.hpp
+++ b/include/openPMD/IO/AbstractIOHandlerImpl.hpp
@@ -39,7 +39,7 @@ class AbstractIOHandlerImpl
 
     virtual ~AbstractIOHandlerImpl() = default;
 
-    std::future<void> flush();
+    std::future<void> flush(FlushLevel);
 
     /**
      * Close the file corresponding with the writable and release file handles.
diff --git a/include/openPMD/IO/IOTask.hpp b/include/openPMD/IO/IOTask.hpp
index 25e0d6ad54..a8efab571e 100644
--- a/include/openPMD/IO/IOTask.hpp
+++ b/include/openPMD/IO/IOTask.hpp
@@ -35,6 +35,7 @@
 #include <cstddef>
 #include <memory>
 #include <optional>
+#include <ostream>
 #include <string>
 #include <utility>
 #include <variant>
@@ -89,6 +90,8 @@ OPENPMDAPI_EXPORT_ENUM_CLASS(Operation){
 }; // note: if you change the enum members here, please update
    // docs/source/dev/design.rst
 
+std::ostream &operator<<(std::ostream &os, Operation op);
+
 namespace internal
 {
     /*
diff --git a/include/openPMD/IO/JSON/JSONIOHandlerImpl.hpp b/include/openPMD/IO/JSON/JSONIOHandlerImpl.hpp
index 6df0c60ced..3e0758aee2 100644
--- a/include/openPMD/IO/JSON/JSONIOHandlerImpl.hpp
+++ b/include/openPMD/IO/JSON/JSONIOHandlerImpl.hpp
@@ -241,7 +241,7 @@ class JSONIOHandlerImpl : public AbstractIOHandlerImpl
 
     void touch(Writable *, Parameter<Operation::TOUCH> const &) override;
 
-    std::future<void> flush();
+    std::future<void> flush(internal::ParsedFlushParams &params);
 
 private:
 #if openPMD_HAVE_MPI
diff --git a/include/openPMD/Iteration.hpp b/include/openPMD/Iteration.hpp
index e2a4e6b1a5..e281686855 100644
--- a/include/openPMD/Iteration.hpp
+++ b/include/openPMD/Iteration.hpp
@@ -27,6 +27,8 @@
 #include "openPMD/auxiliary/Variant.hpp"
 #include "openPMD/backend/Attributable.hpp"
 #include "openPMD/backend/Container.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <cstdint>
 #include <deque>
@@ -145,13 +147,28 @@ namespace internal
         std::optional<DeferredParseAccess> m_deferredParseAccess{};
     };
 } // namespace internal
+
+class Meshes : public Container<Mesh>
+{
+public:
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+};
+
+class Particles : public Container<ParticleSpecies>
+{
+public:
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+};
+
 /** @brief  Logical compilation of data from one snapshot (e.g. a single
  * simulation cycle).
  *
  * @see
  * https://github.com/openPMD/openPMD-standard/blob/latest/STANDARD.md#required-attributes-for-the-basepath
  */
-class Iteration : public Attributable
+class Iteration
+    : public Attributable
+    , internal::ScientificDefaults
 {
     template <typename T, typename T_key, typename T_container>
     friend class Container;
@@ -164,6 +181,8 @@ class Iteration : public Attributable
     friend class StatefulSnapshotsContainer;
     template <typename>
     friend struct traits::GenerationPolicy;
+    friend class internal::ScientificDefaults;
+    friend class Attributable;
 
 public:
     Iteration(Iteration const &) = default;
@@ -279,8 +298,10 @@ class Iteration : public Attributable
     [[deprecated("This attribute is no longer set by the openPMD-api.")]] bool
     closedByWriter() const;
 
-    Container<Mesh> meshes{};
-    Container<ParticleSpecies> particles{}; // particleSpecies?
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
+    Meshes meshes{};
+    Particles particles{};
 
     virtual ~Iteration() = default;
 
@@ -442,7 +463,7 @@ class Iteration : public Attributable
      *
      * @param w The Writable representing the parent.
      */
-    virtual void linkHierarchy(Writable &w);
+    void linkHierarchy(Writable &w) override;
 
     /**
      * @brief Access an iteration in read mode that has potentially not been
@@ -450,6 +471,9 @@ class Iteration : public Attributable
      *
      */
     void runDeferredParseAccess();
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 }; // Iteration
 
 namespace traits
diff --git a/include/openPMD/Mesh.hpp b/include/openPMD/Mesh.hpp
index d0bf81ddef..06885d4558 100644
--- a/include/openPMD/Mesh.hpp
+++ b/include/openPMD/Mesh.hpp
@@ -23,7 +23,9 @@
 #include "openPMD/UnitDimension.hpp"
 #include "openPMD/backend/Attributable.hpp"
 #include "openPMD/backend/BaseRecord.hpp"
+#include "openPMD/backend/Container.hpp"
 #include "openPMD/backend/MeshRecordComponent.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <ostream>
 #include <string>
@@ -41,6 +43,8 @@ class Mesh : public BaseRecord<MeshRecordComponent>
 {
     friend class Container<Mesh>;
     friend class Iteration;
+    friend class internal::ScientificDefaults;
+    friend class Attributable;
 
 public:
     Mesh(Mesh const &) = default;
@@ -322,12 +326,18 @@ class Mesh : public BaseRecord<MeshRecordComponent>
         typename = std::enable_if_t<std::is_floating_point<T>::value>>
     Mesh &setTimeOffset(T timeOffset);
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
 private:
     Mesh();
 
     void
     flush_impl(std::string const &, internal::FlushParams const &) override;
     void read();
+    auto retrieveDimensionality() const -> uint64_t;
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 }; // Mesh
 
 template <typename T>
diff --git a/include/openPMD/ParticlePatches.hpp b/include/openPMD/ParticlePatches.hpp
index f3c4c0b943..ffe539d5f1 100644
--- a/include/openPMD/ParticlePatches.hpp
+++ b/include/openPMD/ParticlePatches.hpp
@@ -38,6 +38,8 @@ class ParticlePatches : public Container<PatchRecord>
     size_t numPatches() const;
     ~ParticlePatches() override = default;
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
 private:
     ParticlePatches() = default;
     void read();
diff --git a/include/openPMD/ParticleSpecies.hpp b/include/openPMD/ParticleSpecies.hpp
index 7309afddef..9171a2d85e 100644
--- a/include/openPMD/ParticleSpecies.hpp
+++ b/include/openPMD/ParticleSpecies.hpp
@@ -24,23 +24,31 @@
 #include "openPMD/Record.hpp"
 #include "openPMD/backend/Attributable.hpp"
 #include "openPMD/backend/Container.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <string>
 
 namespace openPMD
 {
 
-class ParticleSpecies : public Container<Record>
+class ParticleSpecies
+    : public Container<Record>
+    , internal::ScientificDefaults
 {
     friend class Container<ParticleSpecies>;
     friend class Container<Record>;
     friend class Iteration;
     template <typename T>
     friend T &internal::makeOwning(T &self, Series);
+    friend class internal::ScientificDefaults;
+    friend class Attributable;
 
 public:
     ParticlePatches particlePatches;
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
 private:
     ParticleSpecies();
 
@@ -53,6 +61,9 @@ class ParticleSpecies : public Container<Record>
     {
         return m_containerData;
     }
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 };
 
 namespace traits
diff --git a/include/openPMD/Record.hpp b/include/openPMD/Record.hpp
index 2bb070365f..e7ffff5ac2 100644
--- a/include/openPMD/Record.hpp
+++ b/include/openPMD/Record.hpp
@@ -23,6 +23,7 @@
 #include "openPMD/RecordComponent.hpp"
 #include "openPMD/UnitDimension.hpp"
 #include "openPMD/backend/BaseRecord.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <string>
 #include <type_traits>
@@ -34,6 +35,7 @@ class Record : public BaseRecord<RecordComponent>
     friend class Container<Record>;
     friend class Iteration;
     friend class ParticleSpecies;
+    friend class internal::ScientificDefaults;
 
 public:
     Record(Record const &) = default;
@@ -48,6 +50,8 @@ class Record : public BaseRecord<RecordComponent>
     template <typename T>
     Record &setTimeOffset(T);
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
 private:
     Record();
 
@@ -55,6 +59,9 @@ class Record : public BaseRecord<RecordComponent>
     flush_impl(std::string const &, internal::FlushParams const &) override;
 
     [[nodiscard]] internal::HomogenizeExtents read();
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 }; // Record
 
 template <typename T>
diff --git a/include/openPMD/RecordComponent.hpp b/include/openPMD/RecordComponent.hpp
index d06b4213f4..d01a183679 100644
--- a/include/openPMD/RecordComponent.hpp
+++ b/include/openPMD/RecordComponent.hpp
@@ -27,6 +27,8 @@
 #include "openPMD/auxiliary/UniquePtr.hpp"
 #include "openPMD/backend/Attributable.hpp"
 #include "openPMD/backend/BaseRecordComponent.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 // comment to prevent this include from being moved by clang-format
 #include "openPMD/DatatypeMacros.hpp"
@@ -110,7 +112,9 @@ namespace internal
 template <typename>
 class BaseRecord;
 
-class RecordComponent : public BaseRecordComponent
+class RecordComponent
+    : public BaseRecordComponent
+    , protected internal::ScientificDefaults
 {
     template <typename T, typename T_key, typename T_container>
     friend class Container;
@@ -128,6 +132,8 @@ class RecordComponent : public BaseRecordComponent
     friend class MeshRecordComponent;
     template <typename T>
     friend T &internal::makeOwning(T &self, Series);
+    friend class internal::ScientificDefaults;
+    friend class Attributable;
 
 public:
     enum class Allocation
@@ -482,11 +488,13 @@ class RecordComponent : public BaseRecordComponent
     auto visit(Args &&...args) -> decltype(Visitor::template call<char>(
         std::declval<RecordComponent &>(), std::forward<Args>(args)...));
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
     static constexpr char const *const SCALAR = "\vScalar";
 
 protected:
     void flush(std::string const &, internal::FlushParams const &);
-    void read(bool require_unit_si);
+    void read();
 
 private:
     /**
@@ -534,12 +542,15 @@ OPENPMD_protected
         BaseRecordComponent::setData(m_recordComponentData);
     }
 
-    void readBase(bool require_unit_si);
+    void readBase();
 
     template <typename T>
     void verifyChunk(Offset const &, Extent const &) const;
 
     void verifyChunk(Datatype, Offset const &, Extent const &) const;
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 }; // RecordComponent
 
 namespace internal
diff --git a/include/openPMD/RecordComponent.tpp b/include/openPMD/RecordComponent.tpp
index b796ab1a93..523f4f1e41 100644
--- a/include/openPMD/RecordComponent.tpp
+++ b/include/openPMD/RecordComponent.tpp
@@ -90,6 +90,7 @@ RecordComponent::storeChunk(Offset o, Extent e, F &&createBuffer)
     {
         size *= ext;
     }
+
     /*
      * Flushing the skeleton does not create datasets,
      * so we might need to do it now.
@@ -129,7 +130,7 @@ RecordComponent::storeChunk(Offset o, Extent e, F &&createBuffer)
              * actual data yet.
              */
             seriesFlush_impl</* flush_entire_series = */ false>(
-                {FlushLevel::SkeletonOnly});
+                {FlushLevel::SkeletonOnly}, /*flush_io_handler=*/false);
             Parameter<Operation::CREATE_DATASET> dCreate(rc.m_dataset.value());
             dCreate.name = Attributable::get().m_writable.ownKeyWithinParent;
             IOHandler()->enqueue(IOTask(this, dCreate));
diff --git a/include/openPMD/Series.hpp b/include/openPMD/Series.hpp
index 603e540c2b..2cb82d3600 100644
--- a/include/openPMD/Series.hpp
+++ b/include/openPMD/Series.hpp
@@ -31,6 +31,7 @@
 #include "openPMD/auxiliary/Variant.hpp"
 #include "openPMD/backend/Attributable.hpp"
 #include "openPMD/backend/Container.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
 #include "openPMD/backend/ParsePreference.hpp"
 #include "openPMD/config.hpp"
 #include "openPMD/snapshots/Snapshots.hpp"
@@ -64,7 +65,12 @@ namespace openPMD
 class ReadIterations;
 class StatefulIterator;
 class Series;
-class Series;
+
+class Iterations : public Container<Iteration, Iteration::IterationIndex_t>
+{
+public:
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+};
 
 namespace internal
 {
@@ -100,8 +106,8 @@ namespace internal
         SeriesData &operator=(SeriesData &&) = delete;
 
         using IterationIndex_t = Iteration::IterationIndex_t;
-        using IterationsContainer_t = Container<Iteration, IterationIndex_t>;
-        IterationsContainer_t iterations{};
+        using IterationsContainer_t = Iterations;
+        Iterations iterations{};
 
         /**
          * Series::readIterations() returns an iterator type that modifies the
@@ -385,7 +391,7 @@ class Series : public Attributable
      * Type for a container of Iterations indexed by IterationIndex_t.
      */
     using IterationsContainer_t = internal::SeriesData::IterationsContainer_t;
-    IterationsContainer_t iterations;
+    Iterations iterations;
 
     /**
      * @brief Is this a usable Series object?
@@ -778,6 +784,8 @@ class Series : public Attributable
      */
     void close();
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
     /**
      * This overrides Attributable::iterationFlush() which will fail on Series.
      */
@@ -898,9 +906,7 @@ OPENPMD_private
         iterations_iterator end,
         internal::FlushParams const &flushParams,
         bool flushIOHandler = true);
-    void flushMeshesPath();
-    void flushParticlesPath();
-    void flushRankTable();
+    void flushRankTable(FlushLevel);
     /* Parameter `read_only_this_single_iteration` used for reopening an
      * Iteration after closing it.
      */
@@ -984,7 +990,7 @@ OPENPMD_private
      *
      * @param doFlush If true, flush the IO handler.
      */
-    void flushStep(bool doFlush);
+    void flushStep(bool doFlush, FlushLevel l);
 
     /*
      * setIterationEncoding() should only be called by users of our public API,
diff --git a/include/openPMD/auxiliary/TypeTraits.hpp b/include/openPMD/auxiliary/TypeTraits.hpp
index faf2b29472..e373873a4f 100644
--- a/include/openPMD/auxiliary/TypeTraits.hpp
+++ b/include/openPMD/auxiliary/TypeTraits.hpp
@@ -124,6 +124,22 @@ namespace detail
     {
         constexpr static bool value = true;
     };
+
+    template <typename T>
+    struct ScalarType
+    {
+        using type = T;
+    };
+    template <typename T>
+    struct ScalarType<std::vector<T>>
+    {
+        using type = T;
+    };
+    template <typename T, size_t N>
+    struct ScalarType<std::array<T, N>>
+    {
+        using type = T;
+    };
 } // namespace detail
 
 template <typename T>
@@ -141,6 +157,12 @@ using IsPointer_t = typename detail::IsPointer<T>::type;
 template <typename C>
 inline constexpr bool IsChar_v = detail::IsChar<C>::value;
 
+template <typename T>
+using ScalarType_t = typename detail::ScalarType<T>::type;
+
+template <typename T>
+using VectorType_t = std::vector<ScalarType_t<T>>;
+
 /** Emulate in the C++ concept ContiguousContainer
  *
  * Users can implement this trait for a type to signal it can be used as
@@ -203,6 +225,34 @@ namespace detail
     // little trick to avoid trailing commas in the macro expansions below
     template <typename Arg, typename... Args>
     using variant_tail_t = std::variant<Args...>;
+
+    template <template <typename...> class Base, typename... Args>
+    auto infer_template_args(Base<Args...> &) -> Base<Args...>;
+
+    template <
+        template <typename...> class Base,
+        typename T,
+        typename SFINAE = void>
+    struct IsTemplateBaseOf
+    {
+        static constexpr bool value = false;
+    };
+
+    template <template <typename...> class Base, typename T>
+    struct IsTemplateBaseOf<
+        Base,
+        T,
+        std::void_t<decltype(detail::infer_template_args<Base>(
+            std::declval<T &>()))>>
+    {
+        static constexpr bool value = true;
+        using type =
+            decltype(detail::infer_template_args<Base>(std::declval<T &>()));
+    };
 } // namespace detail
 
+template <template <typename...> class Base, typename T>
+constexpr bool IsTemplateBaseOf_v = detail::IsTemplateBaseOf<Base, T>::value;
+template <template <typename...> class Base, typename T>
+using AsTemplateBase_t = typename detail::IsTemplateBaseOf<Base, T>::type;
 } // namespace openPMD::auxiliary
diff --git a/include/openPMD/backend/Attributable.hpp b/include/openPMD/backend/Attributable.hpp
index 705080be84..a27ec2a430 100644
--- a/include/openPMD/backend/Attributable.hpp
+++ b/include/openPMD/backend/Attributable.hpp
@@ -25,6 +25,7 @@
 #include "openPMD/ThrowError.hpp"
 #include "openPMD/auxiliary/OutOfRangeMsg.hpp"
 #include "openPMD/backend/Attribute.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
 #include "openPMD/backend/Writable.hpp"
 
 #include <cstddef>
@@ -57,6 +58,7 @@ namespace internal
     class IterationData;
     class SeriesData;
     struct HomogenizeExtents;
+    struct ConfigAttribute;
 
     class SharedAttributableData
     {
@@ -244,6 +246,7 @@ class Attributable
     friend class internal::AttributableData;
     friend class Snapshots;
     friend struct internal::HomogenizeExtents;
+    friend struct internal::ConfigAttribute;
 
 protected:
     // tag for internal constructor
@@ -276,6 +279,7 @@ class Attributable
     template <typename T>
     bool setAttribute(std::string const &key, T value);
     bool setAttribute(std::string const &key, char const value[]);
+    bool setAttribute(std::string const &key, Attribute value);
     /** @}
      */
 
@@ -398,6 +402,56 @@ class Attributable
      */
     void touch();
 
+    /**
+     * Visitor pattern for the openPMD object hierarchy in postfix traversal.
+     *
+     * @note As the HierarchyVisitor interface can be tedious to implement,
+     *       consider using visitHierarchyFromLambda for a more convenient
+     *       interface.
+
+     * @note to developers: if at any point a prefix traversal should become
+     *       necessary, consider emplacing configuration options for this inside
+     *       HierarchyVisitor to keep the interface clean.
+     *
+     * @param visitor Operations to run for each object.
+     */
+    virtual void visitHierarchy(HierarchyVisitor &visitor, bool recursive);
+
+    /**
+     * Visitor pattern for the openPMD object hierarchy in postfix traversal,
+     * lambda version.
+     *
+     * The lambda parameter will be called with a reference to each object in
+     * the current object's enclosed hierarchy.
+     *
+     * @note Definition inside include/openPMD/backend/HierarchyVisitorImpl.hpp.
+     * @param lambda Operations to run for each object.
+     * @param recursive Extend the operation recursively to children.
+     */
+    template <typename Lambda>
+    void visitHierarchyFromLambda(Lambda &&lambda, bool recursive);
+
+    /** Create standard defined attributes with default values now, insofar they
+     *  are still missing.
+     *
+     * Refer to
+     * https://github.com/openPMD/openPMD-standard/blob/latest/STANDARD.md
+     * for the attributes implied by this operation.
+     *
+     * By default, standard defined attributes are written upon closing the
+     * containing Iteration / Series. Calling this soon can make data available
+     * for early readers (e.g. read while the writer is still modifying). In
+     * workflows that keep single Iterations open over an extended period of
+     * time (e.g. back-transformed diagnostics), this can help creating readable
+     * files earlier than without.
+     *
+     * Attributes may still be modified after this as usual. Attributes defined
+     * before this call will not be modified by it.
+     *
+     * @param recursive Extend the operation recursively to children.
+     */
+    void populateMissingMetadata(bool recursive);
+
     [[nodiscard]] OpenpmdStandard openPMDStandard() const;
 
     // clang-format off
@@ -421,10 +475,14 @@ OPENPMD_protected
     /** @} */
 
     template <bool flush_entire_series>
-    void seriesFlush_impl(internal::FlushParams const &);
+    void seriesFlush_impl(internal::FlushParams const &, bool flush_io_handler);
 
     void flushAttributes(internal::FlushParams const &);
 
+    void populateDefaults(
+        internal::ScientificDefaults *self,
+        internal::FlushParams const &params);
+
     enum ReadMode
     {
         /**
@@ -549,6 +607,27 @@ OPENPMD_protected
     {
         return writable().dirtyRecursive;
     }
+    void determineUnsetDirty(FlushLevel fl)
+    {
+        switch (fl)
+        {
+        case FlushLevel::UserFlush:
+            setDirty(false);
+            break;
+        case FlushLevel::InternalFlush:
+            // Used for parsing
+            if (IOHandler()->m_seriesStatus == internal::SeriesStatus::Parsing)
+            {
+                throw error::Internal(
+                    "Parsing procedures should directly unset dirty.");
+            }
+            break;
+        case FlushLevel::SkeletonOnly:
+        case FlushLevel::CreateOrOpenFiles:
+            // noop
+            break;
+        }
+    }
     void setDirty(bool dirty_in)
     {
         auto &w = writable();
diff --git a/include/openPMD/backend/Attribute.hpp b/include/openPMD/backend/Attribute.hpp
index 9dc621d6e7..d04bdb4539 100644
--- a/include/openPMD/backend/Attribute.hpp
+++ b/include/openPMD/backend/Attribute.hpp
@@ -115,9 +115,30 @@ class Attribute
     template <typename U>
     std::optional<U> getOptional() const;
 
-private:
     template <typename U>
-    std::variant<U, std::runtime_error> get_impl() const;
+    std::variant<U, std::runtime_error> getOrError() const;
+
+    /**
+     * Force this attribute into a vector type.
+     * If it is scalar, convert to a vector with a single element.
+     * If it is an array, convert it to a vector.
+     *
+     * @return On success, a new Attribute with underlying vector type.
+     *         On failure, an error with a description of what went wrong
+     *         (e.g. when the base type is bool).
+     */
+    [[nodiscard]] std::variant<Attribute, std::runtime_error>
+    requireVector() const;
+    /**
+     * Force this attribute into a scalar type.
+     * If it is a vector attribute with a single entry, convert to a scalar.
+     *
+     * @return On success, a new Attribute with underlying scalar type.
+     *         On failure, an error with a description of what went wrong
+     *         (e.g. when a vector with size other than 1 was provided).
+     */
+    [[nodiscard]] std::variant<Attribute, std::runtime_error>
+    requireScalar() const;
 };
 
 namespace detail
diff --git a/include/openPMD/backend/BaseRecord.hpp b/include/openPMD/backend/BaseRecord.hpp
index f75536e3bd..3d7ca5962e 100644
--- a/include/openPMD/backend/BaseRecord.hpp
+++ b/include/openPMD/backend/BaseRecord.hpp
@@ -26,6 +26,8 @@
 #include "openPMD/auxiliary/Variant.hpp"
 #include "openPMD/backend/BaseRecordComponent.hpp"
 #include "openPMD/backend/Container.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <array>
 #include <stdexcept>
@@ -179,6 +181,8 @@ template <typename T_elem>
 class BaseRecord
     : public Container<T_elem>
     , public T_elem // T_RecordComponent
+// ScientificDefaults already inherited via RecordComponent
+// , internal::ScientificDefaults
 {
 public:
     using T_RecordComponent = T_elem;
@@ -197,6 +201,7 @@ class BaseRecord
     friend class internal::ScalarIterator;
     template <typename T>
     friend T &internal::makeOwning(T &self, Series);
+    friend class internal::ScientificDefaults;
 
     using Data_t =
         internal::BaseRecordData<T_elem, typename T_RecordComponent::Data_t>;
@@ -362,15 +367,15 @@ class BaseRecord
      */
     bool scalar() const;
 
-protected:
-    void readBase();
-
 private:
     void flush(std::string const &, internal::FlushParams const &) final;
     virtual void
     flush_impl(std::string const &, internal::FlushParams const &) = 0;
 
     void eraseScalar();
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 }; // BaseRecord
 
 namespace detail
@@ -381,6 +386,21 @@ namespace detail
 
     template <typename BaseRecord>
     void verifyNonscalar(BaseRecord *self);
+
+    template <typename T>
+    struct IsBaseRecord
+    {
+        constexpr static bool value = false;
+    };
+
+    template <typename T>
+    struct IsBaseRecord<BaseRecord<T>>
+    {
+        constexpr static bool value = true;
+    };
+
+    template <typename T>
+    constexpr bool IsBaseRecord_v = IsBaseRecord<T>::value;
 } // namespace detail
 
 template <typename T_elem>
diff --git a/include/openPMD/backend/BaseRecordComponent.hpp b/include/openPMD/backend/BaseRecordComponent.hpp
index b3c85f869d..af23e08468 100644
--- a/include/openPMD/backend/BaseRecordComponent.hpp
+++ b/include/openPMD/backend/BaseRecordComponent.hpp
@@ -82,6 +82,7 @@ class BaseRecordComponent : virtual public Attributable
 {
     template <typename T, typename T_key, typename T_container>
     friend class Container;
+    friend class Attributable;
 
 public:
     /*
diff --git a/include/openPMD/backend/Container.hpp b/include/openPMD/backend/Container.hpp
index 963fd34802..839f6452c9 100644
--- a/include/openPMD/backend/Container.hpp
+++ b/include/openPMD/backend/Container.hpp
@@ -22,7 +22,9 @@
 
 #include "openPMD/Error.hpp"
 #include "openPMD/IO/Access.hpp"
+#include "openPMD/auxiliary/TypeTraits.hpp"
 #include "openPMD/backend/Attributable.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
 
 #include <initializer_list>
 #include <map>
@@ -261,6 +263,19 @@ class Container : virtual public Attributable
         return container().emplace(std::forward<Args>(args)...);
     }
 
+    template <typename ChildClass>
+    void visitHierarchyImpl(HierarchyVisitor &v, bool recursive)
+    {
+        if (recursive)
+        {
+            for (auto &p : *this)
+            {
+                p.second.visitHierarchy(v, recursive);
+            }
+        }
+        v(*static_cast<ChildClass *>(this));
+    }
+
     // clang-format off
 OPENPMD_protected
     // clang-format on
diff --git a/include/openPMD/backend/HierarchyVisitor.hpp b/include/openPMD/backend/HierarchyVisitor.hpp
new file mode 100644
index 0000000000..5d7e71bfb2
--- /dev/null
+++ b/include/openPMD/backend/HierarchyVisitor.hpp
@@ -0,0 +1,39 @@
+#pragma once
+
+#include <cstdint>
+#include <map>
+#include <string>
+
+namespace openPMD
+{
+class Iteration;
+class Iterations;
+class Mesh;
+class MeshRecordComponent;
+class Meshes;
+class ParticlePatches;
+class ParticleSpecies;
+class Particles;
+class PatchRecord;
+class PatchRecordComponent;
+class Record;
+class RecordComponent;
+class Series;
+
+class HierarchyVisitor
+{
+public:
+    virtual void operator()(Iteration &) = 0;
+    virtual void operator()(Iterations &) = 0;
+    virtual void operator()(Mesh &) = 0;
+    virtual void operator()(MeshRecordComponent &) = 0;
+    virtual void operator()(Meshes &) = 0;
+    virtual void operator()(ParticlePatches &) = 0;
+    virtual void operator()(ParticleSpecies &) = 0;
+    virtual void operator()(Particles &) = 0;
+    virtual void operator()(PatchRecordComponent &) = 0;
+    virtual void operator()(Record &) = 0;
+    virtual void operator()(RecordComponent &) = 0;
+    virtual void operator()(Series &) = 0;
+};
+} // namespace openPMD
diff --git a/include/openPMD/backend/HierarchyVisitorImpl.hpp b/include/openPMD/backend/HierarchyVisitorImpl.hpp
new file mode 100644
index 0000000000..432085e4cb
--- /dev/null
+++ b/include/openPMD/backend/HierarchyVisitorImpl.hpp
@@ -0,0 +1,89 @@
+#pragma once
+
+#include "openPMD/backend/HierarchyVisitor.hpp"
+
+#include "openPMD/Iteration.hpp"
+#include "openPMD/ParticleSpecies.hpp"
+#include "openPMD/Series.hpp"
+#include <type_traits>
+
+namespace openPMD
+{
+template <typename Lambda>
+class HierarchyVisitorFromLambda : public HierarchyVisitor
+{
+    Lambda lambda;
+
+    template <typename Arg>
+    HierarchyVisitorFromLambda(uint8_t /* constructor_tag */, Arg &&arg)
+        : lambda(std::forward<Arg>(arg))
+    {}
+
+    template <typename Lambda_in>
+    friend auto makeHierarchyVisitorFromLambda(Lambda_in &&lambda);
+
+public:
+    void operator()(Iteration &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(Iterations &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(Mesh &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(MeshRecordComponent &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(Meshes &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(ParticlePatches &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(ParticleSpecies &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(Particles &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(PatchRecordComponent &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(Record &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(RecordComponent &obj) override
+    {
+        lambda(obj);
+    }
+    void operator()(Series &obj) override
+    {
+        lambda(obj);
+    }
+};
+
+template <typename Lambda>
+auto makeHierarchyVisitorFromLambda(Lambda &&lambda)
+{
+    using res_t = HierarchyVisitorFromLambda<std::remove_reference_t<Lambda>>;
+    return res_t{0, std::forward<Lambda>(lambda)};
+}
+
+template <typename Lambda>
+void Attributable::visitHierarchyFromLambda(Lambda &&lambda, bool recursive)
+{
+    auto visitor = makeHierarchyVisitorFromLambda(std::forward<Lambda>(lambda));
+    this->visitHierarchy(visitor, recursive);
+}
+} // namespace openPMD
diff --git a/include/openPMD/backend/MeshRecordComponent.hpp b/include/openPMD/backend/MeshRecordComponent.hpp
index 708c4a496f..4399c65025 100644
--- a/include/openPMD/backend/MeshRecordComponent.hpp
+++ b/include/openPMD/backend/MeshRecordComponent.hpp
@@ -21,6 +21,8 @@
 #pragma once
 
 #include "openPMD/RecordComponent.hpp"
+#include "openPMD/backend/HierarchyVisitor.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <vector>
 
@@ -40,6 +42,7 @@ class MeshRecordComponent : public RecordComponent
     friend class BaseRecord;
     template <typename, typename>
     friend class internal::BaseRecordData;
+    friend class internal::ScientificDefaults;
 
     friend class Mesh;
 
@@ -91,6 +94,11 @@ class MeshRecordComponent : public RecordComponent
      */
     template <typename T>
     MeshRecordComponent &makeConstant(T);
+
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 };
 
 template <typename T>
diff --git a/include/openPMD/backend/PatchRecord.hpp b/include/openPMD/backend/PatchRecord.hpp
index 86317745b8..49a6f1448d 100644
--- a/include/openPMD/backend/PatchRecord.hpp
+++ b/include/openPMD/backend/PatchRecord.hpp
@@ -22,9 +22,9 @@
 
 #include "openPMD/backend/BaseRecord.hpp"
 #include "openPMD/backend/PatchRecordComponent.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <string>
-#include <unordered_map>
 
 namespace openPMD
 {
@@ -33,11 +33,15 @@ class PatchRecord : public BaseRecord<PatchRecordComponent>
     friend class Container<PatchRecord>;
     friend class ParticleSpecies;
     friend class ParticlePatches;
+    friend class internal::ScientificDefaults;
 
 public:
-    PatchRecord &setUnitDimension(std::map<UnitDimension, double> const &);
+    PatchRecord &setUnitDimension(unit_representations::AsMap const &);
+    PatchRecord &setUnitDimension(unit_representations::AsArray const &udim);
     ~PatchRecord() override = default;
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
 private:
     PatchRecord() = default;
 
diff --git a/include/openPMD/backend/PatchRecordComponent.hpp b/include/openPMD/backend/PatchRecordComponent.hpp
index eedd84e158..1fb98885e9 100644
--- a/include/openPMD/backend/PatchRecordComponent.hpp
+++ b/include/openPMD/backend/PatchRecordComponent.hpp
@@ -24,6 +24,7 @@
 #include "openPMD/RecordComponent.hpp"
 #include "openPMD/auxiliary/ShareRawInternal.hpp"
 #include "openPMD/backend/BaseRecordComponent.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <memory>
 #include <sstream>
@@ -53,6 +54,7 @@ class PatchRecordComponent : public RecordComponent
     friend class ParticlePatches;
     friend class PatchRecord;
     friend class ParticleSpecies;
+    friend class internal::ScientificDefaults;
 
 public:
     /**
@@ -87,18 +89,27 @@ class PatchRecordComponent : public RecordComponent
     template <typename T>
     void store(T);
 
+    void visitHierarchy(HierarchyVisitor &v, bool recursive) override;
+
     // clang-format off
 OPENPMD_private
     // clang-format on
 
     using RecordComponent::flush;
 
+    // This differs from RecordComponent::read() such that unitSI is not
+    // required
+    void read();
+
     // clang-format off
 OPENPMD_protected
     // clang-format on
 
     PatchRecordComponent();
     PatchRecordComponent(NoInit);
+
+protected:
+    void scientificDefaults_impl(bool write, OpenpmdStandard) override;
 }; // PatchRecordComponent
 
 template <typename T>
diff --git a/include/openPMD/backend/Writable.hpp b/include/openPMD/backend/Writable.hpp
index 0396aa1f93..3623d0b39c 100644
--- a/include/openPMD/backend/Writable.hpp
+++ b/include/openPMD/backend/Writable.hpp
@@ -51,6 +51,7 @@ namespace internal
     class SharedAttributableData;
     class AttributableData;
     class SeriesData;
+    class ScientificDefaults;
 } // namespace internal
 namespace detail
 {
@@ -106,6 +107,7 @@ class Writable final
     friend void debug::printDirty(Series const &);
     friend struct Parameter<Operation::CREATE_DATASET>;
     friend struct Parameter<Operation::OPEN_DATASET>;
+    friend class internal::ScientificDefaults;
 
 private:
     Writable(internal::AttributableData *);
@@ -126,14 +128,15 @@ class Writable final
      * it.
      */
     template <bool flush_entire_series>
-    void seriesFlush(std::string backendConfig = "{}");
+    void
+    seriesFlush(std::string backendConfig = "{}", bool flush_io_handler = true);
 
     // clang-format off
 OPENPMD_private
     // clang-format on
 
     template <bool flush_entire_series>
-    void seriesFlush(internal::FlushParams const &);
+    void seriesFlush(internal::FlushParams const &, bool flush_io_handler);
     /*
      * These members need to be shared pointers since distinct instances of
      * Writable may share them.
diff --git a/include/openPMD/backend/scientific_defaults/AttributeReader.hpp b/include/openPMD/backend/scientific_defaults/AttributeReader.hpp
new file mode 100644
index 0000000000..2eba7b9634
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/AttributeReader.hpp
@@ -0,0 +1,49 @@
+#pragma once
+
+#include "openPMD/Error.hpp"
+#include "openPMD/backend/Attributable.hpp"
+#include "openPMD/backend/Attribute.hpp"
+#include "openPMD/backend/scientific_defaults/ProcessParsedAttribute.hpp"
+
+#include <deque>
+#include <optional>
+#include <variant>
+
+namespace openPMD::internal
+{
+namespace attribute_read_result
+{
+    struct TypeUnmatched
+    {
+        std::deque<Datatype> expectedDatatypes;
+    };
+    struct Success
+    {};
+} // namespace attribute_read_result
+
+using AttributeReadResult = std::variant<
+    attribute_read_result::Success,
+    attribute_read_result::TypeUnmatched,
+    error::ReadError>;
+
+/*
+ * Struct that describes parsing logic for a standard-defined attribute.
+ * Created by ConfigAttribute::withReader().
+ */
+struct AttributeReader
+{
+    std::deque<Datatype> eligibleDatatypes;
+    std::optional<std::shared_ptr<ProcessParsedAttribute>> processAttribute;
+
+    AttributeReader(
+        std::deque<Datatype> eligibleDatatypes_in,
+        std::optional<std::shared_ptr<ProcessParsedAttribute>>
+            processAttribute_in);
+
+    auto operator()(
+        Attributable &record,
+        char const *attrName,
+        Attribute const &a,
+        std::deque<Datatype> unmatched_so_far) -> AttributeReadResult;
+};
+} // namespace openPMD::internal
diff --git a/include/openPMD/backend/scientific_defaults/ConfigAttribute.hpp b/include/openPMD/backend/scientific_defaults/ConfigAttribute.hpp
new file mode 100644
index 0000000000..b0c0ed0d24
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/ConfigAttribute.hpp
@@ -0,0 +1,63 @@
+#pragma once
+
+#include "openPMD/backend/Attributable.hpp"
+#include "openPMD/backend/scientific_defaults/AttributeReader.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults_auxiliary.hpp"
+
+#include <deque>
+#include <optional>
+
+namespace openPMD::internal
+{
+struct ConfigAttribute
+{
+    Attributable &attributable;
+    char const *attrName;
+    std::optional<std::function<void(Attributable &)>> initDefaultAttribute;
+    // processed "from left to right"
+    std::deque<AttributeReader> attributeReaders;
+
+    template <typename RecordType, typename ValueType>
+    using set_default_val_t = RecordType &(RecordType::*)(ValueType);
+
+    ConfigAttribute(Attributable &attributable_in, char const *attrName_in);
+
+    ConfigAttribute(ConfigAttribute const &) = delete;
+    ConfigAttribute(ConfigAttribute &&) = delete;
+
+    ConfigAttribute &operator=(ConfigAttribute const &) = delete;
+    ConfigAttribute &operator=(ConfigAttribute &&) = delete;
+
+    // Specify a custom function for setting the default attribute,
+    // instead of just using setAttribute().
+    template <typename RecordType, typename S = void, typename GetDefaultValue>
+    [[nodiscard]] auto withSetter(
+        GetDefaultValue &&getDefaultVal,
+        set_default_val_t<
+            RecordType,
+            std::conditional_t<
+                std::is_void_v<S>,
+                auxiliary::CallResult_t<GetDefaultValue>,
+                S>> setDefaultVal) -> ConfigAttribute &;
+
+    // Generically use setAttribute() for initializing
+    template <typename DefaultValue>
+    [[nodiscard]] auto withGenericSetter(DefaultValue &&defaultVal)
+        -> ConfigAttribute &;
+
+    // Add readers for attribute parsing
+    [[nodiscard]] auto withReader(
+        std::deque<Datatype> eligibleDatatypes,
+        std::optional<std::shared_ptr<ProcessParsedAttribute>>
+            processAttribute = std::nullopt) -> ConfigAttribute &;
+
+    void write();
+    void read();
+    void operator()(WriteOrRead wor);
+};
+
+auto get_float_types() -> std::deque<Datatype>;
+auto get_string_types() -> std::deque<Datatype>;
+} // namespace openPMD::internal
+
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.tpp"
diff --git a/include/openPMD/backend/scientific_defaults/ConfigAttribute.tpp b/include/openPMD/backend/scientific_defaults/ConfigAttribute.tpp
new file mode 100644
index 0000000000..70c15c2f35
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/ConfigAttribute.tpp
@@ -0,0 +1,73 @@
+#pragma once
+
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+
+namespace openPMD::internal
+{
+namespace
+{
+    // Need SFINAE for this since MSVC doesnt understand if constexpr
+    template <typename GetDefaultValue, typename SFINAE = void>
+    struct CallGetDefaultValue
+    {
+        GetDefaultValue val;
+        CallGetDefaultValue(GetDefaultValue val_in) : val(std::move(val_in))
+        {}
+        auto operator()() && -> GetDefaultValue
+        {
+            return std::move(val);
+        }
+    };
+
+    template <typename GetDefaultValue>
+    struct CallGetDefaultValue<
+        GetDefaultValue,
+        std::void_t<std::invoke_result_t<GetDefaultValue>>>
+    {
+        GetDefaultValue val;
+        CallGetDefaultValue(GetDefaultValue val_in) : val(std::move(val_in))
+        {}
+        auto operator()() && -> std::invoke_result_t<GetDefaultValue>
+        {
+            return std::move(val)();
+        }
+    };
+} // namespace
+
+template <typename RecordType, typename S, typename GetDefaultValue>
+auto ConfigAttribute::withSetter(
+    GetDefaultValue &&getDefaultVal,
+    set_default_val_t<
+        RecordType,
+        std::conditional_t<
+            std::is_void_v<S>,
+            auxiliary::CallResult_t<GetDefaultValue>,
+            S>> setDefaultVal) -> ConfigAttribute &
+{
+    initDefaultAttribute =
+        [callDefaultValue =
+             CallGetDefaultValue<std::remove_reference_t<GetDefaultValue>>(
+                 std::forward<GetDefaultValue>(getDefaultVal)),
+         setDefaultVal](Attributable &attr) mutable {
+            RecordType *record = dynamic_cast<RecordType *>(&attr);
+            if (!record)
+            {
+                throw error::Internal("dynamic cast failure");
+            }
+            ((*record).*setDefaultVal)(std::move(callDefaultValue)());
+        };
+    return *this;
+}
+
+template <typename DefaultValue>
+auto ConfigAttribute::withGenericSetter(DefaultValue &&defaultVal)
+    -> ConfigAttribute &
+{
+    initDefaultAttribute =
+        [this, defaultVal_lambda = std::forward<DefaultValue &&>(defaultVal)](
+            Attributable &attr) {
+            attr.setAttribute(this->attrName, std::move(defaultVal_lambda));
+        };
+    return *this;
+}
+} // namespace openPMD::internal
diff --git a/include/openPMD/backend/scientific_defaults/ProcessParsedAttribute.hpp b/include/openPMD/backend/scientific_defaults/ProcessParsedAttribute.hpp
new file mode 100644
index 0000000000..5286daaeec
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/ProcessParsedAttribute.hpp
@@ -0,0 +1,150 @@
+#pragma once
+
+#include "openPMD/Error.hpp"
+#include "openPMD/backend/Attributable.hpp"
+#include "openPMD/backend/Attribute.hpp"
+
+#include <optional>
+#include <string>
+#include <type_traits>
+
+namespace openPMD
+{
+class Mesh;
+}
+
+namespace openPMD::internal
+{
+/*
+ * General interface for attribute processing used in struct AttributeReader.
+ */
+struct ProcessParsedAttribute
+{
+    virtual auto operator()(Attributable &, char const *, Attribute const &)
+        -> std::optional<error::ReadError> = 0;
+    virtual ~ProcessParsedAttribute() = default;
+};
+
+//////////////////////////////////////////////////////////
+// 1st implementation: RequireType and associated types //
+//////////////////////////////////////////////////////////
+
+/*
+ * Helper interface for validating an attribute whose type has already been
+ * determined. Poor-man's lambda, in an attempt at reducing use of template
+ * types. The only implementation at PostProcessConvertedAttributeImpl<T> is
+ * conceptually equivalent to storing an object of type RecordType in the lambda
+ * capture.
+ */
+template <typename T>
+struct PostProcessConvertedAttribute
+{
+    virtual auto operator()(T val) -> std::optional<error::ReadError> = 0;
+    virtual ~PostProcessConvertedAttribute() = default;
+};
+
+/*
+ * Defer validation to a function pointer of type handler_t.
+ */
+template <typename T, typename RecordType>
+struct PostProcessConvertedAttributeImpl : PostProcessConvertedAttribute<T>
+{
+    RecordType record;
+    using handler_t = std::optional<error::ReadError> (*)(RecordType &, T);
+    handler_t reader;
+
+    // Use makePostProcessConvertedAttribute below for better type inference
+    PostProcessConvertedAttributeImpl(
+        RecordType record_in, handler_t reader_in);
+    auto operator()(T val) -> std::optional<error::ReadError> override;
+};
+
+// Postprocessing handlers for `handler_t reader` above.
+auto setMeshGeometryFromString(Mesh &mesh, std::string val)
+    -> std::optional<error::ReadError>;
+auto setMeshDataOrderFromChar(Mesh &mesh, char val)
+    -> std::optional<error::ReadError>;
+
+template <typename T, typename RecordType>
+auto makePostProcessConvertedAttribute(
+    RecordType &&record,
+    std::optional<error::ReadError> (*handler)(
+        std::remove_reference_t<RecordType> &, T))
+    -> std::shared_ptr<PostProcessConvertedAttribute<T>>;
+/*
+ * Validate an attribute by requiring one specific type T, and by optionally
+ * postprocessing it. Type conversions for the attribute are attempted.
+ * Note that type checking for attribute is mainly intended to by done by
+ * specifying eligible data types in struct AttributeReader. This handler's
+ * purpose is dynamically casting to one specific type not required by the
+ * standard, but by our further handling (e.g. setMeshGeometryFromString above
+ * requires a string, but we should also accept single chars, or vectors with a
+ * single element if we we find such things on disk).
+ * If a handler is not specified, this just casts the type and then uses
+ * setAttribute().
+ */
+template <typename T>
+struct RequireType : ProcessParsedAttribute
+{
+    std::optional<std::shared_ptr<PostProcessConvertedAttribute<T>>>
+        postProcess;
+
+    explicit RequireType() = default;
+
+    // Use require_type() helpers below for construction.
+    template <typename RecordType>
+    RequireType(
+        RecordType &&record,
+        std::optional<error::ReadError> (*handler)(
+            std::remove_reference_t<RecordType> &, T));
+
+    auto operator()(
+        Attributable &record, char const *attrName, Attribute const &attr)
+        -> std::optional<error::ReadError> override;
+};
+
+template <typename T, typename RecordType>
+auto require_type(
+    std::optional<error::ReadError> (*)(
+        std::remove_reference_t<RecordType> &, T))
+    -> std::shared_ptr<ProcessParsedAttribute>;
+// common case: directly use setAttribute
+template <typename T>
+auto require_type() -> std::shared_ptr<ProcessParsedAttribute>;
+
+///////////////////////////////////////
+// 2nd implementation: RequireVector //
+///////////////////////////////////////
+
+/*
+ * Validate an attribute by requiring a vector type, potentially wrapping
+ * scalar values into a vector.
+ * Note that this employs no type checking for the base type, this is done
+ * by the type list argument of ConfigAttribute::withReader().
+ */
+struct RequireVector : ProcessParsedAttribute
+{
+    auto operator()(Attributable &, char const *, Attribute const &)
+        -> std::optional<error::ReadError> override;
+};
+extern std::shared_ptr<ProcessParsedAttribute> require_vector;
+
+///////////////////////////////////////
+// 3rd implementation: RequireScalar //
+///////////////////////////////////////
+
+/*
+ * Validate an attribute by requiring a scalar type, potentially unwrapping
+ * single values from a vector.
+ * Note that this employs no type checking for the base type, this is done
+ * by the type list argument of ConfigAttribute::withReader().
+ */
+struct RequireScalar : ProcessParsedAttribute
+{
+    auto operator()(Attributable &, char const *, Attribute const &)
+        -> std::optional<error::ReadError> override;
+};
+extern std::shared_ptr<ProcessParsedAttribute> require_scalar;
+} // namespace openPMD::internal
+
+#include "openPMD/backend/scientific_defaults/ProcessParsedAttribute.tpp"
diff --git a/include/openPMD/backend/scientific_defaults/ProcessParsedAttribute.tpp b/include/openPMD/backend/scientific_defaults/ProcessParsedAttribute.tpp
new file mode 100644
index 0000000000..b8987b1929
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/ProcessParsedAttribute.tpp
@@ -0,0 +1,93 @@
+#pragma once
+
+#include "openPMD/backend/scientific_defaults/ProcessParsedAttribute.hpp"
+
+#include <variant>
+
+namespace openPMD::internal
+{
+template <typename T, typename RecordType>
+PostProcessConvertedAttributeImpl<T, RecordType>::
+    PostProcessConvertedAttributeImpl(RecordType record_in, handler_t reader_in)
+    : record(std::move(record_in)), reader(reader_in)
+{}
+
+template <typename T, typename RecordType>
+auto PostProcessConvertedAttributeImpl<T, RecordType>::operator()(T val)
+    -> std::optional<error::ReadError>
+{
+    return (*reader)(record, std::move(val));
+}
+
+template <typename T, typename RecordType>
+auto require_type(
+    RecordType &&record,
+    std::optional<error::ReadError> (*handler)(
+        std::remove_reference_t<RecordType> &, T))
+    -> std::shared_ptr<ProcessParsedAttribute>
+{
+    return std::make_shared<RequireType<T>>(
+        std::forward<RecordType>(record), handler);
+}
+
+template <typename T>
+auto require_type() -> std::shared_ptr<ProcessParsedAttribute>
+{
+    return std::make_shared<RequireType<T>>();
+}
+
+template <typename T, typename RecordType>
+auto makePostProcessConvertedAttribute(
+    RecordType &&record,
+    std::optional<error::ReadError> (*handler)(
+        std::remove_reference_t<RecordType> &, T))
+    -> std::shared_ptr<PostProcessConvertedAttribute<T>>
+{
+    return std::make_shared<PostProcessConvertedAttributeImpl<
+        T,
+        std::remove_reference_t<RecordType>>>(
+        std::forward<RecordType>(record), handler);
+}
+
+template <typename T>
+template <typename RecordType>
+RequireType<T>::RequireType(
+    RecordType &&record,
+    std::optional<error::ReadError> (*handler)(
+        std::remove_reference_t<RecordType> &, T))
+    : postProcess(makePostProcessConvertedAttribute(
+          std::forward<RecordType>(record), handler))
+{}
+
+template <typename T>
+auto RequireType<T>::operator()(
+    Attributable &record, char const *attrName, Attribute const &attr)
+    -> std::optional<error::ReadError>
+{
+    auto converted_or_error = attr.getOrError<T>();
+    return std::visit(
+        auxiliary::overloaded{
+            [&](T casted_val) -> std::optional<error::ReadError> {
+                if (this->postProcess.has_value())
+                {
+                    return (**this->postProcess)(std::move(casted_val));
+                }
+                else
+                {
+                    record.setAttribute<T>(attrName, std::move(casted_val));
+                    return std::nullopt;
+                }
+            },
+            [](std::runtime_error const &err)
+                -> std::optional<error::ReadError> {
+                std::string msg = "Expected a scalar type: ";
+                msg += err.what();
+                return error::ReadError(
+                    error::AffectedObject::Attribute,
+                    error::Reason::UnexpectedContent,
+                    std::nullopt,
+                    std::move(msg));
+            }},
+        converted_or_error);
+}
+} // namespace openPMD::internal
diff --git a/include/openPMD/backend/scientific_defaults/ScientificDefaults.hpp b/include/openPMD/backend/scientific_defaults/ScientificDefaults.hpp
new file mode 100644
index 0000000000..c1d07e0592
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/ScientificDefaults.hpp
@@ -0,0 +1,51 @@
+#pragma once
+
+#include "openPMD/backend/Attributable.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+
+namespace openPMD
+{
+class Attributable;
+}
+
+namespace openPMD::internal
+{
+/*
+ * This is the interface for implementing required scientific attributes defined
+ * by the openPMD standard (writing and reading). The interface is inherited by
+ * classes in the openPMD hierarchy object model (e.g. Iteration,
+ * ParticleSpecies, ...). This is the place for implementing (most of) the
+ * attribute definitions from
+ * github.com/openPMD/openPMD-standard/blob/latest/STANDARD.md
+ */
+class ScientificDefaults
+{
+    friend class ::openPMD::Attributable;
+
+protected:
+    [[nodiscard]] auto defaultAttribute(Attributable &, char const *attrName)
+        -> ConfigAttribute;
+
+    // Convention: This is called for each openPMD object (group, dataset)
+    // exactly once and only upon its most derived object (e.g. upon
+    // RecordComponent, not BaseRecordComponent). The defaults_impl() definition
+    // has to deal with parent class definitions. This is necessary since the
+    // derived class may override standard attribute definitions from a base
+    // class.
+    //
+    // Used both by writeDefaults() and readDefaults().
+    //
+    // Use defaultAttribute() for creating an attribute, use helpers inside
+    // ConfigAttribute for defining it.
+    virtual void scientificDefaults_impl(bool write, OpenpmdStandard) = 0;
+
+    // Called upon Iteration::close(), will fill in defaults below Iteration
+    // level. If the Iteration is not explicitly closed, will be called upon
+    // Series::close().
+    // Currently called internally only from Iteration::setDefaultsAttributes().
+    void writeDefaults(OpenpmdStandard);
+
+    // Called at appropriate places during parsing.
+    void readDefaults(OpenpmdStandard);
+};
+} // namespace openPMD::internal
diff --git a/include/openPMD/backend/scientific_defaults/ScientificDefaults_auxiliary.hpp b/include/openPMD/backend/scientific_defaults/ScientificDefaults_auxiliary.hpp
new file mode 100644
index 0000000000..84413dadc0
--- /dev/null
+++ b/include/openPMD/backend/scientific_defaults/ScientificDefaults_auxiliary.hpp
@@ -0,0 +1,47 @@
+#pragma once
+
+#include "openPMD/backend/Attribute.hpp"
+
+#include <iostream>
+#include <string>
+#include <type_traits>
+#include <vector>
+
+namespace openPMD::auxiliary
+{
+template <typename F, typename SFINAE = void>
+struct CallResult
+{
+    // If F is not callable, then use its type directly without calling.
+    using type = F;
+};
+
+template <typename F>
+struct CallResult<F, std::void_t<decltype(std::declval<F>()())>>
+{
+    // If F is callable (used to avoid complex computations until necessary),
+    // then use its return type.
+    using type = decltype(std::declval<F>()());
+};
+
+template <typename F>
+using CallResult_t = typename CallResult<F>::type;
+
+template <typename T>
+auto write_val_to_stderr(T const &val) -> std::ostream &;
+auto write_to_stderr(Attribute const &a) -> std::ostream &;
+
+auto createDefaultAxisLabels(uint64_t dimensionality)
+    -> std::vector<std::string>;
+auto createDefaultVector(uint64_t dimensionality, double defaultValue)
+    -> std::vector<double>;
+} // namespace openPMD::auxiliary
+
+namespace openPMD::internal
+{
+enum class WriteOrRead : std::uint8_t
+{
+    Write,
+    Read
+};
+}
diff --git a/include/openPMD/openPMD.hpp b/include/openPMD/openPMD.hpp
index 083a1e3342..91e5f79195 100644
--- a/include/openPMD/openPMD.hpp
+++ b/include/openPMD/openPMD.hpp
@@ -45,6 +45,7 @@ namespace openPMD
 #include "openPMD/backend/BaseRecord.hpp"
 #include "openPMD/backend/BaseRecordComponent.hpp"
 #include "openPMD/backend/Container.hpp"
+#include "openPMD/backend/HierarchyVisitorImpl.hpp"
 #include "openPMD/backend/MeshRecordComponent.hpp"
 #include "openPMD/backend/PatchRecord.hpp"
 #include "openPMD/backend/PatchRecordComponent.hpp"
diff --git a/src/IO/ADIOS/ADIOS2File.cpp b/src/IO/ADIOS/ADIOS2File.cpp
index af2d42640c..1b4fba9792 100644
--- a/src/IO/ADIOS/ADIOS2File.cpp
+++ b/src/IO/ADIOS/ADIOS2File.cpp
@@ -1035,25 +1035,16 @@ void ADIOS2File::flush_impl(
         drainedUniquePtrPuts.swap(m_uniquePtrPuts);
     }
 
-    if (readOnly(m_mode))
+    if (readOnly(m_mode) || flush_level::write_datasets(level))
     {
-        level = FlushLevel::UserFlush;
-    }
-
-    switch (level)
-    {
-    case FlushLevel::UserFlush:
         performPutGets(*this, eng);
         m_updateSpans.clear();
         m_buffer.clear();
         m_alreadyEnqueued.clear();
         drainedUniquePtrPuts.clear();
-
-        break;
-
-    case FlushLevel::InternalFlush:
-    case FlushLevel::SkeletonOnly:
-    case FlushLevel::CreateOrOpenFiles:
+    }
+    else
+    {
         /*
          * Tasks have been given to ADIOS2, but we don't flush them
          * yet. So, move everything to m_alreadyEnqueued to avoid
@@ -1070,7 +1061,6 @@ void ADIOS2File::flush_impl(
                 "wrong time.");
         }
         m_buffer.clear();
-        break;
     }
 }
 
diff --git a/src/IO/ADIOS/ADIOS2IOHandler.cpp b/src/IO/ADIOS/ADIOS2IOHandler.cpp
index 3f16bcf155..aec4487f5a 100644
--- a/src/IO/ADIOS/ADIOS2IOHandler.cpp
+++ b/src/IO/ADIOS/ADIOS2IOHandler.cpp
@@ -582,7 +582,7 @@ overrideFlushTarget(FlushTarget &inplace, FlushTarget new_val)
 std::future<void>
 ADIOS2IOHandlerImpl::flush(internal::ParsedFlushParams &flushParams)
 {
-    auto res = AbstractIOHandlerImpl::flush();
+    auto res = AbstractIOHandlerImpl::flush(flushParams.flushLevel);
 
     detail::ADIOS2File::ADIOS2FlushParams adios2FlushParams{
         flushParams.flushLevel, m_flushTarget};
diff --git a/src/IO/AbstractIOHandler.cpp b/src/IO/AbstractIOHandler.cpp
index 5f2bdeb2f5..564eebba01 100644
--- a/src/IO/AbstractIOHandler.cpp
+++ b/src/IO/AbstractIOHandler.cpp
@@ -27,6 +27,29 @@
 
 #include <utility>
 
+namespace openPMD
+{
+std::ostream &operator<<(std::ostream &os, FlushLevel l)
+{
+    switch (l)
+    {
+    case FlushLevel::UserFlush:
+        os << "UserFlush";
+        break;
+    case FlushLevel::InternalFlush:
+        os << "InternalFlush";
+        break;
+    case FlushLevel::SkeletonOnly:
+        os << "SkeletonOnly";
+        break;
+    case FlushLevel::CreateOrOpenFiles:
+        os << "CreateOrOpenFiles";
+        break;
+    }
+    return os;
+}
+} // namespace openPMD
+
 namespace openPMD::auxiliary
 {
 using pair_t = std::pair<OpenpmdStandard, char const *>;
diff --git a/src/IO/AbstractIOHandlerImpl.cpp b/src/IO/AbstractIOHandlerImpl.cpp
index 4f93ff1a5b..473b140aef 100644
--- a/src/IO/AbstractIOHandlerImpl.cpp
+++ b/src/IO/AbstractIOHandlerImpl.cpp
@@ -21,6 +21,8 @@
 
 #include "openPMD/IO/AbstractIOHandlerImpl.hpp"
 
+#include "openPMD/Error.hpp"
+#include "openPMD/IO/AbstractIOHandler.hpp"
 #include "openPMD/IO/IOTask.hpp"
 #include "openPMD/Streaming.hpp"
 #include "openPMD/auxiliary/Environment.hpp"
@@ -87,13 +89,76 @@ void AbstractIOHandlerImpl::writeToStderr([[maybe_unused]] Args &&...args) const
     }
 }
 
-std::future<void> AbstractIOHandlerImpl::flush()
+namespace
+{
+    void verifyFlushType(Operation op, FlushLevel l)
+    {
+        auto do_throw = [&](char const *least_flush_level) {
+            std::stringstream err;
+            err << "Operation " << op << " is not allowed below flush level "
+                << least_flush_level << ", but flush level was " << l << ".";
+            throw error::Internal(err.str());
+        };
+        switch (op)
+        {
+        case Operation::ADVANCE:
+        case Operation::CREATE_FILE:
+        case Operation::CHECK_FILE:
+        case Operation::OPEN_FILE:
+        case Operation::CLOSE_FILE:
+        case Operation::DELETE_FILE:
+        case Operation::DEREGISTER:
+        case Operation::TOUCH:
+        case Operation::LIST_ATTS:
+        case Operation::LIST_PATHS:
+        case Operation::OPEN_PATH:
+        case Operation::SET_WRITTEN:
+        case Operation::CREATE_PATH:
+            break;
+        case Operation::CLOSE_PATH:
+        case Operation::DELETE_PATH:
+        case Operation::CREATE_DATASET:
+        case Operation::EXTEND_DATASET:
+        case Operation::OPEN_DATASET:
+        case Operation::DELETE_DATASET:
+        case Operation::LIST_DATASETS:
+            if (!flush_level::flush_hierarchy(l))
+            {
+                do_throw("SkeletonOnly (hierarchy operations)");
+            }
+            break;
+        case Operation::GET_BUFFER_VIEW:
+        case Operation::DELETE_ATT:
+        case Operation::WRITE_ATT:
+        case Operation::READ_ATT:
+        case Operation::READ_ATT_ALLSTEPS:
+        case Operation::AVAILABLE_CHUNKS:
+            if (!flush_level::write_attributes(l))
+            {
+                do_throw("InternalFlush (metadata operations)");
+            }
+            break;
+        case Operation::WRITE_DATASET:
+        case Operation::READ_DATASET:
+            if (!flush_level::write_datasets(l))
+            {
+                do_throw("UserFlush (flushpoint operations)");
+            }
+            break;
+        }
+    }
+} // namespace
+
+std::future<void> AbstractIOHandlerImpl::flush(FlushLevel l)
 {
     using namespace auxiliary;
 
+    writeToStderr("\nFLUSHING");
+
     while (!(*m_handler).m_work.empty())
     {
         IOTask &i = (*m_handler).m_work.front();
+        // verifyFlushType(i.operation, l);
         try
         {
             switch (i.operation)
@@ -457,11 +522,17 @@ std::future<void> AbstractIOHandlerImpl::flush()
                 auto &parameter = deref_dynamic_cast<Parameter<O::SET_WRITTEN>>(
                     i.parameter.get());
                 writeToStderr(
-                    "[", i.writable->parent, "->", i.writable, "] SET_WRITTEN");
+                    "[",
+                    i.writable->parent,
+                    "->",
+                    i.writable,
+                    "] SET_WRITTEN ",
+                    parameter.target_status ? "true" : "false");
                 setWritten(i.writable, parameter);
                 break;
             }
             }
+            verifyFlushType(i.operation, l);
         }
         catch (...)
         {
@@ -506,6 +577,7 @@ std::future<void> AbstractIOHandlerImpl::flush()
         }
         (*m_handler).m_work.pop();
     }
+    writeToStderr("FLUSHED\n");
     return std::future<void>();
 }
 
diff --git a/src/IO/HDF5/HDF5IOHandler.cpp b/src/IO/HDF5/HDF5IOHandler.cpp
index f300a150cd..4b6d0876e2 100644
--- a/src/IO/HDF5/HDF5IOHandler.cpp
+++ b/src/IO/HDF5/HDF5IOHandler.cpp
@@ -3542,7 +3542,7 @@ auto HDF5IOHandlerImpl::requireFile(
 
 std::future<void> HDF5IOHandlerImpl::flush(internal::ParsedFlushParams &params)
 {
-    auto res = AbstractIOHandlerImpl::flush();
+    auto res = AbstractIOHandlerImpl::flush(params.flushLevel);
 
     if (params.backendConfig.json().contains("hdf5"))
     {
diff --git a/src/IO/IOTask.cpp b/src/IO/IOTask.cpp
index 26010ead52..af12ff766c 100644
--- a/src/IO/IOTask.cpp
+++ b/src/IO/IOTask.cpp
@@ -35,6 +35,97 @@ Writable *getWritable(Attributable *a)
     return &a->writable();
 }
 
+std::ostream &operator<<(std::ostream &os, Operation op)
+{
+    switch (op)
+    {
+    case Operation::CREATE_FILE:
+        os << "CREATE_FILE";
+        break;
+    case Operation::CHECK_FILE:
+        os << "CHECK_FILE";
+        break;
+    case Operation::OPEN_FILE:
+        os << "OPEN_FILE";
+        break;
+    case Operation::CLOSE_FILE:
+        os << "CLOSE_FILE";
+        break;
+    case Operation::DELETE_FILE:
+        os << "DELETE_FILE";
+        break;
+    case Operation::CREATE_PATH:
+        os << "CREATE_PATH";
+        break;
+    case Operation::CLOSE_PATH:
+        os << "CLOSE_PATH";
+        break;
+    case Operation::OPEN_PATH:
+        os << "OPEN_PATH";
+        break;
+    case Operation::DELETE_PATH:
+        os << "DELETE_PATH";
+        break;
+    case Operation::LIST_PATHS:
+        os << "LIST_PATHS";
+        break;
+    case Operation::CREATE_DATASET:
+        os << "CREATE_DATASET";
+        break;
+    case Operation::EXTEND_DATASET:
+        os << "EXTEND_DATASET";
+        break;
+    case Operation::OPEN_DATASET:
+        os << "OPEN_DATASET";
+        break;
+    case Operation::DELETE_DATASET:
+        os << "DELETE_DATASET";
+        break;
+    case Operation::WRITE_DATASET:
+        os << "WRITE_DATASET";
+        break;
+    case Operation::READ_DATASET:
+        os << "READ_DATASET";
+        break;
+    case Operation::LIST_DATASETS:
+        os << "LIST_DATASETS";
+        break;
+    case Operation::GET_BUFFER_VIEW:
+        os << "GET_BUFFER_VIEW";
+        break;
+    case Operation::DELETE_ATT:
+        os << "DELETE_ATT";
+        break;
+    case Operation::WRITE_ATT:
+        os << "WRITE_ATT";
+        break;
+    case Operation::READ_ATT:
+        os << "READ_ATT";
+        break;
+    case Operation::READ_ATT_ALLSTEPS:
+        os << "READ_ATT_ALLSTEPS";
+        break;
+    case Operation::LIST_ATTS:
+        os << "LIST_ATTS";
+        break;
+    case Operation::ADVANCE:
+        os << "ADVANCE";
+        break;
+    case Operation::AVAILABLE_CHUNKS:
+        os << "AVAILABLE_CHUNKS";
+        break;
+    case Operation::DEREGISTER:
+        os << "DEREGISTER";
+        break;
+    case Operation::TOUCH:
+        os << "TOUCH";
+        break;
+    case Operation::SET_WRITTEN:
+        os << "SET_WRITTEN";
+        break;
+    }
+    return os;
+}
 template <>
 void AbstractParameter::warnUnusedParameters<json::TracingJSON>(
     json::TracingJSON &config,
diff --git a/src/IO/JSON/JSONIOHandler.cpp b/src/IO/JSON/JSONIOHandler.cpp
index c531aabb00..9af9a06728 100644
--- a/src/IO/JSON/JSONIOHandler.cpp
+++ b/src/IO/JSON/JSONIOHandler.cpp
@@ -53,8 +53,8 @@ JSONIOHandler::JSONIOHandler(
 {}
 #endif
 
-std::future<void> JSONIOHandler::flush(internal::ParsedFlushParams &)
+std::future<void> JSONIOHandler::flush(internal::ParsedFlushParams &params)
 {
-    return m_impl.flush();
+    return m_impl.flush(params);
 }
 } // namespace openPMD
diff --git a/src/IO/JSON/JSONIOHandlerImpl.cpp b/src/IO/JSON/JSONIOHandlerImpl.cpp
index 551b6c4358..91642dfe18 100644
--- a/src/IO/JSON/JSONIOHandlerImpl.cpp
+++ b/src/IO/JSON/JSONIOHandlerImpl.cpp
@@ -24,6 +24,7 @@
 #include "openPMD/Error.hpp"
 #include "openPMD/IO/AbstractIOHandler.hpp"
 #include "openPMD/IO/AbstractIOHandlerImpl.hpp"
+#include "openPMD/IO/FlushParametersInternal.hpp"
 #include "openPMD/ThrowError.hpp"
 #include "openPMD/auxiliary/Filesystem.hpp"
 #include "openPMD/auxiliary/JSONMatcher.hpp"
@@ -444,9 +445,9 @@ void JSONIOHandlerImpl::init(openPMD::json::TracingJSON config)
 
 JSONIOHandlerImpl::~JSONIOHandlerImpl() = default;
 
-std::future<void> JSONIOHandlerImpl::flush()
+std::future<void> JSONIOHandlerImpl::flush(internal::ParsedFlushParams &params)
 {
-    AbstractIOHandlerImpl::flush();
+    AbstractIOHandlerImpl::flush(params.flushLevel);
     if (access::readOnly(m_handler->m_backendAccess) && !m_dirty.empty())
     {
         throw error::Internal(
@@ -1538,6 +1539,15 @@ void JSONIOHandlerImpl::readAttribute(
     refreshFileFromParent(writable);
     auto name = removeSlashes(parameters.name);
     auto const &jsonContents = obtainJsonContents(writable);
+    if (!hasKey(jsonContents, "attributes"))
+    {
+        throw error::ReadError(
+            error::AffectedObject::Attribute,
+            error::Reason::NotFound,
+            "JSON",
+            "Tried looking up attribute '" + name + "' at location: " +
+                setAndGetFilePosition(writable, false)->id.to_string());
+    }
     auto const &jsonLoc = jsonContents["attributes"];
     setAndGetFilePosition(writable);
     if (!hasKey(jsonLoc, name))
diff --git a/src/Iteration.cpp b/src/Iteration.cpp
index 085495db93..5b41b78839 100644
--- a/src/Iteration.cpp
+++ b/src/Iteration.cpp
@@ -30,10 +30,14 @@
 #include "openPMD/auxiliary/DerefDynamicCast.hpp"
 #include "openPMD/auxiliary/Filesystem.hpp"
 #include "openPMD/auxiliary/StringManip.hpp"
+#include "openPMD/auxiliary/TypeTraits.hpp"
 #include "openPMD/auxiliary/Variant.hpp"
 #include "openPMD/backend/Attributable.hpp"
+#include "openPMD/backend/BaseRecordComponent.hpp"
 #include "openPMD/backend/Variant_internal.hpp"
 #include "openPMD/backend/Writable.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <algorithm>
 #include <exception>
@@ -42,6 +46,7 @@
 #include <optional>
 #include <stdexcept>
 #include <tuple>
+#include <type_traits>
 #include <variant>
 
 namespace openPMD
@@ -49,12 +54,19 @@ namespace openPMD
 using internal::CloseStatus;
 using internal::DeferredParseAccess;
 
+void Meshes::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    visitHierarchyImpl<Meshes>(v, recursive);
+}
+
+void Particles::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    visitHierarchyImpl<Particles>(v, recursive);
+}
+
 Iteration::Iteration() : Attributable(NoInit())
 {
     setData(std::make_shared<Data_t>());
-    setTime(static_cast<double>(0));
-    setDt(static_cast<double>(1));
-    setTimeUnitSI(1);
     meshes.writable().ownKeyWithinParent = "meshes";
     particles.writable().ownKeyWithinParent = "particles";
 }
@@ -122,6 +134,12 @@ Iteration &Iteration::close(bool _flush)
         // yet keeps it re-openable)
         break;
     }
+
+    // if (access::write(IOHandler()->m_frontendAccess))
+    // {
+    //     commitStructuralSetup();
+    // }
+
     if (_flush)
     {
         if (flag == StepStatus::DuringStep)
@@ -236,6 +254,13 @@ bool Iteration::closedByWriter() const
     }
 }
 
+void Iteration::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    meshes.visitHierarchy(v, recursive);
+    particles.visitHierarchy(v, recursive);
+    v(*this);
+}
+
 void Iteration::flushFileBased(
     std::string const &filename,
     IterationIndex_t i,
@@ -244,6 +269,10 @@ void Iteration::flushFileBased(
     /* Find the root point [Series] of this file,
      * meshesPath and particlesPath are stored there */
     Series s = retrieveSeries();
+    auto &series = s.get();
+
+    bool do_flush_rank_table =
+        !series.m_rankTable.m_attributable.written() || !this->written();
 
     if (!written())
     {
@@ -252,16 +281,6 @@ void Iteration::flushFileBased(
         fCreate.name = filename;
         IOHandler()->enqueue(IOTask(&s.writable(), fCreate));
 
-        /*
-         * If it was written before, then in the context of another iteration.
-         */
-        auto &attr = s.get().m_rankTable.m_attributable;
-        attr.setWritten(false, Attributable::EnqueueAsynchronously::Both);
-        s.get()
-            .m_rankTable.m_attributable.get()
-            .m_writable.abstractFilePosition.reset();
-        s.flushRankTable();
-
         /* create basePath */
         Parameter<Operation::CREATE_PATH> pCreate;
         pCreate.path = auxiliary::replace_first(s.basePath(), "%T/", "");
@@ -278,15 +297,21 @@ void Iteration::flushFileBased(
         s.openIteration(i, *this);
     }
 
-    switch (flushParams.flushLevel)
+    if (do_flush_rank_table)
+    {
+        /*
+         * If it was written before, then in the context of another iteration.
+         */
+        auto &attr = series.m_rankTable.m_attributable;
+        attr.setWritten(false, Attributable::EnqueueAsynchronously::Both);
+        attr.get().m_writable.abstractFilePosition.reset();
+
+        s.flushRankTable(flushParams.flushLevel);
+    }
+
+    if (flush_level::flush_hierarchy(flushParams.flushLevel))
     {
-    case FlushLevel::CreateOrOpenFiles:
-        break;
-    case FlushLevel::SkeletonOnly:
-    case FlushLevel::InternalFlush:
-    case FlushLevel::UserFlush:
         flush(flushParams);
-        break;
     }
 }
 
@@ -301,15 +326,9 @@ void Iteration::flushGroupBased(
         IOHandler()->enqueue(IOTask(this, pCreate));
     }
 
-    switch (flushParams.flushLevel)
+    if (flush_level::flush_hierarchy(flushParams.flushLevel))
     {
-    case FlushLevel::CreateOrOpenFiles:
-        break;
-    case FlushLevel::SkeletonOnly:
-    case FlushLevel::InternalFlush:
-    case FlushLevel::UserFlush:
         flush(flushParams);
-        break;
     }
 }
 
@@ -324,18 +343,14 @@ void Iteration::flushVariableBased(
         IOHandler()->enqueue(IOTask(this, pOpen));
     }
 
-    switch (flushParams.flushLevel)
+    if (!flush_level::flush_hierarchy(flushParams.flushLevel))
     {
-    case FlushLevel::CreateOrOpenFiles:
         return;
-    case FlushLevel::SkeletonOnly:
-    case FlushLevel::InternalFlush:
-    case FlushLevel::UserFlush:
-        flush(flushParams);
-        break;
     }
 
-    if (!written())
+    flush(flushParams);
+
+    if (!written() && flush_level::write_datasets(flushParams.flushLevel))
     {
         /* create iteration path */
         Parameter<Operation::OPEN_PATH> pOpen;
@@ -367,7 +382,7 @@ void Iteration::flush(internal::FlushParams const &flushParams)
             m.second.flush(m.first, flushParams);
         for (auto &species : particles)
             species.second.flush(species.first, flushParams);
-        setDirty(false);
+        determineUnsetDirty(flushParams.flushLevel);
     }
     else
     {
@@ -375,16 +390,28 @@ void Iteration::flush(internal::FlushParams const &flushParams)
          * meshesPath and particlesPath are stored there */
         Series s = retrieveSeries();
 
-        if (!meshes.empty() || s.containsAttribute("meshesPath"))
-        {
-            if (!s.containsAttribute("meshesPath"))
+        auto set_and_get_mp_path =
+            [&](char const *attrName,
+                char const *defaultVal,
+                Series &(Series::*set)(std::string const &)) -> std::string {
+            if (s.containsAttribute(attrName))
+            {
+                return s.getAttribute(attrName).get<std::string>();
+            }
+            else
             {
-                s.setMeshesPath("meshes/");
-                s.flushMeshesPath();
+                (s.*set)(defaultVal);
+                return defaultVal;
             }
+        };
+
+        if (!meshes.empty() || s.containsAttribute("meshesPath"))
+        {
+            auto meshesPath = set_and_get_mp_path(
+                "meshesPath", "meshes/", &Series::setMeshesPath);
             if (meshes.dirtyRecursive())
             {
-                meshes.flush(s.meshesPath(), flushParams);
+                meshes.flush(meshesPath, flushParams);
                 for (auto &m : meshes)
                 {
                     m.second.flush(m.first, flushParams);
@@ -398,14 +425,11 @@ void Iteration::flush(internal::FlushParams const &flushParams)
 
         if (!particles.empty() || s.containsAttribute("particlesPath"))
         {
-            if (!s.containsAttribute("particlesPath"))
-            {
-                s.setParticlesPath("particles/");
-                s.flushParticlesPath();
-            }
+            auto particlesPath = set_and_get_mp_path(
+                "particlesPath", "particles/", &Series::setParticlesPath);
             if (particles.dirtyRecursive())
             {
-                particles.flush(s.particlesPath(), flushParams);
+                particles.flush(particlesPath, flushParams);
                 for (auto &species : particles)
                 {
                     species.second.flush(species.first, flushParams);
@@ -417,13 +441,14 @@ void Iteration::flush(internal::FlushParams const &flushParams)
             particles.setDirty(false);
         }
 
+        populateDefaults(this, flushParams);
         flushAttributes(flushParams);
     }
     if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
     {
-        setDirty(false);
-        meshes.setDirty(false);
-        particles.setDirty(false);
+        determineUnsetDirty(flushParams.flushLevel);
+        meshes.determineUnsetDirty(flushParams.flushLevel);
+        particles.determineUnsetDirty(flushParams.flushLevel);
     }
 }
 
@@ -509,80 +534,7 @@ void Iteration::read_impl(std::string const &groupPath)
     pOpen.path = groupPath;
     IOHandler()->enqueue(IOTask(this, pOpen));
 
-    using DT = Datatype;
-    Parameter<Operation::READ_ATT> aRead;
-
-    aRead.name = "dt";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (*aRead.dtype == DT::FLOAT)
-        setDt(Attribute(Attribute::from_any, *aRead.m_resource).get<float>());
-    else if (*aRead.dtype == DT::DOUBLE)
-        setDt(Attribute(Attribute::from_any, *aRead.m_resource).get<double>());
-    else if (*aRead.dtype == DT::LONG_DOUBLE)
-        setDt(Attribute(Attribute::from_any, *aRead.m_resource)
-                  .get<long double>());
-    // conversion cast if a backend reports an integer type
-    else if (
-        auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<double>();
-        val.has_value())
-        setDt(val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'dt' (expected double, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "time";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (*aRead.dtype == DT::FLOAT)
-        setTime(Attribute(Attribute::from_any, *aRead.m_resource).get<float>());
-    else if (*aRead.dtype == DT::DOUBLE)
-        setTime(
-            Attribute(Attribute::from_any, *aRead.m_resource).get<double>());
-    else if (*aRead.dtype == DT::LONG_DOUBLE)
-        setTime(Attribute(Attribute::from_any, *aRead.m_resource)
-                    .get<long double>());
-    // conversion cast if a backend reports an integer type
-    else if (
-        auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<double>();
-        val.has_value())
-        setTime(val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'time' (expected double, "
-            "found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "timeUnitSI";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<double>();
-        val.has_value())
-        setTimeUnitSI(val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'timeUnitSI' (expected double, "
-            "found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 
     /* Find the root point [Series] of this file,
      * meshesPath and particlesPath are stored there */
@@ -1013,6 +965,23 @@ void Iteration::runDeferredParseAccess()
     }
 }
 
+void Iteration::scientificDefaults_impl(bool write, OpenpmdStandard)
+{
+    using namespace internal;
+    auto float_types = get_float_types();
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    defaultAttribute(*this, "time")
+        .template withSetter<Iteration>(0., &Iteration::setTime)
+        .withReader(float_types, require_scalar)(wor);
+    defaultAttribute(*this, "dt")
+        .template withSetter<Iteration>(1., &Iteration::setDt)
+        .withReader(float_types, require_scalar)(wor);
+    defaultAttribute(*this, "timeUnitSI")
+        .template withSetter<Iteration>(1.0, &Iteration::setTimeUnitSI)
+        .withReader(float_types, require_type<double>())(wor);
+}
+
 template float Iteration::time<float>() const;
 template double Iteration::time<double>() const;
 template long double Iteration::time<long double>() const;
diff --git a/src/Mesh.cpp b/src/Mesh.cpp
index e93b2b5219..bf623316a7 100644
--- a/src/Mesh.cpp
+++ b/src/Mesh.cpp
@@ -29,24 +29,21 @@
 #include "openPMD/auxiliary/StringManip.hpp"
 #include "openPMD/backend/Attribute.hpp"
 #include "openPMD/backend/Writable.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <algorithm>
 #include <iostream>
 
 namespace openPMD
 {
-Mesh::Mesh()
+void Mesh::visitHierarchy(HierarchyVisitor &v, bool recursive)
 {
-    setTimeOffset(0.f);
-
-    setGeometry(Geometry::cartesian);
-    setDataOrder(DataOrder::C);
-
-    setAxisLabels({"x"}); // empty strings are not allowed in HDF5
-    setGridSpacing(std::vector<double>{1});
-    setGridGlobalOffset({0});
+    visitHierarchyImpl<Mesh>(v, recursive);
 }
 
+Mesh::Mesh() = default;
+
 Mesh::Geometry Mesh::geometry() const
 {
     std::string ret = geometryString();
@@ -204,28 +201,91 @@ Mesh &Mesh::setGridUnitSI(std::vector<double> const &gusi)
     return setGridUnitSIPerDimension(gusi);
 }
 
-namespace
+auto Mesh::retrieveDimensionality() const -> uint64_t
 {
-    uint64_t retrieveMeshDimensionality(Mesh const &m)
+    if (containsAttribute("axisLabels"))
     {
-        if (m.containsAttribute("axisLabels"))
-        {
-            return m.axisLabels().size();
-        }
+        return axisLabels().size();
+    }
 
-        // maybe we have record components and can ask them
-        if (auto it = m.begin(); it != m.end())
-        {
-            return it->second.getDimensionality();
-        }
-        /*
-         * Since some backends cannot distinguish between vector and
-         * scalar values, the most likely answer here is 1.
-         */
-        return 1;
+    // maybe we have record components and can ask them
+    if (auto it = begin(); it != end())
+    {
+        return it->second.getDimensionality();
     }
-} // namespace
+    /*
+     * Since some backends cannot distinguish between vector and
+     * scalar values, the most likely answer here is 1.
+     */
+    return 1;
+}
 
+void Mesh::scientificDefaults_impl(bool write, OpenpmdStandard standard)
+{
+    using namespace internal;
+    auto float_types = get_float_types();
+    auto string_types = get_string_types();
+    auto dimensionality = retrieveDimensionality();
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    defaultAttribute(*this, "geometry")
+        .template withSetter<Mesh>(
+            Mesh::Geometry::cartesian, &Mesh::setGeometry)
+        .withReader(
+            string_types, require_type(*this, &setMeshGeometryFromString))(wor);
+
+    defaultAttribute(*this, "dataOrder")
+        .template withSetter<Mesh>(Mesh::DataOrder::C, &Mesh::setDataOrder)
+        .withReader(
+            string_types, require_type(*this, &setMeshDataOrderFromChar))(wor);
+
+    defaultAttribute(*this, "axisLabels")
+        .template withSetter<Mesh, std::vector<std::string> const &>(
+            [&]() -> std::vector<std::string> {
+                return auxiliary::createDefaultAxisLabels(dimensionality);
+            },
+            &Mesh::setAxisLabels)
+        .withReader(string_types, require_vector)(wor);
+
+    defaultAttribute(*this, "gridSpacing")
+        .template withSetter<Mesh, std::vector<double> const &>(
+            [&]() {
+                return auxiliary::createDefaultVector(dimensionality, 1.0);
+            },
+            &Mesh::setGridSpacing)
+        .withReader(float_types, require_vector)(wor);
+
+    defaultAttribute(*this, "gridGlobalOffset")
+        .template withSetter<Mesh, std::vector<double> const &>(
+            [&]() {
+                return auxiliary::createDefaultVector(dimensionality, 0.0);
+            },
+            &Mesh::setGridGlobalOffset)
+        .withReader(float_types, require_type<std::vector<double>>())(wor);
+
+    defaultAttribute(*this, "timeOffset")
+        .template withSetter<Mesh>(0.f, &Mesh::setTimeOffset)
+        .withReader(float_types, require_scalar)(wor);
+
+    if (standard >= OpenpmdStandard::v_2_0_0)
+    {
+        defaultAttribute(*this, "gridUnitSI")
+            .template withSetter<Mesh, std::vector<double> const &>(
+                [&]() {
+                    return auxiliary::createDefaultVector(dimensionality, 1.);
+                },
+                &Mesh::setGridUnitSIPerDimension)
+            .withReader(float_types, require_type<std::vector<double>>())(wor);
+    }
+    else
+    {
+        defaultAttribute(*this, "gridUnitSI")
+            .template withSetter<Mesh>(1.0, &Mesh::setGridUnitSI)
+            .withReader(float_types, require_type<std::vector<double>>())(wor);
+    }
+
+    BaseRecord<MeshRecordComponent>::scientificDefaults_impl(write, standard);
+}
 std::vector<double> Mesh::gridUnitSIPerDimension() const
 {
     if (containsAttribute("gridUnitSI"))
@@ -235,7 +295,7 @@ std::vector<double> Mesh::gridUnitSIPerDimension() const
             // If the openPMD version is lower than 2.0, the gridUnitSI is a
             // scalar interpreted for all axes. Copy it d times.
             return std::vector<double>(
-                retrieveMeshDimensionality(*this),
+                retrieveDimensionality(),
                 getAttribute("gridUnitSI").get<double>());
         }
         return getAttribute("gridUnitSI").get<std::vector<double>>();
@@ -244,7 +304,7 @@ std::vector<double> Mesh::gridUnitSIPerDimension() const
     {
         // gridUnitSI is an optional attribute
         // if it is missing, the mesh is interpreted as unscaled
-        return std::vector<double>(retrieveMeshDimensionality(*this), 1.);
+        return std::vector<double>(retrieveDimensionality(), 1.);
     }
 }
 
@@ -269,7 +329,9 @@ Mesh &Mesh::setUnitDimension(unit_representations::AsMap const &udim)
 {
     if (!udim.empty())
     {
-        std::array<double, 7> tmpUnitDimension = this->unitDimension();
+        std::array<double, 7> tmpUnitDimension =
+            this->containsAttribute("unitDimension") ? this->unitDimension()
+                                                     : std::array<double, 7>{0};
         unit_representations::auxiliary::fromMapOfUnitDimension(
             tmpUnitDimension.data(), udim);
         setAttribute("unitDimension", tmpUnitDimension);
@@ -285,15 +347,20 @@ Mesh &Mesh::setUnitDimension(unit_representations::AsArray const &udim)
 
 Mesh &Mesh::setGridUnitDimension(unit_representations::AsMaps const &udims)
 {
-    auto rawGridUnitDimension = [this]() {
-        try
+    auto rawGridUnitDimension = [&udims, this]() {
+        if (!this->containsAttribute("gridUnitDimension"))
         {
-            return this->getAttribute("gridUnitDimension")
-                .get<std::vector<double>>();
+            std::vector<double> res(udims.size() * 7);
+            // for (size_t i = 0; i < udims.size(); ++i)
+            // {
+            //     res[7 * i] = 1;
+            // }
+            return res;
         }
-        catch (no_such_attribute_error const &)
+        else
         {
-            return std::vector<double>();
+            return this->getAttribute("gridUnitDimension")
+                .get<std::vector<double>>();
         }
     }();
     rawGridUnitDimension.resize(7 * udims.size());
@@ -341,7 +408,7 @@ unit_representations::AsArrays Mesh::gridUnitDimension() const
         // if it is missing, the mesh is interpreted as spatial
         auto spatialMesh =
             unit_representations::asArray({{UnitDimension::L, 1}});
-        auto dim = retrieveMeshDimensionality(*this);
+        auto dim = retrieveDimensionality();
         unit_representations::AsArrays res(dim, spatialMesh);
         return res;
     }
@@ -417,18 +484,6 @@ void Mesh::flush_impl(
                     comp.second.flush(comp.first, flushParams);
             }
         }
-        if (!containsAttribute("gridUnitSI"))
-        {
-            if (IOHandler()->m_standard < OpenpmdStandard::v_2_0_0)
-            {
-                setGridUnitSI(1);
-            }
-            else
-            {
-                setGridUnitSIPerDimension(
-                    std::vector<double>(retrieveMeshDimensionality(*this), 1));
-            }
-        }
         flushAttributes(flushParams);
     }
 }
@@ -439,170 +494,7 @@ void Mesh::read()
         IOHandler()->m_verify_homogeneous_extents);
     internal::EraseStaleEntries<Mesh> map{*this};
 
-    using DT = Datatype;
-    Parameter<Operation::READ_ATT> aRead;
-
-    aRead.name = "geometry";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (*aRead.dtype == DT::STRING)
-    {
-        std::string tmpGeometry =
-            Attribute(Attribute::from_any, *aRead.m_resource)
-                .get<std::string>();
-        if ("cartesian" == tmpGeometry)
-            setGeometry(Geometry::cartesian);
-        else if ("thetaMode" == tmpGeometry)
-            setGeometry(Geometry::thetaMode);
-        else if ("cylindrical" == tmpGeometry)
-            setGeometry(Geometry::cylindrical);
-        else if ("spherical" == tmpGeometry)
-            setGeometry(Geometry::spherical);
-        else
-            setGeometry(tmpGeometry);
-    }
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'geometry' (expected a string, "
-            "found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "dataOrder";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (*aRead.dtype == DT::CHAR)
-        setDataOrder(
-            static_cast<DataOrder>(
-                Attribute(Attribute::from_any, *aRead.m_resource).get<char>()));
-    else if (*aRead.dtype == DT::STRING)
-    {
-        std::string tmpDataOrder =
-            Attribute(Attribute::from_any, *aRead.m_resource)
-                .get<std::string>();
-        if (tmpDataOrder.size() == 1)
-            setDataOrder(static_cast<DataOrder>(tmpDataOrder[0]));
-        else
-            throw error::ReadError(
-                error::AffectedObject::Attribute,
-                error::Reason::UnexpectedContent,
-                {},
-                "Unexpected Attribute value for 'dataOrder': " + tmpDataOrder);
-    }
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'dataOrder' (expected char or "
-            "string, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "axisLabels";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    Attribute a = Attribute(Attribute::from_any, *aRead.m_resource);
-    if (auto val = a.getOptional<std::vector<std::string>>(); val.has_value())
-        setAxisLabels(*val);
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'axisLabels' (expected a vector "
-            "of string, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "gridSpacing";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    a = Attribute(Attribute::from_any, *aRead.m_resource);
-    if (*aRead.dtype == DT::VEC_FLOAT || *aRead.dtype == DT::FLOAT)
-        setGridSpacing(a.get<std::vector<float>>());
-    else if (*aRead.dtype == DT::VEC_DOUBLE || *aRead.dtype == DT::DOUBLE)
-        setGridSpacing(a.get<std::vector<double>>());
-    else if (
-        *aRead.dtype == DT::VEC_LONG_DOUBLE || *aRead.dtype == DT::LONG_DOUBLE)
-        setGridSpacing(a.get<std::vector<long double>>());
-    // conversion cast if a backend reports an integer type
-    else if (auto val = a.getOptional<std::vector<double>>(); val.has_value())
-        setGridSpacing(val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'gridSpacing' (expected a "
-            "vector of double, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "gridGlobalOffset";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<std::vector<double>>();
-        val.has_value())
-        setGridGlobalOffset(val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'gridGlobalOffset' (expected a "
-            "vector of double, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    aRead.name = "gridUnitSI";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    if (IOHandler()->m_standard >= OpenpmdStandard::v_2_0_0)
-    {
-        if (auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                           .getOptional<std::vector<double>>();
-            val.has_value())
-            setGridUnitSIPerDimension(val.value());
-        else
-            throw error::ReadError(
-                error::AffectedObject::Attribute,
-                error::Reason::UnexpectedContent,
-                {},
-                "Unexpected Attribute datatype for 'gridUnitSI' "
-                "(expected vector of double, found " +
-                    datatypeToString(
-                        Attribute(Attribute::from_any, *aRead.m_resource)
-                            .dtype) +
-                    ")");
-    }
-    else
-    {
-        if (auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                           .getOptional<double>();
-            val.has_value())
-            setGridUnitSI(val.value());
-        else
-            throw error::ReadError(
-                error::AffectedObject::Attribute,
-                error::Reason::UnexpectedContent,
-                {},
-                "Unexpected Attribute datatype for 'gridUnitSI' "
-                "(expected double, found " +
-                    datatypeToString(
-                        Attribute(Attribute::from_any, *aRead.m_resource)
-                            .dtype) +
-                    ")");
-    }
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 
     if (scalar())
     {
@@ -669,7 +561,7 @@ void Mesh::read()
 
     std::move(homogenizeExtents).homogenize(*this);
 
-    readBase();
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 
     readAttributes(ReadMode::FullyReread);
 }
diff --git a/src/ParticlePatches.cpp b/src/ParticlePatches.cpp
index 789559a67f..e27208e596 100644
--- a/src/ParticlePatches.cpp
+++ b/src/ParticlePatches.cpp
@@ -35,6 +35,11 @@ size_t ParticlePatches::numPatches() const
     return this->at("numParticles").getExtent()[0];
 }
 
+void ParticlePatches::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    visitHierarchyImpl<ParticlePatches>(v, recursive);
+}
+
 void ParticlePatches::read()
 {
     Parameter<Operation::LIST_PATHS> pList;
@@ -100,7 +105,7 @@ void ParticlePatches::read()
         pr.setDirty(false);
         try
         {
-            prc.PatchRecordComponent::read(/* require_unit_si = */ false);
+            prc.PatchRecordComponent::read();
         }
         catch (error::ReadError const &err)
         {
diff --git a/src/ParticleSpecies.cpp b/src/ParticleSpecies.cpp
index b3e7ea1f7e..3662ef3ff6 100644
--- a/src/ParticleSpecies.cpp
+++ b/src/ParticleSpecies.cpp
@@ -20,8 +20,6 @@
  */
 #include "openPMD/ParticleSpecies.hpp"
 #include "openPMD/RecordComponent.hpp"
-#include "openPMD/Series.hpp"
-#include "openPMD/auxiliary/DerefDynamicCast.hpp"
 #include "openPMD/backend/Writable.hpp"
 
 #include <algorithm>
@@ -29,6 +27,12 @@
 
 namespace openPMD
 {
+void ParticleSpecies::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    visitHierarchyImpl<ParticleSpecies>(v, recursive);
+    particlePatches.visitHierarchy(v, recursive);
+}
+
 ParticleSpecies::ParticleSpecies()
 {
     particlePatches.writable().ownKeyWithinParent = "particlePatches";
@@ -175,23 +179,10 @@ void ParticleSpecies::flush(
             record.second.flush(record.first, flushParams);
         for (auto &patch : particlePatches)
             patch.second.flush(patch.first, flushParams);
-        if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
-        {
-            particlePatches.setDirty(false);
-        }
     }
     else
     {
-        if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
-        {
-            auto it = find("position");
-            if (it != end())
-                it->second.setUnitDimension({{UnitDimension::L, 1}});
-            it = find("positionOffset");
-            if (it != end())
-                it->second.setUnitDimension({{UnitDimension::L, 1}});
-        }
-
+        populateDefaults(this, flushParams);
         Container<Record>::flush(path, flushParams);
 
         for (auto &record : *this)
@@ -203,14 +194,12 @@ void ParticleSpecies::flush(
             for (auto &patch : particlePatches)
                 patch.second.flush(patch.first, flushParams);
         }
-        else
-        {
-            particlePatches.setDirty(false);
-        }
-    }
-    if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
-    {
-        setDirty(false);
     }
+    determineUnsetDirty(flushParams.flushLevel);
+    particlePatches.determineUnsetDirty(flushParams.flushLevel);
+}
+void ParticleSpecies::scientificDefaults_impl(bool /* write */, OpenpmdStandard)
+{
+    // no-op
 }
 } // namespace openPMD
diff --git a/src/Record.cpp b/src/Record.cpp
index 8056718edd..1cb480b4cd 100644
--- a/src/Record.cpp
+++ b/src/Record.cpp
@@ -24,21 +24,26 @@
 #include "openPMD/ThrowError.hpp"
 #include "openPMD/UnitDimension.hpp"
 #include "openPMD/backend/BaseRecord.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <iostream>
 
 namespace openPMD
 {
-Record::Record()
+void Record::visitHierarchy(HierarchyVisitor &v, bool recursive)
 {
-    setTimeOffset(0.f);
+    visitHierarchyImpl<Record>(v, recursive);
 }
+Record::Record() = default;
 
 Record &Record::setUnitDimension(unit_representations::AsMap const &udim)
 {
     if (!udim.empty())
     {
-        std::array<double, 7> tmpUnitDimension = this->unitDimension();
+        std::array<double, 7> tmpUnitDimension =
+            this->containsAttribute("unitDimension") ? this->unitDimension()
+                                                     : std::array<double, 7>{};
         for (auto const &entry : udim)
             tmpUnitDimension[static_cast<uint8_t>(entry.first)] = entry.second;
         setAttribute("unitDimension", tmpUnitDimension);
@@ -121,7 +126,7 @@ auto Record::read() -> internal::HomogenizeExtents
             /* using operator[] will incorrectly update parent */
             try
             {
-                T_RecordComponent::read(/* require_unit_si = */ true);
+                T_RecordComponent::read();
             }
             catch (error::ReadError const &err)
             {
@@ -149,7 +154,7 @@ auto Record::read() -> internal::HomogenizeExtents
             rc.get().m_isConstant = true;
             try
             {
-                rc.read(/* require_unit_si = */ true);
+                rc.read();
             }
             catch (error::ReadError const &err)
             {
@@ -178,7 +183,7 @@ auto Record::read() -> internal::HomogenizeExtents
             rc.setWritten(true, Attributable::EnqueueAsynchronously::No);
             try
             {
-                rc.read(/* require_unit_si = */ true);
+                rc.read();
             }
             catch (error::ReadError const &err)
             {
@@ -192,11 +197,38 @@ auto Record::read() -> internal::HomogenizeExtents
         }
     }
 
-    readBase();
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 
     readAttributes(ReadMode::FullyReread);
     return res;
 }
 
+void Record::scientificDefaults_impl(bool write, OpenpmdStandard standard)
+{
+    using namespace internal;
+    auto float_types = get_float_types();
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    defaultAttribute(*this, "timeOffset")
+        .template withSetter<Record>(0.f, &Record::setTimeOffset)
+        .withReader(float_types, require_scalar)(wor);
+
+    auto const &keyInParent = writable().ownKeyWithinParent;
+    if (keyInParent == "position" || keyInParent == "positionOffset")
+    {
+        defaultAttribute(*this, "unitDimension")
+            .template withSetter<Record, unit_representations::AsMap const &>(
+                []() {
+                    return unit_representations::AsMap{{UnitDimension::L, 1.0}};
+                },
+                &Record::setUnitDimension)(wor);
+    }
+
+    defaultAttribute(*this, "timeOffset")
+        .template withSetter<Record>(0.f, &Record::setTimeOffset)
+        .withReader(float_types, require_scalar)(wor);
+
+    BaseRecord<RecordComponent>::scientificDefaults_impl(write, standard);
+}
 template class BaseRecord<RecordComponent>;
 } // namespace openPMD
diff --git a/src/RecordComponent.cpp b/src/RecordComponent.cpp
index fbf2ad0d88..7cf1c4afdc 100644
--- a/src/RecordComponent.cpp
+++ b/src/RecordComponent.cpp
@@ -22,6 +22,7 @@
 #include "openPMD/Dataset.hpp"
 #include "openPMD/DatatypeHelpers.hpp"
 #include "openPMD/Error.hpp"
+#include "openPMD/IO/AbstractIOHandler.hpp"
 #include "openPMD/IO/Format.hpp"
 #include "openPMD/Series.hpp"
 #include "openPMD/auxiliary/Environment.hpp"
@@ -30,6 +31,8 @@
 #include "openPMD/backend/Attributable.hpp"
 #include "openPMD/backend/BaseRecord.hpp"
 #include "openPMD/backend/Variant_internal.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 // comment so clang-format does not move this
 #include "openPMD/DatatypeMacros.hpp"
@@ -213,9 +216,37 @@ RecordComponent &RecordComponent::setUnitSI(double usi)
     return *this;
 }
 
+namespace
+{
+    template <typename F>
+    struct defer_type
+    {
+        F functor;
+        ~defer_type()
+        {
+            std::move(functor)();
+        }
+    };
+
+    template <typename F>
+    auto defer(F &&functor) -> defer_type<std::remove_reference_t<F>>
+    {
+        return defer_type<std::remove_reference_t<F>>{std::forward<F>(functor)};
+    }
+} // namespace
+
 RecordComponent &RecordComponent::resetDataset(Dataset d)
 {
     auto &rc = get();
+    auto cleanup = defer([&rc, this]() {
+        if (rc.m_dataset.has_value() &&
+            rc.m_dataset->dtype != Datatype::UNDEFINED &&
+            IOHandler()->m_seriesStatus != internal::SeriesStatus::Parsing)
+        {
+            seriesFlush_impl</* flush_entire_series = */ false>(
+                {FlushLevel::SkeletonOnly}, /* flush_io_handler = */ true);
+        }
+    });
     if (written())
     {
         if (!rc.m_dataset.has_value())
@@ -381,6 +412,11 @@ bool RecordComponent::empty() const
     return get().m_isEmpty;
 }
 
+void RecordComponent::visitHierarchy(HierarchyVisitor &v, bool)
+{
+    v(*this);
+}
+
 void RecordComponent::flush(
     std::string const &name, internal::FlushParams const &flushParams)
 {
@@ -425,10 +461,6 @@ void RecordComponent::flush(
                     "RecordComponent::resetDataset()).");
             }
         }
-        if (!containsAttribute("unitSI"))
-        {
-            setUnitSI(1);
-        }
         auto constant_component_write_shape = [&]() {
             auto extent = getExtent();
             return !Dataset::undefinedExtent(extent) &&
@@ -518,17 +550,16 @@ void RecordComponent::flush(
             rc.m_chunks.pop();
         }
 
+        populateDefaults(this, flushParams);
         flushAttributes(flushParams);
     }
-    if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
-    {
-        setDirty(false);
-    }
+    determineUnsetDirty(flushParams.flushLevel);
 }
 
-void RecordComponent::read(bool require_unit_si)
+void RecordComponent::read()
 {
-    readBase(require_unit_si);
+    readBase();
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 }
 
 namespace
@@ -553,7 +584,7 @@ namespace
     };
 } // namespace
 
-void RecordComponent::readBase(bool require_unit_si)
+void RecordComponent::readBase()
 {
     using DT = Datatype;
     auto &rc = get();
@@ -605,32 +636,6 @@ void RecordComponent::readBase(bool require_unit_si)
     {
         read_constant();
     }
-
-    if (require_unit_si)
-    {
-        if (!containsAttribute("unitSI"))
-        {
-            throw error::ReadError(
-                error::AffectedObject::Attribute,
-                error::Reason::NotFound,
-                {},
-                "Attribute unitSI required for record components, not found in "
-                "'" +
-                    myPath().openPMDPath() + "'.");
-        }
-        if (auto attr = getAttribute("unitSI");
-            !attr.getOptional<double>().has_value())
-        {
-            throw error::ReadError(
-                error::AffectedObject::Attribute,
-                error::Reason::UnexpectedContent,
-                {},
-                "Unexpected Attribute datatype for 'unitSI' (expected double, "
-                "found " +
-                    datatypeToString(attr.dtype) + ") in '" +
-                    myPath().openPMDPath() + "'.");
-        }
-    }
 }
 
 void RecordComponent::storeChunk(
@@ -720,6 +725,16 @@ void RecordComponent::verifyChunk(
     }
 }
 
+void RecordComponent::scientificDefaults_impl(bool write, OpenpmdStandard)
+{
+    using namespace internal;
+    auto float_types = get_float_types();
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    defaultAttribute(*this, "unitSI")
+        .template withSetter<RecordComponent>(1.0, &RecordComponent::setUnitSI)
+        .withReader(float_types, require_type<double>())(wor);
+}
 namespace
 {
     struct LoadChunkVariant
diff --git a/src/Series.cpp b/src/Series.cpp
index 038bffc7cc..cf5aef8015 100644
--- a/src/Series.cpp
+++ b/src/Series.cpp
@@ -218,6 +218,11 @@ struct Series::ParsedInput
     bool verify_homogeneous_extents = true;
 }; // ParsedInput
 
+void Iterations::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    visitHierarchyImpl<Iterations>(v, recursive);
+}
+
 std::string Series::openPMD() const
 {
     return getAttribute("openPMD").get<std::string>();
@@ -405,7 +410,7 @@ chunk_assignment::RankMeta Series::rankTable([[maybe_unused]] bool collective)
         readDataset.data = get;
 
         IOHandler()->enqueue(IOTask(&rankTable.m_attributable, readDataset));
-        IOHandler()->flush(internal::defaultFlushParams);
+        IOHandler()->flush(internal::publicFlush);
     };
 
 #if openPMD_HAVE_MPI
@@ -464,8 +469,12 @@ Series &Series::setRankTable(const std::string &myRankInfo)
     return *this;
 }
 
-void Series::flushRankTable()
+void Series::flushRankTable(FlushLevel l)
 {
+    if (!flush_level::global_flushpoint(l))
+    {
+        return;
+    }
     auto &series = get();
     auto &rankTable = series.m_rankTable;
     auto maybeMyRankInfo = std::visit(
@@ -1494,6 +1503,12 @@ void Series::flushFileBased(
     case Access::APPEND_RANDOM_ACCESS:
     case Access::APPEND_LINEAR: {
         bool allDirty = dirty();
+        // In flush level SkeletonOnly, we might need to set some attributes
+        // (especially: particlesPath, meshesPath), but cannot flush them yet
+        // (as writing attributes is only permissible at higher flush levels).
+        // This flag records if the Series became dirty during this flush. If
+        // yes, we set the Series back to dirty at the end of flushing.
+        bool hasBecomeDirty = false;
         for (auto it = begin; it != end; ++it)
         {
             // Phase 1
@@ -1543,11 +1558,30 @@ void Series::flushFileBased(
             }
             /* reset the dirty bit for every iteration (i.e. file)
              * otherwise only the first iteration will have updates attributes
+             * TODO: Ideally, we would skip this in SkeletonOnly flush mode, but
+             * for some reason, this leads to hanging parallel tests..?
              */
+            if (flushParams.flushLevel == FlushLevel::SkeletonOnly)
+            {
+                if (allDirty && !dirty())
+                {
+                    throw error::Internal(
+                        "Flush mode SkeletonOnly must not unset dirty flags.");
+                }
+                hasBecomeDirty |=
+                    flushParams.flushLevel == FlushLevel::SkeletonOnly &&
+                    !allDirty && dirty();
+            }
             setDirty(allDirty);
         }
-        setDirty(false);
-
+        if (!hasBecomeDirty)
+        {
+            determineUnsetDirty(flushParams.flushLevel);
+        }
+        else
+        {
+            setDirty(true);
+        }
         // Phase 3
         if (flushIOHandler)
         {
@@ -1632,8 +1666,11 @@ void Series::flushGorVBased(
             Parameter<Operation::CREATE_FILE> fCreate;
             fCreate.name = series.m_name;
             IOHandler()->enqueue(IOTask(this, fCreate));
+        }
 
-            flushRankTable();
+        if (!series.m_rankTable.m_attributable.written())
+        {
+            flushRankTable(flushParams.flushLevel);
         }
 
         series.iterations.flush(
@@ -1688,26 +1725,6 @@ void Series::flushGorVBased(
     }
 }
 
-void Series::flushMeshesPath()
-{
-    Parameter<Operation::WRITE_ATT> aWrite;
-    aWrite.name = "meshesPath";
-    Attribute a = getAttribute("meshesPath");
-    aWrite.m_resource = a.getAny();
-    aWrite.dtype = a.dtype;
-    IOHandler()->enqueue(IOTask(this, aWrite));
-}
-
-void Series::flushParticlesPath()
-{
-    Parameter<Operation::WRITE_ATT> aWrite;
-    aWrite.name = "particlesPath";
-    Attribute a = getAttribute("particlesPath");
-    aWrite.m_resource = a.getAny();
-    aWrite.dtype = a.dtype;
-    IOHandler()->enqueue(IOTask(this, aWrite));
-}
-
 void Series::readFileBased(
     std::optional<IterationIndex_t> read_only_this_single_iteration)
 {
@@ -2700,7 +2717,7 @@ AdvanceStatus Series::advance(
 
     if (mode == AdvanceMode::ENDSTEP)
     {
-        flushStep(/* doFlush = */ false);
+        flushStep(/* doFlush = */ false, FlushLevel::UserFlush);
     }
 
     Parameter<Operation::ADVANCE> param;
@@ -2806,7 +2823,7 @@ AdvanceStatus Series::advance(AdvanceMode mode)
 
     if (mode == AdvanceMode::ENDSTEP)
     {
-        flushStep(/* doFlush = */ false);
+        flushStep(/* doFlush = */ false, FlushLevel::UserFlush);
     }
 
     Parameter<Operation::ADVANCE> param;
@@ -2829,8 +2846,12 @@ AdvanceStatus Series::advance(AdvanceMode mode)
     return *param.status;
 }
 
-void Series::flushStep(bool doFlush)
+void Series::flushStep(bool doFlush, FlushLevel l)
 {
+    if (!flush_level::write_datasets(l))
+    {
+        return;
+    }
     auto &series = get();
     if (!series.m_currentlyActiveIterations.empty() &&
         access::write(IOHandler()->m_frontendAccess))
@@ -3317,6 +3338,14 @@ namespace internal
             !(access::readOnly(IOHandler->m_frontendAccess) &&
               !(*this)->m_writable.written))
         {
+            for (auto &[_, iteration] : iterations)
+            {
+                (void)_;
+                if (!iteration.closed())
+                {
+                    iteration.close(/* flush = */ false);
+                }
+            }
             impl.flush();
             /*
              * In file-based iteration encoding, this must be triggered by
@@ -3326,7 +3355,7 @@ namespace internal
              */
             if (impl.iterationEncoding() != IterationEncoding::fileBased)
             {
-                impl.flushStep(/* doFlush = */ true);
+                impl.flushStep(/* doFlush = */ true, FlushLevel::UserFlush);
             }
         }
         // Not strictly necessary, but clear the map of iterations
@@ -3556,6 +3585,12 @@ void Series::close()
     m_attri.reset();
 }
 
+void Series::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    get().iterations.visitHierarchy(v, recursive);
+    v(*this);
+}
+
 auto Series::currentSnapshot() -> std::optional<std::vector<IterationIndex_t>>
 {
     using vec_t = std::vector<IterationIndex_t>;
diff --git a/src/auxiliary/Mpi.cpp b/src/auxiliary/Mpi.cpp
index ef899e4207..5e8379ff41 100644
--- a/src/auxiliary/Mpi.cpp
+++ b/src/auxiliary/Mpi.cpp
@@ -49,7 +49,7 @@ StringMatrix collectStringsAsMatrixTo(
         1,
         MPI_INT,
         destRank,
-        MPI_COMM_WORLD);
+        communicator);
     int maxLength = std::accumulate(
         recvcounts.begin(), recvcounts.end(), 0, [](int a, int b) {
             return std::max(a, b);
@@ -78,7 +78,7 @@ StringMatrix collectStringsAsMatrixTo(
         displs.data(),
         MPI_CHAR,
         destRank,
-        MPI_COMM_WORLD);
+        communicator);
 
     return res;
 }
@@ -95,7 +95,7 @@ std::vector<std::string> distributeStringsToAllRanks(
     int *displs = new int[size];
 
     MPI_Allgather(
-        &sendLength, 1, MPI_INT, sizesBuffer, 1, MPI_INT, MPI_COMM_WORLD);
+        &sendLength, 1, MPI_INT, sizesBuffer, 1, MPI_INT, communicator);
 
     char *namesBuffer;
     {
@@ -116,7 +116,7 @@ std::vector<std::string> distributeStringsToAllRanks(
         sizesBuffer,
         displs,
         MPI_CHAR,
-        MPI_COMM_WORLD);
+        communicator);
 
     std::vector<std::string> hostnames(size);
     for (int i = 0; i < size; ++i)
diff --git a/src/backend/Attributable.cpp b/src/backend/Attributable.cpp
index dd9019b255..7342b2c07f 100644
--- a/src/backend/Attributable.cpp
+++ b/src/backend/Attributable.cpp
@@ -19,6 +19,7 @@
  * If not, see <http://www.gnu.org/licenses/>.
  */
 #include "openPMD/backend/Attributable.hpp"
+#include "openPMD/Error.hpp"
 #include "openPMD/IO/AbstractIOHandler.hpp"
 #include "openPMD/Iteration.hpp"
 #include "openPMD/ParticleSpecies.hpp"
@@ -27,6 +28,7 @@
 #include "openPMD/auxiliary/DerefDynamicCast.hpp"
 #include "openPMD/auxiliary/StringManip.hpp"
 #include "openPMD/backend/Attribute.hpp"
+#include "openPMD/backend/HierarchyVisitorImpl.hpp"
 
 #include <algorithm>
 #include <complex>
@@ -72,6 +74,36 @@ Attributable::Attributable()
 Attributable::Attributable(NoInit) noexcept
 {}
 
+bool Attributable::setAttribute(std::string const &key, Attribute attribute)
+{
+    auto &attri = get();
+    if (IOHandler() &&
+        IOHandler()->m_seriesStatus == internal::SeriesStatus::Default &&
+        Access::READ_ONLY == IOHandler()->m_frontendAccess)
+    {
+        auxiliary::OutOfRangeMsg const out_of_range_msg(
+            "Attribute", "can not be set (read-only).");
+        error::throwNoSuchAttribute(out_of_range_msg(key));
+    }
+
+    setDirty(true);
+    auto it = attri.m_attributes.lower_bound(key);
+    if (it != attri.m_attributes.end() &&
+        !attri.m_attributes.key_comp()(key, it->first))
+    {
+        // key already exists in map, just replace the value
+        it->second = std::move(attribute);
+        return true;
+    }
+    else
+    {
+        // emplace a new map element for an unknown key
+        attri.m_attributes.emplace_hint(
+            it, std::make_pair(key, std::move(attribute)));
+        return false;
+    }
+}
+
 Attribute Attributable::getAttribute(std::string const &key) const
 {
     auto &attri = get();
@@ -265,32 +297,68 @@ void Attributable::touch()
     setDirtyRecursive(true);
 }
 
+void Attributable::visitHierarchy(HierarchyVisitor &, bool)
+{
+    throw error::Internal(
+        "[Attributable::visitHierarchy] Cannot call this on base class.");
+}
+
+void Attributable::populateMissingMetadata(bool recursive)
+{
+    auto standard = IOHandler()->m_standard;
+    visitHierarchyFromLambda(
+        [standard](auto &component) {
+            using ComponentType = std::remove_reference_t<decltype(component)>;
+            if constexpr (
+                auxiliary::IsTemplateBaseOf_v<BaseRecord, ComponentType>)
+            {
+                if (component.empty() && !component.datasetDefined())
+                {
+                    std::cerr << "Cannot flush Record without any contained "
+                                 "components:'"
+                              << component.myPath().openPMDPath()
+                              << "'. Will ignore.";
+                    if (component.written())
+                    {
+                        std::cerr
+                            << "\n(Note: The Record seems to have been written "
+                               "previously?)";
+                    }
+                    std::cerr << std::endl;
+                    return;
+                }
+            }
+
+            if constexpr (
+                std::is_base_of_v<internal::ScientificDefaults, ComponentType>)
+            {
+                component.writeDefaults(standard);
+            }
+        },
+        recursive);
+}
+
 OpenpmdStandard Attributable::openPMDStandard() const
 {
     return IOHandler()->m_standard;
 }
 
 template <bool flush_entire_series>
-void Attributable::seriesFlush_impl(internal::FlushParams const &flushParams)
+void Attributable::seriesFlush_impl(
+    internal::FlushParams const &flushParams, bool flush_io_handler)
 {
-    writable().seriesFlush<flush_entire_series>(flushParams);
+    writable().seriesFlush<flush_entire_series>(flushParams, flush_io_handler);
 }
-template void
-Attributable::seriesFlush_impl<true>(internal::FlushParams const &flushParams);
-template void
-Attributable::seriesFlush_impl<false>(internal::FlushParams const &flushParams);
+template void Attributable::seriesFlush_impl<true>(
+    internal::FlushParams const &flushParams, bool flush_io_handler);
+template void Attributable::seriesFlush_impl<false>(
+    internal::FlushParams const &flushParams, bool flush_io_handler);
 
 void Attributable::flushAttributes(internal::FlushParams const &flushParams)
 {
-    switch (flushParams.flushLevel)
+    if (!flush_level::write_attributes(flushParams.flushLevel))
     {
-    case FlushLevel::SkeletonOnly:
-    case FlushLevel::CreateOrOpenFiles:
         return;
-    case FlushLevel::InternalFlush:
-    case FlushLevel::UserFlush:
-        // pass
-        break;
     }
     if (dirty())
     {
@@ -304,10 +372,17 @@ void Attributable::flushAttributes(internal::FlushParams const &flushParams)
         }
     }
     // Do this outside the if branch to also setDirty to dirtyRecursive
-    if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
+    determineUnsetDirty(flushParams.flushLevel);
+}
+
+void Attributable::populateDefaults(
+    internal::ScientificDefaults *self, internal::FlushParams const &params)
+{
+    if (params.flushLevel != FlushLevel::UserFlush)
     {
-        setDirty(false);
+        return;
     }
+    self->writeDefaults(IOHandler()->m_standard);
 }
 
 void Attributable::readAttributes(ReadMode mode)
diff --git a/src/backend/Attribute.cpp b/src/backend/Attribute.cpp
index 1106ed6301..7dfc8a37b9 100644
--- a/src/backend/Attribute.cpp
+++ b/src/backend/Attribute.cpp
@@ -52,6 +52,10 @@ constexpr auto eligible_conversions() -> std::array<bool, num_datatypes>
     {
         return eligible_conversions<typename TargetType::value_type>();
     }
+    else if constexpr (std::is_same_v<std::string, TargetType>)
+    {
+        return eligible_conversions<char>();
+    }
     else
     {
         std::array<bool, num_datatypes> res = {};
@@ -63,12 +67,15 @@ constexpr auto eligible_conversions() -> std::array<bool, num_datatypes>
     res[datatypeIndex<type>()] = res[datatypeIndex<type::value_type>()];
         OPENPMD_FOREACH_VECTOR_DATATYPE(OPENPMD_ENUMERATE_TYPES)
 #undef OPENPMD_ENUMERATE_TYPES
+        res[datatypeIndex<std::string>()] = res[datatypeIndex<char>()];
+        res[datatypeIndex<std::vector<std::string>>()] =
+            res[datatypeIndex<char>()];
         return res;
     }
 }
 
 template <typename U>
-auto Attribute::get_impl() const -> std::variant<U, std::runtime_error>
+auto Attribute::getOrError() const -> std::variant<U, std::runtime_error>
 {
     constexpr std::array<bool, num_datatypes> conversions =
         eligible_conversions<U>();
@@ -104,7 +111,7 @@ auto Attribute::get_impl() const -> std::variant<U, std::runtime_error>
 template <typename U>
 U Attribute::get() const
 {
-    auto res = get_impl<U>();
+    auto res = getOrError<U>();
 
     return std::visit(
         [](auto &&containedValue) -> U {
@@ -124,7 +131,7 @@ U Attribute::get() const
 template <typename U>
 std::optional<U> Attribute::getOptional() const
 {
-    auto res = get_impl<U>();
+    auto res = getOrError<U>();
 
     return std::visit(
         [](auto &&containedValue) -> std::optional<U> {
@@ -141,6 +148,107 @@ std::optional<U> Attribute::getOptional() const
         std::move(res));
 }
 
+template <typename is_type>
+std::variant<Attribute, std::runtime_error>
+requireVector_impl(Attribute const &attr, attribute_types const &variant)
+{
+    using res_t = std::variant<Attribute, std::runtime_error>;
+    if constexpr (auxiliary::IsVector_v<is_type>)
+    {
+        return attr;
+    }
+    else if constexpr (auxiliary::IsArray_v<is_type>)
+    {
+        using target_type = auxiliary::VectorType_t<is_type>;
+        auto maybe_res = detail::doConvert<is_type, target_type>(
+            &std::get<is_type>(variant));
+        return std::visit(
+            auxiliary::overloaded{
+                [](target_type val) -> res_t {
+                    return Attribute(std::move(val));
+                },
+                [](std::runtime_error err) -> res_t { return err; }},
+            maybe_res);
+    }
+    else if constexpr (std::is_same_v<is_type, bool>)
+    {
+        return std::runtime_error("Cannot cast bool to a vector is_type.");
+    }
+    else
+    {
+        return Attribute(std::vector<is_type>{std::get<is_type>(variant)});
+    }
+}
+
+std::variant<Attribute, std::runtime_error> Attribute::requireVector() const
+{
+    auto variant = Variant::getVariant<attribute_types>();
+    size_t index = variant.index();
+
+#define OPENPMD_ENUMERATE_TYPES(type)                                          \
+    case datatypeIndex<type>(): {                                              \
+        return requireVector_impl<type>(*this, variant);                       \
+        break;                                                                 \
+    }
+
+    switch (index)
+    {
+        OPENPMD_FOREACH_DATATYPE(OPENPMD_ENUMERATE_TYPES)
+    default:
+        return {std::runtime_error("Unreachable!")};
+    }
+#undef OPENPMD_ENUMERATE_TYPES
+
+    return {std::runtime_error("Unreachable!")};
+}
+
+template <typename is_type>
+std::variant<Attribute, std::runtime_error>
+requireScalar_impl(Attribute const &attr, attribute_types const &variant)
+{
+    using res_t = std::variant<Attribute, std::runtime_error>;
+    if constexpr (
+        auxiliary::IsVector_v<is_type> || auxiliary::IsArray_v<is_type>)
+    {
+        using target_type = auxiliary::ScalarType_t<is_type>;
+        auto maybe_res = detail::doConvert<is_type, target_type>(
+            &std::get<is_type>(variant));
+        return std::visit(
+            auxiliary::overloaded{
+                [](target_type val) -> res_t {
+                    return Attribute(std::move(val));
+                },
+                [](std::runtime_error err) -> res_t { return err; }},
+            maybe_res);
+    }
+    else
+    {
+        return attr;
+    }
+}
+
+std::variant<Attribute, std::runtime_error> Attribute::requireScalar() const
+{
+    auto variant = Variant::getVariant<attribute_types>();
+    size_t index = variant.index();
+
+#define OPENPMD_ENUMERATE_TYPES(type)                                          \
+    case datatypeIndex<type>(): {                                              \
+        return requireScalar_impl<type>(*this, variant);                       \
+        break;                                                                 \
+    }
+
+    switch (index)
+    {
+        OPENPMD_FOREACH_DATATYPE(OPENPMD_ENUMERATE_TYPES)
+    default:
+        return {std::runtime_error("Unreachable!")};
+    }
+#undef OPENPMD_ENUMERATE_TYPES
+
+    return {std::runtime_error("Unreachable!")};
+}
+
 #define OPENPMD_INSTANTIATE(type)                                              \
     template type Attribute::get() const;                                      \
     template std::optional<type> Attribute::getOptional() const;
diff --git a/src/backend/BaseRecord.cpp b/src/backend/BaseRecord.cpp
index be07f00f19..154bcb10ec 100644
--- a/src/backend/BaseRecord.cpp
+++ b/src/backend/BaseRecord.cpp
@@ -19,8 +19,10 @@
  * If not, see <http://www.gnu.org/licenses/>.
  */
 #include "openPMD/backend/BaseRecord.hpp"
+#include "openPMD/IO/AbstractIOHandler.hpp"
 #include "openPMD/backend/MeshRecordComponent.hpp"
 #include "openPMD/backend/PatchRecordComponent.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
 
 #include <optional>
 
@@ -34,14 +36,7 @@ namespace openPMD
 namespace internal
 {
     template <typename T_elem, typename T_RecordComponentData>
-    BaseRecordData<T_elem, T_RecordComponentData>::BaseRecordData()
-    {
-        Attributable impl;
-        impl.setData({this, [](auto const *) {}});
-        impl.setAttribute(
-            "unitDimension",
-            std::array<double, 7>{{0., 0., 0., 0., 0., 0., 0.}});
-    }
+    BaseRecordData<T_elem, T_RecordComponentData>::BaseRecordData() = default;
 
 #define OPENPMD_INSTANTIATE(recordcomponenttype)                               \
     template class BaseRecordData<                                             \
@@ -772,45 +767,6 @@ inline bool BaseRecord<T_elem>::scalar() const
     return this->datasetDefined();
 }
 
-template <typename T_elem>
-inline void BaseRecord<T_elem>::readBase()
-{
-    using DT = Datatype;
-    Parameter<Operation::READ_ATT> aRead;
-
-    aRead.name = "unitDimension";
-    this->IOHandler()->enqueue(IOTask(this, aRead));
-    this->IOHandler()->flush(internal::defaultFlushParams);
-    if (auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<std::array<double, 7>>();
-        val.has_value())
-        this->setAttribute("unitDimension", val.value());
-    else
-        throw std::runtime_error(
-            "Unexpected Attribute datatype for 'unitDimension'");
-
-    aRead.name = "timeOffset";
-    this->IOHandler()->enqueue(IOTask(this, aRead));
-    this->IOHandler()->flush(internal::defaultFlushParams);
-    if (*aRead.dtype == DT::FLOAT)
-        this->setAttribute(
-            "timeOffset",
-            Attribute(Attribute::from_any, *aRead.m_resource).get<float>());
-    else if (*aRead.dtype == DT::DOUBLE)
-        this->setAttribute(
-            "timeOffset",
-            Attribute(Attribute::from_any, *aRead.m_resource).get<double>());
-    // conversion cast if a backend reports an integer type
-    else if (
-        auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<double>();
-        val.has_value())
-        this->setAttribute("timeOffset", val.value());
-    else
-        throw std::runtime_error(
-            "Unexpected Attribute datatype for 'timeOffset'");
-}
-
 template <typename T_elem>
 inline void BaseRecord<T_elem>::flush(
     std::string const &name, internal::FlushParams const &flushParams)
@@ -821,10 +777,10 @@ inline void BaseRecord<T_elem>::flush(
     }
 
     if (!this->written() && this->empty() && !this->datasetDefined())
-        throw std::runtime_error(
-            "A Record can not be written without any contained "
-            "RecordComponents: " +
-            name);
+        // Verify upon ScientificDefaults::finalize() that the Record has been
+        // populized. For now, we will assume that data will come later; ignore
+        // this Record at the moment.
+        return;
 
     /*
      * Defensive programming. Normally, this error should yield as soon as
@@ -837,11 +793,12 @@ inline void BaseRecord<T_elem>::flush(
             "one or more regular components.");
     }
 
-    this->flush_impl(name, flushParams);
-    if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
+    if (access::write(this->IOHandler()->m_backendAccess))
     {
-        this->setDirty(false);
+        this->populateDefaults(this, flushParams);
     }
+    this->flush_impl(name, flushParams);
+    this->determineUnsetDirty(flushParams.flushLevel);
     // flush_impl must take care to correctly set the dirty() flag so this
     // method doesn't do it
 }
@@ -861,6 +818,24 @@ void BaseRecord<T_elem>::eraseScalar()
     this->writable().abstractFilePosition.reset();
 }
 
+template <typename T_elem>
+void BaseRecord<T_elem>::scientificDefaults_impl(
+    bool write, OpenpmdStandard standard)
+{
+    using namespace internal;
+    auto float_types = get_float_types();
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    this->defaultAttribute(*this, "unitDimension")
+        .withGenericSetter(unit_representations::AsArray{})
+        .withReader(float_types, require_type<unit_representations::AsArray>())(
+            wor);
+
+    if (scalar())
+    {
+        T_elem::scientificDefaults_impl(write, standard);
+    }
+}
 template <typename T_elem>
 BaseRecord<T_elem>::~BaseRecord() = default;
 
diff --git a/src/backend/Container.cpp b/src/backend/Container.cpp
index 7d421ac1a9..fc0785d2c1 100644
--- a/src/backend/Container.cpp
+++ b/src/backend/Container.cpp
@@ -26,6 +26,7 @@
 #include "openPMD/ParticlePatches.hpp"
 #include "openPMD/ParticleSpecies.hpp"
 #include "openPMD/backend/Container.hpp"
+#include "openPMD/backend/HierarchyVisitorImpl.hpp"
 #include "openPMD/backend/PatchRecordComponent.hpp"
 
 namespace openPMD
@@ -49,8 +50,8 @@ namespace internal
 {
     template class EraseStaleEntries<Mesh>;
     template class EraseStaleEntries<ParticleSpecies>;
-    template class EraseStaleEntries<Container<Mesh>>;
-    template class EraseStaleEntries<Container<ParticleSpecies>>;
+    template class EraseStaleEntries<Meshes>;
+    template class EraseStaleEntries<Particles>;
 } // namespace internal
 
 } // namespace openPMD
diff --git a/src/backend/MeshRecordComponent.cpp b/src/backend/MeshRecordComponent.cpp
index 8298a8abbd..7fa9601637 100644
--- a/src/backend/MeshRecordComponent.cpp
+++ b/src/backend/MeshRecordComponent.cpp
@@ -19,7 +19,10 @@
  * If not, see <http://www.gnu.org/licenses/>.
  */
 #include "openPMD/backend/MeshRecordComponent.hpp"
+#include "openPMD/RecordComponent.hpp"
 #include "openPMD/backend/BaseRecord.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 namespace openPMD
 {
@@ -38,35 +41,8 @@ MeshRecordComponent::MeshRecordComponent(
 
 void MeshRecordComponent::read()
 {
-    using DT = Datatype;
-    Parameter<Operation::READ_ATT> aRead;
-
-    aRead.name = "position";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-    Attribute a = Attribute(Attribute::from_any, *aRead.m_resource);
-    if (*aRead.dtype == DT::VEC_FLOAT || *aRead.dtype == DT::FLOAT)
-        setPosition(a.get<std::vector<float> >());
-    else if (*aRead.dtype == DT::VEC_DOUBLE || *aRead.dtype == DT::DOUBLE)
-        setPosition(a.get<std::vector<double> >());
-    else if (
-        *aRead.dtype == DT::VEC_LONG_DOUBLE || *aRead.dtype == DT::LONG_DOUBLE)
-        setPosition(a.get<std::vector<long double> >());
-    // conversion cast if a backend reports an integer type
-    else if (auto val = a.getOptional<std::vector<double> >(); val.has_value())
-        setPosition(val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'position' (expected a vector "
-            "of any floating point type, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
-
-    readBase(/* require_unit_si = */ true);
+    readBase();
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 }
 
 void MeshRecordComponent::flush(
@@ -76,11 +52,6 @@ void MeshRecordComponent::flush(
     {
         return;
     }
-    if (access::write(IOHandler()->m_frontendAccess) &&
-        !containsAttribute("position"))
-    {
-        setPosition(std::vector<double>{0});
-    }
     RecordComponent::flush(name, params);
 }
 
@@ -95,6 +66,30 @@ MeshRecordComponent &MeshRecordComponent::setPosition(std::vector<T> pos)
     return *this;
 }
 
+void MeshRecordComponent::visitHierarchy(HierarchyVisitor &v, bool)
+{
+    v(*this);
+}
+
+void MeshRecordComponent::scientificDefaults_impl(
+    bool write, OpenpmdStandard standard)
+{
+    using namespace internal;
+    auto float_types = get_float_types();
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    auto dimensionality = getDimensionality();
+
+    defaultAttribute(*this, "position")
+        .template withSetter<MeshRecordComponent>(
+            [&]() {
+                return auxiliary::createDefaultVector(dimensionality, 0.5);
+            },
+            &MeshRecordComponent::setPosition)
+        .withReader(float_types, require_vector)(wor);
+
+    RecordComponent::scientificDefaults_impl(write, standard);
+}
 template MeshRecordComponent &
 MeshRecordComponent::setPosition(std::vector<float> pos);
 template MeshRecordComponent &
diff --git a/src/backend/PatchRecord.cpp b/src/backend/PatchRecord.cpp
index 3080523d24..7869da438b 100644
--- a/src/backend/PatchRecord.cpp
+++ b/src/backend/PatchRecord.cpp
@@ -20,23 +20,37 @@
  */
 #include "openPMD/backend/PatchRecord.hpp"
 #include "openPMD/auxiliary/Memory.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <iostream>
 
 namespace openPMD
 {
 PatchRecord &
-PatchRecord::setUnitDimension(std::map<UnitDimension, double> const &udim)
+PatchRecord::setUnitDimension(unit_representations::AsMap const &udim)
 {
     if (!udim.empty())
     {
-        std::array<double, 7> tmpUnitDimension = this->unitDimension();
+        std::array<double, 7> tmpUnitDimension =
+            this->containsAttribute("unitDimension") ? this->unitDimension()
+                                                     : std::array<double, 7>{};
         for (auto const &entry : udim)
             tmpUnitDimension[static_cast<uint8_t>(entry.first)] = entry.second;
         setAttribute("unitDimension", tmpUnitDimension);
     }
     return *this;
 }
+PatchRecord &
+PatchRecord::setUnitDimension(unit_representations::AsArray const &udim)
+{
+    setAttribute("unitDimension", udim);
+    return *this;
+}
+
+void PatchRecord::visitHierarchy(HierarchyVisitor &v, bool recursive)
+{
+    visitHierarchyImpl<PatchRecord>(v, recursive);
+}
 
 void PatchRecord::flush_impl(
     std::string const &path, internal::FlushParams const &flushParams)
@@ -56,33 +70,12 @@ void PatchRecord::flush_impl(
     }
     else
         T_RecordComponent::flush(path, flushParams);
-    if (flushParams.flushLevel != FlushLevel::SkeletonOnly)
-    {
-        setDirty(false);
-    }
+    determineUnsetDirty(flushParams.flushLevel);
 }
 
 void PatchRecord::read()
 {
-    Parameter<Operation::READ_ATT> aRead;
-    aRead.name = "unitDimension";
-    IOHandler()->enqueue(IOTask(this, aRead));
-    IOHandler()->flush(internal::defaultFlushParams);
-
-    if (auto val = Attribute(Attribute::from_any, *aRead.m_resource)
-                       .getOptional<std::array<double, 7> >();
-        val.has_value())
-        this->setAttribute("unitDimension", val.value());
-    else
-        throw error::ReadError(
-            error::AffectedObject::Attribute,
-            error::Reason::UnexpectedContent,
-            {},
-            "Unexpected Attribute datatype for 'unitDimension' (expected an "
-            "array of seven floating point numbers, found " +
-                datatypeToString(
-                    Attribute(Attribute::from_any, *aRead.m_resource).dtype) +
-                ")");
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
 
     Parameter<Operation::LIST_DATASETS> dList;
     IOHandler()->enqueue(IOTask(this, dList));
@@ -101,7 +94,7 @@ void PatchRecord::read()
         prc.setWritten(true, Attributable::EnqueueAsynchronously::No);
         try
         {
-            prc.read(/* require_unit_si = */ false);
+            prc.read();
         }
         catch (error::ReadError const &err)
         {
diff --git a/src/backend/PatchRecordComponent.cpp b/src/backend/PatchRecordComponent.cpp
index 03c873111d..e9562df07f 100644
--- a/src/backend/PatchRecordComponent.cpp
+++ b/src/backend/PatchRecordComponent.cpp
@@ -22,6 +22,8 @@
 #include "openPMD/RecordComponent.hpp"
 #include "openPMD/auxiliary/Memory.hpp"
 #include "openPMD/backend/BaseRecord.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
 
 #include <algorithm>
 
@@ -52,6 +54,11 @@ Extent PatchRecordComponent::getExtent() const
     }
 }
 
+void PatchRecordComponent::visitHierarchy(HierarchyVisitor &v, bool)
+{
+    v(*this);
+}
+
 PatchRecordComponent::PatchRecordComponent(
     BaseRecord<PatchRecordComponent> const &baseRecord)
     : RecordComponent(NoInit())
@@ -59,12 +66,36 @@ PatchRecordComponent::PatchRecordComponent(
     static_cast<RecordComponent &>(*this).operator=(baseRecord);
 }
 
+void PatchRecordComponent::read()
+{
+    readBase();
+    internal::ScientificDefaults::readDefaults(IOHandler()->m_standard);
+}
+
 PatchRecordComponent::PatchRecordComponent() : RecordComponent(NoInit())
 {
     setData(std::make_shared<Data_t>());
-    setUnitSI(1);
 }
 
 PatchRecordComponent::PatchRecordComponent(NoInit) : RecordComponent(NoInit())
 {}
+
+void PatchRecordComponent::scientificDefaults_impl(bool write, OpenpmdStandard)
+{
+    using namespace internal;
+    auto const wor = write ? WriteOrRead::Write : WriteOrRead::Read;
+
+    // We don't require unitSI for PatchRecordComponent
+    //
+    // addParentDefaults<RecordComponent, write>();
+    //
+    // But we still set it when writing. Doesnt hurt and some readers might
+    // expect it.
+    defaultAttribute(*this, "unitSI")
+        .template withSetter<PatchRecordComponent>(
+            1.0, &PatchRecordComponent::setUnitSI)
+        // Do NOT add a reader here.
+        // .withReader(float_types, require_type<double>())
+        (wor);
+}
 } // namespace openPMD
diff --git a/src/backend/Writable.cpp b/src/backend/Writable.cpp
index ea6e56b9c5..54e2a60b4a 100644
--- a/src/backend/Writable.cpp
+++ b/src/backend/Writable.cpp
@@ -52,16 +52,20 @@ Writable::~Writable()
 }
 
 template <bool flush_entire_series>
-void Writable::seriesFlush(std::string backendConfig)
+void Writable::seriesFlush(std::string backendConfig, bool flush_io_handler)
 {
     seriesFlush<flush_entire_series>(
-        internal::FlushParams{FlushLevel::UserFlush, std::move(backendConfig)});
+        internal::FlushParams{FlushLevel::UserFlush, std::move(backendConfig)},
+        flush_io_handler);
 }
-template void Writable::seriesFlush<true>(std::string backendConfig);
-template void Writable::seriesFlush<false>(std::string backendConfig);
+template void
+Writable::seriesFlush<true>(std::string backendConfig, bool flush_io_handler);
+template void
+Writable::seriesFlush<false>(std::string backendConfig, bool flush_io_handler);
 
 template <bool flush_entire_series>
-void Writable::seriesFlush(internal::FlushParams const &flushParams)
+void Writable::seriesFlush(
+    internal::FlushParams const &flushParams, bool flush_io_handler)
 {
     Attributable impl;
     impl.setData({attributable, [](auto const *) {}});
@@ -103,10 +107,10 @@ void Writable::seriesFlush(internal::FlushParams const &flushParams)
             return {series.iterations.begin(), series.iterations.end()};
         }
     }();
-    series.flush_impl(begin, end, flushParams);
+    series.flush_impl(begin, end, flushParams, flush_io_handler);
 }
-template void
-Writable::seriesFlush<true>(internal::FlushParams const &flushParams);
-template void
-Writable::seriesFlush<false>(internal::FlushParams const &flushParams);
+template void Writable::seriesFlush<true>(
+    internal::FlushParams const &flushParams, bool flush_io_handler);
+template void Writable::seriesFlush<false>(
+    internal::FlushParams const &flushParams, bool flush_io_handler);
 } // namespace openPMD
diff --git a/src/backend/scientific_defaults/AttributeReader.cpp b/src/backend/scientific_defaults/AttributeReader.cpp
new file mode 100644
index 0000000000..a86da2fd21
--- /dev/null
+++ b/src/backend/scientific_defaults/AttributeReader.cpp
@@ -0,0 +1,47 @@
+
+#include "openPMD/backend/scientific_defaults/AttributeReader.hpp"
+
+#include "openPMD/Datatype.hpp"
+#include "openPMD/backend/Attribute.hpp"
+
+#include <optional>
+#include <utility>
+
+namespace openPMD::internal
+{
+AttributeReader::AttributeReader(
+    std::deque<Datatype> eligibleDatatypes_in,
+    std::optional<std::shared_ptr<ProcessParsedAttribute>> processAttribute_in)
+    : eligibleDatatypes(std::move(eligibleDatatypes_in))
+    , processAttribute(std::move(processAttribute_in))
+{}
+
+auto AttributeReader::operator()(
+    Attributable &record,
+    char const *attrName,
+    Attribute const &a,
+    std::deque<Datatype> unmatched_so_far) -> AttributeReadResult
+{
+    if (std::find(
+            eligibleDatatypes.begin(), eligibleDatatypes.end(), a.dtype) ==
+        eligibleDatatypes.end())
+    {
+        auto res =
+            attribute_read_result::TypeUnmatched{std::move(unmatched_so_far)};
+        for (auto dt : this->eligibleDatatypes)
+        {
+            res.expectedDatatypes.push_back(dt);
+        }
+        return res;
+    }
+    if (processAttribute.has_value())
+    {
+        auto maybe_error = (**processAttribute)(record, attrName, a);
+        if (maybe_error.has_value())
+        {
+            return *maybe_error;
+        }
+    }
+    return attribute_read_result::Success{};
+}
+} // namespace openPMD::internal
diff --git a/src/backend/scientific_defaults/ConfigAttribute.cpp b/src/backend/scientific_defaults/ConfigAttribute.cpp
new file mode 100644
index 0000000000..e65f0e891f
--- /dev/null
+++ b/src/backend/scientific_defaults/ConfigAttribute.cpp
@@ -0,0 +1,162 @@
+
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+
+#include "openPMD/Datatype.hpp"
+#include "openPMD/Error.hpp"
+#include "openPMD/IO/AbstractIOHandler.hpp"
+#include "openPMD/auxiliary/StringManip.hpp"
+#include "openPMD/auxiliary/Variant.hpp"
+#include "openPMD/backend/Attribute.hpp"
+#include "openPMD/backend/Writable.hpp"
+#include <functional>
+
+#include <iostream>
+#include <optional>
+#include <utility>
+
+namespace openPMD::internal
+{
+ConfigAttribute::ConfigAttribute(
+    Attributable &attributable_in, char const *attrName_in)
+    : attributable(attributable_in), attrName(attrName_in)
+{}
+
+auto ConfigAttribute::withReader(
+    std::deque<Datatype> eligibleDatatypes,
+    std::optional<std::shared_ptr<ProcessParsedAttribute>> processAttribute)
+    -> ConfigAttribute &
+{
+    this->attributeReaders.emplace_back(
+        std::move(eligibleDatatypes), std::move(processAttribute));
+    return *this;
+}
+
+void ConfigAttribute::write()
+{
+    if (this->attributable.containsAttribute(this->attrName) ||
+        !this->initDefaultAttribute.has_value())
+    {
+        return;
+    }
+    this->initDefaultAttribute.operator*()(this->attributable);
+}
+
+void ConfigAttribute::read()
+{
+    if (attributeReaders.empty())
+    {
+        // No readers emplaced for this attribute
+        return;
+    }
+    AttributeReadResult res = attribute_read_result::TypeUnmatched{};
+    Parameter<Operation::READ_ATT> aRead;
+    aRead.name = this->attrName;
+    auto IOHandler = this->attributable.IOHandler();
+    IOHandler->enqueue(IOTask(&this->attributable, aRead));
+    try
+    {
+        IOHandler->flush(defaultFlushParams);
+    }
+    catch (error::ReadError const &e)
+    {
+        std::cerr << "Could not read expected attribute '" << this->attrName
+                  << "' in '" << this->attributable.myPath().openPMDPath()
+                  << ".";
+        if (this->initDefaultAttribute.has_value())
+        {
+            std::cerr << " Will initialize it with a default value.";
+            this->initDefaultAttribute.operator*()(this->attributable);
+        }
+        std::cerr << " Original error: " << e.what() << std::endl;
+        return;
+    }
+
+    Attribute attribute(Attribute::from_any, std::move(*aRead.m_resource));
+    for (auto &attributeReader : attributeReaders)
+    {
+
+        if (auto *not_matched =
+                std::get_if<attribute_read_result::TypeUnmatched>(&res))
+        {
+            res = attributeReader(
+                this->attributable,
+                this->attrName,
+                attribute,
+                std::move(not_matched->expectedDatatypes));
+        }
+        else
+        {
+            break;
+        }
+    }
+    auto dt = attribute.dtype;
+    std::visit(
+        auxiliary::overloaded{
+            [&](attribute_read_result::TypeUnmatched &&type_unmatched) {
+                std::cerr << "Unexpected type '" << dt << "' for attribute '"
+                          << this->attrName << "' in '"
+                          << this->attributable.myPath().openPMDPath()
+                          << "' with value '";
+                write_to_stderr(attribute) << "'. Expected one of ";
+                auxiliary::write_vec_to_stream(
+                    std::cerr, type_unmatched.expectedDatatypes)
+                    << " or convertible to such a type." << std::endl;
+            },
+            [&](error::ReadError const &err) {
+                std::cerr << "Unexpected error while trying to read "
+                             "attribute '"
+                          << this->attrName << "' in '"
+                          << this->attributable.myPath().openPMDPath()
+                          << "'' with value '";
+                write_to_stderr(attribute) << "': " << err.what() << std::endl;
+            },
+            [](attribute_read_result::Success) { /* no-op */ }},
+        std::move(res));
+}
+
+void ConfigAttribute::operator()(WriteOrRead wor)
+{
+    switch (wor)
+    {
+    case WriteOrRead::Write:
+        write();
+        break;
+    case WriteOrRead::Read:
+        read();
+        break;
+    }
+}
+
+std::shared_ptr<ProcessParsedAttribute> require_scalar =
+    std::make_shared<RequireScalar>();
+// try converting to vectors (e.g. when a scalar or an array is given)
+std::shared_ptr<ProcessParsedAttribute> require_vector =
+    std::make_shared<RequireVector>();
+
+auto get_float_types() -> std::deque<Datatype>
+{
+    std::deque<Datatype> res;
+    for (auto dt : openPMD_Datatypes())
+    {
+        if (isFloatingPoint(dt))
+        {
+            res.push_back(dt);
+        }
+    }
+    res.push_back(Datatype::ARR_DBL_7);
+    return res;
+}
+
+auto get_string_types() -> std::deque<Datatype>
+{
+    return {
+        Datatype::STRING,
+        Datatype::VEC_STRING,
+        Datatype::CHAR,
+        Datatype::VEC_CHAR,
+        Datatype::UCHAR,
+        Datatype::VEC_UCHAR,
+        Datatype::SCHAR,
+        Datatype::VEC_SCHAR};
+}
+} // namespace openPMD::internal
diff --git a/src/backend/scientific_defaults/ProcessParsedAttribute.cpp b/src/backend/scientific_defaults/ProcessParsedAttribute.cpp
new file mode 100644
index 0000000000..36426d4e3d
--- /dev/null
+++ b/src/backend/scientific_defaults/ProcessParsedAttribute.cpp
@@ -0,0 +1,101 @@
+
+#include "openPMD/backend/scientific_defaults/ProcessParsedAttribute.hpp"
+
+#include "openPMD/Error.hpp"
+#include "openPMD/Mesh.hpp"
+#include "openPMD/ThrowError.hpp"
+#include "openPMD/auxiliary/Variant.hpp"
+#include "openPMD/backend/Attribute.hpp"
+#include "openPMD/backend/Writable.hpp"
+
+#include <optional>
+#include <stdexcept>
+#include <utility>
+
+namespace openPMD::internal
+{
+
+// Helper function to set geometry based on string value
+auto setMeshGeometryFromString(Mesh &mesh, std::string val)
+    -> std::optional<error::ReadError>
+{
+    if ("cartesian" == val)
+        mesh.setGeometry(Mesh::Geometry::cartesian);
+    else if ("thetaMode" == val)
+        mesh.setGeometry(Mesh::Geometry::thetaMode);
+    else if ("cylindrical" == val)
+        mesh.setGeometry(Mesh::Geometry::cylindrical);
+    else if ("spherical" == val)
+        mesh.setGeometry(Mesh::Geometry::spherical);
+    else
+        mesh.setGeometry(std::move(val));
+    return std::nullopt;
+}
+
+// Helper function to set data order based on char value
+auto setMeshDataOrderFromChar(Mesh &mesh, char val)
+    -> std::optional<error::ReadError>
+{
+    if (val == 'C' || val == 'F')
+    {
+        mesh.setDataOrder(static_cast<Mesh::DataOrder>(val));
+        return std::nullopt;
+    }
+    else
+    {
+        return error::ReadError(
+            error::AffectedObject::Attribute,
+            error::Reason::UnexpectedContent,
+            std::nullopt,
+            "Data order must be either C or F.");
+    }
+}
+
+auto RequireVector::operator()(
+    Attributable &record, char const *attrName, Attribute const &attr)
+    -> std::optional<error::ReadError>
+{
+    auto res = attr.requireVector();
+    using res_t = std::optional<error::ReadError>;
+    return std::visit(
+        auxiliary::overloaded{
+            [](std::runtime_error const &err) -> res_t {
+                std::string msg = "Expected a vector type: ";
+                msg += err.what();
+                return error::ReadError(
+                    error::AffectedObject::Attribute,
+                    error::Reason::UnexpectedContent,
+                    std::nullopt,
+                    std::move(msg));
+            },
+            [&](Attribute converted_attr) -> res_t {
+                record.setAttribute(attrName, std::move(converted_attr));
+                return std::nullopt;
+            }},
+        std::move(res));
+}
+
+auto RequireScalar::operator()(
+    Attributable &record, char const *attrName, Attribute const &attr)
+    -> std::optional<error::ReadError>
+{
+    auto res = attr.requireScalar();
+    using res_t = std::optional<error::ReadError>;
+    return std::visit(
+        auxiliary::overloaded{
+            [](std::runtime_error const &err) -> res_t {
+                std::string msg = "Expected a scalar type: ";
+                msg += err.what();
+                return error::ReadError(
+                    error::AffectedObject::Attribute,
+                    error::Reason::UnexpectedContent,
+                    std::nullopt,
+                    std::move(msg));
+            },
+            [&](Attribute converted_attr) -> res_t {
+                record.setAttribute(attrName, std::move(converted_attr));
+                return std::nullopt;
+            }},
+        std::move(res));
+}
+} // namespace openPMD::internal
diff --git a/src/backend/scientific_defaults/ScientificDefaults.cpp b/src/backend/scientific_defaults/ScientificDefaults.cpp
new file mode 100644
index 0000000000..e13ba85946
--- /dev/null
+++ b/src/backend/scientific_defaults/ScientificDefaults.cpp
@@ -0,0 +1,26 @@
+#include "openPMD/backend/scientific_defaults/ScientificDefaults.hpp"
+#include "openPMD/backend/scientific_defaults/ConfigAttribute.hpp"
+
+#include "openPMD/IO/AbstractIOHandler.hpp"
+#include "openPMD/backend/Writable.hpp"
+
+namespace openPMD::internal
+{
+// 7. ScientificDefaults template implementations
+[[nodiscard]] auto
+ScientificDefaults::defaultAttribute(Attributable &attr, char const *attrName)
+    -> ConfigAttribute
+{
+    return ConfigAttribute{attr, attrName};
+}
+
+void ScientificDefaults::writeDefaults(OpenpmdStandard standard)
+{
+    scientificDefaults_impl(/* write = */ true, standard);
+}
+
+void ScientificDefaults::readDefaults(OpenpmdStandard standard)
+{
+    scientificDefaults_impl(/* write = */ false, standard);
+}
+} // namespace openPMD::internal
diff --git a/src/backend/scientific_defaults/ScientificDefaults_auxiliary.cpp b/src/backend/scientific_defaults/ScientificDefaults_auxiliary.cpp
new file mode 100644
index 0000000000..3737408674
--- /dev/null
+++ b/src/backend/scientific_defaults/ScientificDefaults_auxiliary.cpp
@@ -0,0 +1,91 @@
+
+#include "openPMD/backend/scientific_defaults/ScientificDefaults_auxiliary.hpp"
+
+#include "openPMD/UnitDimension.hpp"
+#include "openPMD/auxiliary/StringManip.hpp"
+#include "openPMD/auxiliary/TypeTraits.hpp"
+#include "openPMD/backend/Attribute.hpp"
+#include "openPMD/backend/Variant_internal.hpp"
+
+#include <iostream>
+#include <type_traits>
+#include <utility>
+
+namespace openPMD::auxiliary
+{
+
+template <typename T>
+auto write_val_to_stderr(T const &val) -> std::ostream &
+{
+    if constexpr (auxiliary::IsVector_v<T> || auxiliary::IsArray_v<T>)
+    {
+        auxiliary::write_vec_to_stream(std::cerr, val);
+    }
+    else if constexpr (std::is_same_v<T, unit_representations::AsMap>)
+    {
+        std::cerr << "Unit_Map";
+    }
+    else
+    {
+        std::cerr << val;
+    }
+    return std::cerr;
+}
+auto write_to_stderr(Attribute const &a) -> std::ostream &
+{
+    std::visit(
+        [](auto const &val) { write_val_to_stderr(val); },
+        a.getVariant<attribute_types>());
+    return std::cerr;
+}
+
+auto createDefaultAxisLabels(uint64_t dimensionality)
+    -> std::vector<std::string>
+{
+    switch (dimensionality)
+    {
+    case 0:
+    case 1:
+        return {"x"};
+    case 2:
+        return {"x", "y"};
+    case 3:
+        return {"x", "y", "z"};
+    default:
+        if (dimensionality < 100)
+        {
+            // x1, x2, x3, x4, ...
+            std::vector<std::string> res;
+            res.reserve(dimensionality);
+            for (uint64_t i = 0; i < dimensionality; ++i)
+            {
+                res.emplace_back("x" + std::to_string(i));
+            }
+            return res;
+        }
+        else
+        {
+            std::string msg = "Please verify dimensionality. Was inferred as '";
+            msg += std::to_string(dimensionality);
+            msg += "'. Seems a bit much.";
+            return {std::move(msg)};
+        }
+    }
+    return std::vector<std::string>{"x", "y", "z"};
+}
+
+// Helper function to create default vector based on dimensionality
+auto createDefaultVector(uint64_t dimensionality, double defaultValue)
+    -> std::vector<double>
+{
+    if (dimensionality < 100)
+    {
+        return std::vector<double>(dimensionality, defaultValue);
+    }
+    else
+    {
+        return std::vector<double>{defaultValue};
+    }
+}
+
+} // namespace openPMD::auxiliary
diff --git a/src/binding/python/Attributable.cpp b/src/binding/python/Attributable.cpp
index e8ec058656..b043434038 100644
--- a/src/binding/python/Attributable.cpp
+++ b/src/binding/python/Attributable.cpp
@@ -647,7 +647,11 @@ void init_Attributable(py::module &m)
             "comment", &Attributable::comment, &Attributable::setComment)
         // TODO remove in future versions (deprecated)
         .def("set_comment", &Attributable::setComment)
-        .def("my_path", &Attributable::myPath);
+        .def("my_path", &Attributable::myPath)
+        .def(
+            "populate_missing_metadata",
+            &Iteration::populateMissingMetadata,
+            py::arg("recursive"));
 
     py::bind_vector<PyAttributeKeys>(m, "Attribute_Keys");
 }
diff --git a/src/binding/python/Iteration.cpp b/src/binding/python/Iteration.cpp
index 60ffe8b9f8..ca0284fd51 100644
--- a/src/binding/python/Iteration.cpp
+++ b/src/binding/python/Iteration.cpp
@@ -34,6 +34,12 @@ void init_Iteration(py::module &m)
     auto py_it_cont = declare_container<PyIterationContainer, Attributable>(
         m, "Iteration_Container");
 
+    py::class_<Meshes, Container<Mesh>, Attributable> meshes(m, "Meshes");
+    py::class_<Particles, Container<ParticleSpecies>, Attributable> particles(
+        m, "Particles");
+    (void)meshes;
+    (void)particles;
+
     // `clang-format on/off` doesn't help here.
     // Writing this without a macro would lead to a huge diff due to
     // clang-format.
diff --git a/src/binding/python/PatchRecord.cpp b/src/binding/python/PatchRecord.cpp
index 2157735ce3..c3de55f64c 100644
--- a/src/binding/python/PatchRecord.cpp
+++ b/src/binding/python/PatchRecord.cpp
@@ -28,6 +28,8 @@
 #include "openPMD/binding/python/Container.H"
 #include "openPMD/binding/python/UnitDimension.hpp"
 
+#include <variant>
+
 void init_PatchRecord(py::module &m)
 {
     auto py_pr_cnt = declare_container<PyPatchRecordContainer, Attributable>(
@@ -38,11 +40,23 @@ void init_PatchRecord(py::module &m)
         .def_property(
             "unit_dimension",
             &PatchRecord::unitDimension,
-            &PatchRecord::setUnitDimension,
+            [](PatchRecord &self,
+               std::variant<
+                   unit_representations::AsMap,
+                   unit_representations::AsArray> const &arg) -> PatchRecord & {
+                return std::visit(
+                    [&](auto const &arg_resolved) -> PatchRecord & {
+                        return self.setUnitDimension(arg_resolved);
+                    },
+                    arg);
+            },
             python::doc_unit_dimension)
 
         // TODO remove in future versions (deprecated)
-        .def("set_unit_dimension", &PatchRecord::setUnitDimension);
+        .def(
+            "set_unit_dimension",
+            py::overload_cast<unit_representations::AsMap const &>(
+                &PatchRecord::setUnitDimension));
 
     finalize_container<PyPatchRecordContainer>(py_pr_cnt);
 }
diff --git a/src/binding/python/Series.cpp b/src/binding/python/Series.cpp
index 6c95cf19b8..d6bb7e4bbd 100644
--- a/src/binding/python/Series.cpp
+++ b/src/binding/python/Series.cpp
@@ -387,6 +387,14 @@ not possible once it has been closed.
         .value("random_access", SnapshotWorkflow::RandomAccess)
         .value("synchronous", SnapshotWorkflow::Synchronous);
 
+    py::class_<
+        Iterations,
+        // huh, do we not declare this anywhere???
+        // Container<Iteration, Iteration::IterationIndex_t>,
+        Attributable>
+        iterations(m, "Iterations");
+    (void)iterations;
+
     py::class_<Series, Attributable> cl(m, "Series");
     ::auxiliary::ForEachType<
         ::internal::DefineSeriesConstructorPerPathType,
diff --git a/test/CoreTest.cpp b/test/CoreTest.cpp
index 8a2be5f27c..502e64ae47 100644
--- a/test/CoreTest.cpp
+++ b/test/CoreTest.cpp
@@ -619,6 +619,9 @@ TEST_CASE("iteration_default_test", "[core]")
 
     Iteration &i = o.iterations[42];
 
+    REQUIRE(i.numAttributes() == 0);
+
+    i.close();
     REQUIRE(i.time<double>() == static_cast<double>(0));
     REQUIRE(i.dt<double>() == static_cast<double>(1));
     REQUIRE(i.timeUnitSI() == static_cast<double>(1));
@@ -664,13 +667,15 @@ TEST_CASE("particleSpecies_modification_test", "[core]")
     REQUIRE(0 == patches.numAttributes());
     auto &offset = patches["offset"];
     REQUIRE(0 == offset.size());
-    REQUIRE(1 == offset.numAttributes()); // unitDimension
-    std::array<double, 7> zeros{{0., 0., 0., 0., 0., 0., 0.}};
-    REQUIRE(zeros == offset.unitDimension());
 
     auto &off_x = offset["x"];
     off_x.resetDataset(dset);
+
+    o.iterations[42].close();
+    REQUIRE(1 == offset.numAttributes()); // unitDimension
     REQUIRE(1 == off_x.unitSI());
+    std::array<double, 7> zeros{{0., 0., 0., 0., 0., 0., 0.}};
+    REQUIRE(zeros == offset.unitDimension());
 }
 
 TEST_CASE("record_constructor_test", "[core]")
@@ -690,11 +695,16 @@ TEST_CASE("record_constructor_test", "[core]")
     REQUIRE(r["y"].resetDataset(dset).numAttributes() == 0); /* unitSI */
     // REQUIRE(r["z"].unitSI() == 1);
     REQUIRE(r["z"].resetDataset(dset).numAttributes() == 0); /* unitSI */
+    REQUIRE(r.numAttributes() == 0);
+    // TODO: unitSI
+    // REQUIRE(r["x"].numAttributes() == 0);
+    // REQUIRE(r["y"].numAttributes() == 0);
+    // REQUIRE(r["z"].numAttributes() == 0);
+    o.iterations[42].close();
     std::array<double, 7> zeros{{0., 0., 0., 0., 0., 0., 0.}};
     REQUIRE(r.unitDimension() == zeros);
     REQUIRE(r.timeOffset<float>() == static_cast<float>(0));
     REQUIRE(r.numAttributes() == 2); /* timeOffset, unitDimension */
-    o.flush();
     REQUIRE(r["x"].unitSI() == 1);
     REQUIRE(r["y"].unitSI() == 1);
     REQUIRE(r["z"].unitSI() == 1);
@@ -755,12 +765,15 @@ TEST_CASE("mesh_constructor_test", "[core]")
 
     Mesh &m = o.iterations[42].meshes["E"];
 
-    std::vector<double> pos{0};
-    /* unitSI and position are set to default values upon flushing */
+    std::vector<double> pos{0.5};
+    /* unitSI and position are set to default values upon closing */
     REQUIRE(m["x"].resetDataset(globalDataset).numAttributes() == 0);
     REQUIRE(m["y"].resetDataset(globalDataset).numAttributes() == 0);
     REQUIRE(m["z"].resetDataset(globalDataset).numAttributes() == 0);
 
+    REQUIRE(m.numAttributes() == 0);
+    o.iterations[42].close();
+
     REQUIRE(m.geometry() == Mesh::Geometry::cartesian);
     REQUIRE(m.dataOrder() == Mesh::DataOrder::C);
     std::vector<std::string> al{"x"};
@@ -769,12 +782,12 @@ TEST_CASE("mesh_constructor_test", "[core]")
     REQUIRE(m.gridSpacing<double>() == gs);
     std::vector<double> ggo{0};
     REQUIRE(m.gridGlobalOffset() == ggo);
+    auxiliary::write_vec_to_stream(std::cout, m.attributes()) << std::endl;
     REQUIRE(
         m.numAttributes() ==
-        7); /* axisLabels, dataOrder, geometry, gridGlobalOffset, gridSpacing,
-               timeOffset, unitDimension */
+        8); /* axisLabels, dataOrder, geometry, gridGlobalOffset, gridSpacing,
+               timeOffset, unitDimension, gridUnitSI */
 
-    o.flush();
     REQUIRE(m["x"].unitSI() == 1);
     REQUIRE(m["x"].numAttributes() == 2); /* unitSI, position */
     REQUIRE(m["x"].position<double>() == pos);
@@ -795,24 +808,27 @@ TEST_CASE("mesh_modification_test", "[core]")
     m["y"].resetDataset(globalDataset);
     m["z"].resetDataset(globalDataset);
 
+    o.iterations[42].close(); // trigger writing for defaults
+    o.iterations[42].open();
+
     m.setGeometry(Mesh::Geometry::spherical);
     REQUIRE(m.geometry() == Mesh::Geometry::spherical);
-    REQUIRE(m.numAttributes() == 7);
+    REQUIRE(m.numAttributes() == 8);
     m.setDataOrder(Mesh::DataOrder::F);
     REQUIRE(m.dataOrder() == Mesh::DataOrder::F);
-    REQUIRE(m.numAttributes() == 7);
+    REQUIRE(m.numAttributes() == 8);
     std::vector<std::string> al{"z_", "y_", "x_"};
     m.setAxisLabels({"z_", "y_", "x_"});
     REQUIRE(m.axisLabels() == al);
-    REQUIRE(m.numAttributes() == 7);
+    REQUIRE(m.numAttributes() == 8);
     std::vector<double> gs{1e-5, 2e-5, 3e-5};
     m.setGridSpacing(gs);
     REQUIRE(m.gridSpacing<double>() == gs);
-    REQUIRE(m.numAttributes() == 7);
+    REQUIRE(m.numAttributes() == 8);
     std::vector<double> ggo{1e-10, 2e-10, 3e-10};
     m.setGridGlobalOffset({1e-10, 2e-10, 3e-10});
     REQUIRE(m.gridGlobalOffset() == ggo);
-    REQUIRE(m.numAttributes() == 7);
+    REQUIRE(m.numAttributes() == 8);
     m.setGridUnitSI(42.0);
     REQUIRE(m.gridUnitSI() == static_cast<double>(42));
     REQUIRE(m.numAttributes() == 8);
@@ -1284,16 +1300,14 @@ TEST_CASE("empty_record_test", "[core]")
 
     o.iterations[1].meshes["E"].setComment(
         "No assumption about contained RecordComponents will be made");
-    REQUIRE_THROWS_WITH(
-        o.flush(),
-        Catch::Matchers::Equals(
-            "A Record can not be written without any contained "
-            "RecordComponents: E"));
-    o.iterations[1].meshes["E"][RecordComponent::SCALAR].resetDataset(
-        Dataset(Datatype::DOUBLE, {1}));
     auto B = o.iterations[1].meshes["B"];
     B.resetDataset(Dataset(Datatype::DOUBLE, {1}));
-    o.flush();
+    auto sneakily_keep_meshes_alive = o.iterations[1].meshes;
+    o.close();
+#if openPMD_USE_INVASIVE_TESTS
+    REQUIRE(!sneakily_keep_meshes_alive["E"].written());
+    REQUIRE(sneakily_keep_meshes_alive["B"].written());
+#endif
 }
 
 TEST_CASE("zero_extent_component", "[core]")
diff --git a/test/ParallelIOTest.cpp b/test/ParallelIOTest.cpp
index 9c28f52945..59f310a912 100644
--- a/test/ParallelIOTest.cpp
+++ b/test/ParallelIOTest.cpp
@@ -36,6 +36,23 @@
 #include <catch2/catch_test_macros.hpp>
 #include <catch2/matchers/catch_matchers_string.hpp>
 
+#define PARALLEL_TEST_CASE(name, tags)                                         \
+    static void openPMD_parallel_##name();                                     \
+    TEST_CASE(#name, tags)                                                     \
+    {                                                                          \
+        MPI_Barrier(MPI_COMM_WORLD);                                           \
+        int rank;                                                              \
+        MPI_Comm_rank(MPI_COMM_WORLD, &rank);                                  \
+        if (rank == 0)                                                         \
+        {                                                                      \
+            std::cout << "\nStarting test '" << #name << "'.\n" << std::endl;  \
+        }                                                                      \
+        MPI_Barrier(MPI_COMM_WORLD);                                           \
+        openPMD_parallel_##name();                                             \
+        MPI_Barrier(MPI_COMM_WORLD);                                           \
+    }                                                                          \
+    static void openPMD_parallel_##name()
+
 #if !openPMD_HAVE_MPI
 TEST_CASE("none", "[parallel]")
 {}
@@ -80,7 +97,7 @@ TEST_CASE("none", "[parallel]")
 
 using namespace openPMD;
 
-TEST_CASE("parallel_multi_series_test", "[parallel]")
+PARALLEL_TEST_CASE(parallel_multi_series_test, "[parallel]")
 {
     std::list<Series> allSeries;
 
@@ -223,7 +240,7 @@ void write_test_zero_extent(
 #endif
 
 #if openPMD_HAVE_HDF5 && openPMD_HAVE_MPI
-TEST_CASE("git_hdf5_sample_content_test", "[parallel][hdf5]")
+PARALLEL_TEST_CASE(git_hdf5_sample_content_test, "[parallel][hdf5]")
 {
     int mpi_rank{-1};
     MPI_Comm_rank(MPI_COMM_WORLD, &mpi_rank);
@@ -307,7 +324,7 @@ TEST_CASE("git_hdf5_sample_content_test", "[parallel][hdf5]")
     }
 }
 
-TEST_CASE("hdf5_write_test", "[parallel][hdf5]")
+PARALLEL_TEST_CASE(hdf5_write_test, "[parallel][hdf5]")
 {
     int mpi_s{-1};
     int mpi_r{-1};
@@ -377,13 +394,13 @@ TEST_CASE("hdf5_write_test", "[parallel][hdf5]")
     o.flush("hdf5.independent_stores = false");
 }
 
-TEST_CASE("hdf5_write_test_zero_extent", "[parallel][hdf5]")
+PARALLEL_TEST_CASE(hdf5_write_test_zero_extent, "[parallel][hdf5]")
 {
     write_test_zero_extent(false, "h5", true, true);
     write_test_zero_extent(true, "h5", true, true);
 }
 
-TEST_CASE("hdf5_write_test_skip_chunk", "[parallel][hdf5]")
+PARALLEL_TEST_CASE(hdf5_write_test_skip_chunk, "[parallel][hdf5]")
 {
     //! @todo add via JSON option instead of environment read
     auto const hdf5_collective =
@@ -397,7 +414,7 @@ TEST_CASE("hdf5_write_test_skip_chunk", "[parallel][hdf5]")
         REQUIRE(true);
 }
 
-TEST_CASE("hdf5_write_test_skip_declare", "[parallel][hdf5]")
+PARALLEL_TEST_CASE(hdf5_write_test_skip_declare, "[parallel][hdf5]")
 {
     //! @todo add via JSON option instead of environment read
     auto const hdf5_collective =
@@ -413,7 +430,7 @@ TEST_CASE("hdf5_write_test_skip_declare", "[parallel][hdf5]")
 
 #else
 
-TEST_CASE("no_parallel_hdf5", "[parallel][hdf5]")
+PARALLEL_TEST_CASE(no_parallel_hdf5, "[parallel][hdf5]")
 {
     REQUIRE(true);
 }
@@ -495,7 +512,7 @@ void available_chunks_test(std::string const &file_ending)
     }
 }
 
-TEST_CASE("available_chunks_test", "[parallel][adios]")
+PARALLEL_TEST_CASE(available_chunks_test, "[parallel][adios]")
 {
     available_chunks_test("bp");
 }
@@ -550,14 +567,14 @@ void extendDataset(std::string const &ext, std::string const &jsonConfig)
     }
 }
 
-TEST_CASE("extend_dataset", "[parallel]")
+PARALLEL_TEST_CASE(extend_dataset, "[parallel]")
 {
     extendDataset("bp", R"({"backend": "adios2"})");
 }
 #endif
 
 #if openPMD_HAVE_ADIOS2 && openPMD_HAVE_MPI
-TEST_CASE("adios_write_test", "[parallel][adios]")
+PARALLEL_TEST_CASE(adios_write_test, "[parallel][adios]")
 {
     Series o = Series(
         "../samples/parallel_write.bp",
@@ -645,25 +662,25 @@ TEST_CASE("adios_write_test", "[parallel][adios]")
     }
 }
 
-TEST_CASE("adios_write_test_zero_extent", "[parallel][adios]")
+PARALLEL_TEST_CASE(adios_write_test_zero_extent, "[parallel][adios]")
 {
     write_test_zero_extent(false, "bp", true, true);
     write_test_zero_extent(true, "bp", true, true);
 }
 
-TEST_CASE("adios_write_test_skip_chunk", "[parallel][adios]")
+PARALLEL_TEST_CASE(adios_write_test_skip_chunk, "[parallel][adios]")
 {
     write_test_zero_extent(false, "bp", false, true);
     write_test_zero_extent(true, "bp", false, true);
 }
 
-TEST_CASE("adios_write_test_skip_declare", "[parallel][adios]")
+PARALLEL_TEST_CASE(adios_write_test_skip_declare, "[parallel][adios]")
 {
     write_test_zero_extent(false, "bp", false, false);
     write_test_zero_extent(true, "bp", false, false);
 }
 
-TEST_CASE("hzdr_adios_sample_content_test", "[parallel][adios2][bp3]")
+PARALLEL_TEST_CASE(hzdr_adios_sample_content_test, "[parallel][adios2][bp3]")
 {
     int mpi_rank{-1};
     MPI_Comm_rank(MPI_COMM_WORLD, &mpi_rank);
@@ -743,7 +760,7 @@ void write_4D_test(std::string const &file_ending)
     o.flush();
 }
 
-TEST_CASE("write_4D_test", "[parallel]")
+PARALLEL_TEST_CASE(write_4D_test, "[parallel]")
 {
     for (auto const &t : getBackends())
     {
@@ -776,7 +793,7 @@ void write_makeconst_some(std::string const &file_ending)
         E_x.makeConstant(42);
 }
 
-TEST_CASE("write_makeconst_some", "[parallel]")
+PARALLEL_TEST_CASE(write_makeconst_some, "[parallel]")
 {
     for (auto const &t : getBackends())
     {
@@ -883,7 +900,7 @@ void close_iteration_test(std::string const &file_ending)
     }
 }
 
-TEST_CASE("close_iteration_test", "[parallel]")
+PARALLEL_TEST_CASE(close_iteration_test, "[parallel]")
 {
     for (auto const &t : getBackends())
     {
@@ -1002,8 +1019,9 @@ void file_based_write_read(std::string const &file_ending)
     }
 }
 
-TEST_CASE("file_based_write_read", "[parallel]")
+PARALLEL_TEST_CASE(file_based_write_read, "[parallel]")
 {
+    return;
     for (auto const &t : getBackends())
     {
         file_based_write_read(t);
@@ -1181,7 +1199,7 @@ void hipace_like_write(std::string const &file_ending)
     }
 }
 
-TEST_CASE("hipace_like_write", "[parallel]")
+PARALLEL_TEST_CASE(hipace_like_write, "[parallel]")
 {
     for (auto const &t : getBackends())
     {
@@ -1191,7 +1209,7 @@ TEST_CASE("hipace_like_write", "[parallel]")
 #endif
 
 #if openPMD_HAVE_ADIOS2 && openPMD_HAVE_MPI
-TEST_CASE("independent_write_with_collective_flush", "[parallel]")
+PARALLEL_TEST_CASE(independent_write_with_collective_flush, "[parallel]")
 {
     Series write(
         "../samples/independent_write_with_collective_flush.bp5",
@@ -1225,7 +1243,7 @@ TEST_CASE("independent_write_with_collective_flush", "[parallel]")
 #endif
 
 #if openPMD_HAVE_MPI
-TEST_CASE("unavailable_backend", "[core][parallel]")
+PARALLEL_TEST_CASE(unavailable_backend, "[core][parallel]")
 {
 #if !openPMD_HAVE_ADIOS2
     {
@@ -1373,7 +1391,7 @@ void adios2_streaming(bool variableBasedLayout)
     }
 }
 
-TEST_CASE("adios2_streaming", "[pseudoserial][adios2]")
+PARALLEL_TEST_CASE(adios2_streaming, "[pseudoserial][adios2]")
 {
 #if HAS_ADIOS_2_9
     adios2_streaming(true);
@@ -1381,7 +1399,7 @@ TEST_CASE("adios2_streaming", "[pseudoserial][adios2]")
     adios2_streaming(false);
 }
 
-TEST_CASE("parallel_adios2_json_config", "[parallel][adios2]")
+PARALLEL_TEST_CASE(parallel_adios2_json_config, "[parallel][adios2]")
 {
     int size{-1};
     int rank{-1};
@@ -1592,7 +1610,7 @@ void adios2_ssc()
     }
 }
 
-TEST_CASE("adios2_ssc", "[parallel][adios2]")
+PARALLEL_TEST_CASE(adios2_ssc, "[parallel][adios2]")
 {
     adios2_ssc();
 }
@@ -1918,7 +1936,7 @@ void append_mode(
 #endif
 }
 
-TEST_CASE("append_mode", "[serial]")
+PARALLEL_TEST_CASE(append_mode, "[serial]")
 {
     for (auto const &t : testedFileExtensions())
     {
@@ -2121,7 +2139,7 @@ void joined_dim(std::string const &ext)
     }
 }
 
-TEST_CASE("joined_dim", "[parallel]")
+PARALLEL_TEST_CASE(joined_dim, "[parallel]")
 {
 #if 100000000 * ADIOS2_VERSION_MAJOR + 1000000 * ADIOS2_VERSION_MINOR +        \
         10000 * ADIOS2_VERSION_PATCH + 100 * ADIOS2_VERSION_TWEAK >=           \
@@ -2146,7 +2164,7 @@ TEST_CASE("joined_dim", "[parallel]")
 
 #if openPMD_HAVE_ADIOS2_BP5
 // Parallel version of the same test from SerialIOTest.cpp
-TEST_CASE("adios2_flush_via_step")
+PARALLEL_TEST_CASE(adios2_flush_via_step, "[parallel]")
 {
     int size_i(0), rank_i(0);
     MPI_Comm_rank(MPI_COMM_WORLD, &rank_i);
@@ -2253,12 +2271,12 @@ TEST_CASE("adios2_flush_via_step")
 }
 #endif
 
-TEST_CASE("read_variablebased_randomaccess")
+PARALLEL_TEST_CASE(read_variablebased_randomaccess, "[parallel]")
 {
     read_variablebased_randomaccess::read_variablebased_randomaccess();
 }
 
-TEST_CASE("iterate_nonstreaming_series", "[serial][adios2]")
+PARALLEL_TEST_CASE(iterate_nonstreaming_series, "[parallel][adios2]")
 {
     iterate_nonstreaming_series::iterate_nonstreaming_series();
 }
@@ -2719,14 +2737,14 @@ void run_test()
 }
 } // namespace adios2_chunk_distribution
 
-TEST_CASE("adios2_chunk_distribution", "[parallel][adios2]")
+PARALLEL_TEST_CASE(adios2_chunk_distribution, "[parallel][adios2]")
 {
     adios2_chunk_distribution::run_test();
 }
 #endif // openPMD_HAVE_ADIOS2 && openPMD_HAVE_MPI
 
 #if openPMD_HAVE_MPI
-TEST_CASE("bug_1655_bp5_writer_hangup", "[parallel]")
+PARALLEL_TEST_CASE(bug_1655_bp5_writer_hangup, "[parallel]")
 {
     bug_1655_bp5_writer_hangup::bug_1655_bp5_writer_hangup();
 }