hdiff output

r33200/key.f90 2017-08-23 17:30:23.726263199 +0100 r33199/key.f90 2017-08-23 17:30:24.746276040 +0100
 18:      &        NTRAPPOW, MAXINTIMAGE, CHECKDID, CHECKREPINTERVAL, INTFREEZEMIN, INTNTRIESMAX, INTIMAGEINCR, & 18:      &        NTRAPPOW, MAXINTIMAGE, CHECKDID, CHECKREPINTERVAL, INTFREEZEMIN, INTNTRIESMAX, INTIMAGEINCR, &
 19:      &        NCONSTRAINTFIX, INTIMAGECHECK, NREPULSIVEFIX, NRANROT, NENDDUP, LOCALPERMNEIGH, & 19:      &        NCONSTRAINTFIX, INTIMAGECHECK, NREPULSIVEFIX, NRANROT, NENDDUP, LOCALPERMNEIGH, &
 20:      &        LOCALPERMMAXSEP, NONEDAPBC, STRUC, QCHEMESNAO, QCHEMESNMO, QCHEMESNZERO, QCHEMESNELEC, PMPATHINR, & 20:      &        LOCALPERMMAXSEP, NONEDAPBC, STRUC, QCHEMESNAO, QCHEMESNMO, QCHEMESNZERO, QCHEMESNELEC, PMPATHINR, &
 21:      &        MULTISUNIT, MULTIFUNIT,NIMAGEINST,NGLJ,ST_TSSTEP,LANSTEP,NONFREEZE, & 21:      &        MULTISUNIT, MULTIFUNIT,NIMAGEINST,NGLJ,ST_TSSTEP,LANSTEP,NONFREEZE, &
 22:      &        MCPATHBINS,MCPATHEQUIL,MCPATHSTEPS,MCPATHPRTFRQ,MCPATHTS,MCPATHSCHECK,RPHSLICES,RPHQBINS, & 22:      &        MCPATHBINS,MCPATHEQUIL,MCPATHSTEPS,MCPATHPRTFRQ,MCPATHTS,MCPATHSCHECK,RPHSLICES,RPHQBINS, &
 23:      &        ITWIST, JTWIST, KTWIST, LTWIST, MCPATHSTART, MCPATHBLOCK, MCPATHOVER, NCPU, MCPATHDOBLOCK, MCMERGES, MCMERGEQ, & 23:      &        ITWIST, JTWIST, KTWIST, LTWIST, MCPATHSTART, MCPATHBLOCK, MCPATHOVER, NCPU, MCPATHDOBLOCK, MCMERGES, MCMERGEQ, &
 24:      &        MCMERGEI,GAUSSIANCHARGE,GAUSSIANMULTI,ITG03, REDOTS, QCIPERMCHECKINT, & 24:      &        MCMERGEI,GAUSSIANCHARGE,GAUSSIANMULTI,ITG03, REDOTS, QCIPERMCHECKINT, &
 25:      &        MLPIN, MLPSTART, MLPOUT, MLPHIDDEN, MLPDATA, NMLP, N_TO_ALIGN, DJWRBID, STM, NHEXAMERS, & 25:      &        MLPIN, MLPSTART, MLPOUT, MLPHIDDEN, MLPDATA, NMLP, N_TO_ALIGN, DJWRBID, STM, NHEXAMERS, &
 26:      &        MLQIN, MLQSTART, MLQOUT, MLQDATA, NMLQ, & 26:      &        MLQIN, MLQSTART, MLQOUT, MLQDATA, NMLQ, &
 27:      &        QCIADDREP, QCIBONDS, QCISECOND, MAXNACTIVE, QCIIMAGE, NADDTARGET, NUMNN, MULTI_COUNT, MULTI_LAST, MULTI_STEP, & 27:      &        QCIADDREP, QCIBONDS, QCISECOND, MAXNACTIVE, QCIIMAGE, NADDTARGET, NUMNN, MULTI_COUNT, MULTI_LAST, MULTI_STEP, &
 28:      &        NDOF, RECCOUNT, MLPPROBPOS, PUSHOPTMAX, MLPNEIGH, QCICYCN 28:      &        NDOF, RECCOUNT, MLPPROBPOS, PUSHOPTMAX, MLPNEIGH
 29:  29: 
 30:       LOGICAL :: DTEST, MASST, RTEST, EFSTEPST, VECTORST, SUMMARYT, DUMPV, DUMPMAG, FREEZE, FREEZERANGE, GRADSQ, & 30:       LOGICAL :: DTEST, MASST, RTEST, EFSTEPST, VECTORST, SUMMARYT, DUMPV, DUMPMAG, FREEZE, FREEZERANGE, GRADSQ, &
 31:      &        PGRAD, VALUEST, ADMT, BFGSMINT, BFGSTST, CHECKINDEX, TOSI, CONTAINER, & 31:      &        PGRAD, VALUEST, ADMT, BFGSMINT, BFGSTST, CHECKINDEX, TOSI, CONTAINER, &
 32:      &        GAUSSIAN, CADPAC, PRESSURE, FTEST, DCHECK, CP2K, DFTP, CPMD, CPMDC, FREEZERES, DF1T, & 32:      &        GAUSSIAN, CADPAC, PRESSURE, FTEST, DCHECK, CP2K, DFTP, CPMD, CPMDC, FREEZERES, DF1T, &
 33:      &        VARIABLES, FIELDT, OHT, IHT, TDT, D5HT, TWOENDS, PV, FRACTIONAL, BLNT, HYBRIDMINT, & 33:      &        VARIABLES, FIELDT, OHT, IHT, TDT, D5HT, TWOENDS, PV, FRACTIONAL, BLNT, HYBRIDMINT, &
 34:      &        INDEXT, LANCZOST, NOSHIFT, GAMESSUS, GAMESSUK, PVTS, RIGIDBODY, CASTEP, ONETEP, QCHEM, QCHEMES, VASP, & 34:      &        INDEXT, LANCZOST, NOSHIFT, GAMESSUS, GAMESSUK, PVTS, RIGIDBODY, CASTEP, ONETEP, QCHEM, QCHEMES, VASP, &
 35:      &        BFGSSTEP, EFOLSTEP, BULKT, HUPDATE, NOHESS, READV, NOIT, THOMSONT, SIO2T, SIO2C6T, BISECTT, BISECTDEBUG, & 35:      &        BFGSSTEP, EFOLSTEP, BULKT, HUPDATE, NOHESS, READV, NOIT, THOMSONT, SIO2T, SIO2C6T, BISECTT, BISECTDEBUG, &
 36:      &        TOSIC6, TOSIPOL, FIXIMAGE, DFTBT, CHECKCONT, CHECKDT, SHIFTED, READSP, DUMPSP, NOFRQS, & 36:      &        TOSIC6, TOSIPOL, FIXIMAGE, DFTBT, CHECKCONT, CHECKDT, SHIFTED, READSP, DUMPSP, NOFRQS, &
 37:      &        ALLSTEPS, ALLVECTORS, MWVECTORS, WELCH, BINARY, READHESS, MOVIE, NORESET, TWOD, & 37:      &        ALLSTEPS, ALLVECTORS, MWVECTORS, WELCH, BINARY, READHESS, MOVIE, NORESET, TWOD, &
 38:      &        DOUBLET, REOPT, PARALLEL, LINEMIN, FIXD, KEEPINDEX, BSMIN, PRINTPTS, RKMIN, REPELTST,& 38:      &        DOUBLET, REOPT, PARALLEL, LINEMIN, FIXD, KEEPINDEX, BSMIN, PRINTPTS, RKMIN, REPELTST,&
 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, QCICYCLEST 61:      &        MBPOLT, MULTIJOB_MACHINET, DUMPDATA_MACHINET, PLUSSIDET, MINUSSIDET, PUSHOPTT, MLPVB3NNT, GAUSSIAN16
 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
