setup2.f

Go to the documentation of this file.

      SUBROUTINE setup2
*
*
*       Initial conditions in astrophysical units (#22 = -1).
*       -----------------------------------------------------
*
      include 'common6.h'
*
*
*       Define GM & PC in cgs units and AU in pc (if preferred).
      gm = 6.67d-08*1.989d+33
      pc = 3.0857d+18
      au = 2.0627e+05
*       Define RBAR (use appropriate value if N-body units skipped in SCALE).
      rbar = 1.0
*
*       Read initial conditions from unit #10.
      zmass = 0.0
      DO 10 i = 1,n
          READ (10,*)  body(i), (x(k,i),k=1,3), (xdot(k,i),k=1,3)
          zmass = zmass + body(i)
   10 CONTINUE
*
*       Form velocity scale in km/sec from mass in Msun & X,Y,Z in pc.
      vstar = 1.0d-05*sqrt(gm/pc)
      vstar = vstar*sqrt(zmass/rbar)
*
*       Convert to N-body units before optional virial theorem scaling.
      DO 20 i = 1,n
          body(i) = body(i)/zmass
          DO 15 k = 1,3
              x(k,i) = x(k,i)/rbar
              xdot(k,i) = xdot(k,i)/vstar
   15     CONTINUE
   20 CONTINUE
*       Note: setting KZ(5) = 3 skips scaling to N-body units (ensure RBAR).
*
      RETURN
*
      END