From d2bf1e62e29ebf40462585121cc3f2939bf8894a Mon Sep 17 00:00:00 2001
From: zhi <zhichen0428@gmail.com>
Date: Mon, 18 Nov 2024 20:36:55 -0500
Subject: [PATCH] update

---
 Source/hydro/Castro_ctu_hydro.cpp | 14 --------------
 Source/hydro/Castro_mol_hydro.cpp | 20 --------------------
 Source/hydro/trace_ppm.cpp        |  2 +-
 3 files changed, 1 insertion(+), 35 deletions(-)

diff --git a/Source/hydro/Castro_ctu_hydro.cpp b/Source/hydro/Castro_ctu_hydro.cpp
index 03ea43e428..8e79739f0b 100644
--- 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
index 1ae6891d7e..69edfb10fe 100644
--- 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
index 6d2150c43f..cd20d147c9 100644
--- 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);