anchor.diff:      name(1:2)(1:1) = [ (StrUpCase(name(j)(1:1)), j=1,2) ]
       Z1 = real(Z(1),dp); Z2 = real(Z(2),dp); A1 = real(A(1),dp); A2 = real(A(2),dp)
   
  +    ! set grid axes
       X1 = linspace(NX,Xmin,Xmax)
       lgrho = linspace(Nrho,lgrho_min,lgrho_max)
       lgT = linspace(NT,lgT_min,lgT_max)
  
positive.code:     type(composition_info_type) :: ionic
      real(dp) :: chi, Xsum
      real(dp), dimension(num_dStar_eos_results) :: res
      character(len=*), parameter :: datadir = '../../data/conductivity'
      type(electron_conductivity_tbl), pointer :: tab
      real(dp), dimension(NX) :: X1
      real(dp), dimension(Nrho) :: lgrho
  
negative.code:         call compute_composition_moments(2,chem_ids,Y,ionic,Xsum, &
              & ncharged, charged_ids, Yion, exclude_neutrons=.TRUE.)
          
          ! table header
          write(output_unit,'(/,2("X(",a,") = ",f3.1,tr4))') &
          &   (trim(name(j)),X(j),j=1,2)
          write(output_unit,'(2(a,f6.2,tr4))') '<Z> = ',ionic%Z,'<A> = ',ionic% A
  