anchor.diff:    elseif(flmatbasedxs) then
       !cross section values
       if(chxstype .ne. 'scatter') &
  -      siga = sigabsave  + sigameanadjust + sqrt(Coabs)  * KL_sum
  +      siga = tsigabsave  + tsigameanadjust + sqrt(Coabs)  * KL_sum
       if(chxstype .ne. 'absorb') &
  -      sigs = sigscatave + sigsmeanadjust + sqrt(Coscat) * KL_sum
  -
  +      sigs = tsigscatave + tsigsmeanadjust + sqrt(Coscat) * KL_sum
       !determine point value
       select case (chxstype)
         case ("total")
  
positive.code: 
  
  
  
    function KLrxi_integral(j,xl,xr,chxstype,tnumEigsin) result(KL_int)
    !This function integrates on KL reconstructed realizations from xl to xr.
    !Integration is always on total cross section, but 'material'- or 'totxs'-
    !cross sections are adjusted with or without meanadjust.
    use genRealzvars, only: lamc, sigave, Coscat, Coabs, scatrat, sigscatave, sigabsave
    use KLvars, only: alpha, Ak, Eig, numEigs, KLrxivals, meanadjust, flmatbasedxs, &
                      sigsmeanadjust, sigameanadjust
  
negative.code: 
  
  
    function KLrxi_point(j,xpos,chxstype,tnumEigsin,orderin) result(KL_point)
    !Evaluates KL reconstructed realizations at a given point.
    !It has options for total, scattering only, absorption only, or scattering ratio.
    !It has an option to solve for less than the available number of eigenvalues.
    !It can solve any derivative order of the KL process with optional argument 'orderin'.
    !It can function when in 'material'-based or 'totxs'-based mode.
    !It can function when adjusting mean or not adjusting mean.
    use genRealzvars, only: lamc, scatrat, Coscat, Coabs
    use KLvars, only: alpha, Ak, Eig, numEigs, KLrxivals, flmatbasedxs
    use utilities, only: Heavi
  
    integer :: j
  