pull get_iatop routine into base Method

Author: wpbonelli

src/Solution/ParticleTracker/Method.f90

@@ -180,21 +180,12 @@ subroutine update(this, particle, celldefn)
 
   end subroutine update
 
-  !> @brief Get the index used to locate top elevation of a cell in the grid.
-  !! This is -1 if the cell is in the top layer and 1 otherwise. (kluge?)
-  function get_iatop(dis, ic) result(iatop)
+  !> @brief Get the index corresponding to top elevation of a cell in the grid.
+  !! This is -1 if the cell is in the top layer and 1 otherwise.
+  function get_iatop(ncpl, icu) result(iatop)
     use BaseDisModule, only: DisBaseType
-    implicit none
-    ! -- dummy
-    class(DisBaseType), intent(in) :: dis
-    integer, intent(in) :: ic
-    ! -- result
+    integer, intent(in) :: ncpl, icu
     integer :: iatop
-    ! -- local
-    integer :: ncpl, icu
-    !
-    ncpl = dis%get_ncpl()
-    icu = dis%get_nodeuser(ic)
     !
     if (icu .le. ncpl) then
       iatop = -1

src/Solution/ParticleTracker/MethodDis.f90

@@ -327,7 +327,8 @@ subroutine load_cellDefn_basic(this, ic, cellDefn)
     cellDefn%npolyverts = 4 ! rectangular cell always has 4 vertices
     !
     ! -- iatop, top, and bot
-    iatop = get_iatop(this%fmi%dis, ic)
+    iatop = get_iatop(this%fmi%dis%get_ncpl(), &
+                      this%fmi%dis%get_nodeuser(ic))
     cellDefn%iatop = iatop
     top = this%fmi%dis%top(ic)
     bot = this%fmi%dis%bot(ic)

src/Solution/ParticleTracker/MethodDisv.f90

@@ -377,7 +377,8 @@ subroutine load_cellDefn_basic(this, ic, cellDefn)
     cellDefn%npolyverts = this%get_npolyverts(ic)
     !
     ! -- iatop, top, and bot
-    iatop = get_iatop(this%fmi%dis, ic)
+    iatop = get_iatop(this%fmi%dis%get_ncpl(), &
+                      this%fmi%dis%get_nodeuser(ic))
     cellDefn%iatop = iatop
     top = this%fmi%dis%top(ic)
     bot = this%fmi%dis%bot(ic)