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