anchor.diff:       &            uu(MAX_STARS_COUNT),
        &            vv(MAX_STARS_COUNT),
        &            ww(MAX_STARS_COUNT),
  -     &            velocities_dispersions(3)
  +     &            velocities_dispersions(3),
  +     &            stellar_galactocentric_distance,
  +     &            xx,
  +     &            sigo2,
  +     &            sigm2,
  +     &            sigr2,
  +     &            dsigr2dr,
  +     &            sigt2,
  +     &            sigr,
  +     &            sigt,
  +     &            vr,
  +     &            vth,
  +     &            vph,
  +     &            delta,
  +     &            sindel,
  +     &            cosdel,
  +     &            vx,
  +     &            vy,
  +     &            vz,
  +     &            sini,
  +     &            cosi
             character(len = 6) :: geometry
             double precision :: coordinate_R(MAX_STARS_COUNT),
        &                        coordinate_Theta(MAX_STARS_COUNT),
  -     &                        coordinate_Zcylindr(MAX_STARS_COUNT)
  +     &                        coordinate_Zcylindr(MAX_STARS_COUNT),
  +     &                        lgac(MAX_STARS_COUNT),
  +     &                        bgac(MAX_STARS_COUNT)
             
             common /vel/ uu, vv, ww
             common /coorcil/ coordinate_R,
  
positive.code:                   velocities_dispersions(1) = 50.0
                    velocities_dispersions(2) = 56.0
                    velocities_dispersions(3) = 34.0
                else
                    write(6, *) "Error in velh:"
                    write(6, *) "    Number of WD's =", numberOfWDs
                    write(6, *) "    index =", i
  
negative.code: C         A, B: Oort constants in km/sKpc (values by Kerr and 
  C              Lynden-Bell 1986)
            real, parameter :: oort_const_A = 14.4,
       &                       oort_const_B = -12.8,
       &                       ro = 10.5,
       &                       l = 5.5,
       &                       sigo = 80.0,
       &                       sigm = 145.0,
       &                       pi = 4.0 * atan(1.0),
       &                       vc = 220.0,
       &                       gamma = 3.4
            double precision, parameter :: solar_galactocentric_distance 
       &                                   = 8.5d0
            integer :: iseed,
  