hdiff output

r33217/key.f90 2017-08-24 15:30:08.389844965 +0100 r33216/key.f90 2017-08-24 15:30:09.293857031 +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, QCICYCLEST, & 61:      &        MBPOLT, MULTIJOB_MACHINET, DUMPDATA_MACHINET, PLUSSIDET, MINUSSIDET, PUSHOPTT, MLPVB3NNT, GAUSSIAN16, QCICYCLEST
 62:      &        QCIDNEBT 
 63:  62: 
 64: ! sy349 > for testing the flatpath after dneb 63: ! sy349 > for testing the flatpath after dneb
 65:       !LOGICAL, ALLOCATABLE :: FLATPATHT(:) 64:       !LOGICAL, ALLOCATABLE :: FLATPATHT(:)
 66:       LOGICAL FLATPATHT 65:       LOGICAL FLATPATHT
 67:  66: 
 68: ! 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)
 69:       LOGICAL :: PORE8T = .FALSE. ! add 8th power cylindrical pore to the potential? 68:       LOGICAL :: PORE8T = .FALSE. ! add 8th power cylindrical pore to the potential?
 70:       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)
 71:       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
 72:       LOGICAL :: HARMPOLYT = .FALSE. ! add harmonic bonds between the beads 71:       LOGICAL :: HARMPOLYT = .FALSE. ! add harmonic bonds between the beads


r33217/keywords.f 2017-08-24 15:30:08.625848116 +0100 r33216/keywords.f 2017-08-24 15:30:09.521860074 +0100
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.707:          QCICYCLEST=.FALSE.
708:          QCICYCDIST=0.0708:          QCICYCDIST=0.0
709:          QCICYCN=100709:          QCICYCN=100
710:          QCIDNEBT=.FALSE. 
711: 710: 
712:          CONPOTT=.FALSE.711:          CONPOTT=.FALSE.
713:          CPCONSTRAINTTOL=0.1D0712:          CPCONSTRAINTTOL=0.1D0
714:          CPCONSTRAINTDEL=1.0D5713:          CPCONSTRAINTDEL=1.0D5
715:          CPCONSTRAINTREP=1.0D0714:          CPCONSTRAINTREP=1.0D0
716:          CPCONSTRAINREPCUT=20.0D0715:          CPCONSTRAINREPCUT=20.0D0
717:          CPCONFRAC=1.0D-4716:          CPCONFRAC=1.0D-4
718:          CPREPSEP=0717:          CPREPSEP=0
719:          CPCONSEP=10000718:          CPCONSEP=10000
720:          CHECKOVERLAPT=.FALSE.719:          CHECKOVERLAPT=.FALSE.
2181:             WRITE(*,'(A,G20.10)') ' keyword> Internal minima terms will be scaled by a factor of ',INTMINFAC2180:             WRITE(*,'(A,G20.10)') ' keyword> Internal minima terms will be scaled by a factor of ',INTMINFAC
2182: !2181: !
2183: ! Maximum active atoms in QCI procedure.2182: ! Maximum active atoms in QCI procedure.
2184: !2183: !
2185:       ELSE IF (WORD.EQ.'QCIMAXACTIVE') THEN2184:       ELSE IF (WORD.EQ.'QCIMAXACTIVE') THEN
2186:          CALL READI(MAXNACTIVE)2185:          CALL READI(MAXNACTIVE)
2187: 2186: 
2188: 2187: 
2189:       ELSE IF (WORD.EQ.'QCICYCLES') THEN2188:       ELSE IF (WORD.EQ.'QCICYCLES') THEN
2190:          QCICYCLEST=.TRUE.2189:          QCICYCLEST=.TRUE.
2191:          QCIDNEBT=.FALSE. 
2192:          CALL READF(QCICYCDIST)2190:          CALL READF(QCICYCDIST)
2193:          CALL READI(QCICYCN)2191:          CALL READI(QCICYCN)
2194:          WRITE(*,'(A,F8.2,A,I8,A)') ' keyword> DNEB for distance <' , QCICYCDIST , ' and after ', QCICYCN,' cycles'   2192:          WRITE(*,'(A,F8.2,A,I8,A)') ' keyword> DNEB for distance <' , QCICYCDIST , ' and after ', QCICYCN,' cycles'   
2195:       2193:       
2196: ! 2194: ! 
2197: ! Absolute distance to allow before turning on constraint potential.2195: ! Absolute distance to allow before turning on constraint potential.
2198: ! 2196: ! 
2199:          ELSE IF (WORD.EQ.'CONCUTABS') THEN2197:          ELSE IF (WORD.EQ.'CONCUTABS') THEN
2200:             CONCUTABST=.TRUE.2198:             CONCUTABST=.TRUE.
2201:             CONCUTFRACT=.FALSE.2199:             CONCUTFRACT=.FALSE.


