hdiff output

r30715/commons.f90 2016-07-08 14:30:16.528757836 +0100 r30714/commons.f90 2016-07-08 14:30:17.388769325 +0100
350:       DOUBLE PRECISION, ALLOCATABLE :: VSITES(:,:) 350:       DOUBLE PRECISION, ALLOCATABLE :: VSITES(:,:) 
351: 351: 
352: !ds656> Stress tensor 352: !ds656> Stress tensor 
353:       LOGICAL :: STRESST353:       LOGICAL :: STRESST
354:       INTEGER :: STRESS_MODE354:       INTEGER :: STRESS_MODE
355:       DOUBLE PRECISION, ALLOCATABLE :: STRESS(:,:,:)355:       DOUBLE PRECISION, ALLOCATABLE :: STRESS(:,:,:)
356: 356: 
357: !ds656> A saw-tooth temperature protocol357: !ds656> A saw-tooth temperature protocol
358:       LOGICAL :: SAWTOOTH358:       LOGICAL :: SAWTOOTH
359:       INTEGER :: SAWTOOTH_NREJMAX359:       INTEGER :: SAWTOOTH_NREJMAX
360:       DOUBLE PRECISION :: SAWTOOTH_TMAX, SAWTOOTH_TFAC, &360:       DOUBLE PRECISION :: SAWTOOTH_TMAX, SAWTOOTH_TFAC
361:            SAWTOOTH_SFAC, SAWTOOTH_SFAC2 
362:       361:       
363: !   allocatables362: !   allocatables
364: 363: 
365:       INTEGER,ALLOCATABLE,DIMENSION(:) :: MOVABLEATOMLIST         ! a list containing the movable atom indices364:       INTEGER,ALLOCATABLE,DIMENSION(:) :: MOVABLEATOMLIST         ! a list containing the movable atom indices
366:       LOGICAL,ALLOCATABLE,DIMENSION(:) :: MOVABLEATOMLISTLOGICAL  ! is atom i movable?365:       LOGICAL,ALLOCATABLE,DIMENSION(:) :: MOVABLEATOMLISTLOGICAL  ! is atom i movable?
367:       INTEGER,ALLOCATABLE,DIMENSION(:) :: ATOMSINBLOCK            ! for BLOCKMOVE, to split movableatoms into separate blocks366:       INTEGER,ALLOCATABLE,DIMENSION(:) :: ATOMSINBLOCK            ! for BLOCKMOVE, to split movableatoms into separate blocks
368:       INTEGER,ALLOCATABLE,DIMENSION(:) :: NSPECIES(:), NSPECIES_INI(:)             ! for multicomponent systems367:       INTEGER,ALLOCATABLE,DIMENSION(:) :: NSPECIES(:), NSPECIES_INI(:)             ! for multicomponent systems
369: 368: 
370:       INTEGER, ALLOCATABLE, DIMENSION(:) :: NT, NQ, BLOCK, JUMPINT, JUMPTO, NCORE, NSURFMOVES  ! dimension will be NPAR369:       INTEGER, ALLOCATABLE, DIMENSION(:) :: NT, NQ, BLOCK, JUMPINT, JUMPTO, NCORE, NSURFMOVES  ! dimension will be NPAR
371:       INTEGER, ALLOCATABLE, DIMENSION(:) :: HISTVALS, LHISTVALS  ! dimension declared HBINS in keyword370:       INTEGER, ALLOCATABLE, DIMENSION(:) :: HISTVALS, LHISTVALS  ! dimension declared HBINS in keyword


