anchor.diff:          &   (trim(name(j)),X(j),j=1,2)
           write(output_unit,'(2(a,f6.2,tr4))') '<Z> = ',ionic%Z,'<A> = ',ionic% A
           write (output_unit, &
  -        &   '(2a6,6a11,a7/,2("======"),6("==========="),"=======")') &
  +        &   '(2a6,7a11,a7/,2("======"),7("==========="),"=======")') &
               & 'lg(r)','lg(T)','Gamma','eta_e', &
  -            & 'K_ee','K_ei','K_tot','K_table','diff'
  +            & 'K_ee','K_ei','K_eQ','K_tot','K_table','diff'
               
  -        temperature: do j = 1, 5
  -            lgT = real(j-1,dp)*0.5_dp + 7.0_dp
  -            T = 10.0_dp**lgT
  -            
  -            density: do i = 1, 13
  -                lgrho = real(i-1,dp)/3.0_dp + 5.0_dp
  -                rho = 10.0_dp**lgrho
  +        temperature: do j = 1, NT
  +            T = 10.0_dp**lgT(j)
  +            density: do i = 1, Nrho
  +                rho = 10.0_dp**lgrho(i)
                   chi = use_default_nuclear_size
                   call eval_crust_eos(eos_handle,rho,T,ionic, &
                       & ncharged, charged_ids, Yion, Tcs,  &
  
positive.code:     real(dp), dimension(num_dStar_eos_results) :: res
      character(len=*), parameter :: datadir = '../../data/conductivity'
      type(electron_conductivity_tbl), pointer :: tab
      real(dp) :: lgrho, lgT, rho, T, eta, Gamma, mu_e, K_e
      type(conductivity_components) :: kappa
      type(crust_eos_component), dimension(num_crust_eos_components) :: &
      &   eos_components
  
negative.code:     use dStar_eos_lib
      use conductivity_lib
      use PPP_electron
      
      implicit none
      integer, parameter :: NX=11,Nrho=13,NT=5
      real(dp), parameter :: Xmin=0.0_dp, Xmax = 1.0_dp
      real(dp), parameter :: lgrho_min=5.0_dp, lgrho_max = 9.0_dp
      real(dp), parameter :: lgT_min=7.0_dp, lgT_max = 9.0_dp
      integer :: eos_handle,cond_handle,ierr,i,j,k
      integer :: Z(2), N(2)
      real(dp) :: A(2), X(2)
  