hdiff output

r29831/keywords.f 2016-07-06 15:38:52.577888203 +0100 r29830/keywords.f 2016-07-06 15:38:53.281897722 +0100
1318: C1318: C
1319:       ELSE IF (WORD.EQ.'ITOL') THEN1319:       ELSE IF (WORD.EQ.'ITOL') THEN
1320:          CALL READF(IDIFFTOL)1320:          CALL READF(IDIFFTOL)
1321: C1321: C
1322: C  Number of jobs to run per node for a distributed memory architecture.1322: C  Number of jobs to run per node for a distributed memory architecture.
1323: C  JPN should always be one for all recent cluster implementations,1323: C  JPN should always be one for all recent cluster implementations,
1324: C  unless you really want to run multiple jobs per core. Note that the1324: C  unless you really want to run multiple jobs per core. Note that the
1325: C  nodes referred to here are really cores.1325: C  nodes referred to here are really cores.
1326: C1326: C
1327:       ELSE IF (WORD.EQ.'JOBSPERNODE') THEN1327:       ELSE IF (WORD.EQ.'JOBSPERNODE') THEN
1328:           WRITE(*,'(A)') 'keywords> ERROR: JOBSPERNODE is deprecated. Please generate a valid nodes.info file and use PBS instead'1328:          PBST=.TRUE.
1329:           STOP1329:          IF (NITEMS.GT.1) CALL READI(JPN)
1330: C         PBST=.TRUE.1330:          CALL GETNODES(NCPU)
1331: C         IF (NITEMS.GT.1) CALL READI(JPN) 
1332: C         CALL GETNODES(NCPU) 
1333:       ELSE IF (WORD.EQ.'KMC') THEN1331:       ELSE IF (WORD.EQ.'KMC') THEN
1334:          KMCT=.TRUE.1332:          KMCT=.TRUE.
1335:          NOPOINTS=.TRUE.1333:          NOPOINTS=.TRUE.
1336:          IF (NITEMS.GT.1) CALL READF(NKMCCYCLES)1334:          IF (NITEMS.GT.1) CALL READF(NKMCCYCLES)
1337:          IF (NITEMS.GT.2) CALL READF(PAIRTHRESH)1335:          IF (NITEMS.GT.2) CALL READF(PAIRTHRESH)
1338:          IF (NITEMS.GT.3) CALL READI(NCONNMIN)1336:          IF (NITEMS.GT.3) CALL READI(NCONNMIN)
1339: C1337: C
1340: C  Set KMC parameters: number of KMC runs for averages and value of PAIRTHRESH.1338: C  Set KMC parameters: number of KMC runs for averages and value of PAIRTHRESH.
1341: C  If the product of branching probabilities p12*p21 > PAIRTHRESH then we1339: C  If the product of branching probabilities p12*p21 > PAIRTHRESH then we
1342: C  renormalise this pair.1340: C  renormalise this pair.


