hdiff output

r29201/QALCS_surface.f90 2015-11-17 23:32:55.920450035 +0000 r29200/QALCS_surface.f90 2015-11-17 23:32:56.112452609 +0000
625:   ! Sort VSITES by weight from biggest to smallest (CHECK!)625:   ! Sort VSITES by weight from biggest to smallest (CHECK!)
626:   CALL SORT3_V2(NVSITES,NVSITES_MAX,VSITE_WEIGHTS,VSITES)626:   CALL SORT3_V2(NVSITES,NVSITES_MAX,VSITE_WEIGHTS,VSITES)
627:   ! weed out vacancies whose avergage coordination is below 627:   ! weed out vacancies whose avergage coordination is below 
628:   !628:   !
629:   !ds656> testing...629:   !ds656> testing...
630:   !WRITE(MYUNIT,'(A)') 'gen_vsites> Vsites sorted by occupancy:'  630:   !WRITE(MYUNIT,'(A)') 'gen_vsites> Vsites sorted by occupancy:'  
631:   !DO I1=1,NVSITES631:   !DO I1=1,NVSITES
632:   !   WRITE(MYUNIT,'(A,I6,I3,3(1X,F10.5))') 'gen_sites>', &632:   !   WRITE(MYUNIT,'(A,I6,I3,3(1X,F10.5))') 'gen_sites>', &
633:   !        I1, VSITE_WEIGHTS(I1), VSITES(I1,1:3)633:   !        I1, VSITE_WEIGHTS(I1), VSITES(I1,1:3)
634:   !ENDDO634:   !ENDDO
635:   ! write(*,'(I5)') natoms+nvsites635:   write(*,'(I5)') natoms+nvsites
636:   ! write(*,'(A)') 'Atoms and all 1st generation vacancies'636:   write(*,'(A)') 'Atoms and all 1st generation vacancies'
637:   ! DO I0=1,NSPECIES(0) !ds656> Should not exceed 10           637:   DO I0=1,NSPECIES(0) !ds656> Should not exceed 10           
638:   !    DO I1=1,ATOMLISTS(I0,1,0)638:      DO I1=1,ATOMLISTS(I0,1,0)
639:   !       J1=ATOMLISTS(I0,1,I1)639:         J1=ATOMLISTS(I0,1,I1)
640:   !       WRITE(*,'(A,I1,3(1X,F20.10))') 'A',I0,&640:         WRITE(*,'(A,I1,3(1X,F20.10))') 'A',I0,&
641:   !            (COORDS(J2,NP), J2=3*(J1-1)+1,3*J1)641:              (COORDS(J2,NP), J2=3*(J1-1)+1,3*J1)
642:   !    ENDDO642:      ENDDO
643:   ! ENDDO643:   ENDDO
644:   ! DO I1=1,NVSITES644:   DO I1=1,NVSITES
645:   !    WRITE(*,'(I2,3(1X,F20.10))') VSITE_WEIGHTS(I1),&645:      WRITE(*,'(I2,3(1X,F20.10))') VSITE_WEIGHTS(I1),&
646:   !         (VSITES(I1, J2), J2=1,3)646:           (VSITES(I1, J2), J2=1,3)
647:   ! ENDDO647:   ENDDO
648:   !<ds656 ...testing648:   !<ds656 ...testing
649:   !649:   !
650:   I2=NVSITES650:   I2=NVSITES
651:   IF(VSITE_WEIGHTS(1) > 1) THEN ! or 1?651:   IF(VSITE_WEIGHTS(1) > 1) THEN ! or 1?
652:      ! Weed out all but the largest WEIGHTS (highest coordination)652:      ! Weed out all but the largest WEIGHTS (highest coordination)
653:      WEED: DO I1=I2,1,-1653:      WEED: DO I1=I2,1,-1
654:         IF(VSITE_WEIGHTS(I1) < VSITE_WEIGHTS(1)) THEN654:         IF(VSITE_WEIGHTS(I1) < VSITE_WEIGHTS(1)) THEN
655:            NVSITES = NVSITES-1655:            NVSITES = NVSITES-1
656:            VSITE_WEIGHTS(I1) = 0656:            VSITE_WEIGHTS(I1) = 0
657:            VSITES(I1,1:3) = 0.0D0657:            VSITES(I1,1:3) = 0.0D0


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0