hdiff output

r29839/disconnectionDPS.f90 2016-01-25 18:30:05.704684105 +0000 r29838/disconnectionDPS.f90 2016-01-25 18:30:05.896686676 +0000
738: 738: 
739:    WRITE (6, '(I6, 3A)') &739:    WRITE (6, '(I6, 3A)') &
740:       & N_TS, ' non-degenerate paths read in from ', TRIM(file_ts), '.'740:       & N_TS, ' non-degenerate paths read in from ', TRIM(file_ts), '.'
741:    WRITE (6, '(A, 2(F18.10))') 'Highest and lowest transition states: ', &741:    WRITE (6, '(A, 2(F18.10))') 'Highest and lowest transition states: ', &
742:       & B_HIGH, B_LOW742:       & B_HIGH, B_LOW
743: 743: 
744:    PRINT *,'n_min,n_ts=',n_min,n_ts744:    PRINT *,'n_min,n_ts=',n_min,n_ts
745:    ALLOCATE (NCONN(N_MIN),PLUS(N_TS),MINUS(N_TS),NDISTA(N_MIN),TSEN(N_TS))745:    ALLOCATE (NCONN(N_MIN),PLUS(N_TS),MINUS(N_TS),NDISTA(N_MIN),TSEN(N_TS))
746:    DO I=1,N_TS746:    DO I=1,N_TS
747:       IF (TS(I)%E.LE.M(TS(I)%MIN1)) THEN747:       IF (TS(I)%E.LE.M(TS(I)%MIN1)) THEN
748:          PRINT '(A,G20.10,A)','WARNING *** energy of ts (= ',TS(I)%E,') is lower than min1:'748:          PRINT '(A,I6,A)','WARNING *** energy of ts ',I,' is lower than min1.'
749:          IF (TS(I)%E-M(TS(I)%MIN1).GE.TSLOWBARTHRESH) THEN749:          IF (TS(I)%E-M(TS(I)%MIN1).GE.TSLOWBARTHRESH) THEN
750:             TS(I)%E=M(TS(I)%MIN1)+DELTA_E/2.0D0750:             TS(I)%E=M(TS(I)%MIN1)+DELTA_E/2.0D0
751:             PRINT '(A)','resetting.'751:             PRINT '(A)','Resetting'
752:          ELSE 752:          ELSE 
753:             PRINT '(A)','excluding this TS; not resetting.'753:             PRINT '(A)','Excluding this TS; not resetting'
754:          ENDIF754:          ENDIF
755:       ENDIF755:       ENDIF
756:       IF (TS(I)%E.LE.M(TS(I)%MIN2)) THEN756:       IF (TS(I)%E.LE.M(TS(I)%MIN2)) THEN
757:          PRINT '(A,G20.10,A)','WARNING *** energy of ts (= ',TS(I)%E,') is lower than min2:'757:          PRINT '(A,I6,A)','WARNING *** energy of ts ',I,' is lower than min2.'
758:          IF (TS(I)%E-M(TS(I)%MIN2).GE.TSLOWBARTHRESH) THEN758:          IF (TS(I)%E-M(TS(I)%MIN2).GE.TSLOWBARTHRESH) THEN
759:             TS(I)%E=M(TS(I)%MIN2)+DELTA_E/2.0D0759:             TS(I)%E=M(TS(I)%MIN2)+DELTA_E/2.0D0
760:             PRINT '(A)','resetting.'760:             PRINT '(A)','Resetting'
761:          ELSE 761:          ELSE 
762:             PRINT '(A)','excluding this TS; not resetting.'762:             PRINT '(A)','Excluding this TS; not resetting'
763:          ENDIF763:          ENDIF
764:       ENDIF764:       ENDIF
765:    ENDDO765:    ENDDO
766:    DO I=1,N_TS766:    DO I=1,N_TS
767:       TSEN(I)=TS(I)%E767:       TSEN(I)=TS(I)%E
768:       PLUS(I)=TS(I)%MIN1768:       PLUS(I)=TS(I)%MIN1
769:       MINUS(I)=TS(I)%MIN2769:       MINUS(I)=TS(I)%MIN2
770:    ENDDO770:    ENDDO
771:    CALL GETNCONN(N_MIN,N_TS,NCONN,PLUS,MINUS,NCONNMIN,NCONNMAX,.FALSE.,TSEN)771:    CALL GETNCONN(N_MIN,N_TS,NCONN,PLUS,MINUS,NCONNMIN,NCONNMAX,.FALSE.,TSEN)
772:    DEADTS(1:N_TS)=.FALSE.772:    DEADTS(1:N_TS)=.FALSE.


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0