fix

MODFLOW-USGS · Feb 5, 2025 · 7b7139a · 7b7139a
1 parent c844f77
commit 7b7139a
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/Model/ParticleTracking/prt.f90 b/src/Model/ParticleTracking/prt.f90
@@ -405,6 +405,7 @@ end subroutine prt_cq
 
   !> @brief Calculate particle mass storage
   subroutine prt_cq_sto(this)
+    use ParticleModule, only: LVL_CELL
     ! modules
     use TdisModule, only: delt
     use PrtPrpModule, only: PrtPrpType
@@ -436,7 +437,7 @@ subroutine prt_cq_sto(this)
           istatus = packobj%particles%istatus(np)
           ! this may need to change if istatus flags change
           if ((istatus > 0) .and. (istatus /= 8)) then
-            n = packobj%particles%itrdomain(np, 2)
+            n = packobj%particles%itrdomain(np, LVL_CELL)
             ! Each particle currently assigned unit mass
             this%masssto(n) = this%masssto(n) + DONE
           end if