anchor.diff:      logical, dimension(num_conductivity_channels), parameter :: &
       &   cond_exclude_neutrons = [ .TRUE., .TRUE., .TRUE., &
       &   .FALSE., .FALSE., .FALSE., .TRUE. ]
  +    logical, dimension(num_conductivity_channels), parameter :: &
  +    &   cond_exclude_radiation = [ .TRUE., .TRUE., .TRUE., &
  +    &   .TRUE., .TRUE., .TRUE., .FALSE. ]
   
       ! flag to control which ee scattering fmla. is used. 
       ! default: Shternin & Yakovlev (2006)
  
positive.code:     logical, dimension(num_conductivity_channels), parameter ::  &
      &   cond_use_all = [ .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE.]
      logical, dimension(num_conductivity_channels), parameter :: &
      &   cond_use_only_conduction = [ .TRUE., .TRUE., .TRUE., &
      &   .FALSE., .FALSE., .FALSE., .FALSE. ]
      logical, dimension(num_conductivity_channels), parameter :: &
      &   cond_use_only_kap = [ .FALSE., .FALSE., .FALSE., &
  
negative.code:     use constants_def, only: dp
  
      ! used for mask array to control which components are included
      ! in regions where a tabulated electron conduction is used, 
      ! only the total can be used. May need to rethink this control.
      !
      integer, parameter :: &
      &       icond_ee    = 1,     &  ! electron conduction, electron scattering
      &       icond_ei    = 2,     &  ! electron conduction, ion scattering
  