anchor.diff:  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, DIMENSION(nnod) :: ss0                   ! Scattering source
  -DOUBLE PRECISION, DIMENSION(nnod) :: ssx1, ssy1, ssz1
  -DOUBLE PRECISION, DIMENSION(nnod) :: ssx2, ssy2, ssz2      ! Scattering source moments
   DOUBLE PRECISION :: f, fc                                  ! new and old integrated fission sources
   DOUBLE PRECISION :: domiR, e1, e2
   INTEGER :: g
  
positive.code:     erro = errn       ! Save old fission source error/difference
      DO g = 1, ng
          !!!Calculate total source
          CALL TSrc2(g, Keo)
          !!!Inner Iteration
          CALL inner4(g)
      END DO
  
negative.code: END SUBROUTINE TSrcT
  
  
  SUBROUTINE TSrcAd(g, Keff)
  !
  ! Purpose:
  !   To update total source
  !
  
  USE sdata, ONLY: nod, nuf, nnod, fs0, fsx1, fsy1, fsz1, fsx2, fsy2, fsz2, &
                   f0, fx1, fy1, fz1, fx2, fy2, fz2, ng, sigs
  
  IMPLICIT NONE
  