r33217/newconnect.f90 2017-08-24 15:30:07.925838777 +0100 r33216/newconnect.f90 2017-08-24 15:30:08.845851054 +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, QCICYCLEST, QCICYCN, QCICYCDIST
 38:   &                      QCIDNEBT 
 39:           USE COMMONS,ONLY : PARAM1, PARAM2, PARAM3, ZSYM 38:           USE COMMONS,ONLY : PARAM1, PARAM2, PARAM3, ZSYM
 40:           USE MODNEB,ONLY :  NEWCONNECTT 39:           USE MODNEB,ONLY :  NEWCONNECTT
 41:           USE GENRIGID 40:           USE GENRIGID
 42:           IMPLICIT NONE 41:           IMPLICIT NONE
 43:  42: 
 44:           INTEGER LNOPT 43:           INTEGER LNOPT
 45:           INTEGER,INTENT(IN)              :: NA 44:           INTEGER,INTENT(IN)              :: NA
 46:           DOUBLE PRECISION           :: ENDPOINTSEP,EII,EFF 45:           DOUBLE PRECISION           :: ENDPOINTSEP,EII,EFF
 47:           DOUBLE PRECISION           QQ(LNOPT),FINFIN(LNOPT) 46:           DOUBLE PRECISION           QQ(LNOPT),FINFIN(LNOPT)
 48:           LOGICAL,INTENT(IN)              :: PTEST 47:           LOGICAL,INTENT(IN)              :: PTEST
193:              ELSE192:              ELSE
194:                 READ(99,*) (MIN2REDO(J1),J1=1,3*NA)193:                 READ(99,*) (MIN2REDO(J1),J1=1,3*NA)
195:              ENDIF194:              ENDIF
196:           ENDIF195:           ENDIF
197: 196: 
198:           DO NCONDONE=1,NCONMAX197:           DO NCONDONE=1,NCONMAX
199:                IF (QCICYCLEST) THEN198:                IF (QCICYCLEST) THEN
200:                   IF (QCICYCN.LT.NCONDONE) THEN199:                   IF (QCICYCN.LT.NCONDONE) THEN
201:                      INTCONSTRAINTT=.FALSE.200:                      INTCONSTRAINTT=.FALSE.
202:                      INTLJT=.FALSE.201:                      INTLJT=.FALSE.
203:                      QCIDNEBT=.TRUE. 
204:                      WRITE(*,'(A)') ' newconnect> turned QCI interpolation off'202:                      WRITE(*,'(A)') ' newconnect> turned QCI interpolation off'
205:                   ENDIF203:                   ENDIF
206:                ENDIF204:                ENDIF
207:                IF (REDOPATH) THEN205:                IF (REDOPATH) THEN
208:                   IF (NCONDONE.EQ.(NREDO-1)/2+1) THEN206:                   IF (NCONDONE.EQ.(NREDO-1)/2+1) THEN
209:                      WRITE(*,'(A)') ' newconnect> no more coordinates in redopoints'207:                      WRITE(*,'(A)') ' newconnect> no more coordinates in redopoints'
210:                      REDOPATH=.FALSE.208:                      REDOPATH=.FALSE.
211:                      CLOSE(99)209:                      CLOSE(99)
212:                   ENDIF210:                   ENDIF
213:                ENDIF211:                ENDIF
558:                         CYCLE556:                         CYCLE
559:                      ENDIF557:                      ENDIF
560: !                    WRITE(*,'(A,I5,A,2I5)') ' newconnect> trying DIJKSTRA pair number ',J1,' minima ',js,jf558: !                    WRITE(*,'(A,I5,A,2I5)') ' newconnect> trying DIJKSTRA pair number ',J1,' minima ',js,jf
561: !559: !
562: ! Do count constrained potential interpolations as connection attempts.560: ! Do count constrained potential interpolations as connection attempts.
563: !561: !
564:                      IF (INTCONSTRAINTT.OR.INTLJT) THEN562:                      IF (INTCONSTRAINTT.OR.INTLJT) THEN
565:                         IF (QCICYCLEST.AND.(DIJPAIRDIST(J1).LT.QCICYCDIST)) THEN563:                         IF (QCICYCLEST.AND.(DIJPAIRDIST(J1).LT.QCICYCDIST)) THEN
566:                            MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1564:                            MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1
567:                            WRITE(*,'(A)') ' newconnect> distance below threshold, use DNEB'565:                            WRITE(*,'(A)') ' newconnect> distance below threshold, use DNEB'
568:                            QCIDNEBT=.TRUE. 
569:                         ELSE  566:                         ELSE  
570:                            MI(JS)%DATA%INTNTRIES(JF)=MI(JS)%DATA%INTNTRIES(JF)+1567:                            MI(JS)%DATA%INTNTRIES(JF)=MI(JS)%DATA%INTNTRIES(JF)+1
571:                            QCIDNEBT=.FALSE. 
572:                         ENDIF568:                         ENDIF
573:                      ELSE569:                      ELSE
574:                         MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1570:                         MI(JS)%DATA%NTRIES(JF)=MI(JS)%DATA%NTRIES(JF)+1
575:                      ENDIF571:                      ENDIF
576: !572: !
577: ! Minimum JF becomes image number one and minimum JS becomes image NIMAGE+2573: ! Minimum JF becomes image number one and minimum JS becomes image NIMAGE+2
578: ! Just to be confusing, tryconnect has the arguments as JS,JF, so the minima574: ! Just to be confusing, tryconnect has the arguments as JS,JF, so the minima
579: ! get swapped around in making the DNEB interpolation. This shouldn't matter,575: ! get swapped around in making the DNEB interpolation. This shouldn't matter,
580: ! except that the reconnection procedure for INTCONSTRAINTT/INTLJT below needs to576: ! except that the reconnection procedure for INTCONSTRAINTT/INTLJT below needs to
581: ! know this order!577: ! know this order!


