anchor.diff:        integer, parameter :: NCOL = 6,
        &                      NCOL2 = 5,
        &                      NROW = 300
  +      real, parameter :: SOLAR_MASS_KG = 1.989e30,
  +     &                   GRAVITATIONAL_CONST_CM_S_KG = 6.67e-5
   
         integer :: i,
        &           model,
  
positive.code:      &        cuv,
       &        zlte,
       &        zlr,
       &        G,
       &        luminosity,
       &        ubvri_u,
       &        ubvri_b,
  
negative.code: 
        effective_temperature = 10.0 ** zlte
  
  C     TODO: find out the meaning of 6.696e10 constant
        surface_gravity_logarithm = log10(GRAVITATIONAL_CONST_CM_S_KG) 
       &                            + log10(wd_mass * SOLAR_MASS_KG) 
       &                            - 2.0 * (zlr + log10(6.696e10))
        ubvri_u = cuv + mv
        ubvri_b = cbv + mv
  