mtrace.f

Go to the documentation of this file.

      SUBROUTINE mtrace(J,DM)
*
*
*       Orbit diagnostics for mass loss.
*       --------------------------------
*
      include 'common6.h'
      LOGICAL  first
      SAVE  first
      DATA  first /.true./
*
*
*       Open unit #14 the first time.
      IF (first) THEN
          OPEN (unit=14,status='NEW',form='FORMATTED',file='MDOT')
          first = .false.
*
*       Print cluster scaling parameters at start of the run.
          IF (nrg.EQ.0) THEN
              WRITE (14,1)  rbar, bodym*zmbar, body1*zmbar, tscale,
     &                      nbin0, nzero
    1         FORMAT (/,6x,'MODEL:    RBAR =',f5.1,'  <M> =',f6.2,
     &                     '  M1 =',f6.1,'  TSCALE =',f6.2,
     &                     '  NB =',i4,'  N0 =',i6,//)
          END IF
*
          WRITE (14,2)
    2     FORMAT ('     TIME  NAME   K*   M    DM     r/Rc   VR',
     &            '     EI        N')
      END IF
*
*       Evaluate potential energy on GRAPE (#I: single star or binary c.m.).
      IF (j.GE.ifirst) THEN
          i = j
          CALL poti(i,potj)
      ELSE
          ks = kvec(j)
          i = n + ks
          CALL poti(i,potj)
      END IF
*
*       Obtain central distance (scaled by core radius) and velocities.
      ri2 = 0.0
      vi2 = 0.0
      vr = 0.0
      DO 10 k = 1,3
          ri2 = ri2 + (x(k,i) - rdens(k))**2
          vi2 = vi2 + xdot(k,i)**2
          vr = vr + (x(k,i) - rdens(k))*xdot(k,i)
   10 CONTINUE
      ri = sqrt(ri2)/rc
      vr = vr/sqrt(ri2)
*
*       Form binding energy per unit mass (note: POTJ < 0).
      ei = 0.5*vi2 + potj
*
*       Include optional external tidal field (note: HT includes mass).
      IF (kz(14).GT.0) THEN
          CALL xtrnlv(i,i)
          ei = ei + ht/(body(i) + 1.0e-20)
      END IF
*
      WRITE (14,20)  ttot, name(j), kstar(j), body(j)*zmbar, dm*zmbar,
     &               ri, vr, ei, n
   20 FORMAT (1x,f8.1,i6,i4,f7.1,f6.2,f7.2,f7.2,f8.3,i6)
      CALL flush(14)
*
      RETURN
*
      END