hdiff output

r32620/fetchz.f 2017-05-26 14:30:10.754428333 +0100 r32619/fetchz.f 2017-05-26 14:30:11.214434327 +0100
268:             CALL SDINIT(NOXYGEN,NCHARGE)268:             CALL SDINIT(NOXYGEN,NCHARGE)
269:             NATOMS=NOPT/3269:             NATOMS=NOPT/3
270:          ELSEIF (RPSYSTEM(1:2).EQ.'TT') THEN270:          ELSEIF (RPSYSTEM(1:2).EQ.'TT') THEN
271:             WRITE(*,'(A,G20.10,A)')  '              1/kT a.u.=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))271:             WRITE(*,'(A,G20.10,A)')  '              1/kT a.u.=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))
272:             RPBETA=RPBETA*0.00159360144367 ! convert to kcal/mol units272:             RPBETA=RPBETA*0.00159360144367 ! convert to kcal/mol units
273:             WRITE(*,'(A,G20.10,A)')  '              1/kT (mol/kcal)=',RPBETA273:             WRITE(*,'(A,G20.10,A)')  '              1/kT (mol/kcal)=',RPBETA
274:             WRITE(*,'(A,I0,A)') '              TTM3-F potential for a cluster of ', RPDOF/9, ' water molecules'274:             WRITE(*,'(A,I0,A)') '              TTM3-F potential for a cluster of ', RPDOF/9, ' water molecules'
275:             NATOMS=NOPT/3275:             NATOMS=NOPT/3
276:          ELSEIF (RPSYSTEM(1:2).EQ.'MB') THEN276:          ELSEIF (RPSYSTEM(1:2).EQ.'MB') THEN
277:             WRITE(*,'(A,G20.10,A)')  '              1/kT a.u.=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))277:             WRITE(*,'(A,G20.10,A)')  '              1/kT a.u.=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))
 278:             RPBETA=RPBETA*0.00159360144367 ! convert to kcal/mol units
278:             WRITE(*,'(A,G20.10,A)')  '              1/kT (mol/kcal)=',RPBETA279:             WRITE(*,'(A,G20.10,A)')  '              1/kT (mol/kcal)=',RPBETA
279:             WRITE(*,'(A,I0,A)') '              MBPOL potential for a cluster of ', RPDOF/9, ' water molecules'280:             WRITE(*,'(A,I0,A)') '              MBPOL potential for a cluster of ', RPDOF/9, ' water molecules'
280:             NATOMS=NOPT/3281:             NATOMS=NOPT/3
281:          ELSEIF (RPSYSTEM(1:3).EQ.'MAL') THEN282:          ELSEIF (RPSYSTEM(1:3).EQ.'MAL') THEN
282:             WRITE(*,'(A,G20.10,A)')  '              1/kT a.u.=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))283:             WRITE(*,'(A,G20.10,A)')  '              1/kT a.u.=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))
283:             WRITE(*,'(A,I0,A)') '              malonaldehyde potential '284:             WRITE(*,'(A,I0,A)') '              malonaldehyde potential '
284:             NATOMS=NOPT/3285:             NATOMS=NOPT/3
285:          ELSEIF (RPSYSTEM(1:3).EQ.'MCY') THEN286:          ELSEIF (RPSYSTEM(1:3).EQ.'MCY') THEN
286:             WRITE(*,'(A,G20.10,A)')  '              1/kT (a.u.)=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))287:             WRITE(*,'(A,G20.10,A)')  '              1/kT (a.u.)=',RPBETA,' potential type=' // TRIM(ADJUSTL(RPSYSTEM))
287:             WRITE(*,'(A,I1,A,I0,A)') '              VRT(MCY-5f) flexible water dimer potential'288:             WRITE(*,'(A,I1,A,I0,A)') '              VRT(MCY-5f) flexible water dimer potential'


