hdiff output

r31591/path.f 2016-12-06 11:30:12.205495235 +0000 r31590/path.f 2016-12-06 11:30:12.461498557 +0000
1893: C        CALL H2OMODES(NATOMS,IPOT,Q,DIAG) ! WCOMMENT1893: C        CALL H2OMODES(NATOMS,IPOT,Q,DIAG) ! WCOMMENT
1894:          CALL H2OMODES(NATOMS/2,IPOT,Q,DIAG)1894:          CALL H2OMODES(NATOMS/2,IPOT,Q,DIAG)
1895: C        WRITE(88,'(3G20.10)') (DIAG(J2),J2=1,6*NATOMS) ! WCOMMENT1895: C        WRITE(88,'(3G20.10)') (DIAG(J2),J2=1,6*NATOMS) ! WCOMMENT
1896:          IF (.NOT.NOFRQS) WRITE(88,'(3G20.10)') (DIAG(J2),J2=1,3*NATOMS)1896:          IF (.NOT.NOFRQS) WRITE(88,'(3G20.10)') (DIAG(J2),J2=1,3*NATOMS)
1897: 1897: 
1898:       ! sn402: For the small number of potentials coded in the fully-rigid form (RBAAT), we can use the NRMLMD subroutine1898:       ! sn402: For the small number of potentials coded in the fully-rigid form (RBAAT), we can use the NRMLMD subroutine
1899:       ! directly to calculate the Hessian and calculate its eigenvalues. We write squared angular frequencies in whatever1899:       ! directly to calculate the Hessian and calculate its eigenvalues. We write squared angular frequencies in whatever
1900:       ! time unit is specified by FRQCONV (see documentation and comments in keywords.f)1900:       ! time unit is specified by FRQCONV (see documentation and comments in keywords.f)
1901: ! hk286 - compute potential for normal modes, notice the toggle1901: ! hk286 - compute potential for normal modes, notice the toggle
1902:       ELSE IF (RBAAT) THEN1902:       ELSE IF (RBAAT) THEN
1903:          CALL POTENTIAL(Q,E,VNEW,.TRUE.,.TRUE.,RMS,.FALSE.,.FALSE.) 
1904:          WRITE(88, '(G20.10)') E 
1905:          ! The following line should definitely not be hard-coded in! 
1906:          WRITE(88, '(A,A)') "1  ", "C1" ! TEMP Solution 
1907:  
1908:          IF (.NOT.NOFRQS) THEN1903:          IF (.NOT.NOFRQS) THEN
1909:             RBAANORMALMODET = .TRUE.           1904:             RBAANORMALMODET = .TRUE.
 1905:             CALL POTENTIAL(Q,E,VNEW,.TRUE.,.TRUE.,RMS,.FALSE.,.FALSE.)
1910:             CALL NRMLMD (Q, DIAG, .FALSE.)1906:             CALL NRMLMD (Q, DIAG, .FALSE.)
1911:             RBAANORMALMODET = .FALSE.1907:             RBAANORMALMODET = .FALSE.
 1908:             WRITE(88, '(G20.10)') E
 1909:             ! The following line should definitely not be hard-coded in!
 1910:             WRITE(88, '(A,A)') "1  ", "C1" ! TEMP Solution
1912:             WRITE(88,'(3G20.10)') (DIAG(J2)*FRQCONV2,J2=1,3*NATOMS)  ! sn402: Added the frequency conversion1911:             WRITE(88,'(3G20.10)') (DIAG(J2)*FRQCONV2,J2=1,3*NATOMS)  ! sn402: Added the frequency conversion
1913:          ENDIF1912:          ENDIF
1914: ! hk2861913: ! hk286
1915: 1914: 
1916:       !!!!! This is the start of the main block to calculate normal mode frequencies. !!!!!1915:       !!!!! This is the start of the main block to calculate normal mode frequencies. !!!!!
1917:       !!!!! A similar procedure is followed for all remaining potentials. !!!!!1916:       !!!!! A similar procedure is followed for all remaining potentials. !!!!!
1918: 1917: 
1919:       ELSE1918:       ELSE
1920: 1919: 
1921:          IF (CHRMMT .OR. AMBERT .OR. AMBER12T .OR. NABT .OR. AMHT) THEN1920:          IF (CHRMMT .OR. AMBERT .OR. AMBER12T .OR. NABT .OR. AMHT) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0