fix(PrtFmiType): accommodate 0s in nodelist

MODFLOW-USGS · Nov 2, 2024 · e31e643 · e31e643
1 parent 3be3c39
commit e31e643
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/Model/ParticleTracking/prt-fmi.f90 b/src/Model/ParticleTracking/prt-fmi.f90
@@ -191,6 +191,7 @@ subroutine accumulate_flows(this)
       end if
       do ib = 1, this%gwfpackages(ip)%nbound
         i = this%gwfpackages(ip)%nodelist(ib)
+        if (i <= 0) cycle
         if (this%ibound(i) <= 0) cycle
         qbnd = this%gwfpackages(ip)%get_flow(ib)
         ! todo, after initial release: default iflowface values for different packages