hdiff output

r33153/dumpit.f 2017-08-10 10:30:12.028099271 +0100 r33152/dumpit.f 2017-08-10 10:30:14.012125655 +0100
 73:       ENDIF 73:       ENDIF
 74:       IF (PV) WRITE(3,'(A,3F20.10)') 'PV ',PRESS,PVCONV,PVTOL 74:       IF (PV) WRITE(3,'(A,3F20.10)') 'PV ',PRESS,PVCONV,PVTOL
 75:       IF (PUSHOFF.NE.0.01D0) WRITE(3,'(A,T15,F20.10)') 'PUSHOFF',PUSHOFF 75:       IF (PUSHOFF.NE.0.01D0) WRITE(3,'(A,T15,F20.10)') 'PUSHOFF',PUSHOFF
 76:       IF (DFTBT) WRITE(3,'(A)') 'DFTB' 76:       IF (DFTBT) WRITE(3,'(A)') 'DFTB'
 77:       IF (CADPAC) WRITE(3,'(A,A,A)') 'CADPAC ',SYS,EDITIT 77:       IF (CADPAC) WRITE(3,'(A,A,A)') 'CADPAC ',SYS,EDITIT
 78:       IF (GAMESSUS) WRITE(3,'(A,A,A)') 'GAMESS-US ',SYS,EDITIT 78:       IF (GAMESSUS) WRITE(3,'(A,A,A)') 'GAMESS-US ',SYS,EDITIT
 79:       IF (GAMESSUK) WRITE(3,'(A,A,A)') 'GAMESS-UK ',SYS,EDITIT 79:       IF (GAMESSUK) WRITE(3,'(A,A,A)') 'GAMESS-UK ',SYS,EDITIT
 80:       IF (GAUSSIAN) WRITE(3,'(A)') 'GAUSSIAN' 80:       IF (GAUSSIAN) WRITE(3,'(A)') 'GAUSSIAN'
 81:       IF (GAUSSIAN03) WRITE(3,'(A)') 'GAUSSIAN03' 81:       IF (GAUSSIAN03) WRITE(3,'(A)') 'GAUSSIAN03'
 82:       IF (GAUSSIAN09) WRITE(3,'(A)') 'GAUSSIAN09' 82:       IF (GAUSSIAN09) WRITE(3,'(A)') 'GAUSSIAN09'
 83:       IF (GAUSSIAN09) WRITE(3,'(A)') 'GAUSSIAN16' 
 84:       IF (ONETEP) THEN 83:       IF (ONETEP) THEN
 85:          IF (DFTP) THEN 84:          IF (DFTP) THEN
 86:             WRITE(3,'(A,A)') 'ONETEP ',ONETEPJOB 85:             WRITE(3,'(A,A)') 'ONETEP ',ONETEPJOB
 87:          ELSE 86:          ELSE
 88:             WRITE(3,'(A,A)') 'ONETEPC ',ONETEPJOB 87:             WRITE(3,'(A,A)') 'ONETEPC ',ONETEPJOB
 89:          ENDIF 88:          ENDIF
 90:       ENDIF 89:       ENDIF
 91:       IF (CASTEP) THEN 90:       IF (CASTEP) THEN
 92:          IF (DFTP) THEN 91:          IF (DFTP) THEN
 93:             WRITE(3,'(A,A)') 'CASTEP ',CASTEPJOB 92:             WRITE(3,'(A,A)') 'CASTEP ',CASTEPJOB


