hdiff output

r28776/disconnectionDPS.f90 2017-01-21 10:41:50.988235139 +0000 r28775/disconnectionDPS.f90 2017-01-21 10:41:51.316247771 +0000
1425:         DO R=1,NMINTR1425:         DO R=1,NMINTR
1426:           DO Q=1, MINRANGE1426:           DO Q=1, MINRANGE
1427:               IF (JOINED(R,Q).GT.0) CYCLE1427:               IF (JOINED(R,Q).GT.0) CYCLE
1428:               K2=BASIN(I2,MINPAIRS(R,Q,2))1428:               K2=BASIN(I2,MINPAIRS(R,Q,2))
1429:               K1=BASIN(I2,MINPAIRS(R,Q,1))1429:               K1=BASIN(I2,MINPAIRS(R,Q,1))
1430:               IF (K2.NE.0) THEN1430:               IF (K2.NE.0) THEN
1431: ! R may already map to another integer.1431: ! R may already map to another integer.
1432: ! Check this first and find the maximum integer for mapping1432: ! Check this first and find the maximum integer for mapping
1433:               RMAX=R 1433:               RMAX=R 
1434:               DO1434:               DO
1435: !               PRINT*, 'Rmax, map(Rmax)= ', RMAX, MAP(RMAX) 1435:                 PRINT*, 'Rmax, map(Rmax)= ', RMAX, MAP(RMAX) 
1436:                 IF (MAP(RMAX).EQ.RMAX) THEN1436:                 IF (MAP(RMAX).EQ.RMAX) THEN
1437:                  EXIT1437:                  EXIT
1438:                 ENDIF1438:                 ENDIF
1439:                 RMAX=MAP(RMAX) 1439:                 RMAX=MAP(RMAX) 
1440:               ENDDO   1440:               ENDDO   
1441:               IF (K1.EQ.K2) THEN1441:               IF (K1.EQ.K2) THEN
1442:      WRITE (20, '(A, I6, A, I6, A, I6, A, I6)')'Nodes ', MINPAIRS(R,Q,2), ' and ' , MINPAIRS(R,Q,1), ' BASIN ', K2, ' joins Colour ', r1442:      WRITE (20, '(A, I6, A, I6, A, I6, A, I6)')'Nodes ', MINPAIRS(R,Q,2), ' and ' , MINPAIRS(R,Q,1), ' BASIN ', K2, ' joins Colour ', r
1443: ! Also need to stop colouring this pair.1443: ! Also need to stop colouring this pair.
1444:                 JOINED(R,Q)=R1444:                 JOINED(R,Q)=R
1445:                 CYCLE1445:                 CYCLE
1446:               ENDIF1446:               ENDIF
1447:               IF ((MARKNODE(I2,K2).NE.0).AND.(MARKNODE(I2,K2).NE.R)) THEN1447:               IF ((MARKNODE(I2,K2).NE.0).AND.(MARKNODE(I2,K2).NE.R)) THEN
1448:                NCHAIN=MARKNODE(I2,K2)1448:                NCHAIN=MARKNODE(I2,K2)
1449:                IF (MARKNODE(I2,K2).GT.R) PRINT*, 'Error- marknode'1449:                IF (MARKNODE(I2,K2).GT.R) PRINT*, 'Error- marknode'
1450: ! Now check mapping for previous value of marknode1450: ! Now check mapping for previous value of marknode
1451:                DO1451:                DO
1452: !                PRINT*, 'nchain, map(nchain)= ', NCHAIN, MAP(NCHAIN) 1452:                 PRINT*, 'nchain, map(nchain)= ', NCHAIN, MAP(NCHAIN) 
1453:                 IF (MAP(NCHAIN).EQ.NCHAIN) THEN1453:                 IF (MAP(NCHAIN).EQ.NCHAIN) THEN
1454:                  EXIT1454:                  EXIT
1455:                 ENDIF 1455:                 ENDIF 
1456:                 NCHAIN=MAP(NCHAIN) 1456:                 NCHAIN=MAP(NCHAIN) 
1457:                ENDDO   1457:                ENDDO   
1458: !              PRINT*, 'K2 Basin:',  K2, ' R= ',R, ' Nchain= ', NCHAIN, 'RMAX= ', RMAX 1458:                PRINT*, 'K2 Basin:',  K2, ' R= ',R, ' Nchain= ', NCHAIN, 'RMAX= ', RMAX 
1459:                IF(NCHAIN.GT.RMAX) THEN1459:                IF(NCHAIN.GT.RMAX) THEN
1460:                 MAP(RMAX)=NCHAIN1460:                 MAP(RMAX)=NCHAIN
1461:                 RMAX=NCHAIN 1461:                 RMAX=NCHAIN 
1462:                ELSE   1462:                ELSE   
1463:                 MAP(NCHAIN)=RMAX1463:                 MAP(NCHAIN)=RMAX
1464:                ENDIF1464:                ENDIF
1465:                MARKNODE(I2,K2)=R1465:                MARKNODE(I2,K2)=R
1466:               ELSE 1466:               ELSE 
1467:                MARKNODE(I2,K2)=R1467:                MARKNODE(I2,K2)=R
1468:               ENDIF 1468:               ENDIF 
1469:               WRITE (20, '(A, I6, A, I6, A, I6)')'Node ', MINPAIRS(R,Q,2), ' BASIN ', K2, ' --> Colour ', r1469:               WRITE (20, '(A, I6, A, I6, A, I6)')'Node ', MINPAIRS(R,Q,2), ' BASIN ', K2, ' --> Colour ', r
1470:               ENDIF 1470:               ENDIF 
1471:               IF (K1.NE.0) THEN1471:               IF (K1.NE.0) THEN
1472: ! R may already map to another integer.1472: ! R may already map to another integer.
1473: ! Check this first and find the maximum integer for mapping1473: ! Check this first and find the maximum integer for mapping
1474:               RMAX=R 1474:               RMAX=R 
1475:               DO1475:               DO
1476: !               PRINT*, 'Rmax, map(Rmax)= ', RMAX, MAP(RMAX) 1476:                 PRINT*, 'Rmax, map(Rmax)= ', RMAX, MAP(RMAX) 
1477:                 IF (MAP(RMAX).EQ.RMAX) THEN1477:                 IF (MAP(RMAX).EQ.RMAX) THEN
1478:                  EXIT1478:                  EXIT
1479:                 ENDIF1479:                 ENDIF
1480:                 RMAX=MAP(RMAX) 1480:                 RMAX=MAP(RMAX) 
1481:               ENDDO   1481:               ENDDO   
1482:               IF ((MARKNODE(I2,K1).NE.0).AND.(MARKNODE(I2,K1).NE.R)) THEN 1482:               IF ((MARKNODE(I2,K1).NE.0).AND.(MARKNODE(I2,K1).NE.R)) THEN 
1483:                NCHAIN=MARKNODE(I2,K1)1483:                NCHAIN=MARKNODE(I2,K1)
1484:                IF (MARKNODE(I2,K1).GT.R) PRINT*, 'Error- marknode'1484:                IF (MARKNODE(I2,K1).GT.R) PRINT*, 'Error- marknode'
1485: ! Now check mapping for previous value of marknode1485: ! Now check mapping for previous value of marknode
1486:                DO1486:                DO
1487: !               PRINT*, 'nchain, map(nchain)= ', NCHAIN, MAP(NCHAIN) 1487:                PRINT*, 'nchain, map(nchain)= ', NCHAIN, MAP(NCHAIN) 
1488:                 IF (MAP(NCHAIN).EQ.NCHAIN) THEN1488:                 IF (MAP(NCHAIN).EQ.NCHAIN) THEN
1489:                  EXIT1489:                  EXIT
1490:                 ENDIF 1490:                 ENDIF 
1491:                 NCHAIN=MAP(NCHAIN) 1491:                 NCHAIN=MAP(NCHAIN) 
1492:                ENDDO   1492:                ENDDO   
1493: !                 PRINT*, 'K1 Basin:',  K1, ' R= ',R, ' Nchain= ', NCHAIN, 'RMAX= ', RMAX 1493:                  PRINT*, 'K1 Basin:',  K1, ' R= ',R, ' Nchain= ', NCHAIN, 'RMAX= ', RMAX 
1494:                IF(NCHAIN.GT.RMAX) THEN1494:                IF(NCHAIN.GT.RMAX) THEN
1495:                 MAP(RMAX)=NCHAIN1495:                 MAP(RMAX)=NCHAIN
1496:                ELSE   1496:                ELSE   
1497:                 MAP(NCHAIN)=RMAX1497:                 MAP(NCHAIN)=RMAX
1498:                ENDIF1498:                ENDIF
1499:                 MARKNODE(I2,K1)=R1499:                 MARKNODE(I2,K1)=R
1500:               ELSE 1500:               ELSE 
1501:                 MARKNODE(I2,K1)=R1501:                 MARKNODE(I2,K1)=R
1502:               ENDIF 1502:               ENDIF 
1503:               WRITE (20, '(A, I6, A, I6, A, I6)')'Node ', MINPAIRS(R,Q,1), ' BASIN ', K1, ' --> Colour ', r 1503:               WRITE (20, '(A, I6, A, I6, A, I6)')'Node ', MINPAIRS(R,Q,1), ' BASIN ', K1, ' --> Colour ', r 
2587: 2587: 
2588:    USE PAGE2588:    USE PAGE
2589:    USE VARS2589:    USE VARS
2590:    IMPLICIT NONE2590:    IMPLICIT NONE
2591: 2591: 
2592:    DOUBLE PRECISION :: X1, X2, Y1, Y2, X_POS, Y_POS2592:    DOUBLE PRECISION :: X1, X2, Y1, Y2, X_POS, Y_POS
2593:    DOUBLE PRECISION, PARAMETER :: TICK_WIDTH=3.0D02593:    DOUBLE PRECISION, PARAMETER :: TICK_WIDTH=3.0D0
2594:    INTEGER :: I2594:    INTEGER :: I
2595: 2595: 
2596: !  Definitions.2596: !  Definitions.
2597:    WRITE (20, '(4(A, /), (A,2I5, /), 18(A, /), A)') '%!PS-Adobe-3.0 EPSF-3.0', &2597:    WRITE (20, '(22(A, /), A)') '%!PS-Adobe-3.0 EPSF-3.0', &
2598:       & '%% Creator: disconnection.f90', &2598:       & '%% Creator: disconnection.f90', &
2599:       & '%% Minima: ' // TRIM(file_min), &2599:       & '%% Minima: ' // TRIM(file_min), &
2600:       & '%% Transition states: ' // TRIM(file_ts), &2600:       & '%% Transition states: ' // TRIM(file_ts)
2601:       & '%%BoundingBox: 0 0 ',PAGE_X,PAGE_Y, &2601:    WRITE (20, '(A,2I5)') '%%BoundingBox: 0 0 ',PAGE_X,PAGE_Y
2602:       & '/ls {lineto stroke} def', &2602:    WRITE (20, '(22(A, /), A)')  '/ls {lineto stroke} def', &
2603:       & '/mt {moveto} def', &2603:       & '/mt {moveto} def', &
2604:       & '0.1 setlinewidth', &2604:       & '0.1 setlinewidth', &
2605:       & '/swidth { 30 mul } def', &2605:       & '/swidth { 30 mul } def', &
2606:       & '/dpi { 30 div } def', &2606:       & '/dpi { 30 div } def', &
2607:       & '/Rez 1000 dpi def % calc lines per point for fill', &2607:       & '/Rez 1000 dpi def % calc lines per point for fill', &
2608:       & '/Step 1 Rez div def', &2608:       & '/Step 1 Rez div def', &
2609:       & '/StartRGB { /B1 exch def /S1 exch def /H1 exch def', &2609:       & '/StartRGB { /B1 exch def /S1 exch def /H1 exch def', &
2610:       & '          } def', &2610:       & '          } def', &
2611:       & '/EndRGB { /B2 exch def /S2 exch def /H2 exch def', &2611:       & '/EndRGB { /B2 exch def /S2 exch def /H2 exch def', &
2612:       & '          /Steps s Step div def', &2612:       & '          /Steps s Step div def', &
2613:       & '          /Hstep H2 H1 sub Steps div def', &2613:       & '          /Hstep H2 H1 sub Steps div def', &
2614:       & '          /Sstep S2 S1 sub Steps div def', &2614:       & '          /Sstep S2 S1 sub Steps div def', &
2615:       & '          /Bstep B2 B1 sub Steps div def', &2615:       & '          /Bstep B2 B1 sub Steps div def', &
2616:       & '        } def', &2616:       & '        } def', &
2617:       & '/StepRGB { /H1 H1 Hstep add def', &2617:       & '/StepRGB { /H1 H1 Hstep add def', &
2618:       & '           /S1 S1 Sstep add def', &2618:       & '           /S1 S1 Sstep add def', &
2619:       & '           /B1 B1 Bstep add def', &2619:       & '           /B1 B1 Bstep add def', &
2620:       & '           H1 S1 B1 setrgbcolor', &2620:       & '           H1 S1 B1 setrgbcolor', &
2621:       & '         } def'2621:       & '         } def' 
 2622: 
2622: 2623: 
2623: 2624: 
2624: !  Vertical energy axis.2625: !  Vertical energy axis.
2625:    X1 = MARGIN_X + SCALE_WIDTH - 10.0D02626:    X1 = MARGIN_X + SCALE_WIDTH - 10.0D0
2626:    Y1 = Y_POS(0.0D0)2627:    Y1 = Y_POS(0.0D0)
2627:    Y2 = Y_POS(1.0D0)2628:    Y2 = Y_POS(1.0D0)
2628:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x1, y2, ' ls'2629:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x1, y2, ' ls'
2629: 2630: 
2630: !  Ticks on axis.2631: !  Ticks on axis.
2631:    X2 = X1 + TICK_WIDTH/2.0D02632:    X2 = X1 + TICK_WIDTH/2.0D0


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0