r29831/printsummary.f90 2016-07-06 15:38:52.921892850 +0100 r29830/printsummary.f90 2016-07-06 15:38:53.617902263 +0100
 21: ! 21: !
 22: ! Printing some summary info about the KTN, to standard out  22: ! Printing some summary info about the KTN, to standard out 
 23:  23: 
 24:       SUBROUTINE PRINTSUMMARY 24:       SUBROUTINE PRINTSUMMARY
 25:       USE COMMONS, ONLY : EMIN, ETS, NMIN, NTS, NCONN, PLUS, MINUS, NCONNMIN, TSTHRESH, MAXBARRIER, MINBARRIER, & 25:       USE COMMONS, ONLY : EMIN, ETS, NMIN, NTS, NCONN, PLUS, MINUS, NCONNMIN, TSTHRESH, MAXBARRIER, MINBARRIER, &
 26:                           KPLUS, KMINUS 26:                           KPLUS, KMINUS
 27:       IMPLICIT NONE 27:       IMPLICIT NONE
 28:  28: 
 29:       INTEGER :: NDEAD, MINTARG, NLEFTTS, NCOUNT, NDONE, J1, NLEFTMIN, MAXTSINDEX, MINTSINDEX, CURRTARG 29:       INTEGER :: NDEAD, MINTARG, NLEFTTS, NCOUNT, NDONE, J1, NLEFTMIN, MAXTSINDEX, MINTSINDEX, CURRTARG
 30:       INTEGER :: NDISTA(NMIN) 30:       INTEGER :: NDISTA(NMIN)
 31:       DOUBLE PRECISION :: MAXTS, MINTS, GMIN, CURRENERGY, LOWBARRIER, HIGHBARRIER 31:       DOUBLE PRECISION :: MAXTS, MINTS, GMIN, CURRENERGY
 32:       DOUBLE PRECISION :: BARRIER_MAX, BARRIER_MIN 
 33:       DOUBLE PRECISION :: CUT_UNDERFLOW=-300.0D0 ! set here but NOT used in the current coding of checkTS 32:       DOUBLE PRECISION :: CUT_UNDERFLOW=-300.0D0 ! set here but NOT used in the current coding of checkTS
 34:       LOGICAL :: CHANGED 33:       LOGICAL :: CHANGED
 35:       LOGICAL :: DEADTS(NTS), INCC(NMIN) 34:       LOGICAL :: DEADTS(NTS), INCC(NMIN)
 36:  35: 
 37:       WRITE(*,'(A)')'*************************************************************************************************' 36:       WRITE(*,'(A)')'*************************************************************************************************'
 38:       WRITE(*,'(A)')'printsummary> Network summary information' 37:       WRITE(*,'(A)')'printsummary> Network summary information'
 39:       WRITE(*,'(A)')'using thresholds:' 38:       WRITE(*,'(A)')'using thresholds:'
 40:       WRITE(*,'(A,I8)')'NCONNMIN=',NCONNMIN 39:       WRITE(*,'(A,I8)')'NCONNMIN=',NCONNMIN
 41:       WRITE(*,'(A,G20.10)')'MAXBARRIER=',MAXBARRIER 40:       WRITE(*,'(A,G20.10)')'MAXBARRIER=',MAXBARRIER
 42:       WRITE(*,'(A,G20.10)')'MINBARRIER=',MINBARRIER 41:       WRITE(*,'(A,G20.10)')'MINBARRIER=',MINBARRIER
 70:                       PLUS(J1),MINUS(J1),.FALSE.,CUT_UNDERFLOW,DEADTS(J1)) 69:                       PLUS(J1),MINUS(J1),.FALSE.,CUT_UNDERFLOW,DEADTS(J1))
 71:          IF (.NOT.DEADTS(J1)) NLEFTTS=NLEFTTS+1 70:          IF (.NOT.DEADTS(J1)) NLEFTTS=NLEFTTS+1
 72:       ENDDO 71:       ENDDO
 73:       WRITE(*,'(A,I8,A,I8)')'There are ',NLEFTTS,' valid transition states out of ',NTS 72:       WRITE(*,'(A,I8,A,I8)')'There are ',NLEFTTS,' valid transition states out of ',NTS
 74:  73: 
 75:       IF ((NMIN-NDEAD.EQ.0).OR.(NLEFTTS.EQ.0)) RETURN ! sanity check 74:       IF ((NMIN-NDEAD.EQ.0).OR.(NLEFTTS.EQ.0)) RETURN ! sanity check
 76:  75: 
 77:       MINTS=HUGE(1.0D0) 76:       MINTS=HUGE(1.0D0)
 78:       MAXTS=-HUGE(1.0D0) 77:       MAXTS=-HUGE(1.0D0)
 79:       DO J1=1,NTS 78:       DO J1=1,NTS
 80:          IF (DEADTS(J1)) CYCLE 79:          IF (ETS(J1).LE.MINTS .AND. (.NOT.DEADTS(J1))) THEN
 81:          IF (ETS(J1).LE.MINTS) THEN 
 82:             MINTS=ETS(J1) 80:             MINTS=ETS(J1)
 83:             MINTSINDEX=J1 81:             MINTSINDEX=J1
 84:          ENDIF 82:          ENDIF
 85:          IF (ETS(J1).GE.MAXTS) THEN 83:          IF (ETS(J1).GE.MAXTS .AND. (.NOT.DEADTS(J1))) THEN
 86:             MAXTS=ETS(J1) 84:             MAXTS=ETS(J1)
 87:             MAXTSINDEX=J1 85:             MAXTSINDEX=J1
 88:          ENDIF 86:          ENDIF
 89:       ENDDO 87:       ENDDO
 90:       WRITE(*,'(A,I8,A,G20.10)')'Highest valid transition state at index ',MAXTSINDEX,' with energy ',MAXTS 88:       WRITE(*,'(A,I8,A,G20.10)')'Highest valid transition state at index ',MAXTSINDEX,' with energy ',MAXTS
 91:       WRITE(*,'(A,I8,A,G20.10)')'Lowest valid transition state at index ',MINTSINDEX,' with energy ',MINTS 89:       WRITE(*,'(A,I8,A,G20.10)')'Lowest valid transition state at index ',MINTSINDEX,' with energy ',MINTS
 92:  90: 
 93:       LOWBARRIER=HUGE(1.0D0) 
 94:       HIGHBARRIER=-HUGE(1.0D0) 
 95:       DO J1=1,NTS 
 96:          IF (DEADTS(J1)) CYCLE 
 97:          print *,ETS(J1)-EMIN(PLUS(J1)) 
 98:          print *,ETS(J1)-EMIN(MINUS(J1)) 
 99:          BARRIER_MAX = MAX(ETS(J1)-EMIN(PLUS(J1)),ETS(J1)-EMIN(MINUS(J1))) 
