change dtheta to rdtheta

AMReX-Astro · Sep 17, 2024 · 0ff40ae · 0ff40ae
1 parent b050905
commit 0ff40ae
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/Source/driver/Castro.cpp b/Source/driver/Castro.cpp
@@ -2884,8 +2884,8 @@ Castro::reflux (int crse_level, int fine_level, bool in_post_timestep)
 
             reg = &getLevel(lev).pres_reg;
 
-            MultiFab dtheta(crse_lev.grids, crse_lev.dmap, 1, 0);
-            dtheta.setVal(crse_lev.geom.CellSize(1));
+            MultiFab rdtheta(crse_lev.grids, crse_lev.dmap, 1, 0);
+            rdtheta.setVal(crse_lev.geom.CellSize(1));
 
             // Need to multiply by r for rdtheta
 
@@ -2904,12 +2904,12 @@ Castro::reflux (int crse_level, int fine_level, bool in_post_timestep)
                 for (OrientationIter fi; fi.isValid(); ++fi)
                 {
                     const FabSet& fs = (*reg)[fi()];
-                    if (fi().coordDir() == 0) {
+                    if (fi().coordDir() == 1) {
                         fs.copyTo(tmp_fluxes, 0, 0, 0, tmp_fluxes.nComp());
                     }
                 }
 
-                MultiFab::Add(crse_lev.P_theta, tmp_fluxes, 0, 0, crse_lev.P_radial.nComp(), 0);
+                MultiFab::Add(crse_lev.P_theta, tmp_fluxes, 0, 0, crse_lev.P_theta.nComp(), 0);
 
             }