const T: include last timestep; do record_all at t=outstart

igfuw · Apr 3, 2024 · 55b400c · 55b400c
1 parent 331ccb7
commit 55b400c
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/libmpdata++/output/detail/output_common.hpp b/libmpdata++/output/detail/output_common.hpp
@@ -158,7 +158,7 @@ namespace libmpdataxx
               record_time = this->time;
               for (int t = 0; t < outwindow; ++t)
               {
-                if ((this->timestep - t) % static_cast<int>(outfreq) == 0 && this->timestep > static_cast<int>(outstart)) record_all();
+                if ((this->timestep - t) % static_cast<int>(outfreq) == 0 && this->timestep >= static_cast<int>(outstart)) record_all();
               }
             }
           }

diff --git a/libmpdata++/output/hdf5.hpp b/libmpdata++/output/hdf5.hpp
@@ -199,7 +199,7 @@ namespace libmpdataxx
             // T
             {
               const hsize_t
-                nt_out = (nt - this->outstart) / this->outfreq + 1; // incl. t=0
+                nt_out = (nt - this->outstart) / this->outfreq + 2; // incl. t=0 and t=outstart
               float dt = this->dt;
 
               blitz::Array<typename solver_t::real_t, 1> coord(nt_out);