anchor.diff:          integer,  intent(in) :: init, mag_length
           real(dp), intent(in) :: teff, logl
           real(dp) :: m_bol
  -        real(dp), parameter  :: m_bol_sun = 4.74
  +        real(dp), parameter    :: m_bol_sun = 4.74
  +        real(dp), dimension(4) :: coeffs
           real(dp), dimension(mag_length), intent(out) :: magnitudes
   
           if (init == 1) then
  
positive.code:         integer :: i
  
          real(dp), intent(in) :: feh, afe
  
          character(len=15), intent(in) :: brand
          character(len=1), dimension(:), intent(in) :: filters
  
negative.code:             write(*, *) 'BC tables were not initialized.'
          else
              ! interpolate using 4-point lagrange interpolation
              m_bol = m_bol_sun - 2.5*logl
  
              ! hunt for Teff index
              do i = 1, n_teffs
  