103:      &        REDOK, REDOFRAC, D1INIT, D2INIT, REDOE1, REDOE2, RPBETA, REPCON, PFORCE, &103:      &        REDOK, REDOFRAC, D1INIT, D2INIT, REDOE1, REDOE2, RPBETA, REPCON, PFORCE, &
104:      &        CPCONSTRAINTTOL, CPCONSTRAINTDEL, CPCONSTRAINTREP, CPCONSTRAINREPCUT, CPCONFRAC, &104:      &        CPCONSTRAINTTOL, CPCONSTRAINTDEL, CPCONSTRAINTREP, CPCONSTRAINREPCUT, CPCONFRAC, &
105:      &        INTLJTOL, INTLJDEL, INTLJEPS, IMSEPMIN, IMSEPMAX, TRAPK, MINOVERLAP, &105:      &        INTLJTOL, INTLJDEL, INTLJEPS, IMSEPMIN, IMSEPMAX, TRAPK, MINOVERLAP, &
106:      &        INTFREEZETOL, LOCALPERMCUT, LOCALPERMCUT2, LOCALPERMCUTINC, CHECKREPCUTOFF, CONCUTABS, &106:      &        INTFREEZETOL, LOCALPERMCUT, LOCALPERMCUT2, LOCALPERMCUTINC, CHECKREPCUTOFF, CONCUTABS, &
107:      &        CONCUTFRAC, ENDNUMHESSDELTA, DNEBEFRAC, QCHEMSCALE, KAA, SIGMAAA, QUIPATOMMASS, TEMPERATURE1, &107:      &        CONCUTFRAC, ENDNUMHESSDELTA, DNEBEFRAC, QCHEMSCALE, KAA, SIGMAAA, QUIPATOMMASS, TEMPERATURE1, &
108:      &        DISTORTINST,DELTAINST,MOLPROSCALE,COVER,STTSRMSCONV,LAN_DIST,LANCONV,LANFACTOR, &108:      &        DISTORTINST,DELTAINST,MOLPROSCALE,COVER,STTSRMSCONV,LAN_DIST,LANCONV,LANFACTOR, &
109:      &        STOCKEXP, JPARAM, MCPATHTEMP, MCPATHDMAX, MCPATHSTEP, MCPATHACCRATIO, BIASFAC, &109:      &        STOCKEXP, JPARAM, MCPATHTEMP, MCPATHDMAX, MCPATHSTEP, MCPATHACCRATIO, BIASFAC, &
110:      &        MCADDDEV, MCPATHQMIN, MCPATHQMAX, RPHQMIN, RPHQMAX, RPHTEMP, TWISTF, TWISTREF, MCPATHADDREF, &110:      &        MCADDDEV, MCPATHQMIN, MCPATHQMAX, RPHQMIN, RPHQMAX, RPHTEMP, TWISTF, TWISTREF, MCPATHADDREF, &
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, PUSHOPTCONV, QCICYCDIST113:      &        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(:)121:       DOUBLE PRECISION, ALLOCATABLE :: SHIFTL(:)
122:       LOGICAL, ALLOCATABLE :: uniaxarray(:)122:       LOGICAL, ALLOCATABLE :: uniaxarray(:)
123:       DOUBLE PRECISION, ALLOCATABLE :: PYADDREP(:,:), PYADDATT(:,:)123:       DOUBLE PRECISION, ALLOCATABLE :: PYADDREP(:,:), PYADDATT(:,:)


