anchor.diff:          kappa% ei = 0.0
           kappa% eQ = 0.0
           kappa% sf = 0.0     
  +        kappa% nQ = 0.0
  +        kappa% np = 0.0
       end subroutine clear_kappa
   end subroutine conductivity
   
  
positive.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
  
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
  