r32620/potential.f 2017-05-26 14:30:10.986431356 +0100 r32619/potential.f 2017-05-26 14:30:11.442437297 +0100
358:                   DO J1=1,RPIMAGES358:                   DO J1=1,RPIMAGES
359:                      ENERGY=ENERGY+BOWMANPOT(COORDS(RPDOF*(J1-1)+1:RPDOF*J1)) ! Hartrees359:                      ENERGY=ENERGY+BOWMANPOT(COORDS(RPDOF*(J1-1)+1:RPDOF*J1)) ! Hartrees
360:                      IF (GTEST) VNEW(RPDOF*(J1-1)+1:RPDOF*J1)=FINDIFGRAD(COORDS(RPDOF*(J1-1)+1:RPDOF*J1), BOWMANPOT, 1.0D-3, GRAD4T)360:                      IF (GTEST) VNEW(RPDOF*(J1-1)+1:RPDOF*J1)=FINDIFGRAD(COORDS(RPDOF*(J1-1)+1:RPDOF*J1), BOWMANPOT, 1.0D-3, GRAD4T)
361:                      IF (STEST) HESS(RPDOF*(J1-1)+1:RPDOF*J1,RPDOF*(J1-1)+1:RPDOF*J1) =361:                      IF (STEST) HESS(RPDOF*(J1-1)+1:RPDOF*J1,RPDOF*(J1-1)+1:RPDOF*J1) =
362:      &               FINDIFHESS_POT(COORDS(RPDOF*(J1-1)+1:RPDOF*J1),BOWMANPOT,1.0D-3)362:      &               FINDIFHESS_POT(COORDS(RPDOF*(J1-1)+1:RPDOF*J1),BOWMANPOT,1.0D-3)
363:                   ENDDO363:                   ENDDO
364:                ELSEIF (TRIM(ADJUSTL(RPSYSTEM)).EQ.'MB') THEN ! MBPOL364:                ELSEIF (TRIM(ADJUSTL(RPSYSTEM)).EQ.'MB') THEN ! MBPOL
365:                   DUMMY2=1.0D0/(RPBETA/RPIMAGES)**2 ! atomic units (hbar = 1)365:                   DUMMY2=1.0D0/(RPBETA/RPIMAGES)**2 ! atomic units (hbar = 1)
366:                   DO J1=1,RPIMAGES366:                   DO J1=1,RPIMAGES
367:                      CALL MBPOL(COORDS(RPDOF*(J1-1)+1:RPDOF*J1),EDUMMY,VNEW(RPDOF*(J1-1)+1:RPDOF*J1),GTEST)367:                      CALL MBPOL(COORDS(RPDOF*(J1-1)+1:RPDOF*J1),EDUMMY,VNEW(RPDOF*(J1-1)+1:RPDOF*J1),GTEST)
368:                      ENERGY=ENERGY+(EDUMMY*0.00159360144367D0)368:                      ENERGY=ENERGY+EDUMMY
369:                      IF (STEST) THEN369:                      IF (STEST) THEN
370:                         PRINT *, 'no hessian for MBPOL'370:                         PRINT *, 'no hessian for MBPOL'
371:                         STOP371:                         STOP
372:                      ENDIF372:                      ENDIF
373:                   ENDDO373:                   ENDDO
374:                   VNEW(:)= VNEW*0.00159360144367D0374: 
375:                ELSEIF (TRIM(ADJUSTL(RPSYSTEM)).EQ.'MAL') THEN ! malonaldehyde375:                ELSEIF (TRIM(ADJUSTL(RPSYSTEM)).EQ.'MAL') THEN ! malonaldehyde
376:                   DUMMY2=1.0D0/(RPBETA/RPIMAGES)**2 ! atomic units (hbar = 1)376:                   DUMMY2=1.0D0/(RPBETA/RPIMAGES)**2 ! atomic units (hbar = 1)
377:                   DO J1=1,RPIMAGES377:                   DO J1=1,RPIMAGES
378:                      CALL MALPES(9, COORDS(RPDOF*(J1-1)+1:RPDOF*J1),  VNEW(RPDOF*(J1-1)+1:RPDOF*J1), DUMMY, GTEST, SSTEST)378:                      CALL MALPES(9, COORDS(RPDOF*(J1-1)+1:RPDOF*J1),  VNEW(RPDOF*(J1-1)+1:RPDOF*J1), DUMMY, GTEST, SSTEST)
379:                      ENERGY=ENERGY+DUMMY ! Hartrees379:                      ENERGY=ENERGY+DUMMY ! Hartrees
380:                      IF (SSTEST) THEN380:                      IF (SSTEST) THEN
381:                         PRINT *, 'no hessian for TTM3-F'381:                         PRINT *, 'no hessian for TTM3-F'
382:                         STOP382:                         STOP
383:                      ENDIF383:                      ENDIF
384:                   ENDDO384:                   ENDDO


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0