sqrt missing in normalization

lanl · Mar 7, 2024 · 6474920 · 6474920
1 parent cae28b4
commit 6474920
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/phoebus_driver.cpp b/src/phoebus_driver.cpp
@@ -1075,7 +1075,8 @@ void UserWorkBeforeOutput(MeshBlock *pmb, ParameterInput *pin) {
         const Real z = coords.Xc<3>(k, j, i);
         const Real sigma = analysis->Param<Real>("sigma");
         const Real pi = 3.14;
-        const Real s0 = s * std::exp(-z * z / sigma / sigma) / pi / sigma; // sigma > 0
+        const Real s0 =
+            s * std::exp(-z * z / sigma / sigma) / std::sqrt(pi) / sigma; // sigma > 0
 
         const Real vp[3] = {v(0, p::velocity(0), k, j, i), v(0, p::velocity(1), k, j, i),
                             v(0, p::velocity(2), k, j, i)};