C CCGDIV    SOURCE    GOUNAND   05/12/21    21:16:26     5281           
      SUBROUTINE CCGDIV(LCOF,
     $     FC,
     $     IMPR,IRET)
      IMPLICIT REAL*8 (A-H,O-Z)
      IMPLICIT INTEGER (I-N)
C***********************************************************************
C NOM         : CCGDIV
C DESCRIPTION : Calcul de la loi de comportement aux points de Gauss : 
C               une divergence
C
C     
C LANGAGE     : ESOPE
C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
C               mél : gounand@semt2.smts.cea.fr
C***********************************************************************
C APPELES          : 
C APPELE PAR       : 
C***********************************************************************
C ENTREES            :
C ENTREES/SORTIES    :
C SORTIES            : -
C TRAVAIL            :
C***********************************************************************
C VERSION    : v1, 11/05/04, version initiale
C HISTORIQUE : v1, 11/05/04, création
C HISTORIQUE :
C HISTORIQUE :
C***********************************************************************
C Prière de PRENDRE LE TEMPS de compléter les commentaires
C en cas de modification de ce sous-programme afin de faciliter
C la maintenance !
C***********************************************************************
-INC CCOPTIO
CBEGININCLUDE SMCHAEL
      SEGMENT MCHAEL
      POINTEUR IMACHE(N1).MELEME
      POINTEUR ICHEVA(N1).MCHEVA
      ENDSEGMENT
      SEGMENT MCHEVA
      REAL*8 VELCHE(NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM)
      ENDSEGMENT
      SEGMENT LCHEVA
      POINTEUR LISCHE(NBCHE).MCHEVA
      ENDSEGMENT
CENDINCLUDE SMCHAEL
      INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM,N1
      POINTEUR FC.MCHEVA
      POINTEUR LCOF.LCHEVA
      POINTEUR MYCOF.MCHEVA
*
      INTEGER IMPR,IRET
*     
* Executable statements
*     
      IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans ccgdiv'
      IF (.NOT.(IDIM.GE.1.AND.IDIM.LE.3)) THEN
         WRITE(IOIMP,*) 'IDIM=',IDIM,' ?'
         GOTO 9999
      ENDIF
      NLFC=FC.VELCHE(/6)
      NPFC=FC.VELCHE(/5)
      DO IIDIM=1,IDIM
         MYCOF=LCOF.LISCHE(IIDIM)
         NLC=MYCOF.VELCHE(/6)
         NPC=MYCOF.VELCHE(/5)
         DO ILFC=1,NLFC
            IF (NLC.EQ.1) THEN
               ILC=1
            ELSE
               ILC=ILFC
            ENDIF
            DO IPFC=1,NPFC
               IF (NPC.EQ.1) THEN
                  IPC=1
               ELSE
                  IPC=IPFC
               ENDIF
               FC.VELCHE(1,1,1,1,IPFC,ILFC)=
     $              FC.VELCHE(1,1,1,1,IPFC,ILFC)+
     $              MYCOF.VELCHE(1,1,1,IIDIM,IPC,ILC)
            ENDDO
         ENDDO
      ENDDO
*     
* Normal termination
*     
      IRET=0
      RETURN
*     
* Format handling
*     
*     
* Error handling
*     
 9999 CONTINUE
      IRET=1
      WRITE(IOIMP,*) 'An error was detected in subroutine ccgdiv'
      RETURN
*     
* End of subroutine CCGDIV
*     
      END