anchor.diff:      eF = sqrt(1.0+xF**2)
       Gamma_e = Gamma/ionic%Z53
       kappa_pre = onethird*pi**2*boltzmann**2*T*ne/Melectron/eF
  -    
  +    kappa_n_pre = 0.0
  +    
  +    ! electrons
       call clear_kappa
       nu_c = 0.0
       nu = 0.0
  
positive.code:     integer, intent(in) :: which_ee, which_eQ
      logical, dimension(num_conductivity_channels), intent(in) :: K_components
      real(dp) :: nn,nion, nu, nu_c, kappa_pre, ne, kF, xF, eF, Gamma_e
      
      ne = rho/amu*ionic%Ye
      kF = (threepisquare*ne)**onethird
  
negative.code:             kappa% eQ = -1.0
          end if
      end if
      
      ! neutrons
      call clear_kappa
      nu_c = 0.0
      nu_n = 0.0  
      if (ionic% Yn > 0.0) then
      nn = rho*ionic% Yn/(1.0-chi)/Mneutron / density_n
      kFn = (threepisquare*nn)**onethird
      xFn = hbar*kFn/Mneutron/clight
      eFn = sqrt(1.0+xFn**2)
      nion = (1.0-ionic%Yn)*rho/Mneutron/ionic% A /density_n    
      kappa_n_pre = onethird*pi**2*boltzmann**2*T*nn/Mneutron/eFn      
      if (K_components(icond_sf)) then
          kappa% sf =  sPh(nn,nion,T,ionic)
      end if
  	if (K_components(icond_nQ)) then
      	nu_c = n_imp(nn,nion,T,ionic)
      	kappa% nQ = kappa_n_pre/nu_c
          nu_n = nu_n + nu_c
      end if    
   	if (K_components(icond_np)) then
          nu_c = n_phonon(nn,nion,T,ionic)
          kappa% np = kappa_n_pre/nu_c
          nu_n = nu_n + nu_c
      end if      
  		kappa% total = kappa_pre/nu + kappa_n_pre/nu_n + kappa% sf
  	else
  		kappa% total = kappa_pre/nu 
  	end if
      
      contains
      subroutine clear_kappa()
  