r33200/keywords.f 2017-08-23 17:30:23.982266421 +0100 r33199/keywords.f 2017-08-23 17:30:24.986279060 +0100
697:          DUMPINTEOS=.FALSE.697:          DUMPINTEOS=.FALSE.
698:          DUMPINTXYZFREQ=100698:          DUMPINTXYZFREQ=100
699:          DUMPINTEOSFREQ=100699:          DUMPINTEOSFREQ=100
700:          KINT=0.0D0700:          KINT=0.0D0
701:          QCIAMBERT=.FALSE.701:          QCIAMBERT=.FALSE.
702:          INTMINT=.FALSE.702:          INTMINT=.FALSE.
703:          INTSPRINGACTIVET=.TRUE.703:          INTSPRINGACTIVET=.TRUE.
704:          INTMINFAC=1.0D0704:          INTMINFAC=1.0D0
705:          QCIRADSHIFTT=.FALSE.705:          QCIRADSHIFTT=.FALSE.
706:          QCIRADSHIFT=1.0D0706:          QCIRADSHIFT=1.0D0
707:          QCICYCLEST=.FALSE. 
708:          QCICYCDIST=0.0 
709:          QCICYCN=100 
710: 707: 
711:          CONPOTT=.FALSE.708:          CONPOTT=.FALSE.
712:          CPCONSTRAINTTOL=0.1D0709:          CPCONSTRAINTTOL=0.1D0
713:          CPCONSTRAINTDEL=1.0D5710:          CPCONSTRAINTDEL=1.0D5
714:          CPCONSTRAINTREP=1.0D0711:          CPCONSTRAINTREP=1.0D0
715:          CPCONSTRAINREPCUT=20.0D0712:          CPCONSTRAINREPCUT=20.0D0
716:          CPCONFRAC=1.0D-4713:          CPCONFRAC=1.0D-4
717:          CPREPSEP=0714:          CPREPSEP=0
718:          CPCONSEP=10000715:          CPCONSEP=10000
719:          CHECKOVERLAPT=.FALSE.716:          CHECKOVERLAPT=.FALSE.
2176: ! 2173: ! 
2177:          ELSE IF ((WORD.EQ.'CONINT').OR.(WORD.EQ.'QCIINT')) THEN2174:          ELSE IF ((WORD.EQ.'CONINT').OR.(WORD.EQ.'QCIINT')) THEN
2178:             CHECKCONINT=.TRUE.2175:             CHECKCONINT=.TRUE.
2179:             IF (NITEMS.GT.1) CALL READF(INTMINFAC)2176:             IF (NITEMS.GT.1) CALL READF(INTMINFAC)
2180:             WRITE(*,'(A,G20.10)') ' keyword> Internal minima terms will be scaled by a factor of ',INTMINFAC2177:             WRITE(*,'(A,G20.10)') ' keyword> Internal minima terms will be scaled by a factor of ',INTMINFAC
2181: !2178: !
2182: ! Maximum active atoms in QCI procedure.2179: ! Maximum active atoms in QCI procedure.
2183: !2180: !
2184:       ELSE IF (WORD.EQ.'QCIMAXACTIVE') THEN2181:       ELSE IF (WORD.EQ.'QCIMAXACTIVE') THEN
2185:          CALL READI(MAXNACTIVE)2182:          CALL READI(MAXNACTIVE)
2186:  
2187:  
2188:       ELSE IF (WORD.EQ.'QCICYCLES') THEN 
2189:          QCICYCLEST=.TRUE. 
2190:          CALL READF(QCICYCDIST) 
2191:          CALL READI(QCICYCN) 
2192:          WRITE(*,'(A,F8.2,A,I8,A)') ' keyword> DNEB for distance <' , QCICYCDIST , ' and after ', QCICYCN,' cycles'    
2193:        
2194: ! 2183: ! 
2195: ! Absolute distance to allow before turning on constraint potential.2184: ! Absolute distance to allow before turning on constraint potential.
2196: ! 2185: ! 
2197:          ELSE IF (WORD.EQ.'CONCUTABS') THEN2186:          ELSE IF (WORD.EQ.'CONCUTABS') THEN
2198:             CONCUTABST=.TRUE.2187:             CONCUTABST=.TRUE.
2199:             CONCUTFRACT=.FALSE.2188:             CONCUTFRACT=.FALSE.
2200:             IF (NITEMS.GT.1) CALL READF(CONCUTABS)2189:             IF (NITEMS.GT.1) CALL READF(CONCUTABS)
2201: ! 2190: ! 
2202: ! Fraction of constraint distance to allow before turning on constraint potential.2191: ! Fraction of constraint distance to allow before turning on constraint potential.
2203: ! 2192: ! 


