hdiff output

r33143/finalq.f 2017-08-08 12:30:17.850420526 +0100 r33142/finalq.f 2017-08-08 12:30:18.086423681 +0100
 67:  67: 
 68: ! csw34> Adding in clarification of what is being tightly quenched 68: ! csw34> Adding in clarification of what is being tightly quenched
 69:       WRITE(MYUNIT,'(A)') 'Tightly converging the SAVE lowest energy minima found' 69:       WRITE(MYUNIT,'(A)') 'Tightly converging the SAVE lowest energy minima found'
 70:       WRITE(MYUNIT,'(A)') 'NOTE: these may NOT match the other output files - see below for a sorted list of Lowest minima' 70:       WRITE(MYUNIT,'(A)') 'NOTE: these may NOT match the other output files - see below for a sorted list of Lowest minima'
 71:       DO J1=1,NSAVE 71:       DO J1=1,NSAVE
 72:          IF (QMIN(J1).LT.1.0D10) THEN 72:          IF (QMIN(J1).LT.1.0D10) THEN
 73:             NATOMS=QMINNATOMS(J1) 73:             NATOMS=QMINNATOMS(J1)
 74:             DO J2=1,3*NATOMS 74:             DO J2=1,3*NATOMS
 75:                COORDS(J2,NP)=QMINP(J1,J2) 75:                COORDS(J2,NP)=QMINP(J1,J2)
 76:             ENDDO 76:             ENDDO
 77:             ! dj337: update box parameters with saved ones 
 78:             if (boxderivt) box_params(1:6) = boxq(j1,1:6) 
 79:              77:             
 80:             !ds656> 78:             !ds656>
 81:             !write(MYUNIT,*) (QMINT(J1,J2),J2=1,NATOMS) 79:             !write(MYUNIT,*) (QMINT(J1,J2),J2=1,NATOMS)
 82:             IF(NSPECIES(0)>1) THEN 80:             IF(NSPECIES(0)>1) THEN
 83:                CALL SET_ATOMLISTS(QMINT(J1,1:NATOMS),1) 81:                CALL SET_ATOMLISTS(QMINT(J1,1:NATOMS),1)
 84:                CALL SET_LABELS(QMINT(J1,1:NATOMS),NP) 82:                CALL SET_LABELS(QMINT(J1,1:NATOMS),NP)
 85:             ENDIF 83:             ENDIF
 86:             !<ds656 84:             !<ds656
 87:              85:             
 88:             NQ(NP)=NQ(NP)+1 86:             NQ(NP)=NQ(NP)+1
134: !               QMIN(J1)=FULLENERGY132: !               QMIN(J1)=FULLENERGY
135: !            ELSE133: !            ELSE
136: !cl457134: !cl457
137:                QMIN(J1)=POTEL135:                QMIN(J1)=POTEL
138: !            ENDIF136: !            ENDIF
139: 137: 
140:             DO J2=1,3*NATOMS138:             DO J2=1,3*NATOMS
141:                QMINP(J1,J2)=COORDS(J2,NP)139:                QMINP(J1,J2)=COORDS(J2,NP)
142:             ENDDO140:             ENDDO
143: 141: 
144:             ! dj337: save box parameters from end of quench 
145:             if (boxderivt) boxq(j1,1:6) = box_params(1:6) 
146:  
147:             IF (CSMT) THEN142:             IF (CSMT) THEN
148:                CSMAV(1:3*NATOMS)=0.0D0143:                CSMAV(1:3*NATOMS)=0.0D0
149:                DO J2=1,CSMGPINDEX144:                DO J2=1,CSMGPINDEX
150: !145: !
151: ! rotate permuted image to best orientation with CSMPMAT146: ! rotate permuted image to best orientation with CSMPMAT
152: ! apply point group operation J2147: ! apply point group operation J2
153: ! 148: ! 
154:                   DO J3=1,NATOMS149:                   DO J3=1,NATOMS
155:                      XTEMP(3*(J3-1)+1)=CSMPMAT(1,1)*CSMIMAGES(3*NATOMS*(J2-1)+3*(J3-1)+1)150:                      XTEMP(3*(J3-1)+1)=CSMPMAT(1,1)*CSMIMAGES(3*NATOMS*(J2-1)+3*(J3-1)+1)
156:      &                                +CSMPMAT(1,2)*CSMIMAGES(3*NATOMS*(J2-1)+3*(J3-1)+2)151:      &                                +CSMPMAT(1,2)*CSMIMAGES(3*NATOMS*(J2-1)+3*(J3-1)+2)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0