hdiff output

r33147/key.f90 2017-08-09 09:30:13.308938959 +0100 r33146/key.f90 2017-08-09 09:30:13.748944809 +0100
 79:       INTEGER :: GTHOMMET, NGTHORI, GTHOMPOT, GTHOMSONBINN 79:       INTEGER :: GTHOMMET, NGTHORI, GTHOMPOT, GTHOMSONBINN
 80:       DOUBLE PRECISION :: GThomsonSigma, GThomsonRho 80:       DOUBLE PRECISION :: GThomsonSigma, GThomsonRho
 81:  81: 
 82:       LOGICAL :: MACHINE=.FALSE., NEBMIND=.TRUE. 82:       LOGICAL :: MACHINE=.FALSE., NEBMIND=.TRUE.
 83:  83: 
 84:       DOUBLE PRECISION :: TRAD, RESIZE, JZ, PUSHOFF, TOLD, TOLE, PCUT, DAMP, RADIUS, TWOEVAL, PVCONV, PVTOL, & 84:       DOUBLE PRECISION :: TRAD, RESIZE, JZ, PUSHOFF, TOLD, TOLE, PCUT, DAMP, RADIUS, TWOEVAL, PVCONV, PVTOL, &
 85:      &        OMEGA, GALPHA, MALPHA1, FOH, FIH, FTD, FD5H, T12FAC, PRESS, MAXERISE, RK_R, RK_THETA, & 85:      &        OMEGA, GALPHA, MALPHA1, FOH, FIH, FTD, FD5H, T12FAC, PRESS, MAXERISE, RK_R, RK_THETA, &
 86:      &        CONVU, CONVR, SYMCUT, MXSTP, GFRACTION, CEIG, REPELPUSH, XMAXERISE, GCCONV, NEBK, & 86:      &        CONVU, CONVR, SYMCUT, MXSTP, GFRACTION, CEIG, REPELPUSH, XMAXERISE, GCCONV, NEBK, &
 87:      &        MFRACTION1, MFRACTION2, MALPHA2, ACCLAN, SHIFTLAN, CUTLAN, GMAX, EVCUT, NEBDGUESS, INTDGUESS,  & 87:      &        MFRACTION1, MFRACTION2, MALPHA2, ACCLAN, SHIFTLAN, CUTLAN, GMAX, EVCUT, NEBDGUESS, INTDGUESS,  &
 88:      &        MINMAX, MAXMAX, C6PP, C6MM, C6PM, PUSHCUT, MAXBFGS, MAXXBFGS, MAXMBFGS, MAXNEBBFGS, MAXINTBFGS, & 88:      &        MINMAX, MAXMAX, C6PP, C6MM, C6PM, PUSHCUT, MAXBFGS, MAXXBFGS, MAXMBFGS, MAXNEBBFGS, MAXINTBFGS, &
 89:      &        APP, AMM, APM, XQP,XQM, ALPHAP, ALPHAM, PHIG, SHIFTV, DGUESS, XDGUESS, TIMELIMIT, TSTART, & 89:      &        APP, AMM, APM, XQP,XQM, ALPHAP, ALPHAM, PHIG, SHIFTV, SHIFTL(6), DGUESS, XDGUESS, TIMELIMIT, TSTART, &
 90:      &        INTEPSILON, FRAMEEDIFF, MORPHERISE, MORPHEMAX, MORPHMXSTP, MAXTSENERGY, MAXBARRIER, GCMXSTP, GCSCALE, & 90:      &        INTEPSILON, FRAMEEDIFF, MORPHERISE, MORPHEMAX, MORPHMXSTP, MAXTSENERGY, MAXBARRIER, GCMXSTP, GCSCALE, &
 91:      &        FREEZETOL, DESMAXEJUMP, DESMAXAVGE, STOCKMU, STOCKLAMBDA, EDIFFTOL, GEOMDIFFTOL, DNEBSWITCH, & 91:      &        FREEZETOL, DESMAXEJUMP, DESMAXAVGE, STOCKMU, STOCKLAMBDA, EDIFFTOL, GEOMDIFFTOL, DNEBSWITCH, &
 92:      &        BHACCREJ, BHSTEPSIZE, BHCONV, BHTEMP, BHDISTTHRESH, BHMAXENERGY, BHK, DIJKSTRALOCAL, KADJUSTFRAC, KADJUSTTOL, & 92:      &        BHACCREJ, BHSTEPSIZE, BHCONV, BHTEMP, BHDISTTHRESH, BHMAXENERGY, BHK, DIJKSTRALOCAL, KADJUSTFRAC, KADJUSTTOL, &
 93:      &        BHSFRAC, HMCEIG, HMMXSTP, HMEVMAX, BBRGAM, BBREPS, BBRSIGMA1, BBRSIGMA2, BBRALPHA, BBRCONV, STOCKZTOL, & 93:      &        BHSFRAC, HMCEIG, HMMXSTP, HMEVMAX, BBRGAM, BBREPS, BBRSIGMA1, BBRSIGMA2, BBRALPHA, BBRCONV, STOCKZTOL, &
 94:      &        GBKAPPA, GBKAPPRM, GBMU, GBNU, GBSIGNOT, GBEPSNOT, GBCHI, GBCHIPRM, & 94:      &        GBKAPPA, GBKAPPRM, GBMU, GBNU, GBSIGNOT, GBEPSNOT, GBCHI, GBCHIPRM, &
 95:      &        PYA1(3), PYA2(3), PYSIGNOT, PYEPSNOT, NEBKINITIAL, NEBKFINAL, NEBFACTOR, & 95:      &        PYA1(3), PYA2(3), PYSIGNOT, PYEPSNOT, NEBKINITIAL, NEBKFINAL, NEBFACTOR, &
 96:      &        DBEPSBB, DBEPSAB, DBSIGBB, DBSIGAB, DBPMU, YEPS, YKAPPA, MAXMAXBARRIER, PAPALP, PAPEPS, PAPS, PAPCD, & 96:      &        DBEPSBB, DBEPSAB, DBSIGBB, DBSIGAB, DBPMU, YEPS, YKAPPA, MAXMAXBARRIER, PAPALP, PAPEPS, PAPS, PAPCD, &
 97:      &        RHOCC0, RHOCC10, RHOCC20,  RHOHH0, RHOHH10, RHOHH20, RHOCH0, RHOC10H, RHOCH10, RHOC20H, RHOCH20, & 97:      &        RHOCC0, RHOCC10, RHOCC20,  RHOHH0, RHOHH10, RHOHH20, RHOCH0, RHOC10H, RHOCH10, RHOC20H, RHOCH20, &
 98:      &        EPS11, EPS22, EPS12, MRHO, MRHO11, MRHO22, MRHO12, REQ11, REQ22, REQ12, & 98:      &        EPS11, EPS22, EPS12, MRHO, MRHO11, MRHO22, MRHO12, REQ11, REQ22, REQ12, &
 99:      &        ALPHACC, ALPHAHH, ALPHACH, DC6CC, DC6HH, DC6CH, KKJ, CCKJ, KINT, & 99:      &        ALPHACC, ALPHAHH, ALPHACH, DC6CC, DC6HH, DC6CH, KKJ, CCKJ, KINT, &
