changes at Field_t

src/field.h

@@ -94,6 +94,7 @@ public:
 
   T getMissValue() const;
   double getMissValueDouble() const;
+  void setMissValue(double value);
 
 private:
   Field(int grid, size_t gridsize, double missval, size_t size);
@@ -177,6 +178,12 @@ inline double Field<T>::getMissValueDouble() const
   return missval;
 }
 
+template<typename T>
+inline void Field<T>::setMissValue(double value)
+{
+  this->missval = static_cast<T>(value);
+}
+
 /************************************************
  ************************************************
  ******************** Field3D *******************
@@ -228,7 +235,8 @@ using _var_ = std::variant<FieldDouble, FieldFloat>;
 class Field_t : public _var_ {
 public:
   // Overwrites of the std::variant<> constructors
-  constexpr Field_t() noexcept;
+  constexpr Field_t() noexcept
+  : _var_() {  }
 
   Field_t(Field_t& other)
   : _var_(other) {  }
@@ -281,6 +289,47 @@ public:
     _var_::operator=(std::move(t));
     return *this;
   }
+
+  void resize(size_t count) {
+    if(std::holds_alternative<FieldFloat>(*this)) {
+      std::get<FieldFloat>(*this).resize(count);
+    } else if(std::holds_alternative<FieldDouble>(*this)) {
+      std::get<FieldDouble>(*this).resize(count);
+    } else {
+      throw std::runtime_error("Not supported Type");
+    }
+  }
+
+  template<typename U>
+  void resize(size_t count, const U value) {
+    if(std::holds_alternative<FieldFloat>(*this)) {
+      std::get<FieldFloat>(*this).resize(count);
+    } else if(std::holds_alternative<FieldDouble>(*this)) {
+      std::get<FieldDouble>(*this).resize(count);
+    } else {
+      throw std::runtime_error("Not supported Type");
+    }
+  }
+
+  void setMissValue(double value) {
+    if(std::holds_alternative<FieldFloat>(*this)) {
+      std::get<FieldFloat>(*this).setMissValue(value);
+    } else if(std::holds_alternative<FieldDouble>(*this)) {
+      std::get<FieldDouble>(*this).setMissValue(value);
+    } else {
+      throw std::runtime_error("Not supported Type");
+    }
+  }
+
+  double getMissValue() const {
+    if(std::holds_alternative<FieldFloat>(*this)) {
+      return std::get<FieldFloat>(*this).getMissValue();
+    } else if(std::holds_alternative<FieldDouble>(*this)) {
+      return  std::get<FieldDouble>(*this).getMissValue();
+    } else {
+      throw std::runtime_error("Not supported Type");
+    }
+  }
 };