Nbody6
 All Files Functions Variables
premsf.f
Go to the documentation of this file.
1  REAL*8 FUNCTION premsf(M,T,RR)
2  REAL*8 rzams,m,t,rr
3  REAL*8 rzamsf
4  EXTERNAL rzamsf
5 *
6  rzams = rzamsf(m)
7 * Coefficients in PREMSR
8  IF(m.LE.1.d0)THEN
9  pre1 = 0d0
10  pre2 = 0d0
11  pre3 = 7.432d-2 - 9.43d-2*m + 7.439d-2*m**2
12  ENDIF
13  IF(m.GT.1d0.AND.m.LT.2d0)THEN
14  pre1 = -4.00772d0 + 4.00772d0*m
15  pre2 = 8.5656d0 - 8.5656d0*m
16  pre3 = -4.50678d0 + 4.56118d0*m
17  ENDIF
18  IF(m.GE.2d0)THEN
19  pre1 = 1.60324d0 + 2.20401d0*m - 0.60433d0*m**2
20  & + 5.172d-2*m**3
21  pre2 = -4.56878d0 - 4.05305d0*m + 1.24575*m**2
22  & - 0.10922d0*m**3
23  pre3 = 3.01153 + 1.85745*m - 0.64290d0*m**2
24  & + 5.759d-2*m**3
25  END IF
26  premsf = rzams*10**((pre1*t**3 + pre2*t**4
27  & + pre3*t**5)/(1.05 - t)) - rr
28  RETURN
29  END