anchor.diff:  REAL, DIMENSION(6) :: bvec, qvec
   
   ! Transverse Leakage Moments(0, Lx1, Ly1, Lz2, Lx2, Ly2, Lz3)
  -REAL, DIMENSION(7) :: Lm
  +REAL, DIMENSION(nnod,7) :: Lm
   
   ! To store old fluxes
   REAL, DIMENSION(nnod) :: flxc
  
positive.code:     flxc = f0(:,g)
      DO n = 1, nnod
  
              IF (ix(n) == ystag(iy(n))%smax) THEN                          ! East (X+) BC
                  CALL bcond(xeast, n, g, 1)
              ELSE
  
negative.code: 
  
              ! Update transverse leakage moments
              CALL TLUpd (n, g, Lm(n,:))
  
        END DO
  
        DO n = nnod,1,-1
              CALL matvec(nod(n,g)%P, nod(n,g)%ji, bvec)
  
              CALL matvec(nod(n,g)%R, nod(n,g)%Q - Lm(n,:), qvec)
  
              ! Update outgoing partial currents
              nod(n,g)%jo = qvec+bvec
  
              ! Update zeroth transverse leakages
              CALL LxyzUpd(n,g)
  
              ! Update flux and flux moments
              CALL FluxUpd4(n, g, Lm(n,:))
      END DO
  
      CALL RelE(f0(:,g), flxc, fer)
  