debugging

Author: wpbonelli

src/Model/ParticleTracking/prt-prp.f90

@@ -429,9 +429,6 @@ subroutine release(this, ip, trelease)
 
     ! Accumulate mass for release point
     this%rptm(ip) = this%rptm(ip) + DONE
-
-    print *, "released particle: ", ip
-
   end subroutine release
 
   subroutine initialize_particle(this, particle, ip, trelease)

src/Solution/ParticleTracker/Method.f90

@@ -188,8 +188,6 @@ subroutine save(this, particle, reason)
     end if
 
     call this%trackctl%save(particle, kper=per, kstp=stp, reason=reason)
-    print *, "saving"
-    ! if (kper == 3) stop
 
   end subroutine save
 
@@ -267,7 +265,6 @@ subroutine check(this, particle, cell_defn)
       ! dry particle
       if (particle%idrymeth == 0) then
         ! drop to water table
-        print *, "dropping"
         particle%z = cell_defn%top
         call this%save(particle, reason=1)
       else if (particle%idrymeth == 1) then

src/Solution/ParticleTracker/MethodCellTernary.f90

@@ -9,7 +9,7 @@ module MethodCellTernaryModule
   use SubcellTriModule, only: SubcellTriType, create_subcell_tri
   use ParticleModule
   use GeomUtilModule, only: area
-  use ConstantsModule, only: DZERO, DONE, DTWO
+  use ConstantsModule, only: DZERO, DONE, DTWO, DSIX
   implicit none
 
   private
@@ -139,18 +139,6 @@ subroutine pass_mct(this, particle)
       inface = this%nverts + 3
     end select
 
-    print *, "--------", isc, particle%iscp, exitFace, particle%iscefp
-    ! if (inface /= exitFace .and. isc == particle%iscp .and. exitFace == particle%iscefp) then
-    !   print *, 'subcell cycle, terminating'
-    !   particle%advancing = .false.
-    !   particle%idomain(3) = particle%iscp
-    !   particle%istatus = 2
-    !   call this%save(particle, reason=3)
-    ! else
-    !   particle%iscp = particle%idomain(3)
-    !   particle%iscefp = exitFace
-    ! end if
-
     if (inface .eq. -1) then
       particle%iboundary(2) = 0
     else if (inface .eq. 0) then
@@ -278,6 +266,9 @@ subroutine load_subcell(this, particle, subcell)
           y1 = this%yvertnext(iv0)
           x2 = this%xctr
           y2 = this%yctr
+
+          print *, "subcell in load_subcell: ", ic, isc, x0, y0, x1, y1, x2, y2
+
           x1rel = x1 - x0
           y1rel = y1 - y0
           x2rel = x2 - x0
@@ -410,12 +401,23 @@ subroutine vertvelo(this)
       ! Cell area
       areacell = area(this%xvert, this%yvert)
 
+      print *, "areacell in vertvelo:, ", areacell
+
       ! Cell centroid (in general, this is NOT the average of the vertex coordinates)
-      sixa = areacell * 6.d0
+      sixa = areacell * DSIX
       wk1 = -(this%xvert * this%yvertnext - this%xvertnext * this%yvert)
+
+      print *, "xvert in vertvelo: ", this%xvert
+      print *, "yvert in vertvelo: ", this%yvert
+      print *, "xvertnext in vertvelo: ", this%xvertnext
+      print *, "yvertnext in vertvelo: ", this%yvertnext
+
       this%xctr = sum((this%xvert + this%xvertnext) * wk1) / sixa
       this%yctr = sum((this%yvert + this%yvertnext) * wk1) / sixa
 
+      print *, "xctr in vertvelo: ", this%xctr
+      print *, "yctr in vertvelo: ", this%yctr
+
       ! Subcell areas
       do i = 1, this%nverts
         xvals(1) = this%xvert(i)

src/Solution/ParticleTracker/MethodSubcellTernary.f90

@@ -62,6 +62,8 @@ end subroutine apply_mst
 
   !> @brief Track a particle across a triangular subcell.
   subroutine track_subcell(this, subcell, particle, tmax)
+    ! modules
+    use TdisModule, only: kper, kstp
     ! dummy
     class(MethodSubcellTernaryType), intent(inout) :: this
     class(SubcellTriType), intent(in) :: subcell
@@ -164,6 +166,8 @@ subroutine track_subcell(this, subcell, particle, tmax)
     vzbot = subcell%vzbot
     vztop = subcell%vztop
 
+    print *, "subcell in track_subcell: ", subcell%icell, subcell%isubcell, x0, y0, x1, y1, x2, y2
+
     ! Transform coordinates to the "canonical" configuration:
     ! barycentric in two dimensions with alpha, beta & gamma
     ! such that at f2 alpha = 0, f0 beta = 0, f1 gamma = 0.
@@ -215,13 +219,11 @@ subroutine track_subcell(this, subcell, particle, tmax)
     ! considering only the velocity component in the direction of
     ! the face. Then compute the particle's exit time.
     if (itrifaceexit /= 0) then
-      print *, "exit thru side"
       ! Exit through lateral subcell face
       exitFace = itrifaceexit
       dtexit = dtexitxy
     else if (dtexitz < dtexitxy) then
       ! Exit through top or bottom
-      print *, "exit thru top or bottom"
       if (itopbotexit == -1) then
         exitFace = 4
       else
@@ -232,16 +234,6 @@ subroutine track_subcell(this, subcell, particle, tmax)
     texit = particle%ttrack + dtexit
     t0 = particle%ttrack
 
-    ! ! don't allow a particle to return to the subcell was just in
-    ! if (particle%idomain(2) == particle%icp .and. particle%idomain(3) == particle%iscp) then
-    !   print *, "isc: ", particle%idomain(3)
-    !   particle%idomain(3) = particle%iscp
-    !   particle%istatus = 5
-    !   particle%advancing = .false.
-    !   call this%save(particle, reason=3)
-    !   return
-    ! end if
-
     ! Select user tracking times to solve. If this is the first time step
     ! of the simulation, include all times before it begins; if it is the
     ! last time step include all times after it ends only if the 'extend'
@@ -302,8 +294,11 @@ subroutine track_subcell(this, subcell, particle, tmax)
     ! save previous subcell id and exit face
     particle%iscp = particle%idomain(3)
     particle%iscefp = particle%iboundary(3)
-    print *, "exit face: ", particle%iboundary(3)
-    print *, "subcell: ", particle%idomain(3)
+    ! print *, "exit face: ", particle%iboundary(3)
+    ! print *, "subcell: ", particle%idomain(3)
+    ! print *, "subcell vertices: ", subcell%x0, subcell%x1, subcell%x2, subcell%y0, subcell%y1, subcell%y2
+    ! print *, "cell vertices: ", this%cell%defn%polyvert
+    if (kper == 52 .and. kstp == 1) stop
 
     call this%save(particle, reason=reason)
   end subroutine track_subcell

src/Timing/tdis.f90

@@ -131,6 +131,7 @@ subroutine tdis_set_counters()
       write (line, fmtvspts) kper, kstp
     case (MNORMAL)
       write (line, fmtspts) cpref, kper, kstp, trim(cend)
+      ! write (line, "(a)") ""
     end select
     if (isim_level >= VALL) &
       call write_message(line)