C BGFCQ4    SOURCE    CHAT      06/03/29    21:15:48     5360
      SUBROUTINE BGFCQ4(IGAU,XEL,SHPTOT,SHP,BGF,DJAC,IARR)
C |====================================================================|
C |                                                                    |
C |  CALCULE LA MATRICE DE GRADIAN EN FLEXION BGF POUR COQ4
C |     CODE SUO X.Z. MARS 87
C |INPUT                                                               |
C |   IGAU=NUMERO DU POINT DE GAUSS                                    |
C |   NBNO=NOMBRE DE NOEUDS                                            |
C |   LRE =NOMBRE DE COLONNES DE LA MATRICE BGR                        |
C |   XEL =COORDONNEES LOCALES DE L ELEMENT                            |
C |   SHPTOT(6,NBNO,NBGAU)=FONCTIONS DE FORMES ET DERIVEES             |
C |ZONE DE TRAVAIL                                                     |
C |   SHP(6,NBNO)=TABLEAU DE TRAVAIL                                   |
C |OUTPUT                                                              |
C |   DJAC=JACOBIEN                                                    |
C |   BGF(9,LRE)=MATRICE BGF                                           |
C |====================================================================|
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      PARAMETER (LRE=24,NBNO=4)
      DIMENSION XEL(3,*),BGF(9,*),SHP(6,*),SHPTOT(6,NBNO,*)
      DATA ZER/0.D0/
      IARR=0
C
      CALL ZERO(BGF,9,LRE)
C
C
      DO 101 NP=1,NBNO
      SHP(1,NP)=SHPTOT(1,NP,IGAU)
      SHP(2,NP)=SHPTOT(2,NP,IGAU)
      SHP(3,NP)=SHPTOT(3,NP,IGAU)
  101 CONTINUE
      CALL JACOBI(XEL,SHP,2,NBNO,DJAC)
      IF(DJAC.LE.ZER) IARR=1
      IF(DJAC.LE.ZER) RETURN
C
C     ====================== CALCUL
C
      K=1
      DO 102 NP=1,NBNO
      BGF(1,K+4)=SHP(2,NP)
      BGF(2,K+4)=SHP(3,NP)
      BGF(3,K+4)=SHP(1,NP)
      BGF(4,K+3)=-SHP(2,NP)
      BGF(5,K+3)=-SHP(3,NP)
      BGF(6,K+3)=-SHP(1,NP)
      K=K+6
  102 CONTINUE
      RETURN
      END