input card: READ(88,*,IOSTAT=iostat) $ IPA,VERTEX(1), VERTEX(2), VERTEX(3), $ PXYZ(1), PXYZ(2), PXYZ(3), $ IB(1),IB(2),IB(3),IB(4),IB(5),IB(6),IB(7), $ IB(8),IB(9),IB(10),IB(11),IB(12) IF(iostat.NE.0)THEN PRINT*,FILEN,': *** Error reading file ***' RETURN ENDIF C decompress the bad va1 integers DO II=1,12 DO JL=1,3 DO IV=1,8 IVA1=(JL-1)*8+IV BADVA1(II,JL,IV)=JBIT( IB(II) , IVA1 ) ENDDO ENDDO ENDDO then look at GPUCOLLECTION