C DEVOLJ    SOURCE    CHAT      05/01/12    22:46:50     5004           
      SUBROUTINE DEVOLJ(XE,XEL,SHP,NBBB,NBNO,IFOU,DJAC)
C=======================================================================
C     
C    CALCULE L'ELEMENT DE VOLUME A PARTIR DU JACOBIEN POUR
C    LE JOINT POREUX
C     
C=======================================================================
C  ENTREES
C     XE(3,NBBB) =COORDONNEES GLOBALES DE L ELEMENT
C     XEL(3,NBBB)=COORDONNEES LOCALES DE L ELEMENT
C     SHP(6,NBNO)=DERIVEES PAR RAPPORT A LA GEOMETRIE DE REFERENCE
C     NBBB       =NOMBRE DE NOEUDS
C     NBNO       =NOMBRE DE FONCTIONS DE FORME
C     IFOU       =IFOUR DE CCOPTIO
C     
C  SORTIE
C     SHP(6,NBNO)=DERIVEES PAR RAPPORT A LA GEOMTRIE DE L ELEMENT
C     DJAC       =JACOBIEN
C     
C=======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C    Include contenant quelques constantes dont XPI :
-INC CCREEL
*     
      DIMENSION XE(3,*),XEL(3,*),SHP(6,*)
*     
      NFAC=(3*NBBB-NBNO)/4
*     
      IF (IFOU.EQ.2) THEN
         CALL JACOBP(XEL,SHP,2,NFAC,NBNO,DJAC)
      ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
         CALL JACOBP(XEL,SHP,1,NFAC,NBNO,DJAC)
      ENDIF
*     
      IF  (IFOU.EQ.0) THEN
         CALL DISTRR(XE,SHP,NFAC,RR)
         DJAC=DJAC*RR*2.D0*XPI
      ENDIF
*     
      RETURN
      END