r33153/fetchz.f 2017-08-10 10:30:12.272102506 +0100 r33152/fetchz.f 2017-08-10 10:30:14.232128584 +0100
1547:          ELSE IF (GAMESSUS) THEN1547:          ELSE IF (GAMESSUS) THEN
1548:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' GAMESS-US atoms'1548:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' GAMESS-US atoms'
1549:             WRITE(*,'(A,A,A,A)') ' fetchz> System name: ',SYS(1:LSYS),', edit command: ',EDITIT1549:             WRITE(*,'(A,A,A,A)') ' fetchz> System name: ',SYS(1:LSYS),', edit command: ',EDITIT
1550:          ELSE IF (GAMESSUK) THEN1550:          ELSE IF (GAMESSUK) THEN
1551:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' GAMESS-UK atoms'1551:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' GAMESS-UK atoms'
1552:             WRITE(*,'(A,A,A,A)') ' fetchz> System name: ',SYS(1:LSYS),', edit command: ',EDITIT1552:             WRITE(*,'(A,A,A,A)') ' fetchz> System name: ',SYS(1:LSYS),', edit command: ',EDITIT
1553:          ELSE IF (GAUSSIAN) THEN1553:          ELSE IF (GAUSSIAN) THEN
1554:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Gaussian atoms'1554:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Gaussian atoms'
1555:          ELSE IF (GAUSSIAN03) THEN1555:          ELSE IF (GAUSSIAN03) THEN
1556:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Gaussian03 atoms'1556:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Gaussian03 atoms'
1557:          ELSE IF (GAUSSIAN09) THEN 
1558:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Gaussian09 atoms' 
1559:          ELSE IF (GAUSSIAN16) THEN 
1560:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Gaussian16 atoms' 
1561:          ELSE IF (DFTBT) THEN1557:          ELSE IF (DFTBT) THEN
1562:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Tiffany TB atoms'1558:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' Tiffany TB atoms'
1563:          ELSE IF (ZSYM(NATOMS).EQ.'SV') THEN1559:          ELSE IF (ZSYM(NATOMS).EQ.'SV') THEN
1564:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' MSEVB atoms '1560:             WRITE(*,'(A,I4,A,I4,A)') ' fetchz> ',NOPT,' Cartesian coordinates will be optimised for ',NATOMS,' MSEVB atoms '
1565: 1561: 
1566: ! Set some values1562: ! Set some values
1567: 1563: 
1568:             num_eig = (NATOMS-1)/31564:             num_eig = (NATOMS-1)/3
1569:             num_hyd = NATOMS - num_eig1565:             num_hyd = NATOMS - num_eig
1570: 1566: 