r30715/keywords.f 2016-07-08 14:30:16.820761741 +0100 r30714/keywords.f 2016-07-08 14:30:17.676773176 +0100
534: 534: 
535: !     ds656> Stress tensor calculation535: !     ds656> Stress tensor calculation
536:       STRESST=.FALSE.536:       STRESST=.FALSE.
537:       STRESS_MODE=1537:       STRESS_MODE=1
538: 538: 
539: !     ds656> Saw-tooth protocol539: !     ds656> Saw-tooth protocol
540:       SAWTOOTH=.FALSE.540:       SAWTOOTH=.FALSE.
541:       SAWTOOTH_NREJMAX=0541:       SAWTOOTH_NREJMAX=0
542:       SAWTOOTH_TMAX=0.3D0542:       SAWTOOTH_TMAX=0.3D0
543:       SAWTOOTH_TFAC=1.0D0543:       SAWTOOTH_TFAC=1.0D0
544:       SAWTOOTH_SFAC=1.0D0 
545:       SAWTOOTH_SFAC2=1.0D0 
546: 544: 
547: !     ds656> Generalised LJ with Yukawa545: !     ds656> Generalised LJ with Yukawa
548:       GLJY = .FALSE.546:       GLJY = .FALSE.
549:       GLJ_EXP = 6547:       GLJ_EXP = 6
550:       YUK_A = 0.0D0548:       YUK_A = 0.0D0
551:       YUK_XI = 1.0D0549:       YUK_XI = 1.0D0
552:       550:       
553:       CHRMMT=.FALSE.551:       CHRMMT=.FALSE.
554:       CHARMMTYPE=1552:       CHARMMTYPE=1
555:       CHARMMDFTBT=.FALSE.553:       CHARMMDFTBT=.FALSE.
1628:          GLJY = .TRUE.1626:          GLJY = .TRUE.
1629:          CALL READI(GLJ_EXP)1627:          CALL READI(GLJ_EXP)
1630:          CALL READF(YUK_A)1628:          CALL READF(YUK_A)
1631:          CALL READF(YUK_XI)1629:          CALL READF(YUK_XI)
1632: ! ds656> Saw-tooth protocol1630: ! ds656> Saw-tooth protocol
1633:       ELSE IF (WORD .EQ. 'SAWTOOTH') THEN1631:       ELSE IF (WORD .EQ. 'SAWTOOTH') THEN
1634:          SAWTOOTH=.TRUE.1632:          SAWTOOTH=.TRUE.
1635:          IF(NITEMS.GT.1) CALL READF(SAWTOOTH_TMAX)1633:          IF(NITEMS.GT.1) CALL READF(SAWTOOTH_TMAX)
1636:          IF(NITEMS.GT.2) CALL READF(SAWTOOTH_TFAC)1634:          IF(NITEMS.GT.2) CALL READF(SAWTOOTH_TFAC)
1637:          IF(NITEMS.GT.3) CALL READI(SAWTOOTH_NREJMAX)1635:          IF(NITEMS.GT.3) CALL READI(SAWTOOTH_NREJMAX)
1638:          IF(NITEMS.GT.4) THEN 
1639:             CALL READF(SAWTOOTH_SFAC) 
1640:             SAWTOOTH_SFAC2=EXP(-LOG(SAWTOOTH_SFAC)/DBLE(SAWTOOTH_NREJMAX)) 
1641:          ENDIF 
1642: ! ds656> Stress tensor calculation1636: ! ds656> Stress tensor calculation
1643:       ELSE IF (WORD .EQ. 'STRESS') THEN1637:       ELSE IF (WORD .EQ. 'STRESS') THEN
1644:          STRESST=.TRUE.1638:          STRESST=.TRUE.
1645:          IF(NITEMS.GT.1) CALL READI(STRESS_MODE)1639:          IF(NITEMS.GT.1) CALL READI(STRESS_MODE)
1646: ! BLN 1640: ! BLN 
1647:       ELSE IF (WORD.EQ.'BLN') THEN1641:       ELSE IF (WORD.EQ.'BLN') THEN
1648:          BLNT=.TRUE.1642:          BLNT=.TRUE.
1649:          CALL READF(RK_R)1643:          CALL READF(RK_R)
1650:          CALL READF(RK_THETA)1644:          CALL READF(RK_THETA)
1651:          ALLOCATE(BEADLETTER(NATOMS),BLNSSTRUCT(NATOMS),1645:          ALLOCATE(BEADLETTER(NATOMS),BLNSSTRUCT(NATOMS),


r30715/mc.F 2016-07-08 14:30:17.104765532 +0100 r30714/mc.F 2016-07-08 14:30:17.964777038 +0100
1458:                      WRITE(MYUNIT,73) JP,J1,POTEL,EPREV(JP),NSUCCESS(JP),NFAIL(JP)1458:                      WRITE(MYUNIT,73) JP,J1,POTEL,EPREV(JP),NSUCCESS(JP),NFAIL(JP)
1459: 73                   FORMAT('JP,J1,POTEL,EPREV,NSUC,NFAIL=',I2,I6,2F15.7,2I6,' TABOO')1459: 73                   FORMAT('JP,J1,POTEL,EPREV,NSUC,NFAIL=',I2,I6,2F15.7,2I6,' TABOO')
1460:                   ENDIF1460:                   ENDIF
1461:                   GOTO 231461:                   GOTO 23
1462:                ENDIF1462:                ENDIF
1463: 1463: 
1464: !ds656> For the saw-tooth temperature protocol, step size is1464: !ds656> For the saw-tooth temperature protocol, step size is
1465: !       dynamically adjusted to help escape from the current 1465: !       dynamically adjusted to help escape from the current 
1466: !       Markov minimum, but not too far...1466: !       Markov minimum, but not too far...
1467:                IF(SAWTOOTH) THEN1467:                IF(SAWTOOTH) THEN
1468:                   ! Energy criterion alone may not be adequate...1468:                   ! Energy criterion may not be enough...
1469:                   IF(ABS(POTEL-EPREV(JP)) < ECONV) THEN1469:                   IF(ABS(POTEL-EPREV(JP)) < ECONV) THEN
1470:                      STEP(JP)=SAWTOOTH_SFAC*STEP(JP)1470:                      STEP(JP)=1.10*STEP(JP)
1471:                   ELSE1471:                   ELSE
1472:                      STEP(JP)=SAWTOOTH_SFAC2*STEP(JP)1472:                      STEP(JP)=0.99*STEP(JP)
1473:                   ENDIF1473:                   ENDIF
1474:                   WRITE(MYUNIT,'(A,F10.5)') 'mc> STEP=', STEP(JP) 
1475:                ENDIF1474:                ENDIF
1476: !ds656> print some energies to test BASWAP1475: !ds656> print some energies to test BASWAP
1477:                IF(BASWAPTEST) THEN1476:                IF(BASWAPTEST) THEN
1478:                   !CALL BGUPTA(COORDS,GRAD,EASWAPQ,.FALSE.)1477:                   !CALL BGUPTA(COORDS,GRAD,EASWAPQ,.FALSE.)
1479:                   IF(NSWAPS == 1) THEN1478:                   IF(NSWAPS == 1) THEN
1480:                      WRITE(MYUNIT, '(2(A,1X,E12.6,1X))') 1479:                      WRITE(MYUNIT, '(2(A,1X,E12.6,1X))') 
1481:      1                    "mc> 1 swap: de_A=",DE1,"and de_B=",DE21480:      1                    "mc> 1 swap: de_A=",DE1,"and de_B=",DE2
1482:                   ENDIF1481:                   ENDIF
1483:                   WRITE(MYUNIT, '(2(A,1X,E12.6,1X),A)') 1482:                   WRITE(MYUNIT, '(2(A,1X,E12.6,1X),A)') 
1484:      2                 "mc> BA swaps: dE=", EASWAP-EPREV(JP), "before and dE=",1483:      2                 "mc> BA swaps: dE=", EASWAP-EPREV(JP), "before and dE=",


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0