anchor.diff:  C     Distances in Kpc
   C     Velocities in Km/s
   C     Proper motions in  arcsec/yr
  +C     More info on conversions at: https://physics.stackexchange.com/questions/88663/converting-between-galactic-and-ecliptic-coordinates
         implicit none
   
         integer, parameter :: MAX_STARS_COUNT = 6000000
  
positive.code:       double precision :: solarGalactocentricDistance,
       &                    ros,
       &                    cos_bgac,
       &                    zz,
       &                    sin_bgac,
       &                    zzx
        double precision :: coordinate_R(MAX_STARS_COUNT),
  
negative.code:      &                   PI = 4.0 * atan(1.0),
       &                   FI = 180.0 / PI,
       &                   PC_PER_KPC = 1e3
  C     THETA - see BK angle at the link above
  C     DELTA_G - declination of NGP 
  C     TODO: find out why we use NGP declination for J1950 and not J2000 
        double precision, parameter :: DELTA_G = 0.478d0, 
       &                               THETA = 2.147d0,
       &                               SIN_DELTA_G = dsin(DELTA_G),
  