r33153/key.f90 2017-08-10 10:30:12.908111296 +0100 r33152/key.f90 2017-08-10 10:30:14.456131565 +0100
 51:      &        NONEBMAX, READMASST, ONEDAPBCT, ONEDPBCT, INVTONEDPBCT, INVTTWODPBCT, TWODAPBCT, TWODPBCT, THREEDAPBCT, & 51:      &        NONEBMAX, READMASST, ONEDAPBCT, ONEDPBCT, INVTONEDPBCT, INVTTWODPBCT, TWODAPBCT, TWODPBCT, THREEDAPBCT, &
 52:      &        THREEDPBCT, FOURDAPBCT, FOURDPBCT, MODEDOWNT, CHEMSHIFT, TTM3T, & 52:      &        THREEDPBCT, FOURDAPBCT, FOURDPBCT, MODEDOWNT, CHEMSHIFT, TTM3T, &
 53:      &        NOINVERSION, INVERTPT, KNOWVECS, PMPATHT, AAORIENTT, MULTIJOBT, QUIPARGSTRT, QUIPPARAMST, HESSDUMPT, & 53:      &        NOINVERSION, INVERTPT, KNOWVECS, PMPATHT, AAORIENTT, MULTIJOBT, QUIPARGSTRT, QUIPPARAMST, HESSDUMPT, &
 54:      &        CLASSICALRATEST, TSPLITTINGT, HESSREADT, INSTANTONOPTT,INSTANTONSTARTDUMPT,VARSTEPOPTT, MOLPRO, REAXFFT, & 54:      &        CLASSICALRATEST, TSPLITTINGT, HESSREADT, INSTANTONOPTT,INSTANTONSTARTDUMPT,VARSTEPOPTT, MOLPRO, REAXFFT, &
 55:      &        EIGENONLY,OVERCONV, GLJT,CLSTRINGT,CLSTRINGTST, PHI4MODT, EX1DT, MCPATHT, MCBIAST, RPHT, TWISTT, MCPATH2T, & 55:      &        EIGENONLY,OVERCONV, GLJT,CLSTRINGT,CLSTRINGTST, PHI4MODT, EX1DT, MCPATHT, MCBIAST, RPHT, TWISTT, MCPATH2T, &
 56:      &        PBST, SSHT, GAUSSIAN03, GAUSSIAN09, CPPNEBT, CUDAT, CUDATIMET, TRUSTMODET,MODELOST, METRICTENSOR, INTSPRINGACTIVET, & 56:      &        PBST, SSHT, GAUSSIAN03, GAUSSIAN09, CPPNEBT, CUDAT, CUDATIMET, TRUSTMODET,MODELOST, METRICTENSOR, INTSPRINGACTIVET, &
 57:      &        PERMGUESS, QCIPERMCHECK, DUMPFRQST, MULTIPOTT, MLP3T, MLPB3T, DUMPBESTPATH, ALIGNRBST, AVOID_COLLISIONS, MLPPROB, & 57:      &        PERMGUESS, QCIPERMCHECK, DUMPFRQST, MULTIPOTT, MLP3T, MLPB3T, DUMPBESTPATH, ALIGNRBST, AVOID_COLLISIONS, MLPPROB, &
 58:      &        MALONALDEHYDE, SIO2PT, MLPNEWREG, DJWRBT, STEALTHYT, STEALTV, LJADDT, MLPB3NEWT, MLPVB3T, & 58:      &        MALONALDEHYDE, SIO2PT, MLPNEWREG, DJWRBT, STEALTHYT, STEALTV, LJADDT, MLPB3NEWT, MLPVB3T, &
 59:      &        QCIPOTT, QCIPOT2T, QCIRADSHIFTT, QCINOREPINT, QCIAMBERT, SLERPT, NOTRANSROTT, MAXGAPT, BULKBOXT, GDSQT, FLATTESTT, & 59:      &        QCIPOTT, QCIPOT2T, QCIRADSHIFTT, QCINOREPINT, QCIAMBERT, SLERPT, NOTRANSROTT, MAXGAPT, BULKBOXT, GDSQT, FLATTESTT, &
 60:      &        MLQT, MLQPROB, LJADD2T, MACROIONT, NOREGBIAS, PYADDT, PYADD2T, SANDBOXT, LJADD3T, LJADD4T, & 60:      &        MLQT, MLQPROB, LJADD2T, MACROIONT, NOREGBIAS, PYADDT, PYADD2T, SANDBOXT, LJADD3T, LJADD4T, &
 61:      &        MBPOLT, MULTIJOB_MACHINET, DUMPDATA_MACHINET, PLUSSIDET, MINUSSIDET, PUSHOPTT, MLPVB3NNT, GAUSSIAN16 61:      &        MBPOLT, MULTIJOB_MACHINET, DUMPDATA_MACHINET, PLUSSIDET, MINUSSIDET, PUSHOPTT, MLPVB3NNT
 62:  62: 
 63: ! sy349 > for testing the flatpath after dneb 63: ! sy349 > for testing the flatpath after dneb
 64:       !LOGICAL, ALLOCATABLE :: FLATPATHT(:) 64:       !LOGICAL, ALLOCATABLE :: FLATPATHT(:)
 65:       LOGICAL FLATPATHT 65:       LOGICAL FLATPATHT
 66:  66: 
 67: ! bf269 > polymer in a pore (non-bonding (LJ) energy from neighbours is not subtracted) 67: ! bf269 > polymer in a pore (non-bonding (LJ) energy from neighbours is not subtracted)
 68:       LOGICAL :: PORE8T = .FALSE. ! add 8th power cylindrical pore to the potential? 68:       LOGICAL :: PORE8T = .FALSE. ! add 8th power cylindrical pore to the potential?
 69:       INTEGER :: PORE8_AXIS = 3 ! principal axis of the cylindric pore (1:x, 2:y, 3:z) 69:       INTEGER :: PORE8_AXIS = 3 ! principal axis of the cylindric pore (1:x, 2:y, 3:z)
 70:       DOUBLE PRECISION :: PORE8_ENERGY = 1.0d1 ! energy of the pore when radius = 1 70:       DOUBLE PRECISION :: PORE8_ENERGY = 1.0d1 ! energy of the pore when radius = 1
 71:       LOGICAL :: HARMPOLYT = .FALSE. ! add harmonic bonds between the beads 71:       LOGICAL :: HARMPOLYT = .FALSE. ! add harmonic bonds between the beads


