hdiff output

r30667/atom_label_swaps.f90 2016-07-06 15:38:15.077381289 +0100 r30666/atom_label_swaps.f90 2016-07-06 15:38:15.725390051 +0100
 69:            WRITE(MYUNIT,'(A,I1,A,I10,A,F20.10,A,I5,A,G12.5,A,F11.1)') & 69:            WRITE(MYUNIT,'(A,I1,A,I10,A,F20.10,A,I5,A,G12.5,A,F11.1)') &
 70:                 '[',NP,']Qu ',NQ(NP),' E=',POTEL,' steps=',ITER, & 70:                 '[',NP,']Qu ',NQ(NP),' E=',POTEL,' steps=',ITER, &
 71:                 ' RMS=',RMS,' t=',TIME 71:                 ' RMS=',RMS,' t=',TIME
 72:         ELSE 72:         ELSE
 73:            WRITE(MYUNIT,'(A,I10,A,F20.10,A,I5,A,G12.5,A,F11.1)') & 73:            WRITE(MYUNIT,'(A,I10,A,F20.10,A,I5,A,G12.5,A,F11.1)') &
 74:                 'Qu ',NQ(NP),' E=',POTEL,' steps=',ITER, & 74:                 'Qu ',NQ(NP),' E=',POTEL,' steps=',ITER, &
 75:                 ' RMS=',RMS,' t=',TIME 75:                 ' RMS=',RMS,' t=',TIME
 76:         ENDIF 76:         ENDIF
 77:      ENDIF 77:      ENDIF
 78:      ! 78:      !
 79:      R=EXP(MIN(2.0D0,(E0-POTEL)/TEMP)) 79:      R=DEXP((E0-POTEL)/TEMP)
 80:      IF(DPRAND() < R) THEN ! Accept 80:      IF(DPRAND() < R) THEN ! Accept
 81:         NACC = NACC + 1 81:         NACC = NACC + 1
 82:         IF(POTEL < EMIN - ECONV) THEN ! Store labels and coordinates 82:         IF(POTEL < EMIN - ECONV) THEN ! Store labels and coordinates
 83:            EMIN=POTEL 83:            EMIN=POTEL
 84:            XMIN(1:3*NATOMS)=COORDS(1:3*NATOMS, NP) 84:            XMIN(1:3*NATOMS)=COORDS(1:3*NATOMS, NP)
 85:            ATOMLISTS_MIN(:,:,:) = ATOMLISTS(:,:,:) 85:            ATOMLISTS_MIN(:,:,:) = ATOMLISTS(:,:,:)
 86:            INVATOMLISTS_MIN(:,:) = INVATOMLISTS(:,:) 86:            INVATOMLISTS_MIN(:,:) = INVATOMLISTS(:,:)
 87:            WRITE(MYUNIT,'(A,F20.10,A,I8,A,F10.8)') & 87:            WRITE(MYUNIT,'(A,F20.10,A,I8,A,F10.8)') &
 88:                 'uswaps> New EMIN= ',EMIN,' on swap', N,& 88:                 'uswaps> New EMIN= ',EMIN,' on swap', N,&
 89:                 ' at T=', TEMP 89:                 ' at T=', TEMP


r30667/msaveit.f 2016-07-06 15:38:15.401385653 +0100 r30666/msaveit.f 2016-07-06 15:38:16.053394485 +0100
 28:       DOUBLE PRECISION VTINT, VTMIN, VTMAX 28:       DOUBLE PRECISION VTINT, VTMIN, VTMAX
 29:       DOUBLE PRECISION IV(3,3), LTEMP(3,3), RM(3,3) 29:       DOUBLE PRECISION IV(3,3), LTEMP(3,3), RM(3,3)
 30:       DOUBLE PRECISION :: ATOMCOORDS(3,3) 30:       DOUBLE PRECISION :: ATOMCOORDS(3,3)
 31:       CHARACTER(LEN=80) ISTR 31:       CHARACTER(LEN=80) ISTR
 32:       CHARACTER(LEN=2) ZSTRING(NATOMS) 32:       CHARACTER(LEN=2) ZSTRING(NATOMS)
 33:       DOUBLE PRECISION ETheta,EPhi,EPsi,LPS(3*NATOMS)  ! sf344> Euler angles and coordinates for ellipsoids 33:       DOUBLE PRECISION ETheta,EPhi,EPsi,LPS(3*NATOMS)  ! sf344> Euler angles and coordinates for ellipsoids
 34:       INTEGER :: NCALLED=0 34:       INTEGER :: NCALLED=0
 35:       SAVE NCALLED 35:       SAVE NCALLED
 36:  36: 
 37:       NCALLED=NCALLED+1 37:       NCALLED=NCALLED+1
 38:       IF (MGUPTAT.AND.SPECLABELST) THEN 38:       IF (.NOT.AMBERT) THEN
 39:          DO J2=1,NSPECIES(0) !ds656> Should not exceed 10 
 40:            DO J3=1,ATOMLISTS(J2,1,0) 
 41:               J4=ATOMLISTS(J2,1,J3) ! Actual atom index 
 42:               ZSTRING(J4)=SPECLABELS(J2) 
 43:            ENDDO 
 44:          ENDDO 
 45:       ELSEIF (.NOT.AMBERT) THEN 
 46:          ZSTRING(1:NATOMS)='X1' 39:          ZSTRING(1:NATOMS)='X1'
 47:          VTMIN=1.0D100 40:          VTMIN=1.0D100
 48:          VTMAX=-1.0D100 41:          VTMAX=-1.0D100
 49:          IF (PYGPERIODICT.OR.PYBINARYT.OR.RIGID) NATOMS=NATOMS/2 ! temporary change! 42:          IF (PYGPERIODICT.OR.PYBINARYT.OR.RIGID) NATOMS=NATOMS/2 ! temporary change!
 50:          DO J1=1,NATOMS 43:          DO J1=1,NATOMS
 51:             IF (VT(J1).GT.VTMAX) VTMAX=VT(J1) 44:             IF (VT(J1).GT.VTMAX) VTMAX=VT(J1)
 52:             IF (VT(J1).LT.VTMIN) VTMIN=VT(J1) 45:             IF (VT(J1).LT.VTMIN) VTMIN=VT(J1)
 53:          ENDDO 46:          ENDDO
 54:          VTINT=VTMAX-VTMIN 47:          VTINT=VTMAX-VTMIN
 55:          IF (VTINT.GT.0.0D0) THEN 48:          IF (VTINT.GT.0.0D0) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0