PROGRAM listpack c Looks through packing.log and lists protein residues that c make contacts CHARACTER*80 line CHARACTER*1 c LOGICAL*1 f(12,600,14)/100800*.FALSE./ c First read and store information OPEN (unit=1, form='FORMATTED', status='OLD', file='packing.log', 1 readonly) DO WHILE (line(1:5) .NE. 'Enter') READ (1, 1000) line 1000 FORMAT (a) endDO READ (1, 1000) DO WHILE (.TRUE.) 1 READ (1, 1001, err=9) n, i, c 1001 FORMAT (i3, 5x, i4, a1, 20x, i4, a1) j = ICHAR(c)-64 IF (i .LE. 600) THEN f(j,i,n) = .TRUE. endIF endDO 9 CLOSE (unit=1) c Now write results DO n = 1,14 nc = 1 DO j = 1,12 c = CHAR(j+64) DO i = 1,600 IF (f(j,i,n)) THEN WRITE (6, 1001) n, i, c, nc nc = nc + 1 endIF endDO endDO endDO end