r33153/keywords.f 2017-08-10 10:30:13.252115551 +0100 r33152/keywords.f 2017-08-10 10:30:14.680134545 +0100
255:          PVSTEPS=100255:          PVSTEPS=100
256:          NBOXTS=1256:          NBOXTS=1
257: 257: 
258:          VARIABLES=.FALSE.258:          VARIABLES=.FALSE.
259:          NZERO=0259:          NZERO=0
260:          EVCUT=0.0D0260:          EVCUT=0.0D0
261: 261: 
262:          GAUSSIAN=.FALSE.262:          GAUSSIAN=.FALSE.
263:          GAUSSIAN03=.FALSE.263:          GAUSSIAN03=.FALSE.
264:          GAUSSIAN09=.FALSE.264:          GAUSSIAN09=.FALSE.
265:          GAUSSIAN16=.FALSE. 
266:          CADPAC=.FALSE.265:          CADPAC=.FALSE.
267:          GAMESSUS=.FALSE.266:          GAMESSUS=.FALSE.
268:          GAMESSUK=.FALSE.267:          GAMESSUK=.FALSE.
269:          CASTEP=.FALSE.268:          CASTEP=.FALSE.
270:          CASTEPJOB=''269:          CASTEPJOB=''
271:          QCHEM=.FALSE.270:          QCHEM=.FALSE.
272:          QCHEMES=.FALSE.271:          QCHEMES=.FALSE.
273:          QCHEMESNAO=0272:          QCHEMESNAO=0
274:          QCHEMESNMO=0273:          QCHEMESNMO=0
275:          QCHEMESNZERO=0274:          QCHEMESNZERO=0
3175:             CALL READI(GAUSSIANMULTI)3174:             CALL READI(GAUSSIANMULTI)
3176:             OPEN(LUNIT,FILE='gaussian.xyz',STATUS='REPLACE')3175:             OPEN(LUNIT,FILE='gaussian.xyz',STATUS='REPLACE')
3177:             CLOSE(LUNIT)3176:             CLOSE(LUNIT)
3178: ! Gaussian09 interface            3177: ! Gaussian09 interface            
3179:          ElSE IF (WORD.EQ.'GAUSSIAN09') THEN3178:          ElSE IF (WORD.EQ.'GAUSSIAN09') THEN
3180:             GAUSSIAN09=.TRUE.3179:             GAUSSIAN09=.TRUE.
3181:             CALL READI(GAUSSIANCHARGE)3180:             CALL READI(GAUSSIANCHARGE)
3182:             CALL READI(GAUSSIANMULTI)3181:             CALL READI(GAUSSIANMULTI)
3183:             OPEN(LUNIT,FILE='gaussian.xyz',STATUS='REPLACE')3182:             OPEN(LUNIT,FILE='gaussian.xyz',STATUS='REPLACE')
3184:             CLOSE(LUNIT)3183:             CLOSE(LUNIT)
3185: ! Gaussian16 interface             
3186:          ElSE IF (WORD.EQ.'GAUSSIAN16') THEN 
3187:             GAUSSIAN16=.TRUE. 
3188:             CALL READI(GAUSSIANCHARGE) 
3189:             CALL READI(GAUSSIANMULTI) 
3190:             OPEN(LUNIT,FILE='gaussian.xyz',STATUS='REPLACE') 
3191:             CLOSE(LUNIT) 
3192: ! DC430 >3184: ! DC430 >
3193: 3185: 
3194:          ELSE IF (WORD .EQ. 'GB') THEN3186:          ELSE IF (WORD .EQ. 'GB') THEN
3195:             GBT   = .TRUE.3187:             GBT   = .TRUE.
3196:             RBAAT = .TRUE.3188:             RBAAT = .TRUE.
3197:             CALL READF(GBKAPPA)3189:             CALL READF(GBKAPPA)
3198:             CALL READF(GBKAPPRM)3190:             CALL READF(GBKAPPRM)
3199:             CALL READF(GBMU)3191:             CALL READF(GBMU)
3200:             CALL READF(GBNU)3192:             CALL READF(GBNU)
3201:             CALL READF(GBSIGNOT)3193:             CALL READF(GBSIGNOT)


