diff --git a/bfps/cpp/full_code/NSVE.cpp b/bfps/cpp/full_code/NSVE.cpp
index d9332589fdcdda8241ae656386cfb3b64d5b8aa2..292d05298a236fa17ff8f30e1fbfb26eb96f3dc2 100644
--- a/bfps/cpp/full_code/NSVE.cpp
+++ b/bfps/cpp/full_code/NSVE.cpp
@@ -364,6 +364,8 @@ int NSVE<rnumber>::do_stats()
                 this->stat_file,
                 "statistics",
                 H5P_DEFAULT);
+    else
+        stat_group = 0;
     fs->compute_velocity(fs->cvorticity);
     *tmp_vec_field = fs->cvelocity->get_cdata();
     tmp_vec_field->compute_stats(