https://github.com/florentrenaud/nbody6tt
Raw File
Tip revision: 8a4716382ead3ece116c48a4ae5c65f8c9534437 authored by Florent on 29 January 2015, 12:19:28 UTC
Nbody6 - 29 January 2015 (added GPU2/Build/.gitkeep)
Tip revision: 8a47163
intide.f
      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
back to top