r33153/potential.f 2017-08-10 10:30:13.516119065 +0100 r33152/potential.f 2017-08-10 10:30:14.912137632 +0100
3213:                   J4=3*(J4-1)+(J2-1)/9+13213:                   J4=3*(J4-1)+(J2-1)/9+1
3214:                   HGAUSS(J4,J3)=HESS(J2,J1)3214:                   HGAUSS(J4,J3)=HESS(J2,J1)
3215:                ENDDO3215:                ENDDO
3216:             ENDDO3216:             ENDDO
3217:             DO J1=1,3*NATOMS3217:             DO J1=1,3*NATOMS
3218:                DO J2=1,3*NATOMS3218:                DO J2=1,3*NATOMS
3219:                   HESS(J2,J1)=HGAUSS(J2,J1)3219:                   HESS(J2,J1)=HGAUSS(J2,J1)
3220:                ENDDO3220:                ENDDO
3221:             ENDDO3221:             ENDDO
3222:             CLOSE(15)3222:             CLOSE(15)
3223:          ELSE IF (GAUSSIAN03.OR.GAUSSIAN09.OR.GAUSSIAN16) THEN3223:          ELSE IF (GAUSSIAN03.OR.GAUSSIAN09) THEN
3224:             IF (GAUSSIAN03) WRITE(*,'(A)') ' potential> Calling GAUSSIAN03'3224:             IF (GAUSSIAN03) WRITE(*,'(A)') ' potential> Calling GAUSSIAN03'
3225:             IF (GAUSSIAN09) WRITE(*,'(A)') ' potential> Calling GAUSSIAN09'3225:             IF (GAUSSIAN09) WRITE(*,'(A)') ' potential> Calling GAUSSIAN09'
3226:             IF (GAUSSIAN16) WRITE(*,'(A)') ' potential> Calling GAUSSIAN16' 
3227:             CALL SYSTEM('cp gaussian.inp.pre gaussian.inp' )3226:             CALL SYSTEM('cp gaussian.inp.pre gaussian.inp' )
3228:             LUNIT=GETUNIT()3227:             LUNIT=GETUNIT()
3229:             OPEN(LUNIT,file="gaussian.inp",STATUS="OLD",POSITION="APPEND")3228:             OPEN(LUNIT,file="gaussian.inp",STATUS="OLD",POSITION="APPEND")
3230:             INQUIRE(FILE='gaussian.chk',EXIST=GAUSSIANTEST) 3229:             INQUIRE(FILE='gaussian.chk',EXIST=GAUSSIANTEST) 
3231:             IF(GAUSSIANTEST)THEN3230:             IF(GAUSSIANTEST)THEN
3232:                WRITE(LUNIT,"(A10)") "GUESS=READ"3231:                WRITE(LUNIT,"(A10)") "GUESS=READ"
3233:             ENDIF3232:             ENDIF
3234:             WRITE(LUNIT,*) ""3233:             WRITE(LUNIT,*) ""
3235:             WRITE(LUNIT,*) "OPTIM INTERFACE"3234:             WRITE(LUNIT,*) "OPTIM INTERFACE"
3236:             WRITE(LUNIT,*) ""3235:             WRITE(LUNIT,*) ""
3253:             IF (SSTEST) THEN3252:             IF (SSTEST) THEN
3254:                CALL SYSTEM('sed "s/Force/Freq=Noraman/Ig" gaussian.inp >gaussian.inp.temp')3253:                CALL SYSTEM('sed "s/Force/Freq=Noraman/Ig" gaussian.inp >gaussian.inp.temp')
3255:                CALL SYSTEM('mv gaussian.inp.temp gaussian.inp')3254:                CALL SYSTEM('mv gaussian.inp.temp gaussian.inp')
3256:             ENDIF3255:             ENDIF
3257:             INQUIRE(FILE='gaussian.out',EXIST=GAUSSIANTEST) 3256:             INQUIRE(FILE='gaussian.out',EXIST=GAUSSIANTEST) 
3258:             IF(GAUSSIANTEST) THEN3257:             IF(GAUSSIANTEST) THEN
3259:                CALL SYSTEM('mv gaussian.out gaussian.out.old')3258:                CALL SYSTEM('mv gaussian.out gaussian.out.old')
3260:             ENDIF3259:             ENDIF
3261:             IF (GAUSSIAN03) CALL SYSTEM('g03 <gaussian.inp> gaussian.out')3260:             IF (GAUSSIAN03) CALL SYSTEM('g03 <gaussian.inp> gaussian.out')
3262:             IF (GAUSSIAN09) CALL SYSTEM('g09 <gaussian.inp> gaussian.out')3261:             IF (GAUSSIAN09) CALL SYSTEM('g09 <gaussian.inp> gaussian.out')
3263:             IF (GAUSSIAN16) CALL SYSTEM('g16 <gaussian.inp> gaussian.out') 
3264:             IF (SSTEST) THEN3262:             IF (SSTEST) THEN
3265:                !WRITE(*,'(A)') ' potential> extract GAUSSIAN second derivatives'3263:                !WRITE(*,'(A)') ' potential> extract GAUSSIAN second derivatives'
3266:                J1=(1+3*NATOMS)*3*NATOMS/23264:                J1=(1+3*NATOMS)*3*NATOMS/2
3267:                IF (.NOT.ALLOCATED(HGAUSS)) ALLOCATE(HGAUSS(J1,1))3265:                IF (.NOT.ALLOCATED(HGAUSS)) ALLOCATE(HGAUSS(J1,1))
3268:                IF(MOD(J1,5)==0) then3266:                IF(MOD(J1,5)==0) then
3269:                   J1=J1/53267:                   J1=J1/5
3270:                ELSE3268:                ELSE
3271:                   J1=J1/5+13269:                   J1=J1/5+1
3272:                ENDIF3270:                ENDIF
3273:                WRITE (STRING, *) J13271:                WRITE (STRING, *) J1


