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% 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
  
negative.code:         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
  
  