hdiff output

r33432/aaextraconstraint.f90 2017-11-02 15:30:16.253488171 +0000 r33431/aaextraconstraint.f90 2017-11-02 15:30:16.505491529 +0000
  1: PROGRAM AAXC  1: PROGRAM AAXC
  2: IMPLICIT NONE  2: IMPLICIT NONE
  3: CHARACTER(LEN=4) DSTRING,ASTRING  3: CHARACTER(LEN=4) DSTRING,ASTRING
  4: CHARACTER(LEN=3) RES  4: CHARACTER(LEN=3) RES
  5: INTEGER NAT, J1, HPOS, OPOS, CPOS, CPOSP  5: INTEGER NAT, J1, HPOS, OPOS, CPOS
  6:   6: 
  7: OPEN(UNIT=9,FILE='coords.inpcrd.pdb',STATUS='UNKNOWN')  7: OPEN(UNIT=9,FILE='coords.inpcrd.pdb',STATUS='UNKNOWN')
  8: OPEN(UNIT=7,FILE='constraintfile',STATUS='UNKNOWN')  8: OPEN(UNIT=7,FILE='constraintfile',STATUS='UNKNOWN')
  9: OPEN(UNIT=8,FILE='aabk',STATUS='UNKNOWN')  9: OPEN(UNIT=8,FILE='aabk',STATUS='UNKNOWN')
 10:  10: 
 11: DO  11: DO 
 12:  12: 
 13: READ(9,'(A4,1X,I6,1X,A4,1X,A3)',END=666) DSTRING,NAT,ASTRING,RES 13: READ(9,'(A4,1X,I6,1X,A4,1X,A3)',END=666) DSTRING,NAT,ASTRING,RES
 14: WRITE(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING,NAT,ASTRING,RES 14: WRITE(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING,NAT,ASTRING,RES
 15: IF (RES.EQ.'PHE') THEN 15: IF (RES.EQ.'PHE') THEN
115: ENDIF115: ENDIF
116: 116: 
117: ENDDO117: ENDDO
118: 118: 
119: 666 CONTINUE119: 666 CONTINUE
120: 120: 
121: REWIND(9)121: REWIND(9)
122: OPOS=0122: OPOS=0
123: HPOS=0123: HPOS=0
124: CPOS=0124: CPOS=0
125: CPOSP=0 
126: 125: 
127: DO126: DO
128: 127: 
129:    READ(9,'(A4,1X,I6,1X,A4,1X,A3)',END=667) DSTRING,NAT,ASTRING,RES128:    READ(9,'(A4,1X,I6,1X,A4,1X,A3)',END=667) DSTRING,NAT,ASTRING,RES
130: 129: 
131:    WRITE(*,'(A3)') ASTRING(2:4)130:    WRITE(*,'(A3)') ASTRING(2:4)
132: 131: 
133:    IF ((ASTRING(2:4).EQ.'N  ').OR. &132:    IF ((ASTRING(2:4).EQ.'N  ').OR. &
134:   &    (ASTRING(2:4).EQ.'H  ').OR. &133:   &    (ASTRING(2:4).EQ.'H  ').OR. &
135:   &    (ASTRING(2:4).EQ.'CA ').OR. &134:   &    (ASTRING(2:4).EQ.'CA ').OR. &
138:   &    (ASTRING(2:4).EQ.'C  ').OR. &137:   &    (ASTRING(2:4).EQ.'C  ').OR. &
139:   &    (ASTRING(2:4).EQ.'O  ').OR. &138:   &    (ASTRING(2:4).EQ.'O  ').OR. &
140:   &    (ASTRING(2:4).EQ.'HA2').OR. &139:   &    (ASTRING(2:4).EQ.'HA2').OR. &
141:   &    (ASTRING(2:4).EQ.'HA3')) THEN140:   &    (ASTRING(2:4).EQ.'HA3')) THEN
142: 141: 
143:       WRITE(8,'(I8)') NAT142:       WRITE(8,'(I8)') NAT
144:       IF (ASTRING(2:4).EQ.'O  ') OPOS=NAT143:       IF (ASTRING(2:4).EQ.'O  ') OPOS=NAT
145:       IF (ASTRING(2:4).EQ.'CA ') THEN144:       IF (ASTRING(2:4).EQ.'CA ') THEN
146:          CPOS=NAT145:          CPOS=NAT
147:          IF ((CPOS.GT.OPOS).AND.(OPOS.GT.0)) WRITE(7,'(2I8)') CPOS,OPOS146:          IF ((CPOS.GT.OPOS).AND.(OPOS.GT.0)) WRITE(7,'(2I8)') CPOS,OPOS
148: ! 
149: ! using this at the moment. 
150: ! 
151:          IF (CPOSP.GT.0) WRITE(7,'(2I8)') CPOSP,CPOS  ! Ca with Ca 
152:          CPOSP=CPOS 
153:       ENDIF147:       ENDIF
154:       IF (ASTRING(2:4).EQ.'H  ') THEN148:       IF (ASTRING(2:4).EQ.'H  ') THEN
155:          HPOS=NAT149:          HPOS=NAT
156:          IF ((HPOS.GT.OPOS).AND.(OPOS.GT.0)) WRITE(7,'(2I8)') OPOS,HPOS150:          IF ((HPOS.GT.OPOS).AND.(OPOS.GT.0)) WRITE(7,'(2I8)') OPOS,HPOS
157:          IF ((HPOS.GT.CPOS).AND.(CPOS.GT.0)) WRITE(7,'(2I8)') CPOS,HPOS151:          IF ((HPOS.GT.CPOS).AND.(CPOS.GT.0)) WRITE(7,'(2I8)') CPOS,HPOS
158:       ENDIF152:       ENDIF
159: 153: 
160:    ENDIF154:    ENDIF
161: 155: 
162: ENDDO156: ENDDO


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0