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:                   velocities_dispersions(1) = 50.0
                    velocities_dispersions(2) = 56.0
                    velocities_dispersions(3) = 34.0
                else if (disk_belonging(i) /= 3) then
                    write(6, *) "Error in velh:"
                    write(6, *) "    Number of WD's =", numberOfWDs
                    write(6, *) "    index =", i
  