Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
c%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
c
c MERCURY.INC (ErikSoft 4 March 2001)
c
c%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
c
c Author: John E. Chambers
c
c Parameters that you may want to alter at some point:
c
c NMAX = maximum number of bodies
c CMAX = maximum number of close-encounter minima monitored simultaneously
c NMESS = maximum number of messages in message.in
c HUGE = an implausibly large number
c NFILES = maximum number of files that can be open at the same time
c
integer NMAX, CMAX, NMESS, NFILES
real*8 HUGE
c
parameter (NMAX = 1000)
parameter (CMAX = 1000)
parameter (NMESS = 200)
parameter (HUGE = 9.9d29)
parameter (NFILES = 90)
c
c------------------------------------------------------------------------------
c
c Constants:
c
c DR = conversion factor from degrees to radians
c K2 = Gaussian gravitational constant squared
c AU = astronomical unit in cm
c MSUN = mass of the Sun in g
c
real*8 PI,TWOPI,PIBY2,DR,K2,AU,MSUN
c
parameter (PI = 3.141592653589793d0)
parameter (TWOPI = PI * 2.d0)
parameter (PIBY2 = PI * .5d0)
parameter (DR = PI / 180.d0)
parameter (K2 = 2.959122082855911d-4)
parameter (AU = 1.4959787d13)
parameter (MSUN = 1.9891d33)