hdiff output

r33428/aaextraconstraint.f90 2017-10-30 15:30:11.531176602 +0000 r33427/aaextraconstraint.f90 2017-10-30 15:30:11.751179523 +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  5: INTEGER NAT, J1
  6:   6: 
  7: OPEN(UNIT=9,FILE='coords.inpcrd.pdb',STATUS='UNKNOWN')  7: OPEN(UNIT=7,FILE='aaxc',STATUS='UNKNOWN')
  8: OPEN(UNIT=7,FILE='constraintfile',STATUS='UNKNOWN') 
  9: OPEN(UNIT=8,FILE='aabk',STATUS='UNKNOWN') 
 10:   8: 
 11: DO   9: DO 
 12:  10: 
 13: READ(9,'(A4,1X,I6,1X,A4,1X,A3)',END=666) DSTRING,NAT,ASTRING,RES 11: READ(*,'(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 12: WRITE(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING,NAT,ASTRING,RES
 15: IF (RES.EQ.'PHE') THEN 13: IF (RES.EQ.'PHE') THEN
 16:  14: 
 17: ! HD2 CE1 15: ! HD2 CE1
 18: ! HE2 CD1 16: ! HE2 CD1
 19: ! HZ  CG 17: ! HZ  CG
 20: ! HE1 CD2 18: ! HE1 CD2
 21: ! HD1 CE2 19: ! HD1 CE2
 22: ! CB CZ 20: ! CB CZ
 23:  21: 
 24:    WRITE(7,'(2I8)') NAT+17,NAT+10 22:    WRITE(7,'(2I8)') NAT+17,NAT+10
 25:    WRITE(7,'(2I8)') NAT+15,NAT+8 23:    WRITE(7,'(2I8)') NAT+15,NAT+8
 26:    WRITE(7,'(2I8)') NAT+13,NAT+7 24:    WRITE(7,'(2I8)') NAT+13,NAT+7
 27:    WRITE(7,'(2I8)') NAT+11,NAT+16 25:    WRITE(7,'(2I8)') NAT+11,NAT+16
 28:    WRITE(7,'(2I8)') NAT+9,NAT+14 26:    WRITE(7,'(2I8)') NAT+9,NAT+14
 29:    WRITE(7,'(2I8)') NAT+4,NAT+12 27:    WRITE(7,'(2I8)') NAT+4,NAT+12
 30:  28: 
 31:    DO J1=1,19 29:    DO J1=1,19
 32:       READ(9,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING 30:       READ(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING
 33:    ENDDO 31:    ENDDO
 34:  32: 
 35: ELSEIF (RES.EQ.'ARG') THEN 33: ELSEIF (RES.EQ.'ARG') THEN
 36:  34: 
 37: ! HH12 HH22 35: ! HH12 HH22
 38: ! HH12 HH21 36: ! HH12 HH21
 39: ! HH11 HH22 37: ! HH11 HH22
 40: ! HH11 HH21 38: ! HH11 HH21
 41: ! HH11 NE 39: ! HH11 NE
 42: ! HH12 NE 40: ! HH12 NE
 46:    WRITE(7,'(2I8)') NAT+18,NAT+21 44:    WRITE(7,'(2I8)') NAT+18,NAT+21
 47:    WRITE(7,'(2I8)') NAT+18,NAT+20 45:    WRITE(7,'(2I8)') NAT+18,NAT+20
 48:    WRITE(7,'(2I8)') NAT+17,NAT+21 46:    WRITE(7,'(2I8)') NAT+17,NAT+21
 49:    WRITE(7,'(2I8)') NAT+17,NAT+20 47:    WRITE(7,'(2I8)') NAT+17,NAT+20
 50:    WRITE(7,'(2I8)') NAT+17,NAT+13 48:    WRITE(7,'(2I8)') NAT+17,NAT+13
 51:    WRITE(7,'(2I8)') NAT+18,NAT+13 49:    WRITE(7,'(2I8)') NAT+18,NAT+13
 52:    WRITE(7,'(2I8)') NAT+20,NAT+13 50:    WRITE(7,'(2I8)') NAT+20,NAT+13
 53:    WRITE(7,'(2I8)') NAT+21,NAT+13 51:    WRITE(7,'(2I8)') NAT+21,NAT+13
 54:  52: 
 55:    DO J1=1,23 53:    DO J1=1,23
 56:       READ(9,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING 54:       READ(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING
 57:    ENDDO 55:    ENDDO
 58:  56: 
 59: ELSEIF (RES.EQ.'TRP') THEN 57: ELSEIF (RES.EQ.'TRP') THEN
 60:  58: 
 61: ! HD1 CD2 59: ! HD1 CD2
 62: ! HE1 CG 60: ! HE1 CG
 63: ! HZ2 CE3 61: ! HZ2 CE3
 64: ! HH2 CD2 62: ! HH2 CD2
 65: ! HZ3 CE2 63: ! HZ3 CE2
 66: ! HE3 CZ2 64: ! HE3 CZ2
 67:  65: 
 68:    WRITE(7,'(2I8)') NAT+9,NAT+21 66:    WRITE(7,'(2I8)') NAT+9,NAT+21
 69:    WRITE(7,'(2I8)') NAT+11,NAT+7 67:    WRITE(7,'(2I8)') NAT+11,NAT+7
 70:    WRITE(7,'(2I8)') NAT+14,NAT+19 68:    WRITE(7,'(2I8)') NAT+14,NAT+19
 71:    WRITE(7,'(2I8)') NAT+16,NAT+21 69:    WRITE(7,'(2I8)') NAT+16,NAT+21
 72:    WRITE(7,'(2I8)') NAT+18,NAT+12 70:    WRITE(7,'(2I8)') NAT+18,NAT+12
 73:    WRITE(7,'(2I8)') NAT+20,NAT+13 71:    WRITE(7,'(2I8)') NAT+20,NAT+13
 74:  72: 
 75:    DO J1=1,23 73:    DO J1=1,23
 76:       READ(9,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING 74:       READ(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING
 77:    ENDDO 75:    ENDDO
 78:  76: 
 79: ELSEIF (RES.EQ.'TYR') THEN 77: ELSEIF (RES.EQ.'TYR') THEN
 80:  78: 
 81: ! HD1 CE2 79: ! HD1 CE2
 82: ! HE1 CD2 80: ! HE1 CD2
 83: ! OH CG 81: ! OH CG
 84: ! HE2 CD1 82: ! HE2 CD1
 85: ! HD2 CE1 83: ! HD2 CE1
 86: ! CB CZ 84: ! CB CZ
 87:  85: 
 88:    WRITE(7,'(2I8)') NAT+9,NAT+15 86:    WRITE(7,'(2I8)') NAT+9,NAT+15
 89:    WRITE(7,'(2I8)') NAT+11,NAT+17 87:    WRITE(7,'(2I8)') NAT+11,NAT+17
 90:    WRITE(7,'(2I8)') NAT+13,NAT+7 88:    WRITE(7,'(2I8)') NAT+13,NAT+7
 91:    WRITE(7,'(2I8)') NAT+16,NAT+8 89:    WRITE(7,'(2I8)') NAT+16,NAT+8
 92:    WRITE(7,'(2I8)') NAT+18,NAT+10 90:    WRITE(7,'(2I8)') NAT+18,NAT+10
 93:    WRITE(7,'(2I8)') NAT+4,NAT+12 91:    WRITE(7,'(2I8)') NAT+4,NAT+12
 94:  92: 
 95:    DO J1=1,20 93:    DO J1=1,20
 96:       READ(9,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING 94:       READ(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING
 97:    ENDDO 95:    ENDDO
 98:  96: 
 99: ELSEIF (RES.EQ.'HIP') THEN 97: ELSEIF (RES.EQ.'HIP') THEN
100:  98: 
101: ! HD1 NE2 99: ! HD1 NE2
102: ! HE1 CD2100: ! HE1 CD2
103: ! HE2 CG101: ! HE2 CG
104: ! HD2 ND1102: ! HD2 ND1
105: 103: 
106:    WRITE(7,'(2I8)') NAT+9,NAT+12104:    WRITE(7,'(2I8)') NAT+9,NAT+12
107:    WRITE(7,'(2I8)') NAT+11,NAT+14105:    WRITE(7,'(2I8)') NAT+11,NAT+14
108:    WRITE(7,'(2I8)') NAT+13,NAT+7106:    WRITE(7,'(2I8)') NAT+13,NAT+7
109:    WRITE(7,'(2I8)') NAT+15,NAT+8107:    WRITE(7,'(2I8)') NAT+15,NAT+8
110: 108: 
111:    DO J1=1,17109:    DO J1=1,17
112:       READ(9,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING110:       READ(*,'(A4,1X,I6,1X,A4,1X,A3)') DSTRING
113:    ENDDO111:    ENDDO
114: 112: 
115: ENDIF113: ENDIF
116: 114: 
117: ENDDO115: ENDDO
118: 116: 
119: 666 CONTINUE117: 666 CONTINUE
120: 118: 
121: REWIND(9) 
122: OPOS=0 
123: HPOS=0 
124:  
125: DO 
126:  
127:    READ(9,'(A4,1X,I6,1X,A4,1X,A3)',END=667) DSTRING,NAT,ASTRING,RES 
128:  
129:    WRITE(*,'(A3)') ASTRING(2:4) 
130:  
131:    IF ((ASTRING(2:4).EQ.'N  ').OR. & 
132:   &    (ASTRING(2:4).EQ.'H  ').OR. & 
133:   &    (ASTRING(2:4).EQ.'CA ').OR. & 
134:   &    (ASTRING(2:4).EQ.'HA ').OR. & 
135:   &    (ASTRING(2:4).EQ.'CB ').OR. & 
136:   &    (ASTRING(2:4).EQ.'C  ').OR. & 
137:   &    (ASTRING(2:4).EQ.'O  ').OR. & 
138:   &    (ASTRING(2:4).EQ.'HA2').OR. & 
139:   &    (ASTRING(2:4).EQ.'HA3')) THEN 
140:  
141:       WRITE(8,'(I8)') NAT 
142:       IF (ASTRING(2:4).EQ.'O  ') OPOS=NAT 
143:       IF (ASTRING(2:4).EQ.'H  ') THEN 
144:          HPOS=NAT 
145:          WRITE(7,'(2I8)') OPOS,HPOS 
146:       ENDIF 
147:  
148:    ENDIF 
149:  
150: ENDDO 
151:  
152: 667 CONTINUE 
153:  
154: CLOSE(7) 
155: CLOSE(8) 
156: CLOSE(9) 
157:  
158:  
159: END PROGRAM AAXC119: END PROGRAM AAXC


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0