hdiff output

r32249/minpermdistrbcom.f90 2017-04-05 10:30:12.407528998 +0100 r32248/minpermdistrbcom.f90 2017-04-05 10:30:12.687531626 +0100
710:          CMX = CMX + T(J-2)710:          CMX = CMX + T(J-2)
711:          CMY = CMY + T(J-1)711:          CMY = CMY + T(J-1)
712:          CMZ = CMZ + T(J)712:          CMZ = CMZ + T(J)
713:       ENDDO713:       ENDDO
714:       CMX = CMX/NATOMS; CMY = CMY/NATOMS; CMZ = CMZ/NATOMS714:       CMX = CMX/NATOMS; CMY = CMY/NATOMS; CMZ = CMZ/NATOMS
715:       DO I = 1, NATOMS715:       DO I = 1, NATOMS
716:          J      = 3*I716:          J      = 3*I
717:          X(J-2) = T(J-2) - CMX717:          X(J-2) = T(J-2) - CMX
718:          X(J-1) = T(J-1) - CMY718:          X(J-1) = T(J-1) - CMY
719: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)719: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)
720:          IF(EFIELDT.AND.MULTISITEPYT) THEN720:          IF(EFIELDT.AND.MULTISITEPYT)
721:           X(J)   = T(J)721:           X(J)   = T(J)
722:          ELSE722:          ELSE
723:           X(J)   = T(J)   - CMZ723:           X(J)   = T(J)   - CMZ
724:          END IF724:          END IF
725:       ENDDO725:       ENDDO
726: 726: 
727: !     extract the rotation matrix corresponding to rotation via Q727: !     extract the rotation matrix corresponding to rotation via Q
728: 728: 
729:       CALL QROTMAT(Q2,RM)729:       CALL QROTMAT(Q2,RM)
730: 730: 
768:          CMXA = CMXA + XA(J2+1)768:          CMXA = CMXA + XA(J2+1)
769:          CMYA = CMYA + XA(J2+2)769:          CMYA = CMYA + XA(J2+2)
770:          CMZA = CMZA + XA(J2+3)770:          CMZA = CMZA + XA(J2+3)
771:       ENDDO771:       ENDDO
772:       CMXA = CMXA/NSIZE; CMYA = CMYA/NSIZE; CMZA = CMZA/NSIZE772:       CMXA = CMXA/NSIZE; CMYA = CMYA/NSIZE; CMZA = CMZA/NSIZE
773:       DO J1 = 1, NSIZE773:       DO J1 = 1, NSIZE
774:          J2 = 3*(J1-1)774:          J2 = 3*(J1-1)
775:          XA(J2+1) = XA(J2+1) - CMXA775:          XA(J2+1) = XA(J2+1) - CMXA
776:          XA(J2+2) = XA(J2+2) - CMYA776:          XA(J2+2) = XA(J2+2) - CMYA
777: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)777: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)
778:          IF(EFIELDT.AND.MULTISITEPYT) THEN778:          IF(EFIELDT.AND.MULTISITEPYT)
779:           XA(J2+3) = XA(J2+3) 779:           XA(J2+3) = XA(J2+3) 
780:          ELSE780:          ELSE
781:           XA(J2+3) = XA(J2+3) - CMZA781:           XA(J2+3) = XA(J2+3) - CMZA
782:          END IF782:          END IF
783:       ENDDO783:       ENDDO
784: 784: 
785:       CMXB = 0.0D0; CMYB = 0.0D0; CMZB = 0.0D0785:       CMXB = 0.0D0; CMYB = 0.0D0; CMZB = 0.0D0
786:       DO J1 = 1, NSIZE786:       DO J1 = 1, NSIZE
787:          J2 = 3*(J1-1)787:          J2 = 3*(J1-1)
788:          CMXB = CMXB + XB(J2+1)788:          CMXB = CMXB + XB(J2+1)
789:          CMYB = CMYB + XB(J2+2)789:          CMYB = CMYB + XB(J2+2)
790:          CMZB = CMZB + XB(J2+3)790:          CMZB = CMZB + XB(J2+3)
791:       ENDDO791:       ENDDO
792:       CMXB = CMXB/NSIZE; CMYB = CMYB/NSIZE; CMZB = CMZB/NSIZE792:       CMXB = CMXB/NSIZE; CMYB = CMYB/NSIZE; CMZB = CMZB/NSIZE
793:       DO J1 = 1, NSIZE793:       DO J1 = 1, NSIZE
794:          J2 = 3*(J1-1)794:          J2 = 3*(J1-1)
795:          XB(J2+1) = XB(J2+1) - CMXB795:          XB(J2+1) = XB(J2+1) - CMXB
796:          XB(J2+2) = XB(J2+2) - CMYB796:          XB(J2+2) = XB(J2+2) - CMYB
797: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)797: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)
798:          IF(EFIELDT.AND.MULTISITEPYT) THEN798:          IF(EFIELDT.AND.MULTISITEPYT)
799:           XB(J2+3) = XB(J2+3)799:           XB(J2+3) = XB(J2+3)
800:          ELSE800:          ELSE
801:           XB(J2+3) = XB(J2+3) - CMZB801:           XB(J2+3) = XB(J2+3) - CMZB
802:          END IF802:          END IF
803:       ENDDO803:       ENDDO
804: 804: 
805:       QMAT(1:4,1:4) = 0.0D0805:       QMAT(1:4,1:4) = 0.0D0
806: 806: 
807:       DO J1 = 1, NSIZE807:       DO J1 = 1, NSIZE
808:          J2 = 3*(J1-1)808:          J2 = 3*(J1-1)
905: 905: 
906:       ENDIF906:       ENDIF
907: 907: 
908:       DIST = SQRT(MINV)908:       DIST = SQRT(MINV)
909: 909: 
910:       DO J1 = 1, NATOMS910:       DO J1 = 1, NATOMS
911:          J2 = 3*(J1-1)911:          J2 = 3*(J1-1)
912:          RB(J2+1) = RB(J2+1) - CMXB912:          RB(J2+1) = RB(J2+1) - CMXB
913:          RB(J2+2) = RB(J2+2) - CMYB913:          RB(J2+2) = RB(J2+2) - CMYB
914: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)914: ! do not move along the z axis in case of a gravity field (currently implemented for MULTISITEPY)
915:          IF(EFIELDT.AND.MULTISITEPYT) THEN915:          IF(EFIELDT.AND.MULTISITEPYT)
916:           RB(J2+3) = RB(J2+3)916:           RB(J2+3) = RB(J2+3)
917:          ELSE917:          ELSE
918:           RB(J2+3) = RB(J2+3) - CMZB918:           RB(J2+3) = RB(J2+3) - CMZB
919:          END IF919:          END IF
920:       ENDDO920:       ENDDO
921: 921: 
922:       CALL NEWROTGEOMA(NATOMS,RB,Q2,RM,CMXA,CMYA,CMZA)922:       CALL NEWROTGEOMA(NATOMS,RB,Q2,RM,CMXA,CMYA,CMZA)
923: 923: 
924:       DEALLOCATE(XA,XB)924:       DEALLOCATE(XA,XB)
925: 925: 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0