hdiff output

r25547/disconnectionDPS.f90 2017-01-21 10:41:47.408096138 +0000 r25546/disconnectionDPS.f90 2017-01-21 10:41:47.648105510 +0000
187: ! <x> and <y> (in postscript units) are optional and overide the default values. - vkd187: ! <x> and <y> (in postscript units) are optional and overide the default values. - vkd
188: !188: !
189: ! TSTHRESH    <threshold> ignore transition states above this threshold.189: ! TSTHRESH    <threshold> ignore transition states above this threshold.
190: ! MAXTSENERGY <threshold> ignore transition states above this threshold.190: ! MAXTSENERGY <threshold> ignore transition states above this threshold.
191: ! MAXTSBARRIER <threshold> ignore transition states with both barriers above this threshold.191: ! MAXTSBARRIER <threshold> ignore transition states with both barriers above this threshold.
192: !192: !
193: ! WEIGHTS <file>193: ! WEIGHTS <file>
194: ! If present, use weights in <file> to scale the horizontal width. The expected 194: ! If present, use weights in <file> to scale the horizontal width. The expected 
195: ! format of <file> is 195: ! format of <file> is 
196: ! bin number  Vmin   Vmax  ln weight196: ! bin number  Vmin   Vmax  ln weight
197: ! 
198: ! ZEROGM 
199: ! Shift all energies so that the global minimum is at 0 
200: 197: 
201: !................................................................................!198: !................................................................................!
202: 199: 
203: !!!200: !!!
204: !!! KEYWORDS LIBRARY201: !!! KEYWORDS LIBRARY
205: !!!202: !!!
206: !!! Subroutines: ma readread_line(unit[,logical])203: !!! Subroutines: ma readread_line(unit[,logical])
207: !!!              get_string(string[,logical])204: !!!              get_string(string[,logical])
208: !!!              get_integer(integer[,logical])205: !!!              get_integer(integer[,logical])
209: !!!              get_dp(dp[,logical])206: !!!              get_dp(dp[,logical])
464:    DOUBLE PRECISION :: DJWSCALE=1.0D0 461:    DOUBLE PRECISION :: DJWSCALE=1.0D0 
465:    DOUBLE PRECISION :: TRMIN=0.0D0462:    DOUBLE PRECISION :: TRMIN=0.0D0
466:    DOUBLE PRECISION :: TRMAX=1.0D0463:    DOUBLE PRECISION :: TRMAX=1.0D0
467:    DOUBLE PRECISION :: BENERGY=1.0D100 464:    DOUBLE PRECISION :: BENERGY=1.0D100 
468:    DOUBLE PRECISION :: TRSCALEX=90.0465:    DOUBLE PRECISION :: TRSCALEX=90.0
469:    DOUBLE PRECISION :: TRSCALEY=50.0466:    DOUBLE PRECISION :: TRSCALEY=50.0
470:    CHARACTER(LEN=120) :: FILE_MIN, FILE_PICK, FILE_TS, FILE_WEIGHTS, &    467:    CHARACTER(LEN=120) :: FILE_MIN, FILE_PICK, FILE_TS, FILE_WEIGHTS, &    
471:       FILE_TRACE468:       FILE_TRACE
472:    CHARACTER(LEN=10) :: LAB_FMT469:    CHARACTER(LEN=10) :: LAB_FMT
473:    LOGICAL :: BARRIERS, CENTRE_GMIN, DUMP_NUMBERS, DUMP_SIZES, EXCLUDEALL, &470:    LOGICAL :: BARRIERS, CENTRE_GMIN, DUMP_NUMBERS, DUMP_SIZES, EXCLUDEALL, &
474:       IDENTIFY, MONOTONIC, SPLIT, WEIGHTS,ZEROGM471:       IDENTIFY, MONOTONIC, SPLIT, WEIGHTS
475:    LOGICAL :: IDENTIFY_NODE=.FALSE.472:    LOGICAL :: IDENTIFY_NODE=.FALSE.
476:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.473:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.
477:    LOGICAL :: IDMINT=.FALSE.474:    LOGICAL :: IDMINT=.FALSE.
478:    LOGICAL :: TRMINT=.FALSE.475:    LOGICAL :: TRMINT=.FALSE.
479:    LOGICAL :: TRPRINT=.FALSE.476:    LOGICAL :: TRPRINT=.FALSE.
480:    LOGICAL :: TRVALRANGET=.FALSE.477:    LOGICAL :: TRVALRANGET=.FALSE.
481:    LOGICAL :: TRVALT=.FALSE.478:    LOGICAL :: TRVALT=.FALSE.
482:    LOGICAL :: BASINT=.FALSE.479:    LOGICAL :: BASINT=.FALSE.
483:    LOGICAL :: WEIGHTLOWER=.FALSE.480:    LOGICAL :: WEIGHTLOWER=.FALSE.
484:    LOGICAL :: TRVALSCALET=.FALSE.481:    LOGICAL :: TRVALSCALET=.FALSE.
582:       WRITE (6, '(/,2A,/)') 'ERROR: Could not open ', TRIM(file_ts)579:       WRITE (6, '(/,2A,/)') 'ERROR: Could not open ', TRIM(file_ts)
583:       STOP580:       STOP
584:    END IF581:    END IF
585:    DO582:    DO
586: !     READ (20, *) ts(i)%n, ts(i)%e, lpafs, pg, h_pg, ts(i)%min1, ts(i)%min2583: !     READ (20, *) ts(i)%n, ts(i)%e, lpafs, pg, h_pg, ts(i)%min1, ts(i)%min2
587:       READ (20, *) TS(I)%E, LPAFS, H_PG, TS(I)%MIN1, TS(I)%MIN2584:       READ (20, *) TS(I)%E, LPAFS, H_PG, TS(I)%MIN1, TS(I)%MIN2
588:       IF (TS(I)%MIN1 /= TS(I)%MIN2) I=I+1585:       IF (TS(I)%MIN1 /= TS(I)%MIN2) I=I+1
589:       IF (I == N_TS+1) EXIT586:       IF (I == N_TS+1) EXIT
590:    END DO587:    END DO
591:    CLOSE (20)588:    CLOSE (20)
592:  
593: !  Set global minimum energy to 0. 
594:    IF (ZEROGM) THEN 
595:       DO I=1,N_MIN 
596:          M(I)=M(I)-E_GMIN 
597:       ENDDO 
598:       DO I=1,N_TS 
599:          TS(I)%E=TS(I)%E-E_GMIN 
600:       ENDDO 
601:    ENDIF 
602:  
603:    WRITE (6, '(I6, 3A)') &589:    WRITE (6, '(I6, 3A)') &
604:       & N_TS, ' non-degenerate paths read in from ', TRIM(file_ts), '.'590:       & N_TS, ' non-degenerate paths read in from ', TRIM(file_ts), '.'
605:    WRITE (6, '(A, 2(F18.10))') 'Highest and lowest transition states: ', &591:    WRITE (6, '(A, 2(F18.10))') 'Highest and lowest transition states: ', &
606:       & B_HIGH, B_LOW592:       & B_HIGH, B_LOW
607: 593: 
608:    PRINT *,'n_min,n_ts=',n_min,n_ts594:    PRINT *,'n_min,n_ts=',n_min,n_ts
609:    ALLOCATE (NCONN(N_MIN),PLUS(N_TS),MINUS(N_TS),NDISTA(N_MIN),TSEN(N_TS))595:    ALLOCATE (NCONN(N_MIN),PLUS(N_TS),MINUS(N_TS),NDISTA(N_MIN),TSEN(N_TS))
610:    DO I=1,N_TS596:    DO I=1,N_TS
611:       TSEN(I)=TS(I)%E597:       TSEN(I)=TS(I)%E
612:       PLUS(I)=TS(I)%MIN1598:       PLUS(I)=TS(I)%MIN1
1655:    DUMP_SIZES = .FALSE.1641:    DUMP_SIZES = .FALSE.
1656:    E_HIGH = 0.0D01642:    E_HIGH = 0.0D0
1657:    EXCLUDEALL = .FALSE.1643:    EXCLUDEALL = .FALSE.
1658:    LOWEST = 01644:    LOWEST = 0
1659:    LAB_FMT = 'F8.2'1645:    LAB_FMT = 'F8.2'
1660:    MONOTONIC = .FALSE.1646:    MONOTONIC = .FALSE.
1661:    N_LEVELS = 01647:    N_LEVELS = 0
1662:    FILE_PICK = ''1648:    FILE_PICK = ''
1663:    FILE_TRACE= ''  1649:    FILE_TRACE= ''  
1664:    SPLIT = .TRUE.1650:    SPLIT = .TRUE.
1665:    ZEROGM = .FALSE. 
1666: 1651: 
1667:     1652:     
1668:    OPEN (UNIT=20, FILE='dinfo', STATUS='OLD', IOSTAT=err)1653:    OPEN (UNIT=20, FILE='dinfo', STATUS='OLD', IOSTAT=err)
1669:    IF (ERR /= 0) THEN1654:    IF (ERR /= 0) THEN
1670:       WRITE (6, '(/,A,/)') 'ERROR: Could not open dinfo file'1655:       WRITE (6, '(/,A,/)') 'ERROR: Could not open dinfo file'
1671:       STOP1656:       STOP
1672:    END IF1657:    END IF
1673:    ALLOCATE(MINIDS(1))1658:    ALLOCATE(MINIDS(1))
1674:    DO1659:    DO
1675:       CALL READ_LINE(20, SUCCESS)1660:       CALL READ_LINE(20, SUCCESS)
1792:          CALL GET_DP(TRMIN)1777:          CALL GET_DP(TRMIN)
1793:          CALL GET_DP(TRMAX)1778:          CALL GET_DP(TRMAX)
1794:       CASE ('TS')1779:       CASE ('TS')
1795:          CALL GET_STRING(FILE_TS)1780:          CALL GET_STRING(FILE_TS)
1796:       CASE ('WEIGHTS')1781:       CASE ('WEIGHTS')
1797:          CALL GET_STRING(FILE_WEIGHTS)1782:          CALL GET_STRING(FILE_WEIGHTS)
1798:          WEIGHTS=.TRUE.1783:          WEIGHTS=.TRUE.
1799:       CASE ('BASINT')1784:       CASE ('BASINT')
1800:          BASINT = .TRUE.1785:          BASINT = .TRUE.
1801:          CALL GET_DP(BENERGY)1786:          CALL GET_DP(BENERGY)
1802:       CASE ('ZEROGM') 
1803:          ZEROGM = .TRUE. 
1804:       CASE DEFAULT1787:       CASE DEFAULT
1805:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)1788:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)
1806:          STOP1789:          STOP
1807:       END SELECT1790:       END SELECT
1808:    END DO1791:    END DO
1809:    CLOSE (20)1792:    CLOSE (20)
1810: 1793: 
1811:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))1794:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))
1812: 1795: 
1813:    IF (DELTA_E <= 0.0D0) THEN1796:    IF (DELTA_E <= 0.0D0) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0