anchor.diff:      real(dp), parameter, dimension(number_layers+1) :: Xn = [ &
       &   0.0, 0.0, 0.0, 0.0, 0.0,  &
       &   0.07, 0.18, 0.29, 0.39, 0.45, 0.50, 0.55, 0.61, 0.70, 0.73, 0.76, 0.80, 0.80]
  -    
  +
  +    
  +    integer, parameter :: number_HZ08_layers = 29
  +    
  +    character(len=iso_name_length), parameter, dimension(number_HZ08_layers+1) :: HZ08_ion_composition = [ &
  +    &   character(len=iso_name_length) ::  &
  +    &   'pd106', 'ru106', 'mo106', 'zr106', 'sr106', 'kr106', 'se106', 'ge106', 'ni92', 'fe86', 'cr80', 'ti74', &
  +    &   'ca68', 'ar62', 's56', 'si50', 'mg42', 'ca72', 'ar66', 's60', 'si54', 'cr92', 'ti86', 'ca80', 'ar74', &
  +    &   's68', 'ni124', 'fe120', 'cr118', 'ti116' ]
  +!     real(dp), parameter, dimension(number_HZ08_layers) :: HZ08_transition_pressures = [ &
  +!     & ]
  +!     real(dp), parameter, dimension(number_HZ08_layers+1) :: HZ08_Xn = [ &
  +!     & ]
   contains
       
       subroutine do_make_crust(lgP, Yion, Xneut, charged_ids, ncharged, ion_info)
  
positive.code:     &   'cr56', &
      &   'fe56' ]
      
      integer, parameter :: number_layers = 17
      
      character(len=iso_name_length), parameter, dimension(number_layers+1) :: ion_composition = [ &
  
negative.code:     real(dp), parameter :: transition_width = 0.02  
      
      integer, parameter :: HZ90_number = 19
      integer, parameter :: HZ08_number = 34
      character(len=iso_name_length), parameter, dimension(HZ90_number) :: HZ90_network = [ character(len=iso_name_length) :: &
      &   'n', &
      &   'mg40', &
  