print

wpbonelli · Dec 23, 2024 · ca581c6 · ca581c6
1 parent 845f47f
commit ca581c6
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src/Solution/ParticleTracker/TernarySolveTrack.f90 b/src/Solution/ParticleTracker/TernarySolveTrack.f90
@@ -297,12 +297,21 @@ subroutine solve_coefs(alpi, beti)
 
   !> @brief Step (evaluate) analytically depending on case
   subroutine step_analytical(t, alp, bet)
+    use TdisModule, only: kper
     ! -- dummy
     real(DP), intent(in) :: t
     real(DP) :: alp
     real(DP) :: bet
 
     if (icase .eq. 1) then
+      if (kper > 86) then
+        print *, "----"
+        print *, ca1
+        print *, ca2
+        print *, ca3
+        print *, waa
+        print *, t
+      end if
       alp = ca1 + ca2 * dexp(waa * t) + ca3 * dexp(wbb * t)
       bet = cb1 + cb2 * dexp(wbb * t)
     else if (icase .eq. -1) then