swift/symba5/symba5_nbodm.f

c*************************************************************************
c                            SYMBA5_NBODM.F
c*************************************************************************
c Returns the location of the last massive body in the list
c
c             Input:
c                 nbod          ==>  number of massive bodies (int scalar)
c                 mass          ==>  mass of bodies (real array)
c                 mtiny         ==>  Small mass  (real array)
c             Output:
c                 nbodm         ==>  location of the last massive body 
c                                    (int scalar)
c
c Remarks:  If all the objects are massive, then nbodm=nbod-1 so that
c           the do loops will have the correct limits.
c Authors:  Hal Levison
c Date:    3/20/97
c Last revision: 1/29/06

      subroutine  symba5_nbodm(nbod,mass,mtiny,nbodm)

      include '../swift.inc'
      include 'symba5.inc'

c...  Inputs Only: 
      integer nbod
      real*8 mass(nbod),mtiny

c...  Outputs only
      integer nbodm

c...  Internals
      integer i

c----
c...  Executable code 

      if(mass(nbod).gt.mtiny) then
         nbodm = nbod - 1
         write(*,*) '    Of ',nbod,' objects, all are massive. '
      else
         nbodm = 1
         do i=1,nbod-1
            if(mass(i).gt.mtiny) then
               nbodm = i
            endif
         enddo
         write(*,*) '    Of ',nbod,' objects, ',nbodm,' are massive. '
      endif

      return
      end             ! symba5_nbodm.f
c--------------------------------------------------------