anchor.diff:            common /index/ flagOfWD,
        &                   numberOfWDs,
        &                   disk_belonging
  +          common /lb/ lgac, bgac
       
   C         Peculiar solar velocities, values from Anguiano et al. 2017
   C         TODO: find out why signs are different from values in Python
  
positive.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
            double precision, parameter :: solar_galactocentric_distance 
       &                                   = 8.5d0
            integer :: iseed,
  
negative.code:      &            uu(MAX_STARS_COUNT),
       &            vv(MAX_STARS_COUNT),
       &            ww(MAX_STARS_COUNT),
       &            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),
       &                        lgac(MAX_STARS_COUNT),
       &                        bgac(MAX_STARS_COUNT)
            
            common /vel/ uu, vv, ww
            common /coorcil/ coordinate_R,
  