anchor.diff:          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 
  +		kappa% total = kappa% total + kappa_n_pre/nu_n + kappa% sf
   	end if
       
       contains
  
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
  