Skip to content

Commit

Permalink
(build) bug fix with line truncation in gfortran v14
Browse files Browse the repository at this point in the history
  • Loading branch information
danieljprice committed May 17, 2024
1 parent cfdbddb commit 70ee53b
Showing 1 changed file with 5 additions and 4 deletions.
9 changes: 5 additions & 4 deletions src/setup/set_hierarchical.f90
Original file line number Diff line number Diff line change
Expand Up @@ -575,13 +575,14 @@ subroutine set_multiple(m1,m2,semimajoraxis,eccentricity, &
q_comp = mtot/m_comp
if (q_comp>1) q_comp=q_comp**(-1)

! Mardling&Aarseth (2001) criterion check

period_ratio = sqrt((a_comp*a_comp*a_comp)/(m_comp+mtot)/(semimajoraxis*semimajoraxis*semimajoraxis)*(mtot)) ! Po/Pi
! Mardling & Aarseth (2001) criterion check
period_ratio = sqrt((a_comp*a_comp*a_comp)/(m_comp+mtot)/&
(semimajoraxis*semimajoraxis*semimajoraxis)*(mtot)) ! Po/Pi
criterion = 4.7*(1-e_comp)**(-1.8)*(1+e_comp)**(0.6)*(1+q_comp)**(0.1)

if (criterion > period_ratio) then
print "(1x,a)",'WARNING: set_multiple: orbital parameters does not satisfy Mardling and Aarseth stability criterion.'
print "(1x,a)",'WARNING: set_multiple: orbital parameters do not satisfy '//&
'Mardling & Aarseth stability criterion.'
endif

q2=m2/m1
Expand Down

0 comments on commit 70ee53b

Please sign in to comment.