anchor.diff:                    velocities_dispersions(1) = 50.0
                     velocities_dispersions(2) = 56.0
                     velocities_dispersions(3) = 34.0
  -              else
  +              else if (disk_belonging(i) /= 3) then
                     write(6, *) "Error in velh:"
                     write(6, *) "    Number of WD's =", numberOfWDs
                     write(6, *) "    index =", i
  
positive.code:      &            uu(MAX_STARS_COUNT),
       &            vv(MAX_STARS_COUNT),
       &            ww(MAX_STARS_COUNT),
       &            velocities_dispersions(3)
            character(len = 6) :: geometry
            double precision :: coordinate_R(MAX_STARS_COUNT),
       &                        coordinate_Theta(MAX_STARS_COUNT),
       &                        coordinate_Zcylindr(MAX_STARS_COUNT)
            
            common /vel/ uu, vv, ww
            common /coorcil/ coordinate_R,
  
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,
  