Fix sign error

ewquon · Jan 27, 2024 · b2d5747 · b2d5747
1 parent 7e294c1
commit b2d5747
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 4 deletions.
diff --git a/Exec/DevTests/TemperatureSource/prob.cpp b/Exec/DevTests/TemperatureSource/prob.cpp
@@ -187,8 +187,7 @@ Problem::update_rhotheta_sources (
         {
             const amrex::Real z_cc = prob_lo[2] + (k+0.5)* dx[2];
             if ((time > 1.0) && (z_cc < 100.)) {
-                // +1 K/hr
-                src[k] = 0.0002777777777777778;
+                src[k] = 1.0; // 1 K/s
             } else {
                 src[k] = 0.0;
             }

diff --git a/Source/TimeIntegration/ERF_slow_rhs_inc.cpp b/Source/TimeIntegration/ERF_slow_rhs_inc.cpp
@@ -648,7 +648,7 @@ void erf_slow_rhs_inc (int /*level*/, int nrk,
             const int n = RhoTheta_comp;
             ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
             {
-                cell_rhs(i, j, k, n) -= dptr_rhotheta_src[k];
+                cell_rhs(i, j, k, n) += dptr_rhotheta_src[k];
             });
         }
 

diff --git a/Source/TimeIntegration/ERF_slow_rhs_pre.cpp b/Source/TimeIntegration/ERF_slow_rhs_pre.cpp
@@ -728,7 +728,7 @@ void erf_slow_rhs_pre (int level, int finest_level,
             const int n = RhoTheta_comp;
             ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
             {
-                cell_rhs(i, j, k, n) -= dptr_rhotheta_src[k];
+                cell_rhs(i, j, k, n) += dptr_rhotheta_src[k];
             });
         }