Sources Esope | Cast3M

calj00
C CALJ00    SOURCE    CHAT      05/01/12    21:46:31     5004
      SUBROUTINE CALJ00(GR,XYZ,ND,NP,NPG,AJ)
C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
C
C     CALCUL DU              JACOBIEN  AJ
C     DANS LES CAS 2D ET 3D
C
C     ND   DIMENSION ESPACE DE DEFINITION DE L'ELEMENT
C     NP   NOMBRE DE NOEUDS DE L'ELEMENT
C     NPG  NOMBRE DE POINTS D'INTEGRATION
C
C     XYZ  COORDONNEES
C     GR   GRADIENT
C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
      REAL*8 GR(ND,NP,1),XYZ(ND,NP)
      REAL*8 AJ(ND,ND,NPG)
C
      DO 20 L=1,NPG
      DO 15 M=1,ND
      DO 15 N=1,ND
      AJ(N,M,L)=0.D0
   15 CONTINUE
C
      DO 10 M=1,ND
      DO 10 N=1,ND
      DO 10 I=1,NP
      AJ(N,M,L)=AJ(N,M,L)+GR(M,I,L)*XYZ(N,I)
   10 CONTINUE
C
   20 CONTINUE
C
      RETURN
 1002 FORMAT(10(1X,1PE11.4))
 1001 FORMAT(20(1X,I5))
      END