hdiff output

r27373/disconnectionDPS.f90 2017-01-21 10:41:52.608297446 +0000 r27372/disconnectionDPS.f90 2017-01-21 10:41:52.848306657 +0000
 51: ! If this keyword is present, the outcome of a basin analysis at the  51: ! If this keyword is present, the outcome of a basin analysis at the 
 52: ! specified energy is written to a file 'basins'. 52: ! specified energy is written to a file 'basins'.
 53: ! 53: !
 54: ! CENTREGMIN 54: ! CENTREGMIN
 55: ! If this keyword is present, then when a node splits into its daughter 55: ! If this keyword is present, then when a node splits into its daughter
 56: ! nodes, the one containing the global minimum is always placed centrally 56: ! nodes, the one containing the global minimum is always placed centrally
 57: ! (even if other nodes carry more minima). This does not guarantee that 57: ! (even if other nodes carry more minima). This does not guarantee that
 58: ! the global minimum is central in the overall diagram because other 58: ! the global minimum is central in the overall diagram because other
 59: ! nodes may push the one containing the global minimum over to one side. 59: ! nodes may push the one containing the global minimum over to one side.
 60: ! 60: !
 61: ! CHOOSECOLOURS 
 62: ! Use in combination with TRMIN to have specific colours for groups of  
 63: ! minima, rather than a rainbow spectrum. 
 64: !  
 65: ! DUMPNUMBERS 61: ! DUMPNUMBERS
 66: ! If present, a file called node_numbers is written, listing the minima 62: ! If present, a file called node_numbers is written, listing the minima
 67: ! associated with each node in each level. Nodes are listed from left to 63: ! associated with each node in each level. Nodes are listed from left to
 68: ! right within each level. 64: ! right within each level.
 69: ! 65: !
 70: ! DUMPSIZES 66: ! DUMPSIZES
 71: ! If present, a file called node_sizes is written, listing how many minima 67: ! If present, a file called node_sizes is written, listing how many minima
 72: ! are represented by each node in each level. Nodes are listed from left to 68: ! are represented by each node in each level. Nodes are listed from left to
 73: ! right in each level. 69: ! right in each level.
 74: ! 70: !
166: ! would cause some branch ends to appear or disappear rather than change node162: ! would cause some branch ends to appear or disappear rather than change node
167: ! if this were not the case). The NOSPLIT option turns this feature off, so that163: ! if this were not the case). The NOSPLIT option turns this feature off, so that
168: ! if two minima are separated by a barrier lower than the level above their own164: ! if two minima are separated by a barrier lower than the level above their own
169: ! energy, they are grouped together. This option should probably never be used.165: ! energy, they are grouped together. This option should probably never be used.
170: !166: !
171: ! TRMIN <n> <max> <file> <file> ...167: ! TRMIN <n> <max> <file> <file> ...
172: ! Label n different sections of the graph in colour as specified by the 168: ! Label n different sections of the graph in colour as specified by the 
173: ! minima in each file, one file for each section.  169: ! minima in each file, one file for each section.  
174: ! Each file is a list of numbers of minima, 170: ! Each file is a list of numbers of minima, 
175: ! one per line as for PICK. max is the total number of minima, not the number 171: ! one per line as for PICK. max is the total number of minima, not the number 
176: ! in the colour files, currently used for array allocation. 172: ! in the colour files. 
 173: ! currently used for array allocation.
