update

zhichen3 · Nov 19, 2024 · d2bf1e6 · d2bf1e6
1 parent 7bf9636
commit d2bf1e6
Show file tree

Hide file tree

Showing 3 changed files with 1 addition and 35 deletions.
diff --git a/Source/hydro/Castro_ctu_hydro.cpp b/Source/hydro/Castro_ctu_hydro.cpp
@@ -1268,20 +1268,6 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)  // NOLINT(readability-co
         scale_rad_flux(nbx, rad_flux_arr, area_arr, dt);
 #endif
 
-#if AMREX_SPACEDIM <= 2
-        // get the scaled radial pressure -- we need to treat this specially
-
-        if (idir == 0 && !mom_flux_has_p(0, 0, coord)) {
-            Array4<Real> pradial_fab = pradial.array();
-
-            amrex::ParallelFor(nbx,
-            [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
-            {
-                pradial_fab(i,j,k) = qex_arr(i,j,k,GDPRES) * dt;
-            });
-        }
-#endif
-
         // Store the fluxes from this advance. For simplified SDC integration we
         // only need to do this on the last iteration.
 

diff --git a/Source/hydro/Castro_mol_hydro.cpp b/Source/hydro/Castro_mol_hydro.cpp
@@ -647,13 +647,6 @@ Castro::construct_mol_hydro_source(Real time, Real dt, MultiFab& A_update)
 #endif
 
         // scale the fluxes
-#if AMREX_SPACEDIM <= 2
-        if (!Geom().IsCartesian()) {
-          pradial.resize(xbx, 1);
-        }
-
-        Array4<Real> pradial_fab = pradial.array();
-#endif
 
         for (int idir = 0; idir < AMREX_SPACEDIM; ++idir) {
 
@@ -664,19 +657,6 @@ Castro::construct_mol_hydro_source(Real time, Real dt, MultiFab& A_update)
 
           scale_flux(nbx, flux_arr, area_arr, dt);
 
-#if AMREX_SPACEDIM <= 2
-          // get the scaled radial pressure -- we need to treat this specially
-
-          if (idir == 0 && !mom_flux_has_p(0, 0, coord)) {
-            Array4<Real> const qex_arr = qe[idir].array();
-
-            amrex::ParallelFor(nbx,
-            [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
-            {
-                pradial_fab(i,j,k) = qex_arr(i,j,k,GDPRES) * dt;
-            });
-#endif
-          }
         }
 
 

diff --git a/Source/hydro/trace_ppm.cpp b/Source/hydro/trace_ppm.cpp
@@ -344,7 +344,7 @@ Castro::trace_ppm(const Box& bx,
 #ifndef AMREX_USE_GPU
     do_trace = do_source_trace[QUN];
 #else
-    if (coord == 2 || (idir == 0 && coord = 1)) {
+    if (coord == 2 || (idir == 0 && coord == 1)) {
         do_trace = 1;
     } else {
         do_trace = check_trace_source(srcQ, idir, i, j, k, QUN);