Esope source files

   1 : C CCGDIV    SOURCE    GOUNAND   05/12/21    21:16:26     5281           
   2 :       SUBROUTINE CCGDIV(LCOF,
   3 :      $     FC,
   4 :      $     IMPR,IRET)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 :       IMPLICIT INTEGER (I-N)
   7 : C***********************************************************************
   8 : C NOM         : CCGDIV
   9 : C DESCRIPTION : Calcul de la loi de comportement aux points de Gauss : 
  10 : C               une divergence
  11 : C
  12 : C     
  13 : C LANGAGE     : ESOPE
  14 : C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  15 : C               mél : gounand@semt2.smts.cea.fr
  16 : C***********************************************************************
  17 : C APPELES          : 
  18 : C APPELE PAR       : 
  19 : C***********************************************************************
  20 : C ENTREES            :
  21 : C ENTREES/SORTIES    :
  22 : C SORTIES            : -
  23 : C TRAVAIL            :
  24 : C***********************************************************************
  25 : C VERSION    : v1, 11/05/04, version initiale
  26 : C HISTORIQUE : v1, 11/05/04, création
  27 : C HISTORIQUE :
  28 : C HISTORIQUE :
  29 : C***********************************************************************
  30 : C Prière de PRENDRE LE TEMPS de compléter les commentaires
  31 : C en cas de modification de ce sous-programme afin de faciliter
  32 : C la maintenance !
  33 : C***********************************************************************
  34 : -INC CCOPTIO
  35 : CBEGININCLUDE SMCHAEL
  36 :       SEGMENT MCHAEL
  37 :       POINTEUR IMACHE(N1).MELEME
  38 :       POINTEUR ICHEVA(N1).MCHEVA
  39 :       ENDSEGMENT
  40 :       SEGMENT MCHEVA
  41 :       REAL*8 VELCHE(NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM)
  42 :       ENDSEGMENT
  43 :       SEGMENT LCHEVA
  44 :       POINTEUR LISCHE(NBCHE).MCHEVA
  45 :       ENDSEGMENT
  46 : CENDINCLUDE SMCHAEL
  47 :       INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM,N1
  48 :       POINTEUR FC.MCHEVA
  49 :       POINTEUR LCOF.LCHEVA
  50 :       POINTEUR MYCOF.MCHEVA
  51 : *
  52 :       INTEGER IMPR,IRET
  53 : *     
  54 : * Executable statements
  55 : *     
  56 :       IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans ccgdiv'
  57 :       IF (.NOT.(IDIM.GE.1.AND.IDIM.LE.3)) THEN
  58 :          WRITE(IOIMP,*) 'IDIM=',IDIM,' ?'
  59 :          GOTO 9999
  60 :       ENDIF
  61 :       NLFC=FC.VELCHE(/6)
  62 :       NPFC=FC.VELCHE(/5)
  63 :       DO IIDIM=1,IDIM
  64 :          MYCOF=LCOF.LISCHE(IIDIM)
  65 :          NLC=MYCOF.VELCHE(/6)
  66 :          NPC=MYCOF.VELCHE(/5)
  67 :          DO ILFC=1,NLFC
  68 :             IF (NLC.EQ.1) THEN
  69 :                ILC=1
  70 :             ELSE
  71 :                ILC=ILFC
  72 :             ENDIF
  73 :             DO IPFC=1,NPFC
  74 :                IF (NPC.EQ.1) THEN
  75 :                   IPC=1
  76 :                ELSE
  77 :                   IPC=IPFC
  78 :                ENDIF
  79 :                FC.VELCHE(1,1,1,1,IPFC,ILFC)=
  80 :      $              FC.VELCHE(1,1,1,1,IPFC,ILFC)+
  81 :      $              MYCOF.VELCHE(1,1,1,IIDIM,IPC,ILC)
  82 :             ENDDO
  83 :          ENDDO
  84 :       ENDDO
  85 : *     
  86 : * Normal termination
  87 : *     
  88 :       IRET=0
  89 :       RETURN
  90 : *     
  91 : * Format handling
  92 : *     
  93 : *     
  94 : * Error handling
  95 : *     
  96 :  9999 CONTINUE
  97 :       IRET=1
  98 :       WRITE(IOIMP,*) 'An error was detected in subroutine ccgdiv'
  99 :       RETURN
 100 : *     
 101 : * End of subroutine CCGDIV
 102 : *     
 103 :       END
 104 : 
 105 :