intide.f

Go to the documentation of this file.

      SUBROUTINE intide
*
*
*       Input & scaling for tidal dissipation.
*       --------------------------------------
*
      include 'common6.h'
*
*
*       Read parameters for tidal capture simulation.
      READ (5,*)  rstar, ims, iev, rms, rev
*
*       Convert radii from S.U. to internal units.
      rstar = rstar/su
*
      WRITE (6,10)  rstar, ims, iev, rms, rev
   10 FORMAT (/,12x,'OLD TIDAL:    RSTAR =',1pe8.1,'  IMS =',0p,i5,
     &              '  IEV =',i4,'  RMS/R* =',f6.2,'  REV/R* =',f6.1)
*
*       Assign individual radii for main-sequence and evolved stars.
      DO 30 i = 1,n
*       Adopt a primitive scheme in case of no stellar evolution.
          IF (i.LE.ims) THEN
              radius(i) = rms*rstar
          ELSE
              radius(i) = rev*rstar
          END IF
   30 CONTINUE
*
*     WRITE (6,40)  BODY(1), BODY(N), RADIUS(1), RADIUS(N)
*  40 FORMAT (/,12X,'SCALED RADII:    M(1) =',F8.4,'  M(N) =',F8.4,
*    &                                '  R(1) =',1PE8.1,'  R(N) =',E8.1)
*
      RETURN
*
      END