111:      &        MCPATHGWS, MCPATHGWQ, MCPATHNEGLECT, MCPATHTOL, FRAMESDIFF,TMRATIO, INTMINFAC, MLPLAMBDA, COLL_TOL, KLIM, SCA, &111:      &        MCPATHGWS, MCPATHGWQ, MCPATHNEGLECT, MCPATHTOL, FRAMESDIFF,TMRATIO, INTMINFAC, MLPLAMBDA, COLL_TOL, KLIM, SCA, &
112:      &        NEBMAXERISE, GDSQ, FLATEDIFF, QCIADDREPCUT, QCIADDREPEPS, QCIRADSHIFT, INTCONCUT, MLQLAMBDA, FRQCONV, FRQCONV2, &112:      &        NEBMAXERISE, GDSQ, FLATEDIFF, QCIADDREPCUT, QCIADDREPEPS, QCIRADSHIFT, INTCONCUT, MLQLAMBDA, FRQCONV, FRQCONV2, &
113:      &        TANHFAC, LJADDCUTOFF,LJADDREFNORM,MAXIMFACTOR, PUSHOPTCONV113:      &        TANHFAC, LJADDCUTOFF,LJADDREFNORM,MAXIMFACTOR, PUSHOPTCONV
114: 114: 
115: !     sf344115: !     sf344
116:       DOUBLE PRECISION :: PCUTOFF,PYA11(3),PYA21(3),PYA12(3),PYA22(3),PEPSILON1(3),PSCALEFAC1(2),PSCALEFAC2(2), &116:       DOUBLE PRECISION :: PCUTOFF,PYA11(3),PYA21(3),PYA12(3),PYA22(3),PEPSILON1(3),PSCALEFAC1(2),PSCALEFAC2(2), &
117:      &                     PEPSILONATTR(2),PSIGMAATTR(2), PYOVERLAPTHRESH, LJSITECOORDS(3), LJGSITESIGMA, LJGSITEEPS, &117:      &                     PEPSILONATTR(2),PSIGMAATTR(2), PYOVERLAPTHRESH, LJSITECOORDS(3), LJGSITESIGMA, LJGSITEEPS, &
118:      &                     PYLOCALSTEP(2),PYCFTHRESH,PYGRAVITYC1,PYGRAVITYC2,PERCCUT118:      &                     PYLOCALSTEP(2),PYCFTHRESH,PYGRAVITYC1,PYGRAVITYC2,PERCCUT
119:  119:  
120:       DOUBLE PRECISION, ALLOCATABLE :: PYADDEPS(:,:)120:       DOUBLE PRECISION, ALLOCATABLE :: PYADDEPS(:,:)
121:       DOUBLE PRECISION, ALLOCATABLE :: SHIFTL(:) 
122:       LOGICAL, ALLOCATABLE :: uniaxarray(:)121:       LOGICAL, ALLOCATABLE :: uniaxarray(:)
123:       DOUBLE PRECISION, ALLOCATABLE :: PYADDREP(:,:), PYADDATT(:,:)122:       DOUBLE PRECISION, ALLOCATABLE :: PYADDREP(:,:), PYADDATT(:,:)
124:       DOUBLE PRECISION, ALLOCATABLE :: POINTSDECA(:), POINTSICOS(:)123:       DOUBLE PRECISION, ALLOCATABLE :: POINTSDECA(:), POINTSICOS(:)
125:       DOUBLE PRECISION, ALLOCATABLE :: VT(:), pya1bin(:,:),pya2bin(:,:),SITE(:,:)124:       DOUBLE PRECISION, ALLOCATABLE :: VT(:), pya1bin(:,:),pya2bin(:,:),SITE(:,:)
126:       LOGICAL          :: LJSITE,BLJSITE,LJSITEATTR,PYBINARYT,PARAMONOVPBCX,PARAMONOVPBCY,PARAMONOVPBCZ,PARAMONOVCUTOFF125:       LOGICAL          :: LJSITE,BLJSITE,LJSITEATTR,PYBINARYT,PARAMONOVPBCX,PARAMONOVPBCY,PARAMONOVPBCZ,PARAMONOVCUTOFF
127:       LOGICAL          :: PYGPERIODICT,ELLIPSOIDT,LJSITECOORDST,REALIGNXYZ,MULTISITEPYT,LJGSITET,NORMALMODET,PYT,PYCOLDFUSION126:       LOGICAL          :: PYGPERIODICT,ELLIPSOIDT,LJSITECOORDST,REALIGNXYZ,MULTISITEPYT,LJGSITET,NORMALMODET,PYT,PYCOLDFUSION
128:       INTEGER          :: MAXINTERACTIONS,PYBINARYTYPE1,MYUNIT,NPYSITE127:       INTEGER          :: MAXINTERACTIONS,PYBINARYTYPE1,MYUNIT,NPYSITE
129:       DOUBLE PRECISION :: GBANISOTROPYR, GBWELLDEPTHR128:       DOUBLE PRECISION :: GBANISOTROPYR, GBWELLDEPTHR
130: !     CHRIS129: !     CHRIS
131: !     MCP130: !     MCP


