PROGRAM editpdb c Changes ILE CD1 to ILE CD and adjust PCD sequence for submission CHARACTER*66 line WRITE (6, 1000) 'Enter input file:' 1000 FORMAT (a) READ (5, 1000) line OPEN (unit=1, file=line, form='FORMATTED', status='OLD', 1 readonly) WRITE (6, 1000) 'Enter output file:' READ (5, 1000) line OPEN (unit=2, file=line, form='FORMATTED', status='NEW', 1 carriagecontrol='LIST') DO WHILE (.TRUE.) READ (1, 1000, end=9) line READ (line(23:26), 1001) nres 1001 FORMAT (i4) IF (nres .EQ. 173) line(18:20) = 'LYS' IF (nres .EQ. 303) line(18:20) = 'GLN' IF (nres .EQ. 532) line(18:20) = 'LYS' IF (line(18:20) .EQ. 'ILE' .AND. line(14:16) .EQ. 'CD1') 1 line(14:16) = 'CD ' WRITE (2, 1000) line endDO 9 CLOSE (unit=1) CLOSE (unit=2) end