Sources Esope | Cast3M

cogauf
C COGAUF    SOURCE    GOUNAND   21/06/02    21:15:29     11022          
      SUBROUTINE COGAUF(JCOEFF,FFPG,SSFACT,NBELEF,
     $     JCOEFG,
     $     IMPR,IRET)
      IMPLICIT REAL*8 (A-H,O-Z)
      IMPLICIT INTEGER (I-N)
C***********************************************************************
C NOM         : COGAUF
C PROJET      : Noyau linéaire NLIN
C DESCRIPTION :
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
C ENTREES/SORTIES    : -
C SORTIES            :
C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
C***********************************************************************
C VERSION    : v2, 03/10/03, refonte complète (modif SMTNLIN)
C VERSION    : v1, 17/01/03, version initiale
C HISTORIQUE : v1, 17/01/03, 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 PPARAM
-INC CCOPTIO
-INC TNLIN      
*-INC SFACTIV
*-INC SMCHAEL
      INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM
      POINTEUR JCOEFF.MCHEVA
      POINTEUR JCOEFG.MCHEVA
* Valeurs des fns d'interpolation du coeff. aux points de Gauss
      POINTEUR FFPG.MCHEVA
*
      INTEGER IMPR,IRET
*
* Executable statements
*
      IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans cogauf'
      SEGACT SSFACT
      NBELFV=SSFACT.LFACTI(/1)
      NBELEV=SSFACT.LFACTI(/2)
      SEGACT JCOEFF
      NDLIG =JCOEFF.WELCHE(/1)
      NDCOL =JCOEFF.WELCHE(/2)
      N2DLIG=JCOEFF.WELCHE(/3)
      N2DCOL=JCOEFF.WELCHE(/4)
      NDNOEU=JCOEFF.WELCHE(/5)
      NBELEM=JCOEFF.WELCHE(/6)
      IF (NDLIG.NE.1.OR.N2DLIG.NE.1.OR.N2DCOL.NE.1
     $     .OR.NDNOEU.NE.1
     $     .OR.(NBELEM.NE.1.AND.NBELEM.NE.NBELEV)) THEN
         WRITE(IOIMP,*) 'Erreur dims JCOEFF'
         GOTO 9999
      ENDIF
      NDDL=NDCOL
      NLVCOF=NBELEM
      SEGACT FFPG
      NDLIG =FFPG.WELCHE(/1)
      NDCOL =FFPG.WELCHE(/2)
      N2DLIG=FFPG.WELCHE(/3)
      N2DCOL=FFPG.WELCHE(/4)
      NDNOEU=FFPG.WELCHE(/5)
      NBELEM=FFPG.WELCHE(/6)
      IF (NDLIG.NE.1.OR.NDCOL.NE.NDDL
     $     .OR.N2DLIG.NE.1.OR.N2DCOL.NE.1
     $     .OR.(NBELEM.NE.1.AND.NBELEM.NE.NBELFV)) THEN
         WRITE(IOIMP,*) 'Erreur dims FFPG'
         GOTO 9999
      ENDIF
      NBPOGO=NDNOEU
      NLFVFF=NBELEM
      IF (NLVCOF.EQ.1.AND.NLFVFF.EQ.1) THEN
         NLFCOG=1
      ELSE
         NLFCOG=NBELEF
      ENDIF
*
* Initialisations...
*
      NBLIG=1
      NBCOL=1
      N2LIG=1
      N2COL=1
      NBPOI=NBPOGO
      NBELM=NLFCOG
      SEGINI JCOEFG
*
* On effectue le calcul du coefficient aux points de Gauss
*
      CALL COGAF1(NDDL,NBPOGO,NBELEV,NBELFV,NBELEF,
     $     NLVCOF,NLFVFF,NLFCOG,
     $     JCOEFF.WELCHE,FFPG.WELCHE,SSFACT.LFACTI,
     $     JCOEFG.WELCHE,
     $     IMPR,IRET)
      IF (IRET.NE.0) GOTO 9999
      SEGDES JCOEFG
      SEGDES SSFACT
      SEGDES JCOEFF
      SEGDES FFPG
*
* Normal termination
*
      IRET=0
      RETURN
*
* Format handling
*
*
* Error handling
*
 9999 CONTINUE
      IRET=1
      WRITE(IOIMP,*) 'An error was detected in subroutine cogauf'
      RETURN
*
* End of subroutine COGAUF
*
      END