*       Include GR shrinkage to eliminate small circularized orbits.
  100 CONTINUE
      IF ((KSTAR(I).EQ.10.OR.KSTAR(I).EQ.0).AND.SEMI*SU.LT.6.0) THEN
          KSPAIR = IPAIR
          CVEL = 3.0D+05/VSTAR
*         ZFAC = BODY(I1)*BODY(I1+1)*BODY(I)*SMU**3
          ZFAC = BODY(I1)*BODY(I1+1)*BODY(I)
          ADOT =  64.0/5.0*ZFAC/(CVEL**5*SEMI**3)
*       Shrink circular orbit according to Peters 1964.
          DT0 = 1000.0*DT0
          SEMI0 = SEMI - ADOT*DT0
          DA = ADOT*DT0/SEMI0
*       Employ contraction for safety scaling of U & UDOT.
          CALL EXPAND2(IPAIR,SEMI0)
