anchor.diff:  		    WRITE(ounit,1020)i, g
   			STOP
   		END IF
  -		IF (xnuf(i,g) > 0.) cnuf = .FALSE.
  -		IF (xsigf(i,g) > 0.) csigf = .FALSE.
  +		IF (xnuf(i,g) > 0.) ccnuf = .FALSE.
  +		IF (xsigf(i,g) > 0.) ccsigf = .FALSE.
   		
   		
   		xD(i,g) = 1.d0/(3.d0*xsigtr(i,g)) 
  
positive.code: INTEGER :: i, j, k, n, ztot
  INTEGER :: ly, lx, ys, xs, yf, xf, lz
  REAL :: summ, vsumm
  REAL, DIMENSION(nzz) :: faxi
  REAL :: totp
  INTEGER :: nfuel
  REAL :: fmax
  
negative.code: 		END DO
  	END DO
  END DO
  
  END SUBROUTINE inp_geom2
  
  
  SUBROUTINE misc ()
  !
  ! Purpose:
  !    To assign material xsec to nodes
  !    To arranges the nodes into 1-D array instead of 3-D array
  !
  
  
  USE sdata, ONLY: nxx, nyy, nzz, ix, iy, iz, xyz, &
  				 nnod, sigtr, siga, nuf, sigf, &
                   chi, sigs, D, sigr, ng, ystag, &
  				 xdel, ydel, zdel, vdel, mode, &
  				 mat, xD, xsigr, xchi, &
  				 bcon, ftem, mtem, cden, bpos
  
  IMPLICIT NONE
  
  INTEGER :: i, j, k, n
  
  ALLOCATE(ix(nnod), iy(nnod), iz(nnod))
  ALLOCATE(xyz(nxx, nyy, nzz))
  