hdiff output

r33408/OPTIM.F 2017-10-23 10:30:12.821198420 +0100 r33407/OPTIM.F 2017-10-23 10:30:13.049201436 +0100
220:       IF ((ZSYMSAVE.EQ.'RM').AND.(PARAM1.NE.0.0D0)) BULKT=.TRUE.220:       IF ((ZSYMSAVE.EQ.'RM').AND.(PARAM1.NE.0.0D0)) BULKT=.TRUE.
221: C     LDUMMY=(NATOMS.EQ.64).AND.(ZSYMSAVE(1:1).EQ.'W') 221: C     LDUMMY=(NATOMS.EQ.64).AND.(ZSYMSAVE(1:1).EQ.'W') 
222:       BULKT=(BULKT.OR.222:       BULKT=(BULKT.OR.
223:      1      (ZSYMSAVE.EQ.'ME').OR.(ZSYMSAVE.EQ.'P6').OR.223:      1      (ZSYMSAVE.EQ.'ME').OR.(ZSYMSAVE.EQ.'P6').OR.
224:      2      (ZSYMSAVE.EQ.'SC').OR.(ZSYMSAVE.EQ.'MS').OR.224:      2      (ZSYMSAVE.EQ.'SC').OR.(ZSYMSAVE.EQ.'MS').OR.
225:      3      (ZSYMSAVE.EQ.'MP').OR.(ZSYMSAVE.EQ.'JM').OR.225:      3      (ZSYMSAVE.EQ.'MP').OR.(ZSYMSAVE.EQ.'JM').OR.
226:      3      (ZSYMSAVE.EQ.'DS').OR.226:      3      (ZSYMSAVE.EQ.'DS').OR.
227:      4      (ZSYMSAVE.EQ.'LP').OR.(ZSYMSAVE.EQ.'LS').OR.(ZSYMSAVE.EQ.'LC').OR.227:      4      (ZSYMSAVE.EQ.'LP').OR.(ZSYMSAVE.EQ.'LS').OR.(ZSYMSAVE.EQ.'LC').OR.
228:      5      (ZSYMSAVE.EQ.'LK'))228:      5      (ZSYMSAVE.EQ.'LK'))
229: 229: 
230:       IF (BULKT.AND.NEWCONNECTT.AND.(.NOT.QUIPT)) THEN230:       IF (BULKT.AND.NEWCONNECTT) THEN
231:          CALL INITIALBOX(Q)231:          CALL INITIALBOX(Q)
232:          CALL INITIALBOX(FIN)232:          CALL INITIALBOX(FIN)
233:       ENDIF233:       ENDIF
234: 234: 
235: C     NOSHIFT=(FIELDT.OR.BFGSMINT.OR.BSMIN.OR.RKMIN.OR.NOHESS.OR.BFGSSTEP)235: C     NOSHIFT=(FIELDT.OR.BFGSMINT.OR.BSMIN.OR.RKMIN.OR.NOHESS.OR.BFGSSTEP)
236: C     IF (INR.GE.0) NOSHIFT=.FALSE. ! changed INR default value in keywords to -1236: C     IF (INR.GE.0) NOSHIFT=.FALSE. ! changed INR default value in keywords to -1
237: C                                     This isn;t good enough - if we are doing a path and237: C                                     This isn;t good enough - if we are doing a path and
238: C                                     EFOL is called once to get the eigenvalues and eigenvectors238: C                                     EFOL is called once to get the eigenvalues and eigenvectors
239: C                                     then we need to SHIFT.239: C                                     then we need to SHIFT.
240: C  Shifting only needs to be turned off if we are doing VARIABLES or BFGSTS without NOIT;240: C  Shifting only needs to be turned off if we are doing VARIABLES or BFGSTS without NOIT;
288: C     DO J1=1,NATOMS288: C     DO J1=1,NATOMS
289: C        IF (IATNUM(J1).NE.0) J2=J2+1289: C        IF (IATNUM(J1).NE.0) J2=J2+1
290: C     ENDDO290: C     ENDDO
291: C     NREAL=J2291: C     NREAL=J2
292: C292: C
293: C     CALL GMETRY(0,VEC,Q)293: C     CALL GMETRY(0,VEC,Q)
294: C294: C
295: C  Don't call symmetry if we're doing Fenske-Hall295: C  Don't call symmetry if we're doing Fenske-Hall
296: C296: C
297:       IF ((ZSYMSAVE.NE.'FH').AND.(.NOT.VARIABLES).AND.(.NOT.AMBERT).AND.(.NOT.CHRMMT).AND.297:       IF ((ZSYMSAVE.NE.'FH').AND.(.NOT.VARIABLES).AND.(.NOT.AMBERT).AND.(.NOT.CHRMMT).AND.
298:      1 (.NOT.NABT).AND.(.NOT.UNRST).AND.(.NOT.RINGPOLYMERT).AND.(.NOT.AMBER12T).AND.(.NOT.BULKT)) THEN298:      1 (.NOT.NABT).AND.(.NOT.UNRST).AND.(.NOT.RINGPOLYMERT).AND.(.NOT.AMBER12T)) THEN
299: C299: C
300: C  For W1/W2/W3/W4 potentials Q contains the centre of mass coordinates300: C  For W1/W2/W3/W4 potentials Q contains the centre of mass coordinates
301: C  followed by the Euler angles and the number of molecules is NATOMS/2301: C  followed by the Euler angles and the number of molecules is NATOMS/2
302: C302: C
303:          IF (ZSYMSAVE(1:1).EQ.'W') THEN303:          IF (ZSYMSAVE(1:1).EQ.'W') THEN
304:             ALLOCATE(QW(9*(NATOMS/2)))304:             ALLOCATE(QW(9*(NATOMS/2)))
305: C           DO J2=1,NOPT+IADD305: C           DO J2=1,NOPT+IADD
306: C              QSAVE(J2)=Q(J2)306: C              QSAVE(J2)=Q(J2)
307: C           ENDDO307: C           ENDDO
308: C           DO J2=1,NATOMS  !  WCOMMENT308: C           DO J2=1,NATOMS  !  WCOMMENT


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0