PROGRAM fixbfpcd c Copies B's to monomers E,F,G,H and set occ = 1 CHARACTER*76 line REAL ab(2000), bb(2000) OPEN (unit=1, form='FORMATTED', status='OLD', file='secondbfpcd.pdb') OPEN (unit=2, form='FORMATTED', status='UNKNOWN', 1 file='secondbfpcdfix.pdb') DO WHILE (.TRUE.) READ (1, 1000, end=9) line 1000 FORMAT (a) IF (line(1:4) .EQ. 'ATOM') THEN IF (line(73:73) .EQ. 'A') THEN line(22:22) = 'A' ELSE IF (line(73:73) .EQ. 'B') THEN line(22:22) = 'A' line(73:73) = 'A' ELSE IF (line(73:73) .EQ. 'C') THEN line(22:22) = 'B' line(73:73) = 'B' ELSE IF (line(73:73) .EQ. 'D') THEN line(22:22) = 'B' line(73:73) = 'B' ELSE IF (line(73:73) .EQ. 'E') THEN line(22:22) = 'C' line(73:73) = 'C' ELSE IF (line(73:73) .EQ. 'F') THEN line(22:22) = 'C' line(73:73) = 'C' ELSE IF (line(73:73) .EQ. 'G') THEN line(22:22) = 'D' line(73:73) = 'D' ELSE IF (line(73:73) .EQ. 'H') THEN line(22:22) = 'D' line(73:73) = 'D' ELSE IF (line(73:73) .EQ. 'I') THEN line(22:22) = 'E' line(73:73) = 'E' ELSE IF (line(73:73) .EQ. 'J') THEN line(22:22) = 'E' line(73:73) = 'E' ELSE IF (line(73:73) .EQ. 'K') THEN line(22:22) = 'F' line(73:73) = 'F' ELSE IF (line(73:73) .EQ. 'L') THEN line(22:22) = 'F' line(73:73) = 'F' endIF line(57:57) = '1' endIF WRITE (2, 1000) line endDO 9 CLOSE (unit=1) CLOSE (unit=2) end