hdiff output

r29342/disconnectionDPS.f90 2015-11-17 23:34:19.725573938 +0000 r29341/disconnectionDPS.f90 2015-11-17 23:34:19.933576721 +0000
 34: ! 34: !
 35: ! FIRST <E1> 35: ! FIRST <E1>
 36: ! Specifies the energy of the highest level on the energy axis. 36: ! Specifies the energy of the highest level on the energy axis.
 37: ! 37: !
 38: ! LEVELS <n> 38: ! LEVELS <n>
 39: ! The number of levels at which to perform the basin analysis. 39: ! The number of levels at which to perform the basin analysis.
 40: ! 40: !
 41: ! MINIMA <file> 41: ! MINIMA <file>
 42: ! Specifies filename for minima info. 42: ! Specifies filename for minima info.
 43: ! 43: !
 44: ! SCALEBAR <ndelta> 
 45: ! specifies a scale bar instead of a vertical axis with height ndelta times delta 
 46: ! 
 47: ! TS <file> 44: ! TS <file>
 48: ! Specifies filename for transition state info. 45: ! Specifies filename for transition state info.
 49: ! 46: !
 50: ! 47: !
 51: ! Optional Keywords 48: ! Optional Keywords
 52: ! ----------------- 49: ! -----------------
 53: ! BASINT <energy> 50: ! BASINT <energy>
 54: ! 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 
 55: ! specified energy is written to a file 'basins'. 52: ! specified energy is written to a file 'basins'.
 56: ! 53: !
