diff --git a/src/outputs/per_block.cpp b/src/outputs/per_block.cpp
index 6d120608..636764ed 100644
--- 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",