r33217/newneb.f90 2017-08-24 15:30:08.145841710 +0100 r33216/newneb.f90 2017-08-24 15:30:09.069854042 +0100
 29:           USE MINIMISER3 29:           USE MINIMISER3
 30:           USE NEBOUTPUT 30:           USE NEBOUTPUT
 31:           USE NEBUTILS 31:           USE NEBUTILS
 32:           USE KEY, ONLY : UNRST, GROWSTRINGT, FREEZENODEST, DESMDEBUG, & 32:           USE KEY, ONLY : UNRST, GROWSTRINGT, FREEZENODEST, DESMDEBUG, &
 33:                & NEBMUPDATE, MUPDATE, BFGSSTEPS, NEBRESEEDT, & 33:                & NEBMUPDATE, MUPDATE, BFGSSTEPS, NEBRESEEDT, &
 34:                & INTCONMAX, ORDERI, ORDERJ, EPSALPHA, REDOBFGSSTEPS, &  34:                & INTCONMAX, ORDERI, ORDERJ, EPSALPHA, REDOBFGSSTEPS, & 
 35:                & NREPMAX, DISTREF, NEBKINITIAL, ADDREPT, REPPOW, REDOTSIM, MIN1REDO, MIN2REDO, PUSHOFF, & 35:                & NREPMAX, DISTREF, NEBKINITIAL, ADDREPT, REPPOW, REDOTSIM, MIN1REDO, MIN2REDO, PUSHOFF, &
 36:                & CONI, CONJ, AMHT, NUMGLY, REPI, REPJ, BULKT, D1INIT, D2INIT, & 36:                & CONI, CONJ, AMHT, NUMGLY, REPI, REPJ, BULKT, D1INIT, D2INIT, &
 37:                & REDOKADD, REDOPATH1, INTCONSTRAINTT, INTNEBIMAGES, & 37:                & REDOKADD, REDOPATH1, INTCONSTRAINTT, INTNEBIMAGES, &
 38:                & REDOPATH2, NREPI, NREPJ, REPCUT, NREPCUT, TWOD, RIGIDBODY, PERMDIST, WHOLEDNEB, & 38:                & REDOPATH2, NREPI, NREPJ, REPCUT, NREPCUT, TWOD, RIGIDBODY, PERMDIST, WHOLEDNEB, &
 39:                & CPPNEBT, VARIABLES, SLERPT, QCIXYZ, QCIDNEBT 39:                & CPPNEBT, VARIABLES, SLERPT, QCIXYZ
 40:           USE GROWSTRINGUTILS, ONLY: GROWSTRING, TOTSTEPS 40:           USE GROWSTRINGUTILS, ONLY: GROWSTRING, TOTSTEPS
 41:           USE GSDATA, ONLY : KEYGSPRINT 41:           USE GSDATA, ONLY : KEYGSPRINT
 42:           USE MODGUESS,ONLY: GUESSPATHT,NINTERP 42:           USE MODGUESS,ONLY: GUESSPATHT,NINTERP
 43:           USE MODMEC,ONLY: MECCANOT           43:           USE MODMEC,ONLY: MECCANOT          
 44:           USE INTCOMMONS, ONLY : DESMINT, INTINTERPT, NINTIM, NDIH, DIHINFO, ALIGNDIR, PREVDIH, NINTC 44:           USE INTCOMMONS, ONLY : DESMINT, INTINTERPT, NINTIM, NDIH, DIHINFO, ALIGNDIR, PREVDIH, NINTC
 45:           USE INTCUTILS, ONLY : INTINTERPOLATE, CART2INT 45:           USE INTCUTILS, ONLY : INTINTERPOLATE, CART2INT
 46:           USE SPFUNCTS, ONLY : DUMPCOORDS 46:           USE SPFUNCTS, ONLY : DUMPCOORDS
 47:           USE NEBTOCONNECT 47:           USE NEBTOCONNECT
 48:           USE AMHGLOBALS, ONLY : NMRES 48:           USE AMHGLOBALS, ONLY : NMRES
 49:           USE COMMONS,ONLY: PARAM1,PARAM2,PARAM3,REDOPATHNEB,ZSYM,DEBUG 49:           USE COMMONS,ONLY: PARAM1,PARAM2,PARAM3,REDOPATHNEB,ZSYM,DEBUG
