Esope source files

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