SUBROUTINE SETBIT(ITYPE,JNUM)
      KTYPE=ITYPE
      IF (IBIT(KTYPE,JNUM).EQ.1) GOTO 20
      JNUM=JNUM-1
      KTYPE=KTYPE+2**JNUM
      ITYPE=KTYPE
      JNUM=JNUM+1
  20  RETURN
      END
      SUBROUTINE REBIT(ITYPE,JNUM)
      KTYPE=ITYPE
      IF (IBIT(KTYPE,JNUM).EQ.0) GOTO 20
      JNUM=JNUM-1
      KTYPE=KTYPE-2**JNUM
      ITYPE=KTYPE
      JNUM=JNUM+1
  20  RETURN
      END
      FUNCTION IBIT (ITYPE, JNUM)
      DIMENSION JBIT(8)
      KTYPE=ITYPE
      DO 20 I=7,0,-1
      IF ((KTYPE/2**I).EQ.1) GOTO 15
      JBIT(I+1)=0
      GO TO 20
  15  JBIT(I+1)=1
      KTYPE =KTYPE-(2**I)
  20  CONTINUE
      IBIT=JBIT(JNUM)
      RETURN
      END