fix

src/Solution/ParticleTracker/MethodDisv.f90

@@ -27,7 +27,6 @@ module MethodDisvModule
     procedure, public :: deallocate !< deallocate arrays and scalars
     procedure, public :: load => load_disv !< load the cell method
     procedure, public :: load_cell_defn !< load cell definition from the grid
-    procedure, public :: load_cell !< load the cell from the grid
     procedure, public :: pass => pass_disv !< pass the particle to the next cell
     procedure :: map_neighbor !< map a location on the cell face to the shared face of a neighbor
     procedure :: update_flowja !< update intercell mass flows
@@ -88,7 +87,9 @@ subroutine load_disv(this, particle, next_level, submethod)
     type is (CellPolyType)
       ic = particle%idomain(next_level)
       call this%load_cell_defn(ic, cell%defn)
-      call this%load_cell(ic, cell)
+      cell%xO = cell%defn%polyvert(1, 1)
+      cell%yO = cell%defn%polyvert(2, 1)
+      cell%zO = cell%defn%bot
       if (this%fmi%ibdgwfsat0(ic) == 0) then
         call method_cell_ptb%init( &
           fmi=this%fmi, &
@@ -333,18 +334,6 @@ subroutine load_cell_defn(this, ic, defn)
     call this%load_flows(defn)
   end subroutine load_cell_defn
 
-  subroutine load_cell(this, ic, cell)
-    ! dummy
-    class(MethodDisvType), intent(inout) :: this
-    integer(I4B), intent(in) :: ic
-    type(CellPolyType), pointer, intent(inout) :: cell
-
-    cell%xO = cell%defn%polyvert(1, 1)
-    cell%yO = cell%defn%polyvert(2, 1)
-    cell%zO = cell%defn%bot
-
-  end subroutine load_cell
-
   !> @brief Loads cell properties to cell definition from the grid.
   subroutine load_properties(this, ic, defn)
     ! dummy