Sources Esope | Cast3M

calj00
C CALJ00    SOURCE    CB215821  26/08/24    21:15:22     12622          
      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 1003 M=1,ND
      DO 15 N=1,ND
      AJ(N,M,L)=0.D0
   15 CONTINUE
 1003 CONTINUE
C
      DO 1005 M=1,ND
      DO 1004 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
 1004 CONTINUE
 1005 CONTINUE
C
   20 CONTINUE
C
      RETURN
 1002 FORMAT(10(1X,1PE11.4))
 1001 FORMAT(20(1X,I5))
      END