r33153/secdiag.f90 2017-08-10 10:30:13.768122415 +0100 r33152/secdiag.f90 2017-08-10 10:30:15.152140824 +0100
 63: !      XRMS = 0.0D0 63: !      XRMS = 0.0D0
 64:       DIFF=1.0D-3 64:       DIFF=1.0D-3
 65:       IF (NIH2LEPST .or. NIMET .or. NIHEAM7T .or. NIHLEPST .or. NIHPAIRONLYT) DIFF=1.0D-5 65:       IF (NIH2LEPST .or. NIMET .or. NIHEAM7T .or. NIHLEPST .or. NIHPAIRONLYT) DIFF=1.0D-5
 66:       IF (PYADD2T) DIFF=1.0D-5 66:       IF (PYADD2T) DIFF=1.0D-5
 67:       IF (CHARMMDFTBT) DIFF=1.0D-2 67:       IF (CHARMMDFTBT) DIFF=1.0D-2
 68:       IF (AMHT) DIFF=1.0D-2 68:       IF (AMHT) DIFF=1.0D-2
 69:       IF (ZSYM(NATOMS).EQ.'GO') DIFF=2.0D-3 69:       IF (ZSYM(NATOMS).EQ.'GO') DIFF=2.0D-3
 70:       IF (GAMESSUK.OR.GAMESSUS) DIFF=1.0D-2 70:       IF (GAMESSUK.OR.GAMESSUS) DIFF=1.0D-2
 71:       IF (CADPAC) DIFF=1.0D-2 71:       IF (CADPAC) DIFF=1.0D-2
 72:       IF (CASTEP) DIFF=0.01D0 72:       IF (CASTEP) DIFF=0.01D0
 73:       IF (GAUSSIAN03.OR.GAUSSIAN09.OR.GAUSSIAN16) DIFF=1.0D-2 73:       IF (GAUSSIAN03.OR.GAUSSIAN09) DIFF=1.0D-2
 74:       IF (QCHEM) DIFF=0.01D0 74:       IF (QCHEM) DIFF=0.01D0
 75:       IF (VASP) DIFF=0.01D0 75:       IF (VASP) DIFF=0.01D0
 76:       IF (ONETEP) DIFF=0.01D0 76:       IF (ONETEP) DIFF=0.01D0
 77:       IF (CP2K) DIFF=0.001D0  77:       IF (CP2K) DIFF=0.001D0 
 78:       IF (CPMD) DIFF=0.04D0 78:       IF (CPMD) DIFF=0.04D0
 79: !     IF (CHRMMT) DIFF=5.0D-2 79: !     IF (CHRMMT) DIFF=5.0D-2
 80:       IF (CHRMMT) DIFF=0.01D0 80:       IF (CHRMMT) DIFF=0.01D0
 81: !     IF (DFTBT) DIFF=1.0D-3 81: !     IF (DFTBT) DIFF=1.0D-3
 82: ! 82: !
 83: !  Must read VEC into LOCALV because we are going to play with the vector in 83: !  Must read VEC into LOCALV because we are going to play with the vector in


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0