hdiff output

r29177/genrigid.f90 2015-11-17 23:33:45.773118611 +0000 r29176/genrigid.f90 2015-11-17 23:33:46.157123760 +0000
127:   DO J1 = 1, NRIGIDBODY127:   DO J1 = 1, NRIGIDBODY
128:      DO J2 = 1, NSITEPERBODY(J1)128:      DO J2 = 1, NSITEPERBODY(J1)
129:         DUMMY=RIGIDGROUPS(J2,J1)129:         DUMMY=RIGIDGROUPS(J2,J1)
130:         SITESRIGIDBODY(J2,:,J1) = INICOORDS(3*DUMMY-2:3*DUMMY)130:         SITESRIGIDBODY(J2,:,J1) = INICOORDS(3*DUMMY-2:3*DUMMY)
131:      ENDDO131:      ENDDO
132:   ENDDO132:   ENDDO
133: 133: 
134:   ! hk286 > determine number of degrees of freedom134:   ! hk286 > determine number of degrees of freedom
135:   DEGFREEDOMS = 0135:   DEGFREEDOMS = 0
136:   DO J1 = 1, NRIGIDBODY136:   DO J1 = 1, NRIGIDBODY
137:   ! Here, DEGFREEDOMS is just being used as a dummy to count the number of atoms which do belong to rigid bodies 
138:      DEGFREEDOMS = DEGFREEDOMS + NSITEPERBODY(J1)137:      DEGFREEDOMS = DEGFREEDOMS + NSITEPERBODY(J1)
139:   ENDDO138:   ENDDO
140:   DEGFREEDOMS = 6 * NRIGIDBODY + 3 * (NATOMS - DEGFREEDOMS)139:   DEGFREEDOMS = 6 * NRIGIDBODY + 3 * (NATOMS - DEGFREEDOMS)
141: 140: 
142: ! hk286 > Allocate further data 141: ! hk286 > Allocate further data 
143:   ALLOCATE (RIGIDSINGLES((DEGFREEDOMS/3 - 2 * NRIGIDBODY)))142:   ALLOCATE (RIGIDSINGLES((DEGFREEDOMS/3 - 2 * NRIGIDBODY)))
144:   ALLOCATE (RIGIDCOORDS(DEGFREEDOMS))143:   ALLOCATE (RIGIDCOORDS(DEGFREEDOMS))
145: 144: 
146:   DUMMY = 0145:   DUMMY = 0
147:   DO J1 = 1, NATOMS146:   DO J1 = 1, NATOMS
148:      SATOMT = .TRUE.147:      SATOMT = .TRUE.
149:      DO J2 = 1, NRIGIDBODY148:      DO J2 = 1, NRIGIDBODY
150:         DO J3 = 1, NSITEPERBODY(J2)149:         DO J3 = 1, NSITEPERBODY(J2)
151:            IF (J1 == RIGIDGROUPS(J3,J2)) SATOMT = .FALSE.150:            IF (J1 == RIGIDGROUPS(J3,J2)) SATOMT = .FALSE.
152:         ENDDO151:         ENDDO
153:      ENDDO152:      ENDDO
154:      IF (SATOMT) THEN153:      IF (SATOMT) THEN
155:         IF (DUMMY.EQ.DEGFREEDOMS/3-2*NRIGIDBODY) THEN 
156:            WRITE(*,*) "genrigid> Error. More free atoms than expected." 
157:            WRITE(*,*) "Likely problem with rbodyconfig." 
158:            STOP 
159:         ENDIF 
160:         DUMMY = DUMMY + 1154:         DUMMY = DUMMY + 1
161:         RIGIDSINGLES(DUMMY) = J1155:         RIGIDSINGLES(DUMMY) = J1
162:      ENDIF156:      ENDIF
163:   ENDDO157:   ENDDO
164:   DO J1 = 1, NRIGIDBODY158:   DO J1 = 1, NRIGIDBODY
165:      XMASS = 0.0D0159:      XMASS = 0.0D0
166:      YMASS = 0.0D0160:      YMASS = 0.0D0
167:      ZMASS = 0.0D0161:      ZMASS = 0.0D0
168:      MASS = 0.0d0162:      MASS = 0.0d0
169:      DO J2 = 1, NSITEPERBODY(J1)163:      DO J2 = 1, NSITEPERBODY(J1)
694: !                WRITE(*,*) "Warning, coordinate transform may have failed."688: !                WRITE(*,*) "Warning, coordinate transform may have failed."
695: !                WRITE(*,*) "Original coordinate:", SAVEXCOORDS(J1)689: !                WRITE(*,*) "Original coordinate:", SAVEXCOORDS(J1)
696: !                WRITE(*,*) "New coordinate:", TEMPXCOORDS(J1)690: !                WRITE(*,*) "New coordinate:", TEMPXCOORDS(J1)
697: !                WRITE(*,*) "Difference:", TEMPXCOORDS(J1)-SAVEXCOORDS(J1)691: !                WRITE(*,*) "Difference:", TEMPXCOORDS(J1)-SAVEXCOORDS(J1)
698: !                FAILCOUNT = FAILCOUNT+1692: !                FAILCOUNT = FAILCOUNT+1
699: !            ENDIF693: !            ENDIF
700: !        ENDIF694: !        ENDIF
701:     ENDDO695:     ENDDO
702:     IF (FAILCOUNT > 0) THEN696:     IF (FAILCOUNT > 0) THEN
703:         WRITE(*,*) "Transformation failed ", FAILCOUNT, "times. Aborting."697:         WRITE(*,*) "Transformation failed ", FAILCOUNT, "times. Aborting."
704:         STOP 99698:         STOP
705:     ENDIF699:     ENDIF
706:   ENDIF700:   ENDIF
707: 701: 
708: END SUBROUTINE TRANSFORMCTORIGID702: END SUBROUTINE TRANSFORMCTORIGID
709: 703: 
710: !-----------------------------------------------------------704: !-----------------------------------------------------------
711: 705: 
712: SUBROUTINE TRANSFORMGRAD (G, XR, GR)706: SUBROUTINE TRANSFORMGRAD (G, XR, GR)
713:   707:   
714:   USE COMMONS, ONLY: NATOMS708:   USE COMMONS, ONLY: NATOMS


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0