Esope source files

   1 : C DCOGA1    SOURCE    GOUNAND   05/12/21    21:17:58     5281           
   2 :       SUBROUTINE DCOGA1(NDDL,IESREL,NBPOGO,NLDF,NLCO,NBELEM,
   3 :      $     JCOEFF,JDFFPG,
   4 :      $     JDCOFG,
   5 :      $     IMPR,IRET)
   6 :       IMPLICIT REAL*8 (A-H,O-Z)
   7 :       IMPLICIT INTEGER (I-N)
   8 : C***********************************************************************
   9 : C NOM         : DCOGA1
  10 : C PROJET      : 
  11 : C DESCRIPTION : Calcul des derivees du coefficient aux points de Gauss
  12 : C      d'un élément 
  13 : C               de référence pour chaque élément réel.
  14 : C             
  15 : C LANGAGE     : Fortran 77 (sauf E/S)
  16 : C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  17 : C               mél : gounand@semt2.smts.cea.fr
  18 : C***********************************************************************
  19 : C APPELES          : -
  20 : C APPELE PAR       : DCOGAU
  21 : C***********************************************************************
  22 : C ENTREES            :
  23 : C ENTREES/SORTIES    :
  24 : C SORTIES            : -
  25 : C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  26 : C***********************************************************************
  27 : C VERSION    : v3.1, 30/07/04, possiblité de travailler 
  28 : C                              dans l'espace de référence
  29 : C VERSION    : v1, 26/08/03, version initiale
  30 : C HISTORIQUE : v1, 26/08/03, création
  31 : C HISTORIQUE :
  32 : C HISTORIQUE :
  33 : C***********************************************************************
  34 : C Prière de PRENDRE LE TEMPS de compléter les commentaires
  35 : C en cas de modification de ce sous-programme afin de faciliter
  36 : C la maintenance !
  37 : C***********************************************************************
  38 : -INC CCOPTIO      
  39 :       INTEGER NDDL,IESREL,NBPOGO,NBELEM
  40 :       REAL*8 JCOEFF(NDDL,NLCO)
  41 :       REAL*8 JDFFPG(NDDL,IESREL,NBPOGO,NLDF)
  42 :       REAL*8 JDCOFG(IESREL,NBPOGO,NBELEM)
  43 : *
  44 :       INTEGER IMPR,IRET
  45 :       INTEGER IELEM,IPOGO,JESREL,IDDL
  46 : *     
  47 : * Executable statements
  48 : *
  49 :       IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans dcoga1'
  50 :       DO IELEM=1,NBELEM
  51 :          IF (NLCO.EQ.1) THEN
  52 :             ILCO=1
  53 :          ELSE
  54 :             ILCO=IELEM
  55 :          ENDIF
  56 :          IF (NLDF.EQ.1) THEN
  57 :             ILDF=1
  58 :          ELSE
  59 :             ILDF=IELEM
  60 :          ENDIF
  61 :          DO IPOGO=1,NBPOGO
  62 :             DO JESREL=1,IESREL
  63 :                DO IDDL=1,NDDL
  64 :                   JDCOFG(JESREL,IPOGO,IELEM)=
  65 :      $                 JDCOFG(JESREL,IPOGO,IELEM)
  66 :      $                 + (JCOEFF(IDDL,ILCO)
  67 :      $                    * JDFFPG(IDDL,JESREL,IPOGO,ILDF))
  68 :                ENDDO
  69 :             ENDDO
  70 :          ENDDO
  71 :       ENDDO
  72 : *
  73 : * Normal termination
  74 : *
  75 :       IRET=0
  76 :       RETURN
  77 : *
  78 : * Format handling
  79 : *
  80 : *
  81 : * Error handling
  82 : *
  83 :  9999 CONTINUE
  84 :       IRET=1
  85 :       WRITE(IOIMP,*) 'An error was detected in subroutine dcoga1'
  86 :       RETURN
  87 : *
  88 : * End of subroutine DCOGA1
  89 : *
  90 :       END
  91 : 
  92 :