refactor(xt3d): accumulate flowja instead of set (#306)

This change is required for transport, which accumulates terms in flowja
MODFLOW-USGS · Feb 12, 2020 · 0f2a497 · 0f2a497
1 parent ef46069
commit 0f2a497
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/Model/ModelUtilities/Xt3dInterface.f90 b/src/Model/ModelUtilities/Xt3dInterface.f90
@@ -870,8 +870,8 @@ subroutine xt3d_flowja(this, hnew, flowja)
         call this%xt3d_qnbrs(nodes, m, n, nnbr1, inbr1, chat1j, hnew, qnbrs)
         qnm = qnm - qnbrs
         ipos = ii01
-        flowja(ipos) = qnm
-        flowja(this%dis%con%isym(ipos)) = -qnm        
+        flowja(ipos) = flowja(ipos) + qnm
+        flowja(this%dis%con%isym(ipos)) = flowja(this%dis%con%isym(ipos)) - qnm        
       enddo
     enddo
     !