Esope source files

   1 : C DEVOLJ    SOURCE    CHAT      05/01/12    22:46:50     5004           
   2 :       SUBROUTINE DEVOLJ(XE,XEL,SHP,NBBB,NBNO,IFOU,DJAC)
   3 : C=======================================================================
   4 : C     
   5 : C    CALCULE L'ELEMENT DE VOLUME A PARTIR DU JACOBIEN POUR
   6 : C    LE JOINT POREUX
   7 : C     
   8 : C=======================================================================
   9 : C  ENTREES
  10 : C     XE(3,NBBB) =COORDONNEES GLOBALES DE L ELEMENT
  11 : C     XEL(3,NBBB)=COORDONNEES LOCALES DE L ELEMENT
  12 : C     SHP(6,NBNO)=DERIVEES PAR RAPPORT A LA GEOMETRIE DE REFERENCE
  13 : C     NBBB       =NOMBRE DE NOEUDS
  14 : C     NBNO       =NOMBRE DE FONCTIONS DE FORME
  15 : C     IFOU       =IFOUR DE CCOPTIO
  16 : C     
  17 : C  SORTIE
  18 : C     SHP(6,NBNO)=DERIVEES PAR RAPPORT A LA GEOMTRIE DE L ELEMENT
  19 : C     DJAC       =JACOBIEN
  20 : C     
  21 : C=======================================================================
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8(A-H,O-Z)
  24 : C    Include contenant quelques constantes dont XPI :
  25 : -INC CCREEL
  26 : *     
  27 :       DIMENSION XE(3,*),XEL(3,*),SHP(6,*)
  28 : *     
  29 :       NFAC=(3*NBBB-NBNO)/4
  30 : *     
  31 :       IF (IFOU.EQ.2) THEN
  31 :          CALL JACOBP(XEL,SHP,2,NFAC,NBNO,DJAC)
  32 :       ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  32 :          CALL JACOBP(XEL,SHP,1,NFAC,NBNO,DJAC)
  33 :       ENDIF
  34 : *     
  35 :       IF  (IFOU.EQ.0) THEN
  35 :          CALL DISTRR(XE,SHP,NFAC,RR)
  36 :          DJAC=DJAC*RR*2.D0*XPI
  37 :       ENDIF
  38 : *     
  39 :       RETURN
  40 :       END
  41 : 
  42 : 
  43 : 
  44 : 
  45 :  
  46 :  
  47 :