anchor.diff:          ! define bolometric magnitude
           m_bol = m_bol_sun - 2.5*logl
   
  -        if (trim(brand) == 'mamjek') then
  +        if (trim(bc_type) == 'mamajek') then
               call bc_mamjek(teff, logl, mag_length, 0, magnitudes)
               return
           end if
  
positive.code:     subroutine bc_mamajek(teff, logl, mag_length, init, magnitudes)
          use interpolate, only: lagrange
  
          integer :: i, j, k
  
          integer,  intent(in) :: init, mag_length
          real(dp), intent(in) :: teff, logl
          real(dp) :: mbol
          real(dp), parameter  :: m_bol_sun = 4.74
          real(dp), dimension(mag_length), intent(out) :: magnitudes
  
  
negative.code:                 call vc_semiemp()
              case ('mamajek')
                  call log_warn('Pecaut and Mamajek only valid for MS, solar metallicity')
                  call bc_mamajek(0.0_dp, 0.0_dp, 1, 1)
              case default
                  call log_warn('invalid bc_type in bc_init: default to marcs08')
                  call marcs(feh, afe)
  