hdiff output

r30963/finalq.f 2016-08-04 16:30:09.327599474 +0100 r30962/finalq.f 2016-08-04 16:30:09.575602770 +0100
 14: C 14: C
 15: C   You should have received a copy of the GNU General Public License 15: C   You should have received a copy of the GNU General Public License
 16: C   along with this program; if not, write to the Free Software 16: C   along with this program; if not, write to the Free Software
 17: C   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA 17: C   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 18: C 18: C
 19:       SUBROUTINE FINALQ 19:       SUBROUTINE FINALQ
 20:       USE COMMONS 20:       USE COMMONS
 21:       USE GENRIGID 21:       USE GENRIGID
 22:       use QMODULE 22:       use QMODULE
 23:       USE MODCHARMM, ONLY: ACESOLV,NCHENCALLS,ACEUPSTEP 23:       USE MODCHARMM, ONLY: ACESOLV,NCHENCALLS,ACEUPSTEP
 24:       USE MODAMBER9, ONLY : SETCHIRAL 
 25:       USE CHIRALITY, ONLY: CIS_TRANS_CHECK, CHIRALITY_CHECK 
 26:       IMPLICIT NONE 24:       IMPLICIT NONE
 27:  25: 
 28:  26: 
 29:       INTEGER J1, J2, ITERATIONS, BRUN,QDONE, J3, NP 27:       INTEGER J1, J2, ITERATIONS, BRUN,QDONE, J3, NP
 30:       DOUBLE PRECISION POTEL, SCREENC(3*NATOMSALLOC), X(3*NATOMS), ENERGY, GRAD(3*NATOMSALLOC), TIME 28:       DOUBLE PRECISION POTEL, SCREENC(3*NATOMSALLOC), X(3*NATOMS), ENERGY, GRAD(3*NATOMSALLOC), TIME
 31:       DOUBLE PRECISION SAVECSMNORM, DUMMY2, DIST2, RMAT(3,3), AVVAL, CSMGRAD(3), XTEMP(1:3*NATOMSALLOC) 29:       DOUBLE PRECISION SAVECSMNORM, DUMMY2, DIST2, RMAT(3,3), AVVAL, CSMGRAD(3), XTEMP(1:3*NATOMSALLOC)
 32:       DOUBLE PRECISION DUMMY(3*NATOMSALLOC), AA(3) 30:       DOUBLE PRECISION DUMMY(3*NATOMSALLOC), AA(3)
 33:  31: 
 34:       COMMON /MYPOT/ POTEL 32:       COMMON /MYPOT/ POTEL
 35:       LOGICAL GUIDECHANGET, GUIDET, CSMDOGUIDET,PASS 33:       LOGICAL GUIDECHANGET, GUIDET, CSMDOGUIDET
 36:       COMMON /GD/ GUIDECHANGET, GUIDET, CSMDOGUIDET 34:       COMMON /GD/ GUIDECHANGET, GUIDET, CSMDOGUIDET
 37:   35:  
 38: C 36: C
 39: C  Make sure the lowest minima are tightly converged and then sort 37: C  Make sure the lowest minima are tightly converged and then sort
 40: C  them just to be on the safe side. 38: C  them just to be on the safe side.
 41: C 39: C
 42:       CSMGUIDET=.FALSE. 40:       CSMGUIDET=.FALSE.
 43: C 41: C
 44: C Turn off guiding potential. 42: C Turn off guiding potential.
 45: C 43: C
 88: !fh301>}}} 86: !fh301>}}}
 89: ! hk286 87: ! hk286
 90:  88: 
 91:             IF (RELAXFQ.OR.HYBRIDMINT) THEN 89:             IF (RELAXFQ.OR.HYBRIDMINT) THEN
 92:                RIGIDINIT = .FALSE. 90:                RIGIDINIT = .FALSE.
 93:             ENDIF 91:             ENDIF
 94: ! hk286 92: ! hk286
 95:             IF(UNFREEZEFINALQ) FROZEN(:)=.FALSE.  ! unfreeze everything before the final quenches 93:             IF(UNFREEZEFINALQ) FROZEN(:)=.FALSE.  ! unfreeze everything before the final quenches
 96:             CALL QUENCH(.TRUE.,NP,ITERATIONS,TIME,BRUN,QDONE,SCREENC) 94:             CALL QUENCH(.TRUE.,NP,ITERATIONS,TIME,BRUN,QDONE,SCREENC)
 97:  95: 
 98:             IF(AMBER12T) THEN 
 99:                IF (CHECKCHIRALITY) THEN 
100:                   PASS=.TRUE. 
101:                   IF (SETCHIRAL) THEN 
102:                      CALL CHIRALITY_CHECK(SCREENC, PASS) 
103:                      IF (.NOT. PASS) THEN 
104:                         WRITE(MYUNIT,*) ' finalq> WARNING: chirality differs from initial structure!' 
105:                      ENDIF 
106:                   ENDIF 
107:                ENDIF 
108:                IF (NOCISTRANS) THEN 
109:                   PASS=.TRUE. 
110:                   CALL CIS_TRANS_CHECK(SCREENC, PASS) 
111:                   IF (.NOT. PASS) THEN 
112:                      WRITE(MYUNIT,*) ' finalq> WARNING: cis/trans differs from initial structure!' 
113:                   END IF 
114:                ENDIF 
115:             ENDIF 
116:  
117: !cl457 96: !cl457
118:             IF (CHEMSHIFT2) THEN 97:             IF (CHEMSHIFT2) THEN
119:                WRITE(MYUNIT,'(A,I6,A,G20.10,A,G20.10,A,I5,A,G15.7,A,F12.2)') 'Final Quench ',NQ(NP),' E=', 98:                WRITE(MYUNIT,'(A,I6,A,G20.10,A,G20.10,A,I5,A,G15.7,A,F12.2)') 'Final Quench ',NQ(NP),' E=',
120:      1            POTEL,' ECAM=', ENERGYCAMSHIFT,' steps=',ITERATIONS,' RMS force=',RMS,' time=',TIME-TSTART 99:      1            POTEL,' ECAM=', ENERGYCAMSHIFT,' steps=',ITERATIONS,' RMS force=',RMS,' time=',TIME-TSTART
121:             ELSE100:             ELSE
122: !cl457101: !cl457
123:                WRITE(MYUNIT,'(A,I6,A,G20.10,A,I5,A,G15.7,A,F12.2)') 'Final Quench ',NQ(NP),' energy=',102:                WRITE(MYUNIT,'(A,I6,A,G20.10,A,I5,A,G15.7,A,F12.2)') 'Final Quench ',NQ(NP),' energy=',
124:      1                POTEL,' steps=',ITERATIONS,' RMS force=',RMS,' time=',TIME-TSTART103:      1                POTEL,' steps=',ITERATIONS,' RMS force=',RMS,' time=',TIME-TSTART
125:             ENDIF104:             ENDIF
126: 105: 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0