Merge branch 'feature/tsuji1/spheral' of ssh://rzgitlab.llnl.gov:7999…

…/ale/rzspheral into ale3d
LLNL · Aug 12, 2024 · 720d9de · 720d9de
2 parents affa3ce + a4bec1d
commit 720d9de
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Physics/GenericHydro.cc b/src/Physics/GenericHydro.cc
@@ -168,7 +168,7 @@ dt(const DataBase<Dimension>& dataBase,
   if (GeometryRegistrar::coords() == CoordinateType::Spherical) {
     Fdiv = +[](const Tensor& DvDxi, const Vector& posi, const Vector& veli) { return DvDxi[0] + 2.0*veli[0]*safeInv(posi[0]); };
   } else if (GeometryRegistrar::coords() == CoordinateType::RZ) {
-    Fdiv = +[](const Tensor& DvDxi, const Vector& posi, const Vector& veli) { return DvDxi.Trace(); };
+    Fdiv = +[](const Tensor& DvDxi, const Vector& posi, const Vector& veli) { return DvDxi.Trace() + veli[1]*safeInv(posi[1]); };
   }
 
   // Loop over every fluid node.