hdiff output

r22707/disconnectionDPS.f90 2015-11-17 23:36:45.787532585 +0000 r22706/disconnectionDPS.f90 2015-11-17 23:36:45.979535156 +0000
450:    INTEGER :: NMINTR=0450:    INTEGER :: NMINTR=0
451:    INTEGER :: CONNECTMIN=0451:    INTEGER :: CONNECTMIN=0
452:    INTEGER :: MINRANGE=10000452:    INTEGER :: MINRANGE=10000
453:    INTEGER, ALLOCATABLE :: MINIDS(:)453:    INTEGER, ALLOCATABLE :: MINIDS(:)
454:    INTEGER, ALLOCATABLE :: MINTRS(:)454:    INTEGER, ALLOCATABLE :: MINTRS(:)
455:    DOUBLE PRECISION, ALLOCATABLE :: DPMINTRS(:)455:    DOUBLE PRECISION, ALLOCATABLE :: DPMINTRS(:)
456:    DOUBLE PRECISION :: DELTA_E, E_HIGH456:    DOUBLE PRECISION :: DELTA_E, E_HIGH
457:    DOUBLE PRECISION :: TSTHRESH=1.0D100457:    DOUBLE PRECISION :: TSTHRESH=1.0D100
458:    DOUBLE PRECISION :: TSBARTHRESH=1.0D100458:    DOUBLE PRECISION :: TSBARTHRESH=1.0D100
459:    DOUBLE PRECISION :: CSCALE=1.0D0 459:    DOUBLE PRECISION :: CSCALE=1.0D0 
460:    DOUBLE PRECISION :: BENERGY=1.0D100  
461:    DOUBLE PRECISION :: TRSCALEX=90.0460:    DOUBLE PRECISION :: TRSCALEX=90.0
462:    DOUBLE PRECISION :: TRSCALEY=50.0461:    DOUBLE PRECISION :: TRSCALEY=50.0
463:    CHARACTER(LEN=120) :: FILE_MIN, FILE_PICK, FILE_TS, FILE_WEIGHTS, &    462:    CHARACTER(LEN=120) :: FILE_MIN, FILE_PICK, FILE_TS, FILE_WEIGHTS, &    
464:       FILE_TRACE463:       FILE_TRACE
465:    CHARACTER(LEN=10) :: LAB_FMT464:    CHARACTER(LEN=10) :: LAB_FMT
466:    LOGICAL :: BARRIERS, CENTRE_GMIN, DUMP_NUMBERS, DUMP_SIZES, EXCLUDEALL, &465:    LOGICAL :: BARRIERS, CENTRE_GMIN, DUMP_NUMBERS, DUMP_SIZES, EXCLUDEALL, &
467:       IDENTIFY, MONOTONIC, SPLIT, WEIGHTS466:       IDENTIFY, MONOTONIC, SPLIT, WEIGHTS
468:    LOGICAL :: IDENTIFY_NODE=.FALSE.467:    LOGICAL :: IDENTIFY_NODE=.FALSE.
469:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.468:    LOGICAL :: IDENTIFY_NODE_SIZE=.FALSE.
470:    LOGICAL :: IDMINT=.FALSE.469:    LOGICAL :: IDMINT=.FALSE.
471:    LOGICAL :: TRMINT=.FALSE.470:    LOGICAL :: TRMINT=.FALSE.
472:    LOGICAL :: TRPRINT=.FALSE.471:    LOGICAL :: TRPRINT=.FALSE.
473:    LOGICAL :: TRVALT=.FALSE.472:    LOGICAL :: TRVALT=.FALSE.
474:    LOGICAL :: BASINT=.FALSE. 
475:    LOGICAL :: WEIGHTLOWER=.FALSE.473:    LOGICAL :: WEIGHTLOWER=.FALSE.
476:    LOGICAL :: TRVALSCALET=.FALSE.474:    LOGICAL :: TRVALSCALET=.FALSE.
477:    LOGICAL :: MATCHMIN=.FALSE.475:    LOGICAL :: MATCHMIN=.FALSE.
478:    INTEGER, ALLOCATABLE :: NCONN(:)476:    INTEGER, ALLOCATABLE :: NCONN(:)
479: 477: 
480: END MODULE VARS478: END MODULE VARS
481: ! }}}479: ! }}}
482: !................................................................................!480: !................................................................................!
483: !MODULE PAGE{{{481: !MODULE PAGE{{{
484: MODULE PAGE482: MODULE PAGE
1149:                  ELSE1147:                  ELSE
1150:                   IF (DPMINTRS(P)==-1.0D0.or.DPMINTRS(P)< DPMARKNODE(I,K)) DPMARKNODE(I,K)=DPMINTRS(P)1148:                   IF (DPMINTRS(P)==-1.0D0.or.DPMINTRS(P)< DPMARKNODE(I,K)) DPMARKNODE(I,K)=DPMINTRS(P)
1151:                  ENDIF  1149:                  ENDIF  
1152: ! DPMINTRS(P)=-1.0D0=initial value. If any minimum is not listed at all the node is not coloured1150: ! DPMINTRS(P)=-1.0D0=initial value. If any minimum is not listed at all the node is not coloured
1153:             END DO 1151:             END DO 
1154: !         END DO1152: !         END DO
1155:        END DO   1153:        END DO   
1156:        WRITE (6, '(A, /)') 'Done.'1154:        WRITE (6, '(A, /)') 'Done.'
1157:    END IF1155:    END IF
1158: 1156: 
1159: ! VKD Print basin information. 
1160: ! For each node at a particular energy level, find all minima for which the node is a parent. 
1161:  
1162:    IF (BASINT) THEN 
1163:        OPEN(UNIT=32, FILE='basins', STATUS='REPLACE') 
1164:        WRITE (6, '(A)') 'Finding section assignments of nodes for printing.' 
1165:        DO I = 1, N_LEVELS 
1166:              IF ((E_HIGH-(I-1)*DELTA_E).LT.BENERGY) THEN 
1167:              DO P = 1, N_MIN 
1168:                IF (NCONN(P).LE.NCONNMIN) CYCLE 
1169: ! BASIN(I,P) is the basin to which minimum P belongs at level I. 
1170:                  K=BASIN(I,P) 
1171:                  IF (K.EQ.0) CYCLE 
1172:                  WRITE(32,*) P, K  
1173:              END DO 
1174:              EXIT 
1175:              ENDIF 
1176:        END DO 
1177:        WRITE (6, '(A, /)') 'Done.' 
1178:        CLOSE(32) 
1179:    END IF 
1180:  
1181:  
1182:  
1183:  
1184: !  Find the parent of each basin.1157: !  Find the parent of each basin.
1185: 1158: 
1186:    WRITE (6, '(A)') 'Finding the parent of each basin.'1159:    WRITE (6, '(A)') 'Finding the parent of each basin.'
1187:    ALLOCATE (PARENT(N_LEVELS, MAXVAL(NODES)))1160:    ALLOCATE (PARENT(N_LEVELS, MAXVAL(NODES)))
1188:    DO I = 2, N_LEVELS1161:    DO I = 2, N_LEVELS
1189:       DO J = 1, N_MIN1162:       DO J = 1, N_MIN
1190:          IF (NCONN(J).LE.NCONNMIN) CYCLE1163:          IF (NCONN(J).LE.NCONNMIN) CYCLE
1191:          IF (BASIN(I, J) > 0) THEN1164:          IF (BASIN(I, J) > 0) THEN
1192:             PARENT(I, BASIN(I, J)) = BASIN(I-1, J)1165:             PARENT(I, BASIN(I, J)) = BASIN(I-1, J)
1193:          END IF1166:          END IF
1331:                IF (DPMARKNODE(I,K).GE.0.0D0) THEN1304:                IF (DPMARKNODE(I,K).GE.0.0D0) THEN
1332:                   R2=DPMARKNODE(I,K)*3.0D0*CSCALE1305:                   R2=DPMARKNODE(I,K)*3.0D0*CSCALE
1333:                   IF (R2.GT.3.0D0) THEN1306:                   IF (R2.GT.3.0D0) THEN
1334:                      WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k1307:                      WRITE(6,*) 'WARNING: value of DPMARKNODE is ',DPMARKNODE(I,K),' for ',i,k
1335:                      R2=3.0D01308:                      R2=3.0D0
1336:                   END IF1309:                   END IF
1337:                    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'1310:                    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'
1338:                    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' 1311:                    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' 
1339:                    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'1312:                    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'
1340:                    IF (R2.LE.3.0d0.AND.R2.GE.2.0d0)  &1313:                    IF (R2.LE.3.0d0.AND.R2.GE.2.0d0)  &
1341:   &         WRITE (20, '(A,2F7.4,A)') '0 ',(3.0d0-R2)/2.0D0, ((2.0D0-R2)-1.67d0)/(-2.67D0), '  setrgbcolor'            1314:   &         WRITE (20, '(A,2F7.4,A)') '0 ',(3.0d0-R2)/2.0D0, ((2.0D0-R2)-1.67d0)/(-2.67D0), '  setrgbcolor'
1342: !                  IF (R2.LE.1.AND.R2.GE.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'1315: !                  IF (R2.LE.1.AND.R2.GE.0)  WRITE (20, '(A,F6.4,A)') '1 0 ',(1-R2)/2,'  setrgbcolor'
1343: !                  IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 1316: !                  IF (R2.LE.3.AND.R2.GT.1)  WRITE (20, '(A,F6.4,A)') '1 ',(R2-1)/2,' 0  setrgbcolor' 
1344: !                  IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'1317: !                  IF (R2.LE.5.AND.R2.GT.3)  WRITE (20, '(F6.4,A)') (5-R2)/2,' 1 0  setrgbcolor'
1345: !                  IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'1318: !                  IF (R2.LE.7.AND.R2.GT.5)  WRITE (20, '(A,F6.4,A)') '0 1 ',(R2-5)/2  ,'  setrgbcolor'
1346: !                  IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'1319: !                  IF (R2.LE.9.AND.R2.GT.7)  WRITE (20, '(A,F6.4,A)') '0 ',(9-R2)/2 ,' 1  setrgbcolor'
1347: !                  IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'1320: !                  IF (R2.LE.10.AND.R2.GT.9)  WRITE (20, '(F6.4,A)') (R2-9)/2 ,' 0 1  setrgbcolor'
1348:                   WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1321:                   WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1349:                   WRITE (20, '(A)') '0 0 0 setrgbcolor'1322:                   WRITE (20, '(A)') '0 0 0 setrgbcolor'
1350:                ELSE1323:                ELSE
1351:                  WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'1324:                  WRITE (20, '(2F7.2, A, 2F7.2, A)') x1, y1, ' mt ', x3, y3, ' ls'
1758:          CLOSE(22)1731:          CLOSE(22)
1759:       CASE ('TRVALSCALE')1732:       CASE ('TRVALSCALE')
1760:          TRVALSCALET = .TRUE.1733:          TRVALSCALET = .TRUE.
1761:          CALL GET_DP(TRSCALEX)1734:          CALL GET_DP(TRSCALEX)
1762:          CALL GET_DP(TRSCALEY)1735:          CALL GET_DP(TRSCALEY)
1763:       CASE ('TS')1736:       CASE ('TS')
1764:          CALL GET_STRING(FILE_TS)1737:          CALL GET_STRING(FILE_TS)
1765:       CASE ('WEIGHTS')1738:       CASE ('WEIGHTS')
1766:          CALL GET_STRING(FILE_WEIGHTS)1739:          CALL GET_STRING(FILE_WEIGHTS)
1767:          WEIGHTS=.TRUE.1740:          WEIGHTS=.TRUE.
1768:       CASE ('BASINT') 
1769:          BASINT = .TRUE. 
1770:          CALL GET_DP(BENERGY) 
1771:       CASE DEFAULT1741:       CASE DEFAULT
1772:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)1742:          WRITE (6, '(2A)') 'Keyword not recognised in dinfo: ', TRIM(keyword)
1773:          STOP1743:          STOP
1774:       END SELECT1744:       END SELECT
1775:    END DO1745:    END DO
1776:    CLOSE (20)1746:    CLOSE (20)
1777: 1747: 
1778:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))1748:    IF (.NOT. ALLOCATED(MINTRS)) ALLOCATE(MINTRS(1))
1779: 1749: 
1780:    IF (DELTA_E <= 0.0D0) THEN1750:    IF (DELTA_E <= 0.0D0) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0