renew.f

Go to the documentation of this file.

      SUBROUTINE renew(I)
*
*
*       Re-activation of dormant KS binary.
*       -----------------------------------
*
      include 'common6.h'
      common/chainc/  xc(3,ncmax),uc(3,ncmax),bodyc(ncmax),ich,
     &                listc(lmax)
*
*
*       Restore masses of binary components and initialize time.
      ipair = i - n
      i1 = 2*ipair - 1
      i2 = i1 + 1
      body(i1) = bodyc(9)
      body(i2) = bodyc(10)
      t0(i1) = time
*
*       Specify semi-major axis and Kepler period.
      semi = -0.5d0*body(i)/h(ipair)
      tk = twopi*semi*sqrt(semi/body(i))
*
*       Estimate consistent unperturbed interval.
      ip = ipair
      CALL tpert(ipair,gmin,dt)
      ipair = ip
*
*       Specify unperturbed motion (subject to DT > 0).
      k = 1 + int(0.5d0*dt/tk)
      k = max(k,1)
      step(i1) = float(k)*tk
*
*       Obtain coordinates & velocities of unperturbed binary.
      CALL resolv(ipair,1)
*
*       Copy global indices of all chain members to JPERT for routine NBPOT.
      DO 1 l = 1,nch
          jpert(l) = jlist(l)
    1 CONTINUE
*
*       Evaluate potential energy with respect to inner components.
      jlist(1) = i1
      jlist(2) = i2
      CALL nbpot(2,nch,pot1)
*
*       Include interaction of inner c.m. & other members to give net effect.
      jlist(1) = i
      CALL nbpot(1,nch,pot2)
      dphi = pot1 - pot2
*
*       Form new KS polynomials if nearest perturber is significant (> GMIN).
      IF (dt.LT.0.0) THEN
          CALL kslist(ipair)
          CALL kspoly(ipair,1)
      END IF
*
*       Subtract binding energy from temporary save in ECOLL.
      eb = body(i1)*body(i2)*h(ipair)/body(i)
      ecoll = ecoll - eb + dphi
*
      IF (kz(30).GT.1) THEN
          WRITE (6,5)  ipair, eb, r(ipair), gamma(ipair), dphi, k
    5     FORMAT (' RENEW:    KS =',i3,'  EB =',f8.2,'  R =',1p,e8.1,
     &            '  G =',e8.1,'  DPHI =',e8.1,'  DT/TK =',0p,i4)
      END IF
*
      RETURN
*
      END