Esope source files

   1 : C CALJAJ    SOURCE    CHAT      05/01/12    21:46:37     5004           
   2 :       SUBROUTINE CALJAJ(GR,XYZ,ND,IES,NP,NPG,AJ)
   3 : C************************************************************************
   4 : C     
   5 : C     CALJAJ  APPELLE CALJ00
   6 : C     
   7 : C     
   8 : C     CALCUL DU JACOBIEN  AJ=1/J
   9 : C     CALCUL DE L'ELEMENT D'AIRE       SQ=DET(J)
  10 : C     DANS LES CAS 1D 2D ET 3D
  11 : C     
  12 : C     ND   DIMENSION ESPACE
  13 : C     NP   NOMBRE DE NOEUDS DE L'ELEMENT
  14 : C     NPG  NOMBRE DE POINTS D'INTEGRATION
  15 : C     
  16 : C     XYZ  COORDONNEES
  17 : C     GR   GRADIENT
  18 : C************************************************************************
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8 (A-H,O-Z)
  21 : C     
  22 :       REAL*8 GR(ND,NP,1),XYZ(IES,NP),XY(20),AJ(3,3,9)
  23 : C     
  24 : C     WRITE(6,*)'SUB CALJAJ ND,IES',ND,IES
  25 :       IF(ND.EQ.IES)GOTO 20
  25 :       CALL CALJQ(XYZ,XY,IES,ND,NP)
  26 : C     WRITE(6,*)'SUB CALJAJ XYZ',XYZ
  27 : C     WRITE(6,*)'SUB CALJAJ XY',XY
  27 :       CALL CALJ00(GR,XY,ND,NP,NPG,AJ)
  28 : C     
  29 :       RETURN
  30 :  20   CONTINUE
  30 :       CALL CALJ00(GR,XYZ,ND,NP,NPG,AJ)
  31 :       RETURN
  32 :  1002 FORMAT(10(1X,1PE11.4))
  33 :  1001 FORMAT(20(1X,I5))
  34 :       END
  35 : 
  36 : 
  37 :  
  38 :