Added back the firstkick logical to manage restarted kicks

MintonGroup · Aug 19, 2021 · 52fd48e · 52fd48e
1 parent c39f6f6
commit 52fd48e
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 3 deletions.
diff --git a/src/setup/setup.f90 b/src/setup/setup.f90
@@ -143,6 +143,8 @@ module subroutine setup_initialize_system(self, param)
       call self%tp%set_mu(self%cb) 
       call self%pl%eucl_index()
       if (.not.param%lrhill_present) call self%pl%set_rhill(self%cb)
+      self%pl%lfirst = param%lfirstkick
+      self%tp%lfirst = param%lfirstkick
       return
    end subroutine setup_initialize_system
 

diff --git a/src/symba/symba_step.f90 b/src/symba/symba_step.f90
@@ -28,10 +28,8 @@ module subroutine symba_step_system(self, param, t, dt)
                call self%reset(param)
                lencounter = pl%encounter_check(self, dt, 0) .or. tp%encounter_check(self, dt, 0)
                if (lencounter) then
-                  tp%lfirst = pl%lfirst
                   call self%interp(param, t, dt)
-                  pl%lfirst = .true.
-                  tp%lfirst = .true.
+                  param%lfirstkick = .true.
                else
                   self%irec = -1
                   call helio_step_system(self, param, t, dt)
@@ -266,6 +264,9 @@ module subroutine symba_step_reset_system(self, param)
 
                   call system%pl_adds%setup(0, param)
                   call system%pl_discards%setup(0, param)
+
+                  tp%lfirst = param%lfirstkick
+                  pl%lfirst = param%lfirstkick
                end associate
             end select
          end select