Sources Esope | Cast3M

caljxy
C CALJXY    SOURCE    CHAT      05/01/12    21:47:03     5004
      SUBROUTINE CALJXY(XYZ,P,IES,XY,ND,NP)
C************************************************************************
C
C     CE SP CALCULE XY = P * XYZ (EN GROS)
C
C
C
C
C************************************************************************
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      REAL*8 XYZ(IES,NP),P(IES,IES),XY(ND,NP)
 
C***
 
 
      DO 43 I=1,NP
      DO 43 N=1,ND
      XY(N,I)=0.
      DO 43 M=1,IES
      XY(N,I)=XY(N,I)+P(N,M)*(XYZ(M,I)-XYZ(M,1))
 43   CONTINUE
 
 
C     WRITE(6,*)'SUB CALJXY XYZ'
C     WRITE(6,1983)XYZ
C1983 FORMAT(3(2X,1PE11.4))
C     WRITE(6,*)'SUB CALJCC XY'
C     WRITE(6,1982)XY
C1982 FORMAT(2(2X,1PE11.4))
 
 
      RETURN
      END