Nbody6
 All Files Functions Variables
intide.f
Go to the documentation of this file.
1  SUBROUTINE intide
2 *
3 *
4 * Input & scaling for tidal dissipation.
5 * --------------------------------------
6 *
7  include 'common6.h'
8 *
9 *
10 * Read parameters for tidal capture simulation.
11  READ (5,*) rstar, ims, iev, rms, rev
12 *
13 * Convert radii from S.U. to internal units.
14  rstar = rstar/su
15 *
16  WRITE (6,10) rstar, ims, iev, rms, rev
17  10 FORMAT (/,12x,'OLD TIDAL: RSTAR =',1pe8.1,' IMS =',0p,i5,
18  & ' IEV =',i4,' RMS/R* =',f6.2,' REV/R* =',f6.1)
19 *
20 * Assign individual radii for main-sequence and evolved stars.
21  DO 30 i = 1,n
22 * Adopt a primitive scheme in case of no stellar evolution.
23  IF (i.LE.ims) THEN
24  radius(i) = rms*rstar
25  ELSE
26  radius(i) = rev*rstar
27  END IF
28  30 CONTINUE
29 *
30 * WRITE (6,40) BODY(1), BODY(N), RADIUS(1), RADIUS(N)
31 * 40 FORMAT (/,12X,'SCALED RADII: M(1) =',F8.4,' M(N) =',F8.4,
32 * & ' R(1) =',1PE8.1,' R(N) =',E8.1)
33 *
34  RETURN
35 *
36  END