r33200/newconnect.f90 2017-08-23 17:30:23.394259019 +0100 r33199/newconnect.f90 2017-08-23 17:30:24.526273269 +0100
 27:           USE KEYCONNECT 27:           USE KEYCONNECT
 28:           USE CONNECTUTILS 28:           USE CONNECTUTILS
 29:           USE MODCHARMM 29:           USE MODCHARMM
 30:           USE DECIDEWHATTOCONNECT 30:           USE DECIDEWHATTOCONNECT
 31:           USE TRYCONNECTMODULE 31:           USE TRYCONNECTMODULE
 32:           USE IDOMODULE 32:           USE IDOMODULE
 33:           USE PORFUNCS 33:           USE PORFUNCS
 34:           USE AMHGLOBALS, ONLY : NMRES 34:           USE AMHGLOBALS, ONLY : NMRES
 35:           USE KEY,ONLY : BHDISTTHRESH, BHINTERPT, BHDEBUG, DIJKSTRALOCAL, DUMPDATAT, REOPTIMISEENDPOINTS, & 35:           USE KEY,ONLY : BHDISTTHRESH, BHINTERPT, BHDEBUG, DIJKSTRALOCAL, DUMPDATAT, REOPTIMISEENDPOINTS, &
 36:   &                      AMHT, SEQ, MIN1REDO, MIN2REDO, REDOE1, REDOE2, BULKT, TWOD, PERMDIST, RIGIDBODY, & 36:   &                      AMHT, SEQ, MIN1REDO, MIN2REDO, REDOE1, REDOE2, BULKT, TWOD, PERMDIST, RIGIDBODY, &
 37:   &                      INTCONSTRAINTT, INTLJT, INTTST, INTNTRIESMAX, QCICYCLEST, QCICYCN, QCICYCDIST 37:   &                      INTCONSTRAINTT, INTLJT, INTTST, INTNTRIESMAX
 38:           USE COMMONS,ONLY : PARAM1, PARAM2, PARAM3, ZSYM 38:           USE COMMONS,ONLY : PARAM1, PARAM2, PARAM3, ZSYM
 39:           USE MODNEB,ONLY :  NEWCONNECTT 39:           USE MODNEB,ONLY :  NEWCONNECTT
 40:           USE GENRIGID 40:           USE GENRIGID
 41:           IMPLICIT NONE 41:           IMPLICIT NONE
 42:  42: 
 43:           INTEGER LNOPT 43:           INTEGER LNOPT
 44:           INTEGER,INTENT(IN)              :: NA 44:           INTEGER,INTENT(IN)              :: NA
 45:           DOUBLE PRECISION           :: ENDPOINTSEP,EII,EFF 45:           DOUBLE PRECISION           :: ENDPOINTSEP,EII,EFF
 46:           DOUBLE PRECISION           QQ(LNOPT),FINFIN(LNOPT) 46:           DOUBLE PRECISION           QQ(LNOPT),FINFIN(LNOPT)
 47:           LOGICAL,INTENT(IN)              :: PTEST 47:           LOGICAL,INTENT(IN)              :: PTEST
