hdiff output

r28691/disconnectionDPS.f90 2017-01-21 10:41:41.415859572 +0000 r28690/disconnectionDPS.f90 2017-01-21 10:41:41.655869144 +0000
689: !     READ (20, *) j, m(i)689: !     READ (20, *) j, m(i)
690: !  Standard PATHSAMPLE.2.0 format690: !  Standard PATHSAMPLE.2.0 format
691:       READ (20, *) M(I)691:       READ (20, *) M(I)
692:       IF (M(I) < E_GMIN) THEN692:       IF (M(I) < E_GMIN) THEN
693:          E_GMIN = M(I)693:          E_GMIN = M(I)
694:          GMIN = I694:          GMIN = I
695:       END IF695:       END IF
696:    END DO696:    END DO
697:    CLOSE (20)697:    CLOSE (20)
698:    WRITE (6, '(I6, 3A)') n_min, ' minima read in from ', TRIM(file_min), '.'698:    WRITE (6, '(I6, 3A)') n_min, ' minima read in from ', TRIM(file_min), '.'
699:    WRITE (6, '(A,I10,F18.10,A)') ' Global minimum at index, energy ',GMIN,E_GMIN,'.'699:    WRITE (6, '(A,2(F18.10),A)') ' Global minimum at index, energy ',GMIN,E_GMIN,'.'
700:    IF ((LOWEST <= 0).OR.(LOWEST > N_MIN)) LOWEST = N_MIN700:    IF ((LOWEST <= 0).OR.(LOWEST > N_MIN)) LOWEST = N_MIN
701: 701: 
702: 702: 
703: !  Procure transition state info.703: !  Procure transition state info.
704:    CALL COUNT_TS(FILE_TS, N_TS, B_HIGH, B_LOW)704:    CALL COUNT_TS(FILE_TS, N_TS, B_HIGH, B_LOW)
705:    ALLOCATE (TS(N_TS),DEADTS(N_TS))705:    ALLOCATE (TS(N_TS),DEADTS(N_TS))
706:    I = 1706:    I = 1
707:    OPEN (UNIT=20, FILE=TRIM(FILE_TS), STATUS='OLD', IOSTAT=err)707:    OPEN (UNIT=20, FILE=TRIM(FILE_TS), STATUS='OLD', IOSTAT=err)
708:    IF (ERR /= 0) THEN708:    IF (ERR /= 0) THEN
709:       WRITE (6, '(/,2A,/)') 'ERROR: Could not open ', TRIM(file_ts)709:       WRITE (6, '(/,2A,/)') 'ERROR: Could not open ', TRIM(file_ts)
836:    NLEFTTS=0836:    NLEFTTS=0
837:    DO J1=1,N_TS837:    DO J1=1,N_TS
838:       IF ((NCONN(PLUS(J1)).LE.NCONNMIN).OR.(NCONN(MINUS(J1)).LE.NCONNMIN)) DEADTS(J1)=.TRUE.838:       IF ((NCONN(PLUS(J1)).LE.NCONNMIN).OR.(NCONN(MINUS(J1)).LE.NCONNMIN)) DEADTS(J1)=.TRUE.
839:       IF (TS(J1)%E.GT.TSTHRESH) DEADTS(J1)=.TRUE.839:       IF (TS(J1)%E.GT.TSTHRESH) DEADTS(J1)=.TRUE.
840:       IF ((TS(J1)%E-M(TS(J1)%MIN1).GT.TSBARTHRESH).AND.(TS(J1)%E-M(TS(J1)%MIN2).GT.TSBARTHRESH)) DEADTS(J1)=.TRUE.840:       IF ((TS(J1)%E-M(TS(J1)%MIN1).GT.TSBARTHRESH).AND.(TS(J1)%E-M(TS(J1)%MIN2).GT.TSBARTHRESH)) DEADTS(J1)=.TRUE.
841:       IF (.NOT.DEADTS(J1)) NLEFTTS=NLEFTTS+1841:       IF (.NOT.DEADTS(J1)) NLEFTTS=NLEFTTS+1
842:    ENDDO842:    ENDDO
843:    PRINT '(A,2I8)','Number of remaining minima and transition states=',NLEFTMIN,NLEFTTS843:    PRINT '(A,2I8)','Number of remaining minima and transition states=',NLEFTMIN,NLEFTTS
844: 844: 
845:    IF (GMIN_CONN /= GMIN) THEN845:    IF (GMIN_CONN /= GMIN) THEN
846:       WRITE (6, '(A,I10,F18.10,A)') ' Lowest minimum satisfying connectivity criteria at index, energy ',GMIN_CONN,E_GMIN_CONN,'.'846:       WRITE (6, '(A,2(F18.10),A)') ' Lowest minimum satisfying connectivity criteria at index, energy ',GMIN_CONN,E_GMIN_CONN,'.'
847:       GMIN = GMIN_CONN847:       GMIN = GMIN_CONN
848:       E_GMIN = E_GMIN_CONN848:       E_GMIN = E_GMIN_CONN
849:       IF (ZEROGM) THEN849:       IF (ZEROGM) THEN
850:          WRITE (6, '(A)') ' Setting the energy of the lowest minimum on the tree to zero.'850:          WRITE (6, '(A)') ' Setting the energy of the lowest minimum on the tree to zero.'
851:          DO I=1,N_MIN851:          DO I=1,N_MIN
852:             M(I)=M(I)-E_GMIN_CONN852:             M(I)=M(I)-E_GMIN_CONN
853:          ENDDO853:          ENDDO
854:          DO I=1,N_TS854:          DO I=1,N_TS
855:             TS(I)%E=TS(I)%E-E_GMIN_CONN855:             TS(I)%E=TS(I)%E-E_GMIN_CONN
856:          ENDDO856:          ENDDO


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0