hdiff output

r29910/genrigid.f90 2016-02-05 18:30:10.120364013 +0000 r29909/genrigid.f90 2016-02-05 18:30:10.312366514 +0000
2594: DOUBLE PRECISION, INTENT(OUT) :: SUBSET_BAND(:)2594: DOUBLE PRECISION, INTENT(OUT) :: SUBSET_BAND(:)
2595: INTEGER :: J1, J2, J3, NEW_NATOMS, OFFSET2595: INTEGER :: J1, J2, J3, NEW_NATOMS, OFFSET
2596: 2596: 
2597: ! May eventually change this so we can specify a list of RBs rather than just using the ones which come first in the list2597: ! May eventually change this so we can specify a list of RBs rather than just using the ones which come first in the list
2598: DO J1 = 1, MAX_RB2598: DO J1 = 1, MAX_RB
2599:     ! Add up how many atoms are in the subset2599:     ! Add up how many atoms are in the subset
2600:     NEW_NATOMS = NEW_NATOMS + NSITEPERBODY(J1)2600:     NEW_NATOMS = NEW_NATOMS + NSITEPERBODY(J1)
2601: ENDDO2601: ENDDO
2602: 2602: 
2603: IF (ATOMRIGIDCOORDT) THEN2603: IF (ATOMRIGIDCOORDT) THEN
 2604:     CALL FLUSH()
2604:     DO J1 = 1, NCONFIGS2605:     DO J1 = 1, NCONFIGS
2605:         OFFSET = 02606:         OFFSET = 0
2606:         DO J2 = 1, MAX_RB2607:         DO J2 = 1, MAX_RB
2607:             DO J3 = 1, NSITEPERBODY(J2)2608:             DO J3 = 1, NSITEPERBODY(J2)
2608:                 SUBSET_BAND(3*NEW_NATOMS*(J1-1)+OFFSET+3*J3-2:3*NEW_NATOMS*J1+OFFSET+3*J3) = &2609:                 SUBSET_BAND(3*NEW_NATOMS*(J1-1)+OFFSET+3*J3-2:3*NEW_NATOMS*J1+OFFSET+3*J3) = &
2609:                     BAND(NOPT*(J1-1)+3*RIGIDGROUPS(J3,J2)-2:NOPT*(J1-1)+3*RIGIDGROUPS(J3,J2))2610:                     BAND(NOPT*(J1-1)+3*RIGIDGROUPS(J3,J2)-2:NOPT*(J1-1)+3*RIGIDGROUPS(J3,J2))
2610:             ENDDO2611:             ENDDO
2611:             OFFSET = OFFSET + NSITEPERBODY(J2)  ! Keep a record of how many atoms we've already added in for this configuration2612:             OFFSET = OFFSET + NSITEPERBODY(J2)  ! Keep a record of how many atoms we've already added in for this configuration
2612:         ENDDO2613:         ENDDO
2613:     ENDDO2614:     ENDDO
2614: 2615: 
2615: ELSE2616: ELSE
2616:     SUBSET_BAND(:) = 0.0D02617:     SUBSET_BAND(:) = 0.0D0
 2618:     CALL FLUSH()
2617:     DO J1 = 1, NCONFIGS2619:     DO J1 = 1, NCONFIGS
2618:         DO J2 = 1, MAX_RB2620:         DO J2 = 1, MAX_RB
2619:             SUBSET_BAND(3*NEW_NATOMS*(J1-1)+3*J2-2:3*NEW_NATOMS*(J1-1)+3*J2) = & ! centre-of-mass coordinates2621:             SUBSET_BAND(3*NEW_NATOMS*(J1-1)+3*J2-2:3*NEW_NATOMS*(J1-1)+3*J2) = & ! centre-of-mass coordinates
2620:                     BAND(NOPT*(J1-1)+3*J2-2:NOPT*(J1-1)+3*J2)2622:                     BAND(NOPT*(J1-1)+3*J2-2:NOPT*(J1-1)+3*J2)
2621:             SUBSET_BAND(3*NEW_NATOMS*(J1-1)+3*MAX_RB+3*J2-2:3*NEW_NATOMS*(J1-1)+3*MAX_RB+3*J2) = & 2623:             SUBSET_BAND(3*NEW_NATOMS*(J1-1)+3*MAX_RB+3*J2-2:3*NEW_NATOMS*(J1-1)+3*MAX_RB+3*J2) = & 
2622:                     BAND(NOPT*(J1-1)+3*NRIGIDBODY+3*J2-2:NOPT*(J1-1)+3*NRIGIDBODY+3*J2)  ! AA coordinates2624:                     BAND(NOPT*(J1-1)+3*NRIGIDBODY+3*J2-2:NOPT*(J1-1)+3*NRIGIDBODY+3*J2)  ! AA coordinates
2623:         ENDDO2625:         ENDDO
2624:     ENDDO2626:     ENDDO
2625: ENDIF2627: ENDIF
2626: 2628: 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0