anchor.diff:                  &   chi,Gamma,eta,mu_e,ionic,Tcs(neutron_1S0),kappa)
                   call eval_PPP_electron_table(rho,T,sqrt(ionic% Z2),K_e,ierr)
                   diff(i,j,k) = (K_e - kappa% electron_total)/kappa% electron_total
  -                write (output_unit, '(3f6.2,6es11.3,f7.3)') &
  -                    & lgrho,lgT,ionic%Z, &
  +                write (output_unit, '(2f6.2,6es11.3,f7.3)') &
  +                    & lgrho,lgT, &
                       & Gamma,mu_e*mev_to_ergs/boltzmann/T, &
                       & kappa%ee,kappa%ei,kappa%electron_total, K_e, diff(i,j,k)
               end do density
  
positive.code:         
          write(output_unit,'(/,2("X(",a,") = ",f3.1,tr4))') &
          &   (trim(name(j)),X(j),j=1,2)
          write (output_unit, &
          &   '(3a6,6a11,a7/,3("======"),6("==========="),"=======")') &
              & 'lg(r)','lg(T)','<Z>','Gamma','eta_e', &
              & 'K_ee','K_ei','K_tot','K_table','diff'
              
          temperature: do j = 1, 5
  
negative.code:     do j = 1,2
          name(j)(1:1) = StrUpCase(name(j)(1:1))
      end do
  
      composition: do k = 1, 11
          X(1) = 0.1_dp*real(k-1,dp)
          X(2) = 1.0_dp -X(1)
  