561:    DOUBLE PRECISION :: TRMIN=0.0D0558:    DOUBLE PRECISION :: TRMIN=0.0D0
562:    DOUBLE PRECISION :: TRMAX=1.0D0559:    DOUBLE PRECISION :: TRMAX=1.0D0
563:    DOUBLE PRECISION :: BENERGY=1.0D100 560:    DOUBLE PRECISION :: BENERGY=1.0D100 
564:    DOUBLE PRECISION :: TRSCALEX=90.0561:    DOUBLE PRECISION :: TRSCALEX=90.0
565:    DOUBLE PRECISION :: TRSCALEY=50.0562:    DOUBLE PRECISION :: TRSCALEY=50.0
566:    CHARACTER(LEN=120) :: FILE_MIN, FILE_PICK, FILE_TS, FILE_WEIGHTS, &    563:    CHARACTER(LEN=120) :: FILE_MIN, FILE_PICK, FILE_TS, FILE_WEIGHTS, &    
567:       FILE_TRACE, FILE_IDMIN, FILE_PAIRS564:       FILE_TRACE, FILE_IDMIN, FILE_PAIRS
568:    CHARACTER(LEN=10) :: LAB_FMT565:    CHARACTER(LEN=10) :: LAB_FMT
569:    LOGICAL :: BARRIERS, CENTRE_GMIN, DUMP_NUMBERS, DUMP_SIZES, EXCLUDEALL, &566:    LOGICAL :: BARRIERS, CENTRE_GMIN, DUMP_NUMBERS, DUMP_SIZES, EXCLUDEALL, &
570:       IDENTIFY, MONOTONIC, SPLIT, WEIGHTS,ZEROGM567:       IDENTIFY, MONOTONIC, SPLIT, WEIGHTS,ZEROGM
571:    INTEGER :: NDELTASB 
572:    LOGICAL :: SCALEBAR 
573:    LOGICAL :: ORDER_BY_SIZE568:    LOGICAL :: ORDER_BY_SIZE
574:    LOGICAL :: ORDER_BY_ENERGY569:    LOGICAL :: ORDER_BY_ENERGY
575:    LOGICAL :: IDENTIFY_NODE=.FALSE.570:    LOGICAL :: IDENTIFY_NODE=.FALSE.
576:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.571:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.
577:    LOGICAL :: IDMINT=.FALSE.572:    LOGICAL :: IDMINT=.FALSE.
578:    LOGICAL :: TRMINT=.FALSE.573:    LOGICAL :: TRMINT=.FALSE.
579:    LOGICAL :: TRPRINT=.FALSE.574:    LOGICAL :: TRPRINT=.FALSE.
580:    LOGICAL :: TRVALRANGET=.FALSE.575:    LOGICAL :: TRVALRANGET=.FALSE.
581:    LOGICAL :: TRVALT=.FALSE.576:    LOGICAL :: TRVALT=.FALSE.
582:    LOGICAL :: PAIRST=.FALSE.577:    LOGICAL :: PAIRST=.FALSE.
2211:      & 'stack', &2206:      & 'stack', &
2212:      & '0 Step s {', &2207:      & '0 Step s {', &
2213:      & 'StepRGB', &2208:      & 'StepRGB', &
2214:      & 'newpath', &2209:      & 'newpath', &
2215:      & '0.0 swidth exch moveto % position left hand side', &2210:      & '0.0 swidth exch moveto % position left hand side', &
2216:      & '1.0 swidth 0 rlineto % draw a line', &2211:      & '1.0 swidth 0 rlineto % draw a line', &
2217:      & 'stroke', &2212:      & 'stroke', &
2218:      & '} for', &2213:      & '} for', &
2219:      & '%', & 2214:      & '%', & 
2220:      & '%Draw Outline to Box', &2215:      & '%Draw Outline to Box', &
2221:      & '2.0 setlinewidth', &2216:      & '1.0 setlinewidth', &
2222:      & '0.0 0.0 0.0 setrgbcolor', &2217:      & '0.0 0.0 0.0 setrgbcolor', &
2223:      & '0.0 swidth -2.0 swidth translate % origin to bottom right', &2218:      & '0.0 swidth -2.0 swidth translate % origin to bottom right', &
2224:      & '1 swidth 3 swidth moveto', &2219:      & '1 swidth 3 swidth moveto', &
2225:      & '1 swidth 0 swidth lineto', &2220:      & '1 swidth 0 swidth lineto', &
2226:      & '0 swidth 0 swidth lineto', &2221:      & '0 swidth 0 swidth lineto', &
2227:      & '0 swidth 3 swidth lineto', &2222:      & '0 swidth 3 swidth lineto', &
2228:      & 'closepath', &2223:      & 'closepath', &
2229:      & 'stroke', &2224:      & 'stroke', &
2230:      & '0.7 swidth 0.75 swidth moveto', &2225:      & '0.7 swidth 0.75 swidth moveto', &
2231:      & '1.0 swidth 0.75 swidth lineto', &2226:      & '1.0 swidth 0.75 swidth lineto', &
2413:       CASE ('PICK')2408:       CASE ('PICK')
2414:          CALL GET_STRING(FILE_PICK)2409:          CALL GET_STRING(FILE_PICK)
2415:       CASE ('PS_PAGE_X')2410:       CASE ('PS_PAGE_X')
2416:          CALL GET_INTEGER(PAGE_X)2411:          CALL GET_INTEGER(PAGE_X)
2417:       CASE ('PS_PAGE_Y')2412:       CASE ('PS_PAGE_Y')
2418:          CALL GET_INTEGER(PAGE_Y)2413:          CALL GET_INTEGER(PAGE_Y)
2419:       CASE ('ORDER_BY_ENERGY')2414:       CASE ('ORDER_BY_ENERGY')
2420:          ORDER_BY_ENERGY = .TRUE.2415:          ORDER_BY_ENERGY = .TRUE.
2421:       CASE ('ORDER_BY_SIZE')2416:       CASE ('ORDER_BY_SIZE')
2422:          ORDER_BY_SIZE = .TRUE.2417:          ORDER_BY_SIZE = .TRUE.
2423:       CASE ('SCALEBAR') 
2424:          SCALEBAR = .TRUE. 
2425:          CALL GET_INTEGER(NDELTASB) 
2426:       CASE ('TRMIN')2418:       CASE ('TRMIN')
2427:          TRMINT = .TRUE.2419:          TRMINT = .TRUE.
2428:          CALL GET_INTEGER(NDUMMY)2420:          CALL GET_INTEGER(NDUMMY)
2429:          CALL GET_INTEGER(MINRANGE)2421:          CALL GET_INTEGER(MINRANGE)
2430:          IF (NDUMMY.eq.1) THEN2422:          IF (NDUMMY.eq.1) THEN
2431:          WRITE (6, '(I6, A)') ndummy, ' set of minima will be traced in colour'2423:          WRITE (6, '(I6, A)') ndummy, ' set of minima will be traced in colour'
2432:          ELSE2424:          ELSE
2433:          WRITE (6, '(I6, A)') ndummy, ' sets of minima will be traced in colour'2425:          WRITE (6, '(I6, A)') ndummy, ' sets of minima will be traced in colour'
2434:          ENDIF2426:          ENDIF
2435:          ALLOCATE(MINTRS(1:MINRANGE))2427:          ALLOCATE(MINTRS(1:MINRANGE))
2604:    INTEGER :: I2596:    INTEGER :: I
2605: 2597: 
2606: !  Definitions.2598: !  Definitions.
2607:    WRITE (20, '(4(A, /), (A,2I5, /), 18(A, /), A)') '%!PS-Adobe-3.0 EPSF-3.0', &2599:    WRITE (20, '(4(A, /), (A,2I5, /), 18(A, /), A)') '%!PS-Adobe-3.0 EPSF-3.0', &
2608:       & '%% Creator: disconnection.f90', &2600:       & '%% Creator: disconnection.f90', &
2609:       & '%% Minima: ' // TRIM(file_min), &2601:       & '%% Minima: ' // TRIM(file_min), &
2610:       & '%% Transition states: ' // TRIM(file_ts), &2602:       & '%% Transition states: ' // TRIM(file_ts), &
2611:       & '%%BoundingBox: 0 0 ',PAGE_X,PAGE_Y, &2603:       & '%%BoundingBox: 0 0 ',PAGE_X,PAGE_Y, &
2612:       & '/ls {lineto stroke} def', &2604:       & '/ls {lineto stroke} def', &
2613:       & '/mt {moveto} def', &2605:       & '/mt {moveto} def', &
2614:       & '1.0 setlinewidth', &2606:       & '0.1 setlinewidth', &
2615:       & '/swidth { 30 mul } def', &2607:       & '/swidth { 30 mul } def', &
2616:       & '/dpi { 30 div } def', &2608:       & '/dpi { 30 div } def', &
2617:       & '/Rez 1000 dpi def % calc lines per point for fill', &2609:       & '/Rez 1000 dpi def % calc lines per point for fill', &
2618:       & '/Step 1 Rez div def', &2610:       & '/Step 1 Rez div def', &
2619:       & '/StartRGB { /B1 exch def /S1 exch def /H1 exch def', &2611:       & '/StartRGB { /B1 exch def /S1 exch def /H1 exch def', &
2620:       & '          } def', &2612:       & '          } def', &
2621:       & '/EndRGB { /B2 exch def /S2 exch def /H2 exch def', &2613:       & '/EndRGB { /B2 exch def /S2 exch def /H2 exch def', &
2622:       & '          /Steps s Step div def', &2614:       & '          /Steps s Step div def', &
2623:       & '          /Hstep H2 H1 sub Steps div def', &2615:       & '          /Hstep H2 H1 sub Steps div def', &
2624:       & '          /Sstep S2 S1 sub Steps div def', &2616:       & '          /Sstep S2 S1 sub Steps div def', &
2625:       & '          /Bstep B2 B1 sub Steps div def', &2617:       & '          /Bstep B2 B1 sub Steps div def', &
2626:       & '        } def', &2618:       & '        } def', &
2627:       & '/StepRGB { /H1 H1 Hstep add def', &2619:       & '/StepRGB { /H1 H1 Hstep add def', &
2628:       & '           /S1 S1 Sstep add def', &2620:       & '           /S1 S1 Sstep add def', &
2629:       & '           /B1 B1 Bstep add def', &2621:       & '           /B1 B1 Bstep add def', &
2630:       & '           H1 S1 B1 setrgbcolor', &2622:       & '           H1 S1 B1 setrgbcolor', &
2631:       & '         } def'2623:       & '         } def'
2632: 2624: 
2633: 2625: 
2634: IF (SCALEBAR) THEN 
2635:    X1 = MARGIN_X + SCALE_WIDTH - 10.0D0 
2636:    Y1 = Y_POS(1.0D0*(N_LEVELS-NDELTASB-1)/(N_LEVELS-1)) 
2637:    Y2 = Y_POS(1.0D0) 
2638:    WRITE (20, '(A)') '1.0 0.0 0.0 setrgbcolor' 
2639:    WRITE (20, '(A)') '4.0 setlinewidth' 
2640:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x1, y2, ' ls' 
2641:    X2 = X1 + TICK_WIDTH*2 
2642:    X1 = X1 - TICK_WIDTH*2 
2643:    Y1 = Y_POS(1.0D0) 
2644:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y1, ' ls' 
2645:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y1, ' ls' 
2646:    Y1 = Y_POS(1.0D0*(N_LEVELS-NDELTASB-1)/(N_LEVELS-1)) 
2647:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y1, ' ls' 
2648:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x2, y1, ' ls' 
2649:    WRITE (20, '(A)') '0.0 0.0 0.0 setrgbcolor' 
2650:    WRITE (20, '(A)') '1.0 setlinewidth' 
2651:    WRITE (20, '(A, I2, A)') & 
2652:       & '/Times-Roman findfont ', font_size, ' scalefont setfont' 
2653:    X1 = X1 + TICK_WIDTH*4 
2654:    Y1 = Y_POS(1.0D0-1.0D0*NDELTASB/(2*N_LEVELS-2)) - 5.0D0 
2655:    WRITE (20, '(2F7.2, A)') & 
2656:       & X1, Y1, ' mt (epsilon) show' 
2657: ELSE 
2658: !  Vertical energy axis.2626: !  Vertical energy axis.
2659:    X1 = MARGIN_X + SCALE_WIDTH - 10.0D02627:    X1 = MARGIN_X + SCALE_WIDTH - 10.0D0
2660:    Y1 = Y_POS(0.0D0)2628:    Y1 = Y_POS(0.0D0)
2661:    Y2 = Y_POS(1.0D0)2629:    Y2 = Y_POS(1.0D0)
2662:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x1, y2, ' ls'2630:    WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x1, y2, ' ls'
2663: 2631: 
2664: !  Ticks on axis.2632: !  Ticks on axis.
2665:    X2 = X1 + TICK_WIDTH/2.0D02633:    X2 = X1 + TICK_WIDTH/2.0D0
2666:    X1 = X1 - TICK_WIDTH/2.0D02634:    X1 = X1 - TICK_WIDTH/2.0D0
2667:    DO I=1, N_LEVELS2635:    DO I=1, N_LEVELS
2672: !  Tick labels.2640: !  Tick labels.
2673:    WRITE (20, '(A, I2, A)') &2641:    WRITE (20, '(A, I2, A)') &
2674:       & '/Times-Roman findfont ', font_size, ' scalefont setfont'2642:       & '/Times-Roman findfont ', font_size, ' scalefont setfont'
2675:    DO I=1, N_LEVELS2643:    DO I=1, N_LEVELS
2676:       X1 = 15.02644:       X1 = 15.0
2677:       Y1 = Y_POS(1.0D0*(N_LEVELS-I)/(N_LEVELS-1)) - 5.0D02645:       Y1 = Y_POS(1.0D0*(N_LEVELS-I)/(N_LEVELS-1)) - 5.0D0
2678:       WRITE (20, '(2F7.2, A, '//TRIM(lab_fmt)//', A)') &2646:       WRITE (20, '(2F7.2, A, '//TRIM(lab_fmt)//', A)') &
2679: !     WRITE (20, '(2F7.2, A, F8.3, A)') &2647: !     WRITE (20, '(2F7.2, A, F8.3, A)') &
2680:          & X1, Y1, ' mt (', e_high - (n_levels-i)*delta_e, ') show'2648:          & X1, Y1, ' mt (', e_high - (n_levels-i)*delta_e, ') show'
2681:    END DO2649:    END DO
2682: ENDIF 
2683: 2650: 
2684: !  Landmark for manipulate.f90.2651: !  Landmark for manipulate.f90.
2685:    WRITE (20, '(/, A, /)') '% Manipulate landmark.'2652:    WRITE (20, '(/, A, /)') '% Manipulate landmark.'
2686: 2653: 
2687: END SUBROUTINE HEADER2654: END SUBROUTINE HEADER
2688: !................................................................................!2655: !................................................................................!
2689: !DP X_POS(DP X) 2656: !DP X_POS(DP X) 
2690: FUNCTION X_POS(X)2657: FUNCTION X_POS(X)
2691: 2658: 
2692:    USE PAGE2659:    USE PAGE


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0