anchor.diff:      fs0c  = fs0       ! Save old fission source
       erro = errn       ! Save old fission source error/difference
       DO g = 1, ng
  -        !!!Calculate Scattering source
  -        CALL SSrc(g, ss0, ssx1, ssy1, ssz1, ssx2, ssy2, ssz2)
           !!!Calculate total source
           CALL TSrcT(g, ht)
           !!!Inner Iteration
  
positive.code:     nod(n,g)%Q(7) = chi(mat(n),g) * fsz2(n) + sz2(n)
  END DO
  
  END SUBROUTINE TSrcFx2
  
  
  SUBROUTINE TSrcT(g, ht)
  
negative.code: DOUBLE PRECISION :: Keo                                    !Old Multiplication factor (Keff)
  DOUBLE PRECISION, DIMENSION(nnod) :: fs0c                  !Old fission source
  DOUBLE PRECISION, DIMENSION(nnod,ng) :: f0c                !Old flux
  DOUBLE PRECISION :: f, fc                                  ! new and old integrated fission sources
  DOUBLE PRECISION :: domiR, e1, e2
  INTEGER :: g
  