anchor.diff:                      nu_c = ee_SY06(ne,T)
                   end select
                   nu = nu + nu_c
  -                kappa% ee = kappa_e_pre/nu_c
  +                if (nu_c > tiny(1.0_dp)) kappa% ee = kappa_e_pre/nu_c
               end if
               if (K_components(icond_ei)) then
                   nu_c = eion( &
  
positive.code:                 v = sqrt(1.0-tau)*(1.456-0.157/sqrt(tau) + 1.764/tau)
                  R = (0.4186+sqrt(1.007**2+(0.5010*v)**2))**2.5 * &
                  &   exp(1.456-sqrt(1.456**2+v**2))
                  if (R < 1.0E-10) R = 0.0
              end if
              
              if (K_components(icond_nQ)) then
  
negative.code:         real(dp) :: nn, nion, ne, kF, xF, eF, Gamma_e, kappa_e_pre
          real(dp) :: kappa_n_pre, kfn, mnstar, tau, v, R
          real(dp) :: nu, nu_c, K_opacity
          real(dp), parameter :: eQ_threshold = 1.0e-8_dp
          real(dp), parameter :: sf_reduction_threshold = 1.0e-10_dp
          integer :: ierr
  
          call clear_kappa
  