324:                    DEALLOCATE(DELTAX,VNEW,LCOORDS)324:                    DEALLOCATE(DELTAX,VNEW,LCOORDS)
325:                    KNOWE=.FALSE.; KNOWG=.FALSE.; KNOWH=.FALSE.325:                    KNOWE=.FALSE.; KNOWG=.FALSE.; KNOWH=.FALSE.
326: !                  REDOKADD=.FALSE.326: !                  REDOKADD=.FALSE.
327:                    REDOPATH2=.FALSE.327:                    REDOPATH2=.FALSE.
328:                 ELSEIF (REDOPATH) THEN328:                 ELSEIF (REDOPATH) THEN
329:                    XYZ(NOPT+1:NOPT*2) = TSREDO(1:NOPT)329:                    XYZ(NOPT+1:NOPT*2) = TSREDO(1:NOPT)
330:                    PRINT '(A)','newneb> Setting tangent vector and hence initial eigenvector guess to difference between minima'330:                    PRINT '(A)','newneb> Setting tangent vector and hence initial eigenvector guess to difference between minima'
331:                    TANVEC(1:NOPT,1)=XYZ(1:NOPT)-XYZ(NOPT*2+1:NOPT*3)331:                    TANVEC(1:NOPT,1)=XYZ(1:NOPT)-XYZ(NOPT*2+1:NOPT*3)
332:                    EEE(2)=1.0D100332:                    EEE(2)=1.0D100
333:                 ELSE333:                 ELSE
334:                    IF (INTCONSTRAINTT.AND.(.NOT.(QCIDNEBT))) THEN334:                    IF (INTCONSTRAINTT) THEN
335:                       ! QCI method to construct initial band.335:                       ! QCI method to construct initial band.
336: 336: 
337: !                     XYZ(1:NOPT)=QQ(1:NOPT)337: !                     XYZ(1:NOPT)=QQ(1:NOPT)
338: !                     XYZ(NOPT+1:NOPT*(NIMAGE+1))=INTNEBIMAGES(1:NOPT*NIMAGE)338: !                     XYZ(NOPT+1:NOPT*(NIMAGE+1))=INTNEBIMAGES(1:NOPT*NIMAGE)
339:                       IF (.NOT.ALLOCATED(QCIXYZ)) THEN339:                       IF (.NOT.ALLOCATED(QCIXYZ)) THEN
340:                          PRINT *,'newneb> ERROR *** QCIXYZ is not allocated'340:                          PRINT *,'newneb> ERROR *** QCIXYZ is not allocated'
341:                          STOP341:                          STOP
342:                       ENDIF342:                       ENDIF
343:                       IF (MOREPRINTING) CALL KEYNEBPRINT343:                       IF (MOREPRINTING) CALL KEYNEBPRINT
344:                       XYZ(1:NOPT*(NIMAGE+2))=QCIXYZ(1:NOPT*(NIMAGE+2))344:                       XYZ(1:NOPT*(NIMAGE+2))=QCIXYZ(1:NOPT*(NIMAGE+2))


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0