anchor.diff:              kappa% eQ = -1.0
           end if
       end if
  +    kappa% total = kappa_pre/nu 
       
       ! 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
  
positive.code:         kappa% ei = 0.0
          kappa% eQ = 0.0
          kappa% sf = 0.0     
      end subroutine clear_kappa
  end subroutine conductivity
  
  
negative.code:     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
  