diff --git a/src/boundaries.F90 b/src/boundaries.F90
index a598dd3..3be0c59 100644
--- a/src/boundaries.F90
+++ b/src/boundaries.F90
@@ -5246,6 +5246,9 @@ module boundaries
 
             qn(1:nv,it,jl:ju,kl:ku) =   qn(1:nv,is,jl:ju,kl:ku)
             qn(ivx ,it,jl:ju,kl:ku) = - qn(ivx ,is,jl:ju,kl:ku)
+            if (ibx > 0) then
+              qn(ibx ,it,jl:ju,kl:ku) = - qn(ibx ,is,jl:ju,kl:ku)
+            end if
           end do
         else
           do i = 1, ng
@@ -5254,6 +5257,9 @@ module boundaries
 
             qn(1:nv,it,jl:ju,kl:ku) =   qn(1:nv,is,jl:ju,kl:ku)
             qn(ivx ,it,jl:ju,kl:ku) = - qn(ivx ,is,jl:ju,kl:ku)
+            if (ibx > 0) then
+              qn(ibx ,it,jl:ju,kl:ku) = - qn(ibx ,is,jl:ju,kl:ku)
+            end if
           end do
         end if
 
@@ -5476,6 +5482,9 @@ module boundaries
 
             qn(1:nv,il:iu,jt,kl:ku) =   qn(1:nv,il:iu,js,kl:ku)
             qn(ivy ,il:iu,jt,kl:ku) = - qn(ivy ,il:iu,js,kl:ku)
+            if (iby > 0) then
+              qn(iby ,il:iu,jt,kl:ku) = - qn(iby ,il:iu,js,kl:ku)
+            end if
           end do
         else
           do j = 1, ng
@@ -5484,6 +5493,9 @@ module boundaries
 
             qn(1:nv,il:iu,jt,kl:ku) =   qn(1:nv,il:iu,js,kl:ku)
             qn(ivy ,il:iu,jt,kl:ku) = - qn(ivy ,il:iu,js,kl:ku)
+            if (iby > 0) then
+              qn(iby ,il:iu,jt,kl:ku) = - qn(iby ,il:iu,js,kl:ku)
+            end if
           end do
         end if
 
@@ -5566,6 +5578,9 @@ module boundaries
 
             qn(1:nv,il:iu,jl:ju,kt) =   qn(1:nv,il:iu,jl:ju,ks)
             qn(ivz ,il:iu,jl:ju,kt) = - qn(ivz ,il:iu,jl:ju,ks)
+            if (ibz > 0) then
+              qn(ibz ,il:iu,jl:ju,kt) = - qn(ibz ,il:iu,jl:ju,ks)
+            end if
           end do
         else
           do k = 1, ng
@@ -5574,6 +5589,9 @@ module boundaries
 
             qn(1:nv,il:iu,jl:ju,kt) =   qn(1:nv,il:iu,jl:ju,ks)
             qn(ivz ,il:iu,jl:ju,kt) = - qn(ivz ,il:iu,jl:ju,ks)
+            if (ibz > 0) then
+              qn(ibz ,il:iu,jl:ju,kt) = - qn(ibz ,il:iu,jl:ju,ks)
+            end if
           end do
         end if