hdiff output

r27359/disconnectionDPS.f90 2017-01-21 10:41:48.404134980 +0000 r27358/disconnectionDPS.f90 2017-01-21 10:41:48.652144632 +0000
204: ! MAXTSBARRIER <threshold> ignore transition states with both barriers above this threshold.204: ! MAXTSBARRIER <threshold> ignore transition states with both barriers above this threshold.
205: !205: !
206: ! WEIGHTS <file>206: ! WEIGHTS <file>
207: ! If present, use weights in <file> to scale the horizontal width. The expected 207: ! If present, use weights in <file> to scale the horizontal width. The expected 
208: ! format of <file> is 208: ! format of <file> is 
209: ! bin number  Vmin   Vmax  ln weight209: ! bin number  Vmin   Vmax  ln weight
210: !210: !
211: ! ZEROGM211: ! ZEROGM
212: ! Shift all energies so that the global minimum is at 0212: ! Shift all energies so that the global minimum is at 0
213: 213: 
214: ! COLOURFRONT 
215: ! put the coloured brances in the front of default black branches. 
216:  
217: !................................................................................!214: !................................................................................!
218: 215: 
219: !!!216: !!!
220: !!! KEYWORDS LIBRARY217: !!! KEYWORDS LIBRARY
221: !!!218: !!!
222: !!! Subroutines: ma readread_line(unit[,logical])219: !!! Subroutines: ma readread_line(unit[,logical])
223: !!!              get_string(string[,logical])220: !!!              get_string(string[,logical])
224: !!!              get_integer(integer[,logical])221: !!!              get_integer(integer[,logical])
225: !!!              get_dp(dp[,logical])222: !!!              get_dp(dp[,logical])
226: !!!              get_logical(lgcl[,logical])223: !!!              get_logical(lgcl[,logical])
505:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.502:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.
506:    LOGICAL :: IDMINT=.FALSE.503:    LOGICAL :: IDMINT=.FALSE.
507:    LOGICAL :: TRMINT=.FALSE.504:    LOGICAL :: TRMINT=.FALSE.
508:    LOGICAL :: TRPRINT=.FALSE.505:    LOGICAL :: TRPRINT=.FALSE.
509:    LOGICAL :: TRVALRANGET=.FALSE.506:    LOGICAL :: TRVALRANGET=.FALSE.
510:    LOGICAL :: TRVALT=.FALSE.507:    LOGICAL :: TRVALT=.FALSE.
511:    LOGICAL :: BASINT=.FALSE.508:    LOGICAL :: BASINT=.FALSE.
512:    LOGICAL :: WEIGHTLOWER=.FALSE.509:    LOGICAL :: WEIGHTLOWER=.FALSE.
513:    LOGICAL :: TRVALSCALET=.FALSE.510:    LOGICAL :: TRVALSCALET=.FALSE.
514:    LOGICAL :: MATCHMIN=.FALSE.511:    LOGICAL :: MATCHMIN=.FALSE.
515:    LOGICAL :: COLOURFRONTT 
516:    INTEGER, ALLOCATABLE :: NCONN(:)512:    INTEGER, ALLOCATABLE :: NCONN(:)
517: 513: 
518: END MODULE VARS514: END MODULE VARS
519: !................................................................................!515: !................................................................................!
520: !MODULE PAGE516: !MODULE PAGE
521: MODULE PAGE517: MODULE PAGE
522: 518: 
523:    IMPLICIT NONE519:    IMPLICIT NONE
524:    INTEGER :: PAGE_X=595,  PAGE_Y=842520:    INTEGER :: PAGE_X=595,  PAGE_Y=842
525:    INTEGER :: MARGIN_X=40, MARGIN_Y=50521:    INTEGER :: MARGIN_X=40, MARGIN_Y=50
1549:                   WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y2, ' ls'1545:                   WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y2, ' ls'
1550:                ENDIF1546:                ENDIF
1551:             ELSE1547:             ELSE
1552:                WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y2, ' ls'1548:                WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y2, ' ls'
1553:             ENDIF1549:             ENDIF
1554:          ENDIF1550:          ENDIF
1555: 1234  ENDDO1551: 1234  ENDDO
1556:    ENDDO1552:    ENDDO
1557:    write(*,*) nn1553:    write(*,*) nn
1558: 1554: 
1559:    IF (COLOURFRONTT) THEN 
1560:      DO I = 2, N_LEVELS 
1561:         CALL ENDPOINTS(I, END_E, END_M) 
1562:         WRITE (20, '(/, A, I2, A, I2, /)') '% Level: ', i-1, ' to ', i 
1563:         IF (NODES(I-1) == 0) EXIT 
1564:         ENERGY = E_HIGH - (I-2)*DELTA_E 
1565:         Y1 = Y_POS(1.0D0*(I-2)/(N_LEVELS-1)) 
1566:         Y2 = Y_POS(1.0D0*(I-1)/(N_LEVELS-1)) 
1567:         DO J = 1, NODES(I) 
1568:            K = ORDER(I, J) 
1569:            X1 = (COL_0(I-1, PARENT(I,K)) + NODE_SIZE(I-1, PARENT(I,K))/2.0 - 1.0) / N_COLS 
1570:            X2 = (COL_0(I, K) + NODE_SIZE(I, K)/2.0 - 1.0) / N_COLS 
1571:            IF (WEIGHTS) THEN 
1572:               X1 = (X1-CENTRESPAN(I-1)) * (LEVELWEIGHTS(I-1)/LEVELWEIGHTS(1)) & 
1573:                  * (1.0D0*COLSPAN(1)/COLSPAN(I-1)) + 0.5 
1574:               X2 = (X2-CENTRESPAN(I)) * (LEVELWEIGHTS(I)/LEVELWEIGHTS(1)) & 
1575:                  * (1.0D0*COLSPAN(1)/COLSPAN(I)) + 0.5 
1576:            ENDIF 
1577: !          PRINT '(A,2I6,4F15.5)','i,j,x1,x2,x_pos1,x_pos2: ',i,j,x1,x2,x_pos(x1),x_pos(x2) 
1578:            X1 = X_POS(X1) 
1579:            X2 = X_POS(X2) 
1580:  
1581: !          IF ( BRANCHES(I,K) > 1) THEN 
1582:            IF ( ALLOCATED(BRANCH_XY) ) THEN 
1583:               PRINT '(A,2I8)','I,K=',I,K 
1584:               PRINT '(A,G20.10)','branch',BRANCH_XY(I,K,1) 
1585:               BRANCH_XY(I,K,1) = X1 
1586:               BRANCH_XY(I,K,2) = Y1 
1587:               BRANCH_XY(I,K,3) = X2 
1588:               BRANCH_XY(I,K,4) = Y2 
1589:            ENDIF 
1590:  
1591:            NN = NN + 1 
1592:         
1593:            IF (BRANCHES(I,K) == 0) THEN 
1594:  
1595:               FRAC = (ENERGY-END_E(K))/DELTA_E 
1596:               Y3 = FRAC*(Y2-Y1) + Y1 
1597:               X3 = FRAC*(X2-X1) + X1 
1598:               IF (Y3 < Y2) Y3=Y2 
1599:               END_X(END_M(K)) = X3 
1600:               END_Y(END_M(K)) = Y3 
1601:               IF (TRMINT) THEN ! otherwise MARKNODE is not allocated 
1602:                  IF (MARKNODE(I,K).GT.0) THEN 
1603:  
1604: !   COLOURMARKER - This is where the colours are specified for minima. 
1605: !   To choose your own colours, comment the following section and use the  
1606: !   IF (MARKNODE(I,K) ==#) statements below, where # is the section number. 
1607: !   Colours should be the same for the nodes below. 
1608:  
1609: !                   R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10       
1610: !                   IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor' 
1611: !                   IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor'  
1612: !                   IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor' 
1613: !                   IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor' 
1614: !                   IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor' 
1615: !                   IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor' 
1616:                    IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red 
1617:                    IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green 
1618:                    IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue 
1619:                    IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magneta 
1620:                    IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange 
1621:                    IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor' !scarlet 
1622:                    IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black 
1623:                    IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'   
1624:                    IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'   
1625:                    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls' 
1626:                    WRITE (20, '(A)') '0 0 0 setrgbcolor' 
1627:                  ENDIF 
1628:               ELSE IF (TRVALT) THEN ! otherwise DPMARKNODE is not allocated 
1629:                  IF (DPMARKNODE(I,K).GE.0.0D0) THEN 
1630:                     R2=(CSHIFT+DPMARKNODE(I,K)*CSCALE)*3.0D0 
1631:                     IF (R2.GT.3.0D0) THEN 
1632:                        WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k 
1633:                        R2=3.0D0 
1634:                     END IF 
1635:                      IF (R2.LT.1.0d0.AND.R2.GE.0.0d0)  WRITE (20, '(2F7.4,A)') (2.67d0-R2)/2.67D0 ,R2/2.0D0, ' 0  setrgbcolor' 
1636:                      IF (R2.LT.1.5d0.AND.R2.GE.1.0d0)  WRITE (20, '(F7.4,A)')  0.625D0-(R2-1.0D0)*1.25D0, ' 0.50  0  setrgbcolor'  
1637:                      IF (R2.LT.2.0d0.AND.R2.GE.1.5d0)  WRITE (20, '(A,F7.4,A)') '0 0.5 ',0.625D0+(R2-2.0d0)*1.25D0,'  setrgbcolor' 
1638:                      IF (R2.LE.3.0d0.AND.R2.GE.2.0d0)  & 
1639:     &         WRITE (20, '(A,2F7.4,A)') '0 ',(3.0d0-R2)/2.0D0, ((2.0D0-R2)-1.67d0)/(-2.67D0), '  setrgbcolor'             
1640: !                    IF (R2.LE.1.AND.R2.GE.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor' 
1641: !                    IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor'  
1642: !                    IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor' 
1643: !                    IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor' 
1644: !                    IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor' 
1645: !                    IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor' 
1646:                     WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls' 
1647:                     WRITE (20, '(A)') '0 0 0 setrgbcolor' 
1648:                  ENDIF 
1649:               ENDIF 
1650:            ELSE 
1651:  
1652:               IF (TRMINT) THEN 
1653:                  IF (MARKNODE(I,K).GT.0) THEN 
1654:  
1655: !   COLOURMARKER - This is where the colours are specified for nodes. 
1656: !   To choose your own colours, comment the following section and use the  
1657: !   IF (MARKNODE(I,K) ==#) statements below, where # is the section number. 
1658: !   Colours should be the same for the minima above. 
1659:       
1660: !                    R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10 
1661: !                    IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor' 
1662: !                    IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor'  
1663: !                    IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor' 
1664: !                    IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor' 
1665: !                    IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor' 
1666: !                    IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor' 
1667:                    IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red 
1668:                    IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green 
1669:                    IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue 
1670:                    IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magneta 
1671:                    IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange 
1672:                    IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor' !terracotta 
1673:                    IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black 
1674:                    IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'   
1675:                    IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'   
1676:      
1677: !                   IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' 
1678: !                   IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' 
1679: !                   IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '1 1 0 setrgbcolor' 
1680: !                   IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '0.5 1 0 setrgbcolor' 
1681: !                   IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '0 1 0 setrgbcolor' 
1682: !                   IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '0 1 0.5 setrgbcolor' 
1683: !                   IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 1 1 setrgbcolor' 
1684: !                   IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '0 0.5 1 setrgbcolor' 
1685: !                   IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '0 0 1 setrgbcolor' 
1686: !                   IF (MARKNODE(I,K) ==10)  WRITE (20, '(A)') '0.5 0 1 setrgbcolor' 
1687:                     WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y2, ' ls' 
1688:                     WRITE (20, '(A)') '0 0 0 setrgbcolor' 
1689:                  ENDIF 
1690:               ELSE IF (TRVALT) THEN 
1691:                  IF (DPMARKNODE(I,K).GE.0.0D0) THEN 
1692:                      R2=(CSHIFT+DPMARKNODE(I,K)*CSCALE)*3.0D0 
1693:                      IF (R2.GT.3.0D0) THEN 
1694:                         WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k 
1695:                         R2=3.0D0 
1696:                      END IF 
1697:                      IF (R2.LT.1.0d0.AND.R2.GE.0.0d0)  WRITE (20, '(2F7.4,A)') (2.67d0-R2)/2.67D0 ,R2/2.0D0, ' 0  setrgbcolor' 
1698:                      IF (R2.LT.1.5d0.AND.R2.GE.1.0d0)  WRITE (20, '(F7.4,A)')  0.625D0-(R2-1.0D0)*1.25D0, ' 0.50  0  setrgbcolor'  
1699:                      IF (R2.LT.2.0d0.AND.R2.GE.1.5d0)  WRITE (20, '(A,F7.4,A)') '0 0.5 ',0.625D0+(R2-2.0d0)*1.25D0,'  setrgbcolor' 
1700:                      IF (R2.LE.3.0d0.AND.R2.GE.2.0d0)  & 
1701:     &                   WRITE (20, '(A,2F7.4,A)') & 
1702:     &                   '0 ',(3.0d0-R2)/2.0D0, ((2.0D0-R2)-1.67d0)/(-2.67D0), '  setrgbcolor' 
1703: !                    IF (R2.LT.1.0d0.AND.R2.GE.0.0d0)  WRITE (20, '(A,F6.4,A)') '1 ',R2,' 0  setrgbcolor' 
1704: !                    IF (R2.LT.2.0d0.AND.R2.GE.1.0d0)  WRITE (20, '(F6.4,A)')  2.0d0-R2,' 1 0  setrgbcolor'  
1705: !                    IF (R2.LT.3.0d0.AND.R2.GE.2.0d0)  WRITE (20, '(A,F6.4,A)') '0 1 ',R2-2.0d0,'  setrgbcolor' 
1706: !                    IF (R2.LE.4.0d0.AND.R2.GE.3.0d0)  WRITE (20, '(A,F6.4,A)') '0 ',4.0d0-R2,' 1  setrgbcolor' 
1707: !                    IF (R2.LE.1.AND.R2.GE.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor' 
1708: !                    IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor'  
1709: !                    IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor' 
1710: !                    IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor' 
1711: !                    IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor' 
1712: !                    IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor' 
1713:                      WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y2, ' ls' 
1714:                      WRITE (20, '(A)') '0 0 0 setrgbcolor' 
1715:                  ENDIF 
1716:               ENDIF 
1717:            ENDIF 
1718:         ENDDO 
1719:      ENDDO 
1720:    ENDIF 
1721:    IF (IDENTIFY) THEN1555:    IF (IDENTIFY) THEN
1722:       WRITE (20, '(A)') '1 0 0 setrgbcolor'1556:       WRITE (20, '(A)') '1 0 0 setrgbcolor'
1723:       WRITE (20, '(A, I2, A)') &1557:       WRITE (20, '(A, I2, A)') &
1724:          & '/Times-Roman findfont ', LABEL_SIZE, ' scalefont setfont'1558:          & '/Times-Roman findfont ', LABEL_SIZE, ' scalefont setfont'
1725:       DO I = 1, N_MIN1559:       DO I = 1, N_MIN
1726:          IF (NCONN(I).LE.NCONNMIN) CYCLE1560:          IF (NCONN(I).LE.NCONNMIN) CYCLE
1727:          IF ( (END_X(I) > 0.0D0).AND.(END_Y(I) > 0.0D0) ) THEN1561:          IF ( (END_X(I) > 0.0D0).AND.(END_Y(I) > 0.0D0) ) THEN
1728:             WRITE (MIN_TRIM, '(I8)') I1562:             WRITE (MIN_TRIM, '(I8)') I
1729:             WRITE (20, '(2F7.2, 3A)') END_X(I), END_Y(I), &1563:             WRITE (20, '(2F7.2, 3A)') END_X(I), END_Y(I), &
1730:                & ' mt (', TRIM(ADJUSTL(MIN_TRIM)), ') show'1564:                & ' mt (', TRIM(ADJUSTL(MIN_TRIM)), ') show'
2068:       CASE ('TS')1902:       CASE ('TS')
2069:          CALL GET_STRING(FILE_TS)1903:          CALL GET_STRING(FILE_TS)
2070:       CASE ('WEIGHTS')1904:       CASE ('WEIGHTS')
2071:          CALL GET_STRING(FILE_WEIGHTS)1905:          CALL GET_STRING(FILE_WEIGHTS)
2072:          WEIGHTS=.TRUE.1906:          WEIGHTS=.TRUE.
2073:       CASE ('BASINT')1907:       CASE ('BASINT')
2074:          BASINT = .TRUE.1908:          BASINT = .TRUE.
2075:          CALL GET_DP(BENERGY)1909:          CALL GET_DP(BENERGY)
2076:       CASE ('ZEROGM')1910:       CASE ('ZEROGM')
2077:          ZEROGM = .TRUE.1911:          ZEROGM = .TRUE.
2078:       CASE ('COLOURFRONT') 
2079:          COLOURFRONTT=.TRUE. 
2080:       CASE DEFAULT1912:       CASE DEFAULT
2081:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)1913:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)
2082:          STOP1914:          STOP
2083:       END SELECT1915:       END SELECT
2084:    END DO1916:    END DO
2085:    CLOSE (20)1917:    CLOSE (20)
2086: 1918: 
2087:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))1919:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))
2088: 1920: 
2089:    IF (DELTA_E <= 0.0D0) THEN1921:    IF (DELTA_E <= 0.0D0) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0