Allowing non zero problo for sponge (erf-model#1578)

Co-authored-by: Mahesh Natarajan <[email protected]>
ewquon · Apr 12, 2024 · 7ef6c86 · 7ef6c86
1 parent 19de6a6
commit 7ef6c86
Showing 1 changed file with 12 additions and 12 deletions.
diff --git a/Source/SourceTerms/ERF_ApplySpongeZoneBCs.cpp b/Source/SourceTerms/ERF_ApplySpongeZoneBCs.cpp
@@ -59,9 +59,9 @@ ApplySpongeZoneBCsForCC (
         int jj = amrex::min(amrex::max(j, domlo_y), domhi_y);
         int kk = amrex::min(amrex::max(k, domlo_z), domhi_z);
 
-        Real x = (ii+0.5) * dx[0];
-        Real y = (jj+0.5) * dx[1];
-        Real z = (kk+0.5) * dx[2];
+        Real x = ProbLoArr[0] + (ii+0.5) * dx[0];
+        Real y = ProbLoArr[1] + (jj+0.5) * dx[1];
+        Real z = ProbLoArr[2] + (kk+0.5) * dx[2];
 
         // x left sponge
         if(use_xlo_sponge_damping){
@@ -171,9 +171,9 @@ ApplySpongeZoneBCsForMom (
         int jj = amrex::min(amrex::max(j, domlo_y), domhi_y);
         int kk = amrex::min(amrex::max(k, domlo_z), domhi_z);
 
-        Real x = ii * dx[0];
-        Real y = (jj+0.5) * dx[1];
-        Real z = (kk+0.5) * dx[2];
+        Real x = ProbLoArr[0] + ii * dx[0];
+        Real y = ProbLoArr[1] + (jj+0.5) * dx[1];
+        Real z = ProbLoArr[2] + (kk+0.5) * dx[2];
 
         // x lo sponge
         if(use_xlo_sponge_damping){
@@ -230,9 +230,9 @@ ApplySpongeZoneBCsForMom (
         int jj = amrex::min(amrex::max(j, domlo_y), domhi_y);
         int kk = amrex::min(amrex::max(k, domlo_z), domhi_z);
 
-        Real x = (ii+0.5) * dx[0];
-        Real y = jj * dx[1];
-        Real z = (kk+0.5) * dx[2];
+        Real x = ProbLoArr[0] + (ii+0.5) * dx[0];
+        Real y = ProbLoArr[1] + jj * dx[1];
+        Real z = ProbLoArr[2] + (kk+0.5) * dx[2];
 
         // x lo sponge
         if(use_xlo_sponge_damping){
@@ -289,9 +289,9 @@ ApplySpongeZoneBCsForMom (
         int jj = amrex::min(amrex::max(j, domlo_y), domhi_y);
         int kk = amrex::min(amrex::max(k, domlo_z), domhi_z);
 
-        Real x = (ii+0.5) * dx[0];
-        Real y = (jj+0.5) * dx[1];
-        Real z = kk * dx[2];
+        Real x = ProbLoArr[0] + (ii+0.5) * dx[0];
+        Real y = ProbLoArr[1] + (jj+0.5) * dx[1];
+        Real z = ProbLoArr[2] + kk * dx[2];
 
         // x left sponge
         if(use_xlo_sponge_damping){