Add vort_mag to stats

parthenon-hpc-lab · Nov 22, 2023 · 07a8ee9 · 07a8ee9
1 parent 71004d5
commit 07a8ee9
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/outputs/per_block.cpp b/src/outputs/per_block.cpp
@@ -411,6 +411,10 @@ void UserOutput::WriteOutputFile(Mesh *pm, ParameterInput *pin, SimTime *tm,
     stats.emplace_back(Stats("vel_mag_mw", "prim", {IV1, IV2, IV3}, Weight::Mass));
     stats.emplace_back(Stats("rho", "prim", IDN));
     stats.emplace_back(Stats("pressure", "prim", IPR));
+    if (pin->GetOrAddBoolean("problem/turbulence", "calc_vorticity_mag", false)) {
+      stats.emplace_back(Stats("vort_mag", "vorticity_mag", 0));
+      stats.emplace_back(Stats("vort_mag_mw", "vorticity_mag", 0), Weight::Mass);
+    }
 
     const std::vector<std::string> stat_types = {
         "min", "max",    "absmin", "absmax", "mean",