100:          BARRIER_MIN = MIN(ETS(J1)-EMIN(PLUS(J1)),ETS(J1)-EMIN(MINUS(J1))) 
101:          IF (BARRIER_MIN.LE.LOWBARRIER) THEN 
102:             LOWBARRIER=BARRIER_MIN 
103:             MINTSINDEX=J1 
104:          ENDIF 
105:          IF (BARRIER_MAX.GE.HIGHBARRIER) THEN 
106:             HIGHBARRIER=BARRIER_MAX 
107:             MAXTSINDEX=J1 
108:          ENDIF 
109:       ENDDO 
110:       WRITE(*,'(A,G20.10,A,I8)')'Highest valid barrier ',HIGHBARRIER,' for TS index ',MAXTSINDEX 
111:       WRITE(*,'(A,G20.10,A,I8)')'Lowest valid barrier ',LOWBARRIER,' for TS index ',MINTSINDEX 
112:  
113:       WRITE(*,'(A)')'Connected components for valid minima and TSs:' 91:       WRITE(*,'(A)')'Connected components for valid minima and TSs:'
114:       WRITE(*,'(A)')'CC #        # minima    index of lowest min    energy of lowest min' 92:       WRITE(*,'(A)')'CC #        # minima    index of lowest min    energy of lowest min'
115:       NCOUNT=0 93:       NCOUNT=0
116:       NDONE=0 94:       NDONE=0
117:       CURRTARG=MINTARG ! start with the valid global minimum  95:       CURRTARG=MINTARG ! start with the valid global minimum 
118:       CURRENERGY=GMIN 96:       CURRENERGY=GMIN
119:       INCC(1:NMIN)=.FALSE. 97:       INCC(1:NMIN)=.FALSE.
120:       DO  98:       DO 
121:          NCOUNT=NCOUNT+1 99:          NCOUNT=NCOUNT+1
122:          NDISTA(1:NMIN)=1000000100:          NDISTA(1:NMIN)=1000000


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0