fix incorrect calculation

zhichen3 · Sep 11, 2024 · ca04f9e · ca04f9e
1 parent 79562a8
commit ca04f9e
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Source/hydro/Castro_ctu.cpp b/Source/hydro/Castro_ctu.cpp
@@ -74,8 +74,8 @@ Castro::consup_hydro(const Box& bx,
     } else if (n == UMY && !mom_flux_has_p(1, 1, geomdata.Coord())) {
         // Add gradp term to polar(theta) momentum equation for Spherical 2D geometry
 
-        Real r = geomdata.ProbLoArray()[0] + (static_cast<Real>(i) + 0.5_rt) * geomdata.Coord()[0];
-        U_new(i,j,k,UMY) += - dt * (qy(i,j+1,k,GDPRES) - qy(i,j,k,GDPRES)) / (r * geomdata.CellSize()[0]);
+        Real r = geomdata.ProbLoArray()[0] + (static_cast<Real>(i) + 0.5_rt) * geomdata.Cellsize()[0];
+        U_new(i,j,k,UMY) += - dt * (qy(i,j+1,k,GDPRES) - qy(i,j,k,GDPRES)) / (r * geomdata.CellSize()[1]);
 
     }
   });

diff --git a/Source/hydro/Castro_mol.cpp b/Source/hydro/Castro_mol.cpp
@@ -270,7 +270,7 @@ Castro::mol_consup(const Box& bx,  // NOLINT(readability-convert-member-function
         } else if (n == UMY && !mom_flux_has_p(1, 1, coord)) {
             // Add gradp term to polar(theta) momentum equation for Spherical 2D geometry
 
-            Real r = prob_lo[0] + (static_cast<Real>(i) + 0.5_rt)*dx[0];
+            Real r = prob_lo[0] + (static_cast<Real>(i) + 0.5_rt) * dx[0];
             update(i,j,k,UMY) -= (q1(i,j+1,k,GDPRES) - q1(i,j,k,GDPRES)) / (r * dx[1]);
         }
     }