diff --git a/sn3d.cc b/sn3d.cc
index 1e2552db4..c2054a78c 100644
--- a/sn3d.cc
+++ b/sn3d.cc
@@ -872,7 +872,7 @@ auto main(int argc, char *argv[]) -> int {
   // initialise or read in virtual packet spectra
   vpkt_init(globals::timestep, my_rank, globals::simulation_continued_from_saved);
 
-  while (globals::timestep < std::min(globals::timestep_finish, globals::ntimesteps - 1) && !terminate_early) {
+  while (globals::timestep < std::min(globals::timestep_finish, globals::ntimesteps) && !terminate_early) {
     MPI_Barrier(MPI_COMM_WORLD);
 
     // titer example: Do 3 iterations on timestep 0-6