fix white space

erf-model · Nov 23, 2024 · 2ba3848 · 2ba3848
1 parent b770b0c
commit 2ba3848
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 7 deletions.
diff --git a/Source/LinearSolvers/ERF_TerrainPoisson.cpp b/Source/LinearSolvers/ERF_TerrainPoisson.cpp
@@ -30,22 +30,22 @@ void TerrainPoisson::apply(MultiFab& lhs, MultiFab const& rhs)
     xx.FillBoundary(m_geom.periodicity());
 
     if (!m_geom.isPeriodic(0)) {
-        for (MFIter mfi(xx,true); mfi.isValid(); ++mfi) 
+        for (MFIter mfi(xx,true); mfi.isValid(); ++mfi)
         {
             Box bx = mfi.tilebox();
             const Array4<Real>& x_arr = xx.array(mfi);
             ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k)
             {
                 if (i == domlo.x) {
                     x_arr(i-1,j,k) = x_arr(i,j,k);
-                } else if (i == domhi.x) {
-                    x_arr(i+1,j,k) = x_arr(i,j,k);
+                } else if (i == domhi.x) { // OUTFLOW
+                    x_arr(i+1,j,k) = -x_arr(i,j,k);
                 }
             });
         }
     }
     if (!m_geom.isPeriodic(1)) {
-        for (MFIter mfi(xx,true); mfi.isValid(); ++mfi) 
+        for (MFIter mfi(xx,true); mfi.isValid(); ++mfi)
         {
             Box bx = mfi.tilebox();
             Box bx2(bx); bx2.grow(0,1);
@@ -61,7 +61,7 @@ void TerrainPoisson::apply(MultiFab& lhs, MultiFab const& rhs)
         } // mfi
     }
 
-    for (MFIter mfi(xx,true); mfi.isValid(); ++mfi) 
+    for (MFIter mfi(xx,true); mfi.isValid(); ++mfi)
     {
         Box bx = mfi.tilebox();
         Box bx2(bx); bx2.grow(0,1);

diff --git a/Source/LinearSolvers/ERF_TerrainPoisson_3D_K.H b/Source/LinearSolvers/ERF_TerrainPoisson_3D_K.H
@@ -40,7 +40,7 @@ T terrpoisson_flux_x (int i, int j, int k,
     px_lo -= Real(0.5) * ( pz_lo_md_hi + pz_lo_md_lo );
 
     return -px_lo;
-} 
+}
 template <typename T>
 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
 T terrpoisson_flux_y (int i, int j, int k,

diff --git a/Source/LinearSolvers/ERF_compute_divergence.cpp b/Source/LinearSolvers/ERF_compute_divergence.cpp
@@ -74,7 +74,7 @@ void ERF::compute_divergence (int lev, MultiFab& rhs, Vector<MultiFab>& mom_mf,
             });
         } // mfi
 
-    } else { // no terrain 
+    } else { // no terrain
         computeDivergence(rhs, rho0_u_const, geom_at_lev);
 
     }