*     Include file for N-body problem in HW 02 Question 3.
*
* PARAMETER Declarations

      real*8 G_univ  ! Universal gravity constant
      parameter ( G_univ        = 6.67428D-20          )  ! km**2/(kg sec**2)

      integer*4 max_body  ! Maximum number of bodies allowed
      parameter ( max_body = 1000 ) 

****  Declare the variables in the common block for this program.

      integer*4 num_body ! Number of bodies
     .,         close_body(2) ! Closest pair of bodies

      real*8 dur   ! Duation of integration in days
     .,      tol   ! Tolerance for intergation (factional)
     .,      closest  ! Closest approach (km) 

      real*8 mass(max_body)   ! Masses of bodies
     .,      posi(2,max_body) ! Initial XY coordinates (km)
     .,      veli(2,max_body) ! INitial XY velocities (km/s)


      character*128 infile  ! Name of input file with mass pos vel of each body

      common / NbodyI / num_body, close_body
      common / NbodyR / mass, posi, veli, dur, tol, closest
      common / NbodyC / infile