hdiff output

r33102/mymylbfgs.f 2017-08-01 15:30:10.381439796 +0100 r33101/mymylbfgs.f 2017-08-01 15:30:10.609442825 +0100
180:          ENDIF180:          ENDIF
181:       ENDIF181:       ENDIF
182:       CALL POTENTIAL(XCOORDS,GRAD,ENERGY,.TRUE.,.FALSE.)182:       CALL POTENTIAL(XCOORDS,GRAD,ENERGY,.TRUE.,.FALSE.)
183: !183: !
184: !  Catch cold fusion for ionic potentials and discard.184: !  Catch cold fusion for ionic potentials and discard.
185: !185: !
186: !  Changed EREAL for cold fusion to 1.0D6 rather than 0.0D0, which could result in steps being accepted186: !  Changed EREAL for cold fusion to 1.0D6 rather than 0.0D0, which could result in steps being accepted
187: !  for systems with positive energies. - khs26 26/11/09187: !  for systems with positive energies. - khs26 26/11/09
188: !188: !
189:       IF ((TOSI.OR.WELCH.OR.RGCL2.OR.AMBER.OR.ARNO.OR.PACHECO.OR.TIP.OR.CHRMMT.OR.AMBERT.OR.SILANET 189:       IF ((TOSI.OR.WELCH.OR.RGCL2.OR.AMBER.OR.ARNO.OR.PACHECO.OR.TIP.OR.CHRMMT.OR.AMBERT.OR.SILANET 
190:      &   .OR.PYGPERIODICT.OR.PYBINARYT.OR.PYT.OR.JMT.OR.DMACRYST.OR.TTM3T.OR.AMBER12T.OR.MBPOLT 190:      &   .OR.PYGPERIODICT.OR.PYBINARYT.OR.PYT.OR.JMT.OR.DMACRYST.OR.TTM3T.OR.AMBER12T.OR.MBPOLT)
191:      &   .OR.CLATHRATET) 
192:      &   .AND.((ENERGY.LT.COLDFUSIONLIMIT).OR.PYCOLDFUSION)) THEN191:      &   .AND.((ENERGY.LT.COLDFUSIONLIMIT).OR.PYCOLDFUSION)) THEN
193:          WRITE(MYUNIT,'(A,G20.10)') 'ENERGY=',ENERGY192:          WRITE(MYUNIT,'(A,G20.10)') 'ENERGY=',ENERGY
194:          WRITE(MYUNIT,'(A,2G20.10)') ' Cold fusion diagnosed - step discarded; energy and threshold=',ENERGY,COLDFUSIONLIMIT193:          WRITE(MYUNIT,'(A,2G20.10)') ' Cold fusion diagnosed - step discarded; energy and threshold=',ENERGY,COLDFUSIONLIMIT
195:          ENERGY=1.0D6194:          ENERGY=1.0D6
196:          POTEL=1.0D6195:          POTEL=1.0D6
197:          RMS=1.0D0196:          RMS=1.0D0
198: !     csw34> set COLDFUSION=.TRUE. so that ATEST=.FALSE. in MC197: !     csw34> set COLDFUSION=.TRUE. so that ATEST=.FALSE. in MC
199:          COLDFUSION=.TRUE.198:          COLDFUSION=.TRUE.
200: !        IF (QUENCHDOS) DEALLOCATE(FRAMES, PE, MODGRAD)199: !        IF (QUENCHDOS) DEALLOCATE(FRAMES, PE, MODGRAD)
201:          RETURN200:          RETURN
753: !         IF (QUENCHDOS) DEALLOCATE(FRAMES, PE, MODGRAD)752: !         IF (QUENCHDOS) DEALLOCATE(FRAMES, PE, MODGRAD)
754:           RETURN753:           RETURN
755:       ENDIF754:       ENDIF
756: !755: !
757: !  Catch cold fusion for ionic potentials and discard.756: !  Catch cold fusion for ionic potentials and discard.
758: !757: !
759: !  Changed EREAL for cold fusion to 1.0D6 rather than 0.0D0, which could result in steps being accepted758: !  Changed EREAL for cold fusion to 1.0D6 rather than 0.0D0, which could result in steps being accepted
760: !  for systems with positive energies. - khs26 26/11/09759: !  for systems with positive energies. - khs26 26/11/09
761: !760: !
762:       IF ((TOSI.OR.WELCH.OR.RGCL2.OR.AMBER.OR.ARNO.OR.PACHECO.OR.TIP.OR.CHRMMT.OR.AMBERT.OR.SILANET761:       IF ((TOSI.OR.WELCH.OR.RGCL2.OR.AMBER.OR.ARNO.OR.PACHECO.OR.TIP.OR.CHRMMT.OR.AMBERT.OR.SILANET
763:      &   .OR.PYGPERIODICT.OR.PYBINARYT.OR.PYT.OR.JMT.OR.DMACRYST.OR.TTM3T.OR.AMBER12T.OR.MBPOLT 762:      &   .OR.PYGPERIODICT.OR.PYBINARYT.OR.PYT.OR.JMT.OR.DMACRYST.OR.TTM3T.OR.AMBER12T.OR.MBPOLT)
764:      &   .OR.CLATHRATET) 
765:      &   .AND.((ENEW.LT.COLDFUSIONLIMIT).OR.PYCOLDFUSION)) THEN763:      &   .AND.((ENEW.LT.COLDFUSIONLIMIT).OR.PYCOLDFUSION)) THEN
766:          WRITE(MYUNIT,'(A,2G20.10)') ' Cold fusion diagnosed - step discarded; energy and threshold=',ENEW,COLDFUSIONLIMIT764:          WRITE(MYUNIT,'(A,2G20.10)') ' Cold fusion diagnosed - step discarded; energy and threshold=',ENEW,COLDFUSIONLIMIT
767:          ENERGY=1.0D6765:          ENERGY=1.0D6
768:          ENEW=1.0D6766:          ENEW=1.0D6
769:          POTEL=1.0D6767:          POTEL=1.0D6
770:          RMS=1.0D0768:          RMS=1.0D0
771: !     csw34> set COLDFUSION=.TRUE. so that ATEST=.FALSE. in MC769: !     csw34> set COLDFUSION=.TRUE. so that ATEST=.FALSE. in MC
772:          COLDFUSION=.TRUE.770:          COLDFUSION=.TRUE.
773: !        IF (QUENCHDOS) DEALLOCATE(FRAMES, PE, MODGRAD)771: !        IF (QUENCHDOS) DEALLOCATE(FRAMES, PE, MODGRAD)
774:          RETURN772:          RETURN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0