hdiff output

r33518/congrad.f90 2017-11-25 23:30:08.376290984 +0000 r33517/congrad.f90 2017-11-25 23:30:08.596293902 +0000
451:          ENDDO451:          ENDDO
452: !     ENDIF452: !     ENDIF
453:    ENDDO453:    ENDDO
454: ENDIF454: ENDIF
455: SEPARATION=SUM(DVEC(1:INTIMAGE+1))455: SEPARATION=SUM(DVEC(1:INTIMAGE+1))
456: DEVIATION(1:INTIMAGE+1)=ABS(100*((INTIMAGE+1)*DVEC(1:INTIMAGE+1)/SEPARATION-1.0D0))456: DEVIATION(1:INTIMAGE+1)=ABS(100*((INTIMAGE+1)*DVEC(1:INTIMAGE+1)/SEPARATION-1.0D0))
457: QCIAVDEV=SUM(DEVIATION)/(INTIMAGE+1)457: QCIAVDEV=SUM(DEVIATION)/(INTIMAGE+1)
458: WRITE(*,'(A,I6,A,I6,A,2I6)') ' congrad> Highest spring  contribution for any image in image ',IMAX458: WRITE(*,'(A,I6,A,I6,A,2I6)') ' congrad> Highest spring  contribution for any image in image ',IMAX
459: IF (DEBUG) THEN459: IF (DEBUG) THEN
460:    WRITE(*, '(A,3G20.10)') ' congrad> ECON,EREP,ESPRING=',ECON,EREP,ESPRING460:    WRITE(*, '(A,3G20.10)') ' congrad> ECON,EREP,ESPRING=',ECON,EREP,ESPRING
461: !  WRITE(*,'(A)') '   edge         gap                deviation      '461:    WRITE(*,'(A)') '   edge         gap                deviation      '
462: !  WRITE(*,'(I6,3X,G20.10,G20.10)') (J1,DVEC(J1),DEVIATION(J1),J1=1,INTIMAGE+1)462:    WRITE(*,'(I6,3X,G20.10,G20.10)') (J1,DVEC(J1),DEVIATION(J1),J1=1,INTIMAGE+1)
463:    WRITE(*, '(A,2G20.10)') ' congrad> mean gap and mean deviation=',SEPARATION/(INTIMAGE+1),QCIAVDEV463:    WRITE(*, '(A,G20.10)') ' congrad> mean deviation=',QCIAVDEV
464: ENDIF464: ENDIF
465: !465: !
466: ! Set gradients on frozen atoms to zero.466: ! Set gradients on frozen atoms to zero.
467: !467: !
468: IF (FREEZE) THEN468: IF (FREEZE) THEN
469:    DO J1=2,INTIMAGE+1  469:    DO J1=2,INTIMAGE+1  
470:       DO J2=1,NATOMS470:       DO J2=1,NATOMS
471:          IF (FROZEN(J2)) THEN471:          IF (FROZEN(J2)) THEN
472:             GGG((3*NATOMS)*(J1-1)+3*(J2-1)+1)=0.0D0472:             GGG((3*NATOMS)*(J1-1)+3*(J2-1)+1)=0.0D0
473:             GGG((3*NATOMS)*(J1-1)+3*(J2-1)+2)=0.0D0473:             GGG((3*NATOMS)*(J1-1)+3*(J2-1)+2)=0.0D0


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0