188:                       NCOUNT=NCOUNT+1188:                       NCOUNT=NCOUNT+1
189:                       READ(99,*) MIN2REDO(3*(NCOUNT-1)+1),MIN2REDO(3*(NCOUNT-1)+2),MIN2REDO(3*(NCOUNT-1)+3)189:                       READ(99,*) MIN2REDO(3*(NCOUNT-1)+1),MIN2REDO(3*(NCOUNT-1)+2),MIN2REDO(3*(NCOUNT-1)+3)
190:                    ENDIF190:                    ENDIF
191:                 ENDDO191:                 ENDDO
192:              ELSE192:              ELSE
193:                 READ(99,*) (MIN2REDO(J1),J1=1,3*NA)193:                 READ(99,*) (MIN2REDO(J1),J1=1,3*NA)
194:              ENDIF194:              ENDIF
195:           ENDIF195:           ENDIF
196: 196: 
197:           DO NCONDONE=1,NCONMAX197:           DO NCONDONE=1,NCONMAX
198:                IF (QCICYCLEST) THEN 
199:                   IF (QCICYCN.LT.NCONDONE) THEN 
200:                      INTCONSTRAINTT=.FALSE. 
201:                      INTLJT=.FALSE. 
202:                      WRITE(*,'(A)') ' newconnect> turned QCI interpolation off' 
203:                   ENDIF 
204:                ENDIF 
205:                IF (REDOPATH) THEN198:                IF (REDOPATH) THEN
206:                   IF (NCONDONE.EQ.(NREDO-1)/2+1) THEN199:                   IF (NCONDONE.EQ.(NREDO-1)/2+1) THEN
207:                      WRITE(*,'(A)') ' newconnect> no more coordinates in redopoints'200:                      WRITE(*,'(A)') ' newconnect> no more coordinates in redopoints'
208:                      REDOPATH=.FALSE.201:                      REDOPATH=.FALSE.
209:                      CLOSE(99)202:                      CLOSE(99)
210:                   ENDIF203:                   ENDIF
211:                ENDIF204:                ENDIF
212:                WRITE(CHR,'(i5)') NConDone205:                WRITE(CHR,'(i5)') NConDone
213:                WRITE(STR,'(3a)') '>>>>>>>>>>>>>>>>>>>>> CONNECT CYCLE ',trim(adjustl(chr)),' >>>>>>>>>>>>>>>>>>>>>'206:                WRITE(STR,'(3a)') '>>>>>>>>>>>>>>>>>>>>> CONNECT CYCLE ',trim(adjustl(chr)),' >>>>>>>>>>>>>>>>>>>>>'
214:                WRITE(CHR,'(i5)') Nmin207:                WRITE(CHR,'(i5)') Nmin
553:                      IF (MI(JS)%DATA%D(JF).EQ.0.0D0) THEN546:                      IF (MI(JS)%DATA%D(JF).EQ.0.0D0) THEN
554:                         PRINT '(A,I8,A,2I8,A,G20.10)',' newconnect> skipping DIJKSTRA pair number ',J1,' minima ', &547:                         PRINT '(A,I8,A,2I8,A,G20.10)',' newconnect> skipping DIJKSTRA pair number ',J1,' minima ', &
555:   &                            JS,JF,' distance is now ',MI(JS)%DATA%D(JF)548:   &                            JS,JF,' distance is now ',MI(JS)%DATA%D(JF)
556:                         CYCLE549:                         CYCLE
557:                      ENDIF550:                      ENDIF
558: !                    WRITE(*,'(A,I5,A,2I5)') ' newconnect> trying DIJKSTRA pair number ',J1,' minima ',js,jf551: !                    WRITE(*,'(A,I5,A,2I5)') ' newconnect> trying DIJKSTRA pair number ',J1,' minima ',js,jf
559: !552: !
560: ! Do count constrained potential interpolations as connection attempts.553: ! Do count constrained potential interpolations as connection attempts.
561: !554: !
562:                      IF (INTCONSTRAINTT.OR.INTLJT) THEN555:                      IF (INTCONSTRAINTT.OR.INTLJT) THEN
563:                         IF (QCICYCLEST.AND.(DIJPAIRDIST(J1).LT.QCICYCDIST)) THEN556:                         MI(JS)%DATA%INTNTRIES(JF)=MI(JS)%DATA%INTNTRIES(JF)+1
564:                            MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1 
565:                            WRITE(*,'(A)') ' newconnect> distance below threshold, use DNEB' 
566:                         ELSE   
567:                            MI(JS)%DATA%INTNTRIES(JF)=MI(JS)%DATA%INTNTRIES(JF)+1 
568:                         ENDIF 
569:                      ELSE557:                      ELSE
570:                         MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1558:                         MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1
571:                      ENDIF559:                      ENDIF
572: !560: !
573: ! Minimum JF becomes image number one and minimum JS becomes image NIMAGE+2561: ! Minimum JF becomes image number one and minimum JS becomes image NIMAGE+2
574: ! Just to be confusing, tryconnect has the arguments as JS,JF, so the minima562: ! Just to be confusing, tryconnect has the arguments as JS,JF, so the minima
575: ! get swapped around in making the DNEB interpolation. This shouldn't matter,563: ! get swapped around in making the DNEB interpolation. This shouldn't matter,
576: ! except that the reconnection procedure for INTCONSTRAINTT/INTLJT below needs to564: ! except that the reconnection procedure for INTCONSTRAINTT/INTLJT below needs to
577: ! know this order!565: ! know this order!
578: !566: !
579:                      IF ((INTCONSTRAINTT.OR.INTLJT).AND.(QCICYCLEST.AND.(DIJPAIRDIST(J1).LT.QCICYCDIST))) THEN567:                      CALL TRYCONNECT(JF,JS,TSREDO,REDOPATH,REDOPATHXYZ,INTCONSTRAINTT,INTLJT)
580:                         CALL TRYCONNECT(JF,JS,TSREDO,REDOPATH,REDOPATHXYZ,.FALSE.,.FALSE.) 
581:                      ELSE 
582:                         CALL TRYCONNECT(JF,JS,TSREDO,REDOPATH,REDOPATHXYZ,INTCONSTRAINTT,INTLJT) 
583:                      ENDIF 
584: !                    IF (INTCONSTRAINTT.OR.INTLJT) THEN568: !                    IF (INTCONSTRAINTT.OR.INTLJT) THEN
585: !                       IF ((MI(JS)%DATA%INTNTRIES(JF).GE.INTNTRIESMAX).AND.(MI(JS)%DATA%D(JF).GT.1.0D-5)) THEN569: !                       IF ((MI(JS)%DATA%INTNTRIES(JF).GE.INTNTRIESMAX).AND.(MI(JS)%DATA%D(JF).GT.1.0D-5)) THEN
586: !                          MI(JS)%DATA%D(JF)=HUGE(1.0D0)570: !                          MI(JS)%DATA%D(JF)=HUGE(1.0D0)
587: !                          MI(JS)%DATA%INTERP(JF)=HUGE(1.0D0)571: !                          MI(JS)%DATA%INTERP(JF)=HUGE(1.0D0)
588: !                       ENDIF572: !                       ENDIF
589: !                    ELSE573: !                    ELSE
590:                         IF ((MI(JS)%DATA%NTRIES(JF).GE.NTRIESMAX).AND.(MI(JS)%DATA%D(JF).GT.1.0D-5)) THEN574:                         IF ((MI(JS)%DATA%NTRIES(JF).GE.NTRIESMAX).AND.(MI(JS)%DATA%D(JF).GT.1.0D-5)) THEN
591:                            MI(JS)%DATA%D(JF)=HUGE(1.0D0)575:                            MI(JS)%DATA%D(JF)=HUGE(1.0D0)
592:                            IF (INTERPCOSTFUNCTION) MI(JS)%DATA%INTERP(JF)=HUGE(1.0D0)576:                            IF (INTERPCOSTFUNCTION) MI(JS)%DATA%INTERP(JF)=HUGE(1.0D0)
593:                         ENDIF577:                         ENDIF


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0