Fixed problem where bad angular momentum conservation was not flagged as a failed try

MintonGroup · May 23, 2021 · 7abe593 · 7abe593
1 parent 243ffcd
commit 7abe593
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
@@ -178,6 +178,7 @@ subroutine symba_frag_pos(param, symba_plA, family, x, v, L_spin, Ip, mass, radi
    call define_pre_collisional_family()
 
    try = 1
+   lmerge = .false.
    do while (nfrag >= NFRAG_MIN)
       !write(*,        "(' -------------------------------------------------------------------------------------')")
       !write(*,*) "Try number: ",try, ' of ',ntry
@@ -214,6 +215,7 @@ subroutine symba_frag_pos(param, symba_plA, family, x, v, L_spin, Ip, mass, radi
             lmerge = .true.
          else if (abs(dLmag) > Ltol) then
             write(*,*) 'Failed try ',try,': Angular momentum error too big: ',dLmag
+            lmerge = .true.
          end if
 
          !write(*,        "(' -------------------------------------------------------------------------------------')")
@@ -238,6 +240,7 @@ subroutine symba_frag_pos(param, symba_plA, family, x, v, L_spin, Ip, mass, radi
       if (.not.lmerge) exit
       call restructure_failed_fragments()
       try = try + 1
+      lmerge = .false.
    end do
    !write(*,        "(' -------------------------------------------------------------------------------------')")
    !write(*,        "('  Final diagnostic')")