CNEWSF fixed for multispecies

MODFLOW-USGS · Jan 10, 2018 · d916dd4 · d916dd4
1 parent 0a5995e
commit d916dd4
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/src/gnt1.f b/src/gnt1.f
@@ -33,9 +33,10 @@ SUBROUTINE GNT1FM(ICOMP)
      &                    numactive,level, ierr)
           ENDIF
           iter = Mxiterxmd
-          CNEWSFTMP=CNEWSF
-          CALL xmdsolv(AMATSF, RHSSF, CNEWSF, hclosexmd, rrctol, IASF, 
-     &                 JASF,NJASF,numactive, north, iter, iacl, ierr)
+          CNEWSFTMP(:,ICOMP)=CNEWSF(:,ICOMP)
+          CALL xmdsolv(AMATSF, RHSSF, CNEWSF(:,ICOMP), hclosexmd, 
+     &                 rrctol, IASF, JASF,NJASF,numactive, north,
+     &                 iter, iacl, ierr)
           n_iter = iter
 C         
           IF(IDSCALE.EQ.1)THEN