anchor.diff:  END IF
   
   1201 FORMAT(2X, 'NUMBER OF CONTROL ROD BANK  :', I3)
  -1216 FORMAT(2X, 'MAX. NUMBER OF STEPS  :', I4)
  +1216 FORMAT(2X, 'MAX. NUMBER OF STEPS        :', I4)
   1202 FORMAT(2X, 'FULLY INSERTED POSITION (cm): ', F4.1, ' (FROM BOTTOM OF THE CORE)')
   1203 FORMAT(2X, 'STEP SIZE (cm)              : ', F8.4)
   1204 FORMAT(2X, 10(:, 2X, 'Bank ', I2))
  
positive.code: !    To read control rod position
  
  USE sdata, ONLY: nx, ny, nmat, ng, xdiv, ydiv, &
                   nxx, nyy, bpos, nb, &
                   dsigtr, dsiga, dnuf, dsigf, dsigs
  
  IMPLICIT NONE
  
negative.code: INTEGER, INTENT(IN) :: xbunit
  
  INTEGER :: ln   !Line number
  INTEGER :: ios, iosc  ! IOSTAT status
  
  INTEGER :: i, j, g, h
  INTEGER, DIMENSION(nx, ny) :: bmap       ! Radial control rod bank map (assembly wise)
  