177: ! Colours are chosen automatically to spread over a rainbow spectrum  174: ! Colours are chosen automatically to spread over a rainbow spectrum  
178: ! (from red to purple) in the order the files are specified, unless175: ! (from red to purple) in the order the files are specified but colours can 
179: ! CHOOSECOLOURS is also set, in which case colours are 176: ! be specified individually at both COLOURMARKER in this file. - vkd 
180: ! specified individually (search for "COLOURMARKER" in this file to change the 
181: ! particular colours). - vkd  
182: !177: !
183: ! TRVAL <max> <filename> {<cscale> <weightlower> <cshift>}178: ! TRVAL <max> <filename> {<cscale> <weightlower> <cshift>}
184: ! Colour the graph according to an order parameter value for each 179: ! Colour the graph according to an order parameter value for each 
185: ! minimum. The order parameters are read in from the named file, which should contain 180: ! minimum. The order parameters are read in from the named file, which should contain 
186: ! one line per minimum.  The expected range of the order parameters is [0,1] inclusive.181: ! one line per minimum.  The expected range of the order parameters is [0,1] inclusive.
187: ! max is the total number of minima.182: ! max is the total number of minima.
188: ! Colours are chosen automatically to be evenly distributed with order parameter value 183: ! Colours are chosen automatically to be evenly distributed with order parameter value 
189: ! along the edges of the RBG colour cube: red -> yellow -> green -> cyan ->184: ! along the edges of the RBG colour cube: red -> yellow -> green -> cyan ->
190: ! blue.185: ! blue.
191: ! The optional arguments <cscale> and <cshift> scale and shift order parameters that do186: ! The optional arguments <cscale> and <cshift> scale and shift order parameters that do
510:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.505:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.
511:    LOGICAL :: IDMINT=.FALSE.506:    LOGICAL :: IDMINT=.FALSE.
512:    LOGICAL :: TRMINT=.FALSE.507:    LOGICAL :: TRMINT=.FALSE.
513:    LOGICAL :: TRPRINT=.FALSE.508:    LOGICAL :: TRPRINT=.FALSE.
514:    LOGICAL :: TRVALRANGET=.FALSE.509:    LOGICAL :: TRVALRANGET=.FALSE.
515:    LOGICAL :: TRVALT=.FALSE.510:    LOGICAL :: TRVALT=.FALSE.
516:    LOGICAL :: BASINT=.FALSE.511:    LOGICAL :: BASINT=.FALSE.
517:    LOGICAL :: WEIGHTLOWER=.FALSE.512:    LOGICAL :: WEIGHTLOWER=.FALSE.
518:    LOGICAL :: TRVALSCALET=.FALSE.513:    LOGICAL :: TRVALSCALET=.FALSE.
519:    LOGICAL :: MATCHMIN=.FALSE.514:    LOGICAL :: MATCHMIN=.FALSE.
520:    LOGICAL :: COLOURFRONTT=.FALSE.515:    LOGICAL :: COLOURFRONTT
521:    LOGICAL :: CHOOSECOLOURST=.FALSE. 
522:    INTEGER, ALLOCATABLE :: NCONN(:)516:    INTEGER, ALLOCATABLE :: NCONN(:)
523: 517: 
524: END MODULE VARS518: END MODULE VARS
525: !................................................................................!519: !................................................................................!
526: !MODULE PAGE520: !MODULE PAGE
527: MODULE PAGE521: MODULE PAGE
528: 522: 
529:    IMPLICIT NONE523:    IMPLICIT NONE
530:    INTEGER :: PAGE_X=595,  PAGE_Y=842524:    INTEGER :: PAGE_X=595,  PAGE_Y=842
531:    INTEGER :: MARGIN_X=40, MARGIN_Y=50525:    INTEGER :: MARGIN_X=40, MARGIN_Y=50
1423: 1417: 
1424:             FRAC = (ENERGY-END_E(K))/DELTA_E1418:             FRAC = (ENERGY-END_E(K))/DELTA_E
1425:             Y3 = FRAC*(Y2-Y1) + Y11419:             Y3 = FRAC*(Y2-Y1) + Y1
1426:             X3 = FRAC*(X2-X1) + X11420:             X3 = FRAC*(X2-X1) + X1
1427:             IF (Y3 < Y2) Y3=Y21421:             IF (Y3 < Y2) Y3=Y2
1428:             END_X(END_M(K)) = X31422:             END_X(END_M(K)) = X3
1429:             END_Y(END_M(K)) = Y31423:             END_Y(END_M(K)) = Y3
1430:             IF (TRMINT) THEN ! otherwise MARKNODE is not allocated1424:             IF (TRMINT) THEN ! otherwise MARKNODE is not allocated
1431:                IF (MARKNODE(I,K).GT.0) THEN1425:                IF (MARKNODE(I,K).GT.0) THEN
1432: 1426: 
1433: ! COLOURMARKER - This is one of the places where the colours are specified for minima.1427: ! COLOURMARKER - This is where the colours are specified for minima.
1434: ! To choose your own colours, use the 1428: ! To choose your own colours, comment the following section and use the 
1435: ! IF (MARKNODE(I,K) ==#) statements below, where # is the section number.1429: ! IF (MARKNODE(I,K) ==#) statements below, where # is the section number.
1436: ! Colours should be the same for the nodes below.1430: ! Colours should be the same for the nodes below.
1437: 1431: 
1438:                  IF (.NOT. CHOOSECOLOURST) THEN1432: !                 R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10      
1439:                   R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10      1433: !                 IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'
1440:                   IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'1434: !                 IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 
1441:                   IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 1435: !                 IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'
1442:                   IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'1436: !                 IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'
1443:                   IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'1437: !                 IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'
1444:                   IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'1438: !                 IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'
1445:                   IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'1439:                  IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red
1446:                  ELSE1440:                  IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green
1447:                   IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red1441:                  IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue
1448:                   IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green1442:                  IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magneta
1449:                   IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue1443:                  IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange
1450:                   IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magenta1444:                  IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor' !scarlet
1451:                   IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange1445:                  IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black
1452:                   IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor'1446:                  IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'  
1453:                   IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black1447:                  IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'  
1454:                   IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'   
1455:                   IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'   
1456:                  ENDIF 
1457:                  WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1448:                  WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1458:                  WRITE (20, '(A)') '0 0 0 setrgbcolor'1449:                  WRITE (20, '(A)') '0 0 0 setrgbcolor'
1459:                ELSE1450:                ELSE
1460:                  WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1451:                  WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1461:                ENDIF1452:                ENDIF
1462:             ELSE IF (TRVALT) THEN ! otherwise DPMARKNODE is not allocated1453:             ELSE IF (TRVALT) THEN ! otherwise DPMARKNODE is not allocated
1463:                IF (DPMARKNODE(I,K).GE.0.0D0) THEN1454:                IF (DPMARKNODE(I,K).GE.0.0D0) THEN
1464:                   R2=(CSHIFT+DPMARKNODE(I,K)*CSCALE)*3.0D01455:                   R2=(CSHIFT+DPMARKNODE(I,K)*CSCALE)*3.0D0
1465:                   IF (R2.GT.3.0D0) THEN1456:                   IF (R2.GT.3.0D0) THEN
1466:                      WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k1457:                      WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k
1484:                   WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1475:                   WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1485:                ENDIF1476:                ENDIF
1486:             ELSE1477:             ELSE
1487:                WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1478:                WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1488:             ENDIF1479:             ENDIF
1489:          ELSE1480:          ELSE
1490: 1481: 
1491:             IF (TRMINT) THEN1482:             IF (TRMINT) THEN
1492:                IF (MARKNODE(I,K).GT.0) THEN1483:                IF (MARKNODE(I,K).GT.0) THEN
1493: 1484: 
1494: ! COLOURMARKER - This is one of the places where the colours are specified for nodes.1485: ! COLOURMARKER - This is where the colours are specified for nodes.
1495: ! To choose your own colours, use the 1486: ! To choose your own colours, comment the following section and use the 
1496: ! IF (MARKNODE(I,K) ==#) statements below, where # is the section number.1487: ! IF (MARKNODE(I,K) ==#) statements below, where # is the section number.
1497: ! Colours should be the same for the minima above.1488: ! Colours should be the same for the minima above.
1498:    1489:    
1499:                   IF (.NOT. CHOOSECOLOURST) THEN1490: !                  R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10
1500:                    R2=(MARKNODE(I,K)/FLOAT(NMINTR))*101491: !                  IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'
1501:                    IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'1492: !                  IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 
1502:                    IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 1493: !                  IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'
1503:                    IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'1494: !                  IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'
1504:                    IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'1495: !                  IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'
1505:                    IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'1496: !                  IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'
1506:                    IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'1497:                  IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red
1507:                   ELSE1498:                  IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green
1508:                    IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red1499:                  IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue
1509:                    IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green1500:                  IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magneta
1510:                    IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue1501:                  IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange
1511:                    IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magenta1502:                  IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor' !terracotta
1512:                    IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange1503:                  IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black
1513:                    IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor'1504:                  IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'  
1514:                    IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black1505:                  IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'  
1515:                    IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'   
1516:                    IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'   
1517:                   ENDIF 
1518:  1506:  
1519: !                 IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor'1507: !                 IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor'
1520: !                 IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor'1508: !                 IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor'
1521: !                 IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '1 1 0 setrgbcolor'1509: !                 IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '1 1 0 setrgbcolor'
1522: !                 IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '0.5 1 0 setrgbcolor'1510: !                 IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '0.5 1 0 setrgbcolor'
1523: !                 IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '0 1 0 setrgbcolor'1511: !                 IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '0 1 0 setrgbcolor'
1524: !                 IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '0 1 0.5 setrgbcolor'1512: !                 IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '0 1 0.5 setrgbcolor'
1525: !                 IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 1 1 setrgbcolor'1513: !                 IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 1 1 setrgbcolor'
1526: !                 IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '0 0.5 1 setrgbcolor'1514: !                 IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '0 0.5 1 setrgbcolor'
1527: !                 IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '0 0 1 setrgbcolor'1515: !                 IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '0 0 1 setrgbcolor'
1608: 1596: 
1609:               FRAC = (ENERGY-END_E(K))/DELTA_E1597:               FRAC = (ENERGY-END_E(K))/DELTA_E
1610:               Y3 = FRAC*(Y2-Y1) + Y11598:               Y3 = FRAC*(Y2-Y1) + Y1
1611:               X3 = FRAC*(X2-X1) + X11599:               X3 = FRAC*(X2-X1) + X1
1612:               IF (Y3 < Y2) Y3=Y21600:               IF (Y3 < Y2) Y3=Y2
1613:               END_X(END_M(K)) = X31601:               END_X(END_M(K)) = X3
1614:               END_Y(END_M(K)) = Y31602:               END_Y(END_M(K)) = Y3
1615:               IF (TRMINT) THEN ! otherwise MARKNODE is not allocated1603:               IF (TRMINT) THEN ! otherwise MARKNODE is not allocated
1616:                  IF (MARKNODE(I,K).GT.0) THEN1604:                  IF (MARKNODE(I,K).GT.0) THEN
1617: 1605: 
1618: !   COLOURMARKER - This is one of the places where the colours are specified for minima.1606: !   COLOURMARKER - This is where the colours are specified for minima.
1619: !   To choose your own colours, use the 1607: !   To choose your own colours, comment the following section and use the 
1620: !   IF (MARKNODE(I,K) ==#) statements below, where # is the section number.1608: !   IF (MARKNODE(I,K) ==#) statements below, where # is the section number.
1621: !   Colours should be the same for the nodes below.1609: !   Colours should be the same for the nodes below.
1622: 1610: 
1623:                    IF (.NOT. CHOOSECOLOURST) THEN1611: !                   R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10      
1624:                     R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10      1612: !                   IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'
1625:                     IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'1613: !                   IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 
1626:                     IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 1614: !                   IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'
1627:                     IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'1615: !                   IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'
1628:                     IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'1616: !                   IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'
1629:                     IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'1617: !                   IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'
1630:                     IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'1618:                    IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red
1631:                    ELSE1619:                    IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green
1632:                     IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red1620:                    IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue
1633:                     IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green1621:                    IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magneta
1634:                     IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue1622:                    IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange
1635:                     IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magenta1623:                    IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor' !scarlet
1636:                     IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange1624:                    IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black
1637:                     IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor'1625:                    IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'  
1638:                     IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black1626:                    IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'  
1639:                     IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'   
1640:                     IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'   
1641:                    ENDIF 
1642:                    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1627:                    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1643:                  ENDIF1628:                  ENDIF
1644:               ELSE IF (TRVALT) THEN ! otherwise DPMARKNODE is not allocated1629:               ELSE IF (TRVALT) THEN ! otherwise DPMARKNODE is not allocated
1645:                  IF (DPMARKNODE(I,K).GE.0.0D0) THEN1630:                  IF (DPMARKNODE(I,K).GE.0.0D0) THEN
1646:                     R2=(CSHIFT+DPMARKNODE(I,K)*CSCALE)*3.0D01631:                     R2=(CSHIFT+DPMARKNODE(I,K)*CSCALE)*3.0D0
1647:                     IF (R2.GT.3.0D0) THEN1632:                     IF (R2.GT.3.0D0) THEN
1648:                        WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k1633:                        WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k
1649:                        R2=3.0D01634:                        R2=3.0D0
1650:                     END IF1635:                     END IF
1651:                      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.0d0.AND.R2.GE.0.0d0)  WRITE (20, '(2F7.4,A)') (2.67d0-R2)/2.67D0 ,R2/2.0D0, ' 0  setrgbcolor'
1660: !                    IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'1645: !                    IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'
1661: !                    IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'1646: !                    IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'
1662:                     WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1647:                     WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1663:                  ENDIF1648:                  ENDIF
1664:               ENDIF1649:               ENDIF
1665:            ELSE1650:            ELSE
1666: 1651: 
1667:               IF (TRMINT) THEN1652:               IF (TRMINT) THEN
1668:                  IF (MARKNODE(I,K).GT.0) THEN1653:                  IF (MARKNODE(I,K).GT.0) THEN
1669: 1654: 
1670: !   COLOURMARKER - This is one of the places where the colours are specified for nodes.1655: !   COLOURMARKER - This is where the colours are specified for nodes.
1671: !   To choose your own colours, use the 1656: !   To choose your own colours, comment the following section and use the 
1672: !   IF (MARKNODE(I,K) ==#) statements below, where # is the section number.1657: !   IF (MARKNODE(I,K) ==#) statements below, where # is the section number.
1673: !   Colours should be the same for the minima above.1658: !   Colours should be the same for the minima above.
1674:      1659:      
1675:                    IF (.NOT. CHOOSECOLOURST) THEN1660: !                    R2=(MARKNODE(I,K)/FLOAT(NMINTR))*10
1676:                     R2=(MARKNODE(I,K)/FLOAT(NMINTR))*101661: !                    IF (R2.LE.1.AND.R2.GT.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'
1677:                     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' 
1678:                     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'
1679:                     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'
1680:                     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'
1681:                     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'
1682:                     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
1683:                    ELSE1668:                    IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green
1684:                     IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor' !red1669:                    IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue
1685:                     IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '0 1 0 setrgbcolor'  !green1670:                    IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magneta
1686:                     IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '0 0 1 setrgbcolor'  !blue1671:                    IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange
1687:                     IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '1 0 1 setrgbcolor' !magenta1672:                    IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor' !terracotta
1688:                     IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor' !orange1673:                    IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black
1689:                     IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '1.0 0.2 0.2 setrgbcolor'1674:                    IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'  
1690:                     IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 0 0 setrgbcolor'    !black1675:                    IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'  
1691:                     IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '1 0.5 1 setrgbcolor'   
1692:                     IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '1 0 0.5 setrgbcolor'   
1693:                    ENDIF 
1694:     1676:     
1695: !                   IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor'1677: !                   IF (MARKNODE(I,K) ==1)  WRITE (20, '(A)') '1 0 0  setrgbcolor'
1696: !                   IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor'1678: !                   IF (MARKNODE(I,K) ==2)  WRITE (20, '(A)') '1 0.5 0 setrgbcolor'
1697: !                   IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '1 1 0 setrgbcolor'1679: !                   IF (MARKNODE(I,K) ==3)  WRITE (20, '(A)') '1 1 0 setrgbcolor'
1698: !                   IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '0.5 1 0 setrgbcolor'1680: !                   IF (MARKNODE(I,K) ==4)  WRITE (20, '(A)') '0.5 1 0 setrgbcolor'
1699: !                   IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '0 1 0 setrgbcolor'1681: !                   IF (MARKNODE(I,K) ==5)  WRITE (20, '(A)') '0 1 0 setrgbcolor'
1700: !                   IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '0 1 0.5 setrgbcolor'1682: !                   IF (MARKNODE(I,K) ==6)  WRITE (20, '(A)') '0 1 0.5 setrgbcolor'
1701: !                   IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 1 1 setrgbcolor'1683: !                   IF (MARKNODE(I,K) ==7)  WRITE (20, '(A)') '0 1 1 setrgbcolor'
1702: !                   IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '0 0.5 1 setrgbcolor'1684: !                   IF (MARKNODE(I,K) ==8)  WRITE (20, '(A)') '0 0.5 1 setrgbcolor'
1703: !                   IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '0 0 1 setrgbcolor'1685: !                   IF (MARKNODE(I,K) ==9)  WRITE (20, '(A)') '0 0 1 setrgbcolor'
2086:       CASE ('WEIGHTS')2068:       CASE ('WEIGHTS')
2087:          CALL GET_STRING(FILE_WEIGHTS)2069:          CALL GET_STRING(FILE_WEIGHTS)
2088:          WEIGHTS=.TRUE.2070:          WEIGHTS=.TRUE.
2089:       CASE ('BASINT')2071:       CASE ('BASINT')
2090:          BASINT = .TRUE.2072:          BASINT = .TRUE.
2091:          CALL GET_DP(BENERGY)2073:          CALL GET_DP(BENERGY)
2092:       CASE ('ZEROGM')2074:       CASE ('ZEROGM')
2093:          ZEROGM = .TRUE.2075:          ZEROGM = .TRUE.
2094:       CASE ('COLOURFRONT')2076:       CASE ('COLOURFRONT')
2095:          COLOURFRONTT=.TRUE.2077:          COLOURFRONTT=.TRUE.
2096:       CASE ('CHOOSECOLOURS') 
2097:          CHOOSECOLOURST=.TRUE. 
2098:       CASE DEFAULT2078:       CASE DEFAULT
2099:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)2079:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)
2100:          STOP2080:          STOP
2101:       END SELECT2081:       END SELECT
2102:    END DO2082:    END DO
2103:    CLOSE (20)2083:    CLOSE (20)
2104: 2084: 
2105:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))2085:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))
2106: 2086: 
2107:    IF (DELTA_E <= 0.0D0) THEN2087:    IF (DELTA_E <= 0.0D0) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0