Implement IceModelVec3::allocate_copy()

pism · Aug 17, 2021 · 031c5a5 · 031c5a5
1 parent 4eb7027
commit 031c5a5
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 0 deletions.
diff --git a/src/util/iceModelVec.hh b/src/util/iceModelVec.hh
@@ -495,6 +495,9 @@ public:
   void sum_columns(double A, double B, IceModelVec2S &output) const;
 
   void copy_from(const IceModelVec3 &input);
+
+  std::shared_ptr<IceModelVec3> allocate_copy() const;
+
 protected:
   bool legal_level(double z) const;
 };

diff --git a/src/util/iceModelVec3.cc b/src/util/iceModelVec3.cc
@@ -305,5 +305,13 @@ void IceModelVec3::copy_from(const IceModelVec3 &input) {
   inc_state_counter();
 }
 
+std::shared_ptr<IceModelVec3> IceModelVec3::allocate_copy() const {
+  auto name = get_name();
+  auto z_name = metadata().get_z().get_name();
+  auto z_attrs = metadata().get_z().get_all_strings();
+
+  return std::make_shared<IceModelVec3>(grid(), name, z_name, levels(), z_attrs);
+}
+
 
 } // end of namespace pism