Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
simplified encounter saves
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed Sep 17, 2021
1 parent 2268173 commit 85b296b
Showing 1 changed file with 4 additions and 6 deletions.
10 changes: 4 additions & 6 deletions src/symba/symba_encounter_check.f90
Original file line number Diff line number Diff line change
Expand Up @@ -175,8 +175,8 @@ module function symba_encounter_check_pl(self, param, system, dt, irec) result(l
call symba_encounter_check_all_triangular(npl, nplm, pl%xh, pl%vh, pl%rhill, dt, irec, lvdotr, index1, index2, nenc)
lany_encounter = nenc > 0
if (lany_encounter) then
call move_alloc(lvdotr, plplenc_list%lvdotr)
call move_alloc(index1, plplenc_list%index1)
call move_alloc(lvdotr, plplenc_list%lvdotr)
call move_alloc(index1, plplenc_list%index1)
call move_alloc(index2, plplenc_list%index2)
end if
end if
Expand All @@ -185,10 +185,8 @@ module function symba_encounter_check_pl(self, param, system, dt, irec) result(l
do k = 1, nenc
i = plplenc_list%index1(k)
j = plplenc_list%index2(k)
call util_index_eucl_ij_to_k(npl, i, j, kenc)
plplenc_list%id1(k) = pl%id(plplenc_list%index1(k))
plplenc_list%id2(k) = pl%id(plplenc_list%index2(k))
plplenc_list%status(k) = ACTIVE
plplenc_list%id1(k) = pl%id(i)
plplenc_list%id2(k) = pl%id(j)
plplenc_list%level(k) = irec
pl%lencounter(i) = .true.
pl%lencounter(j) = .true.
Expand Down

0 comments on commit 85b296b

Please sign in to comment.