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

Commit

Permalink
Commented out failure messages in fragment solver
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed Aug 11, 2021
1 parent 46f333e commit 3774b40
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions src/fragmentation/fragmentation.f90
Original file line number Diff line number Diff line change
Expand Up @@ -104,23 +104,23 @@ module subroutine fragmentation_initialize(system, param, family, x, v, L_spin,
!write(*,*) 'Trying new arrangement'
end do
ke_avg_deficit = ke_avg_deficit / subtry
if (lfailure) write(*,*) 'Failed to find tangential velocities'
!if (lfailure) write(*,*) 'Failed to find tangential velocities'

if (.not.lfailure) then
call calculate_system_energy(linclude_fragments=.true.)
ke_radial = -dEtot - Qloss
call set_fragment_radial_velocities(lfailure)
if (lfailure) write(*,*) 'Failed to find radial velocities'
! if (lfailure) write(*,*) 'Failed to find radial velocities'
if (.not.lfailure) then
call calculate_system_energy(linclude_fragments=.true.)
! write(*,*) 'Qloss : ',Qloss
! write(*,*) '-dEtot: ',-dEtot
! write(*,*) 'delta : ',abs((dEtot + Qloss))
if ((abs(dEtot + Qloss) > Etol) .or. (dEtot > 0.0_DP)) then
write(*,*) 'Failed due to high energy error: ',dEtot, abs(dEtot + Qloss) / Etol
!write(*,*) 'Failed due to high energy error: ',dEtot, abs(dEtot + Qloss) / Etol
lfailure = .true.
else if (abs(dLmag) / Lmag_before > Ltol) then
write(*,*) 'Failed due to high angular momentum error: ', dLmag / Lmag_before
!write(*,*) 'Failed due to high angular momentum error: ', dLmag / Lmag_before
lfailure = .true.
end if
end if
Expand Down

0 comments on commit 3774b40

Please sign in to comment.