Sources Esope | Cast3M

tadve
C TADVE     SOURCE    CB215821  22/04/01    14:03:15     11326          
C=========================================================
C=========================================================
C TADVE    SOURCE    JEANV     96/10/21    21:21:56     2334
      SUBROUTINE TADVE
************************************************************************
*
*                             T A D V E
*                             ---------
*
*     SE SOUS-PROGRAMME EST ASSOCIE A L'OPERATEUR "ADVE"
*
* FONCTION:
* ---------
*     CREATION DE LA MATRICE DE ADVECTION
*     DISPONIBLE POUR DES ELEMENTS MASSIFS A INTEGRATION NUMERIQUE
*
* PHRASE D'APPEL (EN GIBIANE):
* ----------------------------
*
*     ADV1 = ADVE MILIEU1 CHMAT1 ;
*
* OPERANDES ET RESULTAT:
* ----------------------
*
*     CND1     'RIGIDITE'  MATRICE DE ADVECTION
*     MILIEU1  'MODELE'    STRUCTURE MODELISEE
*     CHMAT1   'CHAMELEM'  DE SOUS-TYPE 'CARACTERISTIQUES' AVEC DEUX COMPOSANTES
*
* MODULES UTILISES:
* -----------------
*
      IMPLICIT INTEGER(I-N)
 
-INC PPARAM
-INC CCOPTIO
-INC SMCHAML
-INC SMCOORD
*
*
* VARIABLES:
* ----------
*
      CHARACTER*8 LETYPE
      CHARACTER*4 LISM(1)
      DATA LISM/'SYMM'/
*
* AUTEUR,DATE DE CREATION:
* ------------------------
*
*     MARINO ARROYO, 18 MAI 1999
*
* LANGAGE:
* --------
*
*       ESOPE + FORTRAN77
*
************************************************************************
*
      MOTERR(1:8)=' MODELE '
      CALL MESLIR(-137)
      CALL LIROBJ('MMODEL',IPMODE,1,IRET1)
      IF (IERR.NE.0) RETURN
      CALL ACTOBJ('MMODEL  ',IPMODE,1)
*
      CALL MESLIR(-135)
      CALL LIROBJ('MCHAML',IPIN,1,IRET1)
      IF (IERR.NE.0) RETURN
      CALL ACTOBJ('MCHAML  ',IPIN,1)
      CALL REDUAF(IPIN,IPMODE,MCHELM,0,IR,KER)
      IF(IR   .NE. 1) CALL ERREUR(KER)
      IF(IERR .NE. 0) RETURN
*
*     TEST SUR LE SOUS-TYPE DE LE CHAMELEM
*
      LETYPE = TITCHE
      IF (LETYPE.NE.'CARACTER') THEN
         MOTERR='CARACTERISTIQUES'
         CALL ERREUR(291)
         RETURN
      ENDIF
*
*     ON LIT SI LE MOT CLE DE SYMMETRIE EST PRESENT
*
      CALL LIRMOT(LISM,1,KIKI,0)
      IF (KIKI.EQ.1) THEN
         ISYMM = 1
      ELSE
         ISYMM = 0
      ENDIF
 
      SEGACT,MCOORD
      CALL TADVE1(IPMODE,MCHELM,IPRIGI,ISYMM)
      SEGDES,MCOORD
      IF (IERR.NE.0) RETURN
 
      CALL ECROBJ('RIGIDITE',IPRIGI)
      IF (IERR.NE.0) RETURN
 
      END