Esope source files

   1 : C DETJAC    SOURCE    CHAT      05/01/12    22:44:08     5004           
   2 : C           
   3 :       SUBROUTINE DETJAC(XNOEU,XY,ITYP,IDIM,NOEL1,IFORMU,SHP,DET)
   4 : C*******************************************************************
   5 : C     CALCUL DU JACOBIEN POUR LE PASSAGE D'UN SCALAIRE  DE 
   6 : C            L'ELEMENT DE REFERENCE A L'ELEMENT REEL
   7 : C********************************************************************           
   8 : C   ENTREES
   9 : C                                  
  10 : C     XY COORDONNEES DE REFERENCES DU POINT DE CALCUL                           
  11 : C     ITYP TYPE DE L ELEMENT                                                    
  12 : C     XNOEU COORDONNEES REELLES DES NOEUDS DE L ELEMENT                         
  13 : C     NOEL1  NOMBRE DE NOEUDS DE L ELEMENT POSITION                             
  14 : C                   
  15 : C           COORDONNEES DE REFERENCE AU POINT DE CALCUL                         
  16 : C
  17 : C     SHP(1,I) FONCTIONS DE FORMES AU POINT DE CALCUL
  18 : C     SHP(2,I) DERIVEE DES FONCTIONS DE FORME / XSI AU POINT DE CALCUL
  19 : C     SHP(3,I)    "     "      "     "    "   / ETA "    "   "    "
  20 : C     SHP(4,I)    "     "      "      "   "   /DZETA "   "    "   "
  21 : C
  22 : C   SORTIE
  23 : C     DET VALEUR DU JACOBIEN
  24 : C ********************************************************************          
  25 : C                                                                               
  26 : C                                                                               
  27 :       IMPLICIT INTEGER(I-N)
  28 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  29 : C
  30 : C                                                                               
  31 : C                     
  32 :       DIMENSION XNOEU(3,*),SHP(6,*),XY(*)
  33 : C                                                                               
  34 : C                                                                               
  35 : C***  CALCUL DES XI ET DES DXI                                                                                                                     
  36 : C   
  36 :       CALL SHAPE(XY(1),XY(2),XY(3),ITYP,SHP,IRET)                              
  37 : C                                                                               
  38 : C***  COEFFICIENTS DE DERIVATION                                                
  39 : C                      
  40 :       DET=0.D0
  41 :       IF(IDIM.EQ.2)THEN
  42 : C                                                                               
  43 : C***  BIDIM                                                                     
  44 : C                                                                               
  45 :             A=0.D0
  46 :             B=0.D0
  47 :             C=0.D0
  48 :             D=0.D0
  49 :             DO 63 I=1,NOEL1                                                           
  50 :                  A=A+SHP(2,I)*XNOEU(1,I)
  51 :                  B=B+SHP(3,I)*XNOEU(1,I)
  52 :                  C=C+SHP(2,I)*XNOEU(2,I)
  53 :                  D=D+SHP(3,I)*XNOEU(2,I)
  54 :   63        CONTINUE 
  55 :             DET=A*D-B*C
  56 : C                                                                               
  57 : C***  TRIDIM                                                                    
  58 : C                                                                               
  59 :       ELSE
  60 :             A1=0.D0
  61 :             A2=0.D0
  62 :             A3=0.D0
  63 :             B1=0.D0
  64 :             B2=0.D0
  65 :             B3=0.D0
  66 :             C1=0.D0
  67 :             C2=0.D0
  68 :             C3=0.D0
  69 :             DO 75 I=1,NOEL1                                                           
  70 :                  A1=A1+SHP(2,I)*XNOEU(1,I)
  71 :                  A2=A2+SHP(3,I)*XNOEU(1,I)
  72 :                  A3=A3+SHP(4,I)*XNOEU(1,I)
  73 :                  B1=B1+SHP(2,I)*XNOEU(2,I)
  74 :                  B2=B2+SHP(3,I)*XNOEU(2,I)
  75 :                  B3=B3+SHP(4,I)*XNOEU(2,I)
  76 :                  C1=C1+SHP(2,I)*XNOEU(3,I)
  77 :                  C2=C2+SHP(3,I)*XNOEU(3,I)
  78 :                  C3=C3+SHP(4,I)*XNOEU(3,I)
  79 :   75        CONTINUE
  80 :             DET=A1*(B2*C3-C2*B3)+A2*(C1*B3-B1*C3)+A3*(B1*C2-B2*C1)                    
  81 : C                                                                               
  82 :       ENDIF
  83 :       RETURN
  84 :       END          
  85 : 
  86 : 
  87 :  
  88 :