r33147/OPTIM.F 2017-08-09 09:30:13.088936031 +0100 r33146/OPTIM.F 2017-08-09 09:30:13.528941882 +0100
101:       RBATOMSMAX=10101:       RBATOMSMAX=10
102: 102: 
103:       CALL KEYWORDS(Q)103:       CALL KEYWORDS(Q)
104:       ALLOCATE (STPMAX(NOPT))104:       ALLOCATE (STPMAX(NOPT))
105:       IF (UNRST.AND.(RKMIN.OR.BSMIN.OR.(INR.GT.-1))) THEN105:       IF (UNRST.AND.(RKMIN.OR.BSMIN.OR.(INR.GT.-1))) THEN
106:          PRINT '(A)','UNRES not coded for requested optimisation option'106:          PRINT '(A)','UNRES not coded for requested optimisation option'
107:          CALL FLUSH(6)107:          CALL FLUSH(6)
108:          STOP108:          STOP
109:       ENDIF109:       ENDIF
110: 110: 
111:       ALLOCATE(SHIFTL(NOPT)) 
112:       IF (PYGPERIODICT.OR.PYBINARYT) CALL INITIALISEPYGPERIODIC111:       IF (PYGPERIODICT.OR.PYBINARYT) CALL INITIALISEPYGPERIODIC
113:       IF (MULTISITEPYT) CALL DEFINEPYMULTISITES112:       IF (MULTISITEPYT) CALL DEFINEPYMULTISITES
114:       IF (LJGSITET) CALL DEFINELJMULTISITES113:       IF (LJGSITET) CALL DEFINELJMULTISITES
115:       IF (LOCALPERMDIST) THEN114:       IF (LOCALPERMDIST) THEN
116:          ALLOCATE(RBGROUP(RBATOMSMAX))115:          ALLOCATE(RBGROUP(RBATOMSMAX))
117:          ALLOCATE(RBNINGROUP(NATOMS))116:          ALLOCATE(RBNINGROUP(NATOMS))
118:       ENDIF117:       ENDIF
119: 118: 
120:       IF (CHEMSHIFT) CALL CAMSHIFTREAD(20)119:       IF (CHEMSHIFT) CALL CAMSHIFTREAD(20)
121:       IF (MIEFT) CALL MIEF_INI()120:       IF (MIEFT) CALL MIEF_INI()
370:             CALL SYMMETRY(HORDER,.TRUE.,QW,INERTIA)369:             CALL SYMMETRY(HORDER,.TRUE.,QW,INERTIA)
371:             NATOMS=NATOMSSAVE370:             NATOMS=NATOMSSAVE
372:          ELSE371:          ELSE
373:             CALL SYMMETRY(HORDER,.TRUE.,Q,INERTIA)372:             CALL SYMMETRY(HORDER,.TRUE.,Q,INERTIA)
374:          ENDIF373:          ENDIF
375:       ENDIF374:       ENDIF
376:       IF (ADMT) CALL ADM(Q)375:       IF (ADMT) CALL ADM(Q)
377: C376: C
378: C  Eigenvalue shifting377: C  Eigenvalue shifting
379: C378: C
380:       DO J1=1,NOPT379:       DO J1=1,6
381:          SHIFTL(J1)=0.0D0380:          SHIFTL(J1)=0.0D0
382:       ENDDO381:       ENDDO
383: !382: !
384: ! Set the  SHIFTL elements even if NOHESS is .TRUE. 383: ! Set the  SHIFTL elements even if NOHESS is .TRUE. 
385: ! This should do nothing if we never use a Hessian, but allows for a second-order384: ! This should do nothing if we never use a Hessian, but allows for a second-order
386: ! pathway search even if NOHESS is used for the transition state search.385: ! pathway search even if NOHESS is used for the transition state search.
387: !386: !
388:       IF (ZSYM(1).EQ.'CK') THEN387:       IF (ZSYM(1).EQ.'CK') THEN
389:          WRITE(*,'(A,G20.10)') ' OPTIM> Using z rotational ev shift=',SHIFTV388:          WRITE(*,'(A,G20.10)') ' OPTIM> Using z rotational ev shift=',SHIFTV
390:          SHIFTL(6)=SHIFTV389:          SHIFTL(6)=SHIFTV
1033:       IF (ALLOCATED(FIN)) DEALLOCATE(FIN)1032:       IF (ALLOCATED(FIN)) DEALLOCATE(FIN)
1034:       IF (ALLOCATED(START)) DEALLOCATE(START)1033:       IF (ALLOCATED(START)) DEALLOCATE(START)
1035:       IF (ALLOCATED(QW)) DEALLOCATE(QW)1034:       IF (ALLOCATED(QW)) DEALLOCATE(QW)
1036:       IF (ALLOCATED(MIEF_SIG)) DEALLOCATE(MIEF_SIG)1035:       IF (ALLOCATED(MIEF_SIG)) DEALLOCATE(MIEF_SIG)
1037:       IF (ALLOCATED(MIEF_EPS)) DEALLOCATE(MIEF_EPS)1036:       IF (ALLOCATED(MIEF_EPS)) DEALLOCATE(MIEF_EPS)
1038:       IF (ALLOCATED(MIEF_SITES)) DEALLOCATE(MIEF_SITES)1037:       IF (ALLOCATED(MIEF_SITES)) DEALLOCATE(MIEF_SITES)
1039:       IF (ALLOCATED(MIEF_U_RCUT)) DEALLOCATE(MIEF_U_RCUT)1038:       IF (ALLOCATED(MIEF_U_RCUT)) DEALLOCATE(MIEF_U_RCUT)
1040:       IF (ALLOCATED(MIEF_DUDR_RCUT)) DEALLOCATE(MIEF_DUDR_RCUT)1039:       IF (ALLOCATED(MIEF_DUDR_RCUT)) DEALLOCATE(MIEF_DUDR_RCUT)
1041:       IF (ALLOCATED(MLPDAT)) DEALLOCATE(MLPDAT) 1040:       IF (ALLOCATED(MLPDAT)) DEALLOCATE(MLPDAT) 
1042:       IF (ALLOCATED(MLPOUTCOME)) DEALLOCATE(MLPOUTCOME) 1041:       IF (ALLOCATED(MLPOUTCOME)) DEALLOCATE(MLPOUTCOME) 
1043:       IF (ALLOCATED(SHIFTL)) DEALLOCATE(SHIFTL)  
1044:       1042:       
1045:       IF (DUMPDATAT .OR. DUMPDATA_MACHINET) CLOSE(881)1043:       IF (DUMPDATAT .OR. DUMPDATA_MACHINET) CLOSE(881)
1046:       IF (DUMPDATA_MACHINET) CLOSE(882)     1044:       IF (DUMPDATA_MACHINET) CLOSE(882)     
1047:       IF (DUMPMAG) CLOSE(321)1045:       IF (DUMPMAG) CLOSE(321)
1048:       CALL FLUSH(6) 1046:       CALL FLUSH(6) 
1049:       1047:       
1050:       RETURN1048:       RETURN
1051: 1049: 
1052:       END1050:       END
1053: 1051: 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0