Télécharger idmatr.eso

Retour à la liste

idmatr
C IDMATR    SOURCE    OF166741  26/02/23    21:15:15     12480          
*--------------------------------------------------------------------*
*       Noms de composantes de materiaux                             *
*--------------------------------------------------------------------*
*                                                                    *
*  ENTREES:                                                          *
*     MFR      Numero de formulation                                 *
*     IPMODL   objet modele elementaire ( segment actif )            *
*                                                                    *
*  SORTIES:                                                          *
*     IPNOMC   pointeur sur les listes de noms de composantes        *
*              obligatoires et facultatives                          *
*     NBROBL   nombre de composantes obligatoires                    *
*     NBRFAC   nombre de composantes facultatives                    *
*                                                                    *
*                                                                    *
*  Remarque : Voir INOMID ou CCOPTIO pour signification IFOUR        *
*--------------------------------------------------------------------*
C
      SUBROUTINE IDMATR(MFR,IPMODL,IPNOMC,NBROBL,NBRFAC)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
-INC PPARAM
-INC CCOPTIO
C==DEB= FORMULATION HHO == Include specifique ==========================
-INC CCHHOPA
C==FIN= FORMULATION HHO ================================================
C
-INC SMLMOTS
-INC SMMODEL
      POINTEUR NOMID1.NOMID
 
******Modif SELLIER ****************************************************
*      augmentation nombre de caracteristiques      
      PARAMETER (ITA=500)
      CHARACTER*16 MOMODL(500)
c     nombre de variables et de parametres par variables
c      include './nombre_helmholtz.h' 
-INC HNBRHEL   
      CHARACTER*1 motnl1,motnlv1,motnld1,motnla1
      CHARACTER*2 motnl2
      CHARACTER*3 motnl3 
      CHARACTER*4 motnl4    
******fin modif SELLIER ************************************************
 
      CHARACTER*(LOCOMP) TABOBL(ITA),TABFAC(ITA)
      LOGICAL bDARCY
 
      IMODEL=IPMODL
      NOMID =IMODEL.LNOMID(6)
C
C     Ne pas recreer le NOMID si deja present
      IF (NOMID.NE.0) THEN
        NBROBL=LESOBL(/2)
        NBRFAC=LESFAC(/2)
        IPNOMC=NOMID
        RETURN
      ENDIF
C
      JGOBL=0
      JGFAC=0
      IPNOMC=0
      IRET = 1
C
      NMAT=MATMOD(/2)
      NFOR=FORMOD(/2)
      MELE=NEFMOD
      MFR2=NUMFOR(IMODEL)
C
C     =================================================================
C     FORMULATION THERMOHYDRIQUE
C     =================================================================
      IF (MFR2.EQ.65) THEN
C
        IF (CMATEE.EQ.'SCHREFLE') THEN
          JGOBL = 20
          TABOBL(1)='KGG'
          TABOBL(2)='KGC'
          TABOBL(3)='KGT'
          TABOBL(4)='KCG'
          TABOBL(5)='KCC'
          TABOBL(6)='KCT'
          TABOBL(7)='KTG'
          TABOBL(8)='KTC'
          TABOBL(9)='KTT'
          TABOBL(10)='CGG'
          TABOBL(11)='CGC'
          TABOBL(12)='CGT'
          TABOBL(13)='CCG'
          TABOBL(14)='CCC'
          TABOBL(15)='CCT'
          TABOBL(16)='CTG'
          TABOBL(17)='CTC'
          TABOBL(18)='CTT'
          TABOBL(19)='KTGG'
          TABOBL(20)='KTCG'
        ELSEIF (CMATEE.EQ.'BETON_TH') THEN
          JGOBL = 44
          TABOBL(1)='KGG'
          TABOBL(2)='KGC'
          TABOBL(3)='KGT'
          TABOBL(4)='KCG'
          TABOBL(5)='KCC'
          TABOBL(6)='KCT'
          TABOBL(7)='KTG'
          TABOBL(8)='KTC'
          TABOBL(9)='KTT'
          TABOBL(10)='CGG'
          TABOBL(11)='CGC'
          TABOBL(12)='CGT'
          TABOBL(13)='CCG'
          TABOBL(14)='CCC'
          TABOBL(15)='CCT'
          TABOBL(16)='CTG'
          TABOBL(17)='CTC'
          TABOBL(18)='CTT'
          TABOBL(19)='KTGG'
          TABOBL(20)='KTCG'
          TABOBL(21)='CS'
          TABOBL(22)='KS'
          TABOBL(23)='POR1'
          TABOBL(24)='APOR'
          TABOBL(25)='VGA'
          TABOBL(26)='VGB'
          TABOBL(27)='VGC'
          TABOBL(28)='GAM0'
          TABOBL(29)='K1'
          TABOBL(30)='AK'                  
          TABOBL(31)='BK'
          TABOBL(32)='KLI1'
          TABOBL(33)='AKLI'                
          TABOBL(34)='BKLI'                
          TABOBL(35)='EAR'
          TABOBL(36)='QLAT'
          TABOBL(37)='AFF1'
          TABOBL(38)='AFF2'
          TABOBL(39)='AFF3'
          TABOBL(40)='AFF4'
          TABOBL(41)='WINF'   
          TABOBL(42)='DHYD'
          TABOBL(43)='DDCH'
          TABOBL(44)='DSWE'
        ENDIF
C     =================================================================
C     FORMULATION THERMIQUE/DARCY
C     =================================================================
      ELSE IF ((MFR2.EQ.29) .OR. (MFR2.EQ.34)) THEN
C
        bDARCY=(MFR2.EQ.34)
C
C       MATERIAU THERMIQUE CONDUCTION OU MATERIAU DARCY
C       ===============================================
        CALL PLACE(MATMOD,NMAT,ICONDU,'CONDUCTION')
        IF (ICONDU.NE.0.OR.bDARCY) THEN
          IF (MFR.EQ.75) THEN
            JGFAC    = 4
            TABFAC(1)='M'
            TABFAC(2)='C'
            TABFAC(3)='TINI'
          ELSE
            JGFAC    = 4
            TABFAC(1)='RHO '
            TABFAC(2)='C   '
            TABFAC(3)='H   '
            TABFAC(4)='TINI'
          ENDIF
C
C         Isotrope
          IF (IMATEE.EQ.1) THEN
            JGOBL = 1
            IF(MFR .EQ. 75)THEN
              TABOBL(1)='KT'
            ELSE
              TABOBL(1)='K   '
            ENDIF
C
C         Orthotrope
          ELSEIF (IMATEE.EQ.2) THEN
            IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
              JGOBL=5
              TABOBL(1)='K1  '
              TABOBL(2)='K2  '
              TABOBL(3)='K3  '
              TABOBL(4)='V1X '
              TABOBL(5)='V1Y '
            ELSEIF (MFR.EQ.1.OR.MFR.EQ.45) THEN
              IF (IDIM.EQ.2) THEN
                IF (IFOMOD.EQ.1) THEN
                  JGOBL=5
                  TABOBL(1)='K1  '
                  TABOBL(2)='K2  '
                  TABOBL(3)='V1X '
                  TABOBL(4)='V1Y '
                  TABOBL(5)='K3  '
                ELSE
                  JGOBL=4
                  TABOBL(1)='K1  '
                  TABOBL(2)='K2  '
                  TABOBL(3)='V1X '
                  TABOBL(4)='V1Y '
                ENDIF
              ELSEIF (IDIM.EQ.3) THEN
                JGOBL=9
                TABOBL(1)='K1  '
                TABOBL(2)='K2  '
                TABOBL(3)='K3  '
                TABOBL(4)='V1X '
                TABOBL(5)='V1Y '
                TABOBL(6)='V1Z '
                TABOBL(7)='V2X '
                TABOBL(8)='V2Y '
                TABOBL(9)='V2Z '
              ENDIF
            ENDIF
C
C         Anisotrope
          ELSEIF (IMATEE.EQ.3) THEN
            IF (MFR.EQ.1.OR.MFR.EQ.45) THEN
              IF (IDIM.EQ.2) THEN
                IF (IFOMOD.EQ.1) THEN
                  JGOBL=6
                  TABOBL(1)='K11 '
                  TABOBL(2)='K22 '
                  TABOBL(3)='K21 '
                  TABOBL(4)='V1X '
                  TABOBL(5)='V1Y '
                  TABOBL(6)='K33 '
                ELSE
                  JGOBL=5
                  TABOBL(1)='K11 '
                  TABOBL(2)='K22 '
                  TABOBL(3)='K21 '
                  TABOBL(4)='V1X '
                  TABOBL(5)='V1Y '
                ENDIF
              ELSEIF (IDIM.EQ.3) THEN
                JGOBL=12
                TABOBL(1)='K11 '
                TABOBL(2)='K22 '
                TABOBL(3)='K33 '
                TABOBL(4)='K21 '
                TABOBL(5)='K31 '
                TABOBL(6)='K32 '
                TABOBL(7)='V1X '
                TABOBL(8)='V1Y '
                TABOBL(9)='V1Z '
                TABOBL(10)='V2X '
                TABOBL(11)='V2Y '
                TABOBL(12)='V2Z '
              ENDIF
            ENDIF
          ENDIF
        ENDIF
C
        IF (bDARCY) GOTO 9999
C
C       ===============================================================
C       MATERIAU THERMIQUE PHASE
C       ===============================================================
        CALL PLACE(MATMOD,NMAT,IPHA,'PHASE')
        IF (IPHA.NE.0) THEN
          JGOBL0 = JGOBL
          JGOBL  = JGOBL0+4
          TABOBL(JGOBL0+1)='RHO '
          TABOBL(JGOBL0+2)='C   '
          TABOBL(JGOBL0+3)='QLAT'
          TABOBL(JGOBL0+4)='TPHA'
 
          JGFAC=1
          TABFAC(1)='H   '
        ENDIF
C       ===============================================================
C       MATERIAU THERMIQUE CONVECTION
C       ===============================================================
        IF (CMATEE.EQ.'CONVECTI') THEN
C         IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
CC          Coques thermiques
C           JGOBL    =2
C           TABOBL(1)='HINF'
C           TABOBL(2)='HSUP'
C           JGFAC    =2
C           TABFAC(1)='TCINF'
C           TABFAC(2)='TCSUP'
C         else
C           Elements massifs
            JGOBL = 1
            TABOBL(1)='H'
            JGFAC=1
            TABFAC(1)='TC'
C         endif
        ENDIF
C       ===============================================================
C       MATERIAU THERMIQUE RAYONNEMENT
C       ===============================================================
        IF (CMATEE.EQ.'RAYONNEM') THEN
C
C         RAYONNEMENT en CAVITE
          IF (IMATEE.EQ.4) THEN
            IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
              JGOBL=2
              TABOBL(1)='EINF'
              TABOBL(2)='ESUP'
            ELSE
              JGOBL=1
              TABOBL(1)='EMIS'
            ENDIF
 
            JGFAC=3
            TABFAC(1)='CABS'
            TABFAC(2)='TABS'
            TABFAC(3)='H'
 
C         RAYONNEMENT FACE_A_FACE
          ELSEIF (IMATEE.EQ.2) THEN
            JGOBL=1
            TABOBL(1)='EMIS'
            JGFAC = 1
            TABFAC(1)='H'
 
C         RAYONNEMENT a l'INFINI
          ELSEIF (IMATEE.EQ.3) THEN
            IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
              JGOBL=3
              TABOBL(1)='EINF'
              TABOBL(2)='ESUP'
              TABOBL(3)='T_IN'
            ELSE
              JGOBL=2
              TABOBL(1)='EMIS'
              TABOBL(2)='T_IN'
            ENDIF
 
            JGFAC = 2
            TABFAC(1)='E_IN'
            TABFAC(2)='H'
          ELSE
            IRET = 0
            CALL ERREUR(5)
          ENDIF
        ENDIF
C       ===============================================================
C       MATERIAU THERMIQUE ADVECTION
C       ===============================================================
        CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
        IF (IADVE.NE.0) THEN
C         Cas des Tuyaux 1D (MFR=79)
          IF (MFR .EQ. 79) THEN
            JGFAC   = 3
            TABFAC(1)='RHO '
            TABFAC(2)='C   '
            TABFAC(3)='VITE'
 
C         Cas Massif (MFR=1)
          ELSEIF (MFR.EQ.1) THEN
            IF (IFOUR.EQ.-3) THEN
              JGFAC   = 5
              TABFAC(1)='RHO '
              TABFAC(2)='C   '
              TABFAC(3)='VITX'
              TABFAC(4)='VITY'
              TABFAC(5)='VITZ'
            ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
              JGFAC   = 4
              TABFAC(1)='RHO '
              TABFAC(2)='C   '
              TABFAC(3)='VITX'
              TABFAC(4)='VITY'
            ELSEIF (IFOUR .EQ. 0) THEN
              JGFAC   = 4
              TABFAC(1)='RHO '
              TABFAC(2)='C   '
              TABFAC(3)='VITR'
              TABFAC(4)='VITZ'
            ELSEIF (IFOUR.EQ.1) THEN
              CALL ERREUR(21)
              RETURN
C             JGFAC   = 5C
C             TABFAC(1)='RHO '
C             TABFAC(2)='C   '
C             TABFAC(3)='VITR'
C             TABFAC(4)='VITZ'
C             TABFAC(5)='VITT'
            ELSEIF (IFOUR.EQ.2) THEN
              JGFAC   = 5
              TABFAC(1)='RHO '
              TABFAC(2)='C   '
              TABFAC(3)='VITX'
              TABFAC(4)='VITY'
              TABFAC(5)='VITZ'
            ELSEIF (IFOUR.GE.3.AND.IFOUR.LE.15) THEN
              IF (IFOUR.LE.6) THEN
                JGFAC   = 3
                TABFAC(1)='RHO '
                TABFAC(2)='C   '
                TABFAC(3)='VITX'
              ELSEIF (IFOUR.EQ.7.OR.IFOUR.EQ.8) THEN
                JGFAC   = 4
C               Verifier l'utilite des composantes au dela de 4
                TABFAC(1)='RHO '
                TABFAC(2)='C   '
                TABFAC(3)='VITX'
                TABFAC(4)='VITY'
              ELSEIF (IFOUR.EQ.9.OR.IFOUR.EQ.10) THEN
                JGFAC   = 4
C               Verifier l'utilite des composantes au dela de 4
                TABFAC(1)='RHO '
                TABFAC(2)='C   '
                TABFAC(3)='VITX'
                TABFAC(4)='VITZ'
              ELSEIF (IFOUR.EQ.11) THEN
                JGFAC   = 5
C               Verifier l'utilite des composantes au dela de 4
                TABFAC(1)='RHO '
                TABFAC(2)='C   '
                TABFAC(3)='VITX'
                TABFAC(4)='VITY'
                TABFAC(5)='VITZ'
             ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN
                JGFAC   = 3
                TABFAC(1)='RHO '
                TABFAC(2)='C   '
                TABFAC(3)='VITR'
              ELSEIF (IFOUR.EQ.14) THEN
                JGFAC   = 4
C               Verifier l'utilite des composantes au dela de 4
                TABFAC(1)='RHO '
                TABFAC(2)='C   '
                TABFAC(3)='VITR'
                TABFAC(4)='VITZ'
              ENDIF
            ENDIF
          ENDIF
        ENDIF
C       ===============================================================
C       MATERIAU THERMIQUE SOURCE
C       ===============================================================
        IF (CMATEE.EQ.'SOURCE') THEN
          IF (INATUU.EQ.1.AND.IMATEE.EQ.1) THEN
            JGOBL = 1
            TABOBL(1)='QVOL'
            JGFAC = 2
            TABFAC(1)='QINF'
            TABFAC(2)='QSUP'
          ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.1) THEN
            JGOBL = 3
            TABOBL(1)='QTOT'
            TABOBL(2)='ORIG'
            TABOBL(3)='RGAU'
          ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.2) THEN
            JGOBL = 5
            TABOBL(1)='QTOT'
            TABOBL(2)='ORIG'
            TABOBL(3)='RGAU'
            TABOBL(4)='DIRE'
            TABOBL(5)='ZGAU'
          ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.3) THEN
            JGOBL = 7
            TABOBL(1)='QTOT'
            TABOBL(2)='ORIG'
            TABOBL(3)='RGAU'
            TABOBL(4)='DIRE'
            TABOBL(5)='ZGAU'
            TABOBL(6)='DIRL'
            TABOBL(7)='LGAU'
          ELSE
            WRITE(6,*) ' Dans IDMATR : numero IMATEE non prevu'
            CALL ERREUR(21)
            RETURN
          ENDIF
        ENDIF
C     =================================================================
C     FORMULATION LIQUIDE/MECANIQUE+LIQUIDE
C     =================================================================
      ELSE IF ((MFR2.EQ.11) .OR. (MFR2.EQ.44)) THEN
C
        IF (MFR.NE.41.AND.MFR.NE.43) THEN
          JGOBL = 6
          TABOBL(1)='RHO '
          TABOBL(2)='CSON'
          TABOBL(3)='RORF'
          TABOBL(4)='CREF'
          TABOBL(5)='LCAR'
          TABOBL(6)='G   '
        ELSEIF (MFR.EQ.41) THEN
          JGOBL = 5
          TABOBL(1)='RHO '
          TABOBL(2)='CSON'
          TABOBL(3)='RORF'
          TABOBL(4)='CREF'
          TABOBL(5)='LCAR'
        ELSE
          JGOBL = 3
          TABOBL(1)='RHO '
          TABOBL(2)='LCAR'
          TABOBL(3)='RORF'
        ENDIF
C     =================================================================
C     FORMULATION MECANIQUE
C     =================================================================
      ELSE IF (MFR2.EQ.2) THEN
 
C== FORMULATION HHO == MECANIQUE = idem qu'en MASSIF ==> MFR = 1 =======
C== Les composantes specifiques seront ajoutees dans HHOIDC via INOMID.
        IF (MFR.EQ.HHO_MFR_ELEMENT) MFR = 1
C== FORMULATION HHO ====================================================
C
        CALL MODLIN(MOMODL,NMOD)
        CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
        IF (IRET.EQ.0) GOTO 9999
C
        IF (NMAT.LT.2) GOTO 13
C
        CALL MODELA(MOMODL,NMOD)
        CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
C
        INMAT=3
        IF (IPLAC.EQ.0) INMAT=2
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE IMPEDANCE
C       =============================================================
        CALL PLACE(MATMOD,NMAT,IIMPE,'IMPEDANCE')
        IF (IIMPE.EQ.0) GOTO 13
C
        INMAT = INMAT + 1
        JGOBL = 0
        JGFAC = 0
 
        IF(CMATEE.EQ.'IMPELAST') THEN
          JGOBL = 1
          TABOBL(1)='RAID'
 
          JGFAC = 7
          TABFAC(1)='AMOR'
          TABFAC(2)='CPLE'
          TABFAC(3)='INER'
          TABFAC(4)='AROT'
          TABFAC(5)='ZNU'
          TABFAC(6)='MASS'
          TABFAC(7)='ALPH'
        ELSEIF (CMATEE.EQ.'IMPVOIGT'.or.CMATEE.EQ.'IMPREUSS') THEN
          JGOBL = 2
          TABOBL(1)='RAID'
          TABOBL(2)='VISC'
 
          JGFAC = 5
          TABFAC(1)='MASS'
          TABFAC(2)='AMOR'
          TABFAC(3)='CPLE'
          TABFAC(4)='INER'
          TABFAC(5)='AROT'
        ELSEIF (CMATEE.EQ.'IMPCOMPL') then
          JGOBL = 1
          TABOBL(1)='RAID'
 
          JGFAC = 7
          TABFAC(1)='MOCO'
          TABFAC(2)='VISC'
          TABFAC(3)='MASS'
          TABFAC(4)='AMOR'
          TABFAC(5)='CPLE'
          TABFAC(6)='INER'
          TABFAC(7)='AROT'
        ENDIF
        GOTO 9999
C       =============================================================
 13     CONTINUE
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE ISOTROPE (PAR DEFAUT)
C       =============================================================
        IF (IMATEE.EQ.1) THEN
          IF (MFR.EQ.35) THEN
            JGOBL=2
            TABOBL(1)='KS  '
            TABOBL(2)='KN  '
 
            JGFAC=2
            TABFAC(1)='RHO '
            TABFAC(2)='ALPN'
          ELSEIF (MFR.EQ.78) THEN
            JGOBL=2
            TABOBL(1)='KS  '
            TABOBL(2)='KN  '
          ELSEIF (MFR.EQ.53) THEN
            JGOBL=1
            TABOBL(1)='KS  '
*
            JGFAC=2
            TABFAC(1)='RHO '
            TABFAC(2)='ALPN'
          ELSE
            JGOBL=2
            TABOBL(1)='YOUN'
            TABOBL(2)='NU  '
*
            IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.63).AND.
     &       IFOUR.EQ.-2) THEN
              JGFAC=4
              TABFAC(3)='DIM3'
              TABFAC(4)='VISQ'
            ELSE
              JGFAC=3
              TABFAC(3)='VISQ'
            ENDIF
            TABFAC(1)='RHO '
            TABFAC(2)='ALPH'
 
* AM  CAS FORMULATION NON-LOCALE
 
            INLOC = 0
            MN3 = INFMOD(/1)
            IF (MN3.GE.13) INLOC=-1*INFMOD(13)
            IF (INLOC.GT.0) THEN
              JGOBL0 = JGFAC
*  moyenne
              IF(INLOC.EQ.1) THEN
                JGFAC=JGOBL0+1
                TABFAC(JGOBL0+1)='LCAR'
*  stress-based
              ELSEIF(INLOC.EQ.2) THEN
                JGFAC=JGOBL0+2
                TABFAC(JGOBL0+1)='LCAR'
                TABFAC(JGOBL0+2)='SBFT'
 
*  helmholtz
              ELSE IF(INLOC.EQ.3) THEN
*  sellier 03/07/20
*                JGFAC=JGOBL0+6
*                TABFAC(JGOBL0+1)='LCAR'
*                TABFAC(JGOBL0+2)='LCF1'
*                TABFAC(JGOBL0+3)='LCF2'
*                TABFAC(JGOBL0+4)='LCF3'
*                TABFAC(JGOBL0+5)='LCF4'
*                TABFAC(JGOBL0+6)='LCF5'       
                JGFAC=JGOBL0+NB_PARA_HELM
c               variables de Helmholtz a diffusions ORThotropes
c               modif liste para heilmholtz sellier 10//11/22
c               print*,'Dans IDMATR ', NB_HELM, 'formulationss Helmholtz'
                do I=1,NB_HELM
c                   numero de la variable de Helmholtz                 
                    write (motnl1,'(I1)') I
c                   Capacité (infinie si BH=1)           
                    motnl3='CAP' 
                    motnl4=motnl3//motnl1 
                    TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+1)=motnl4
c                   print*,motnl4
c                   Blocage de Dirichlet si 1 sinon 0           
                    motnl3='BLO' 
                    motnl4=motnl3//motnl1          
                    TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+2)=motnl4
c                   print*,motnl4
c                   Valeur imposee a la zone de blocage de Dirichlet             
                    motnl3='DEP' 
                    motnl4=motnl3//motnl1          
                    TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+3)=motnl4
c                   print*,motnl4
c                   Valeur initiale de la variable de Helmholtz
                    motnl3='INI' 
                    motnl4=motnl3//motnl1          
                    TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+4)=motnl4
c                   print*,motnl4
c                   indicateur logique de lineratite pour ne pas passer par istep=3 si 1
                    motnl3='LIN' 
                    motnl4=motnl3//motnl1          
                    TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+5)=motnl4
c                   print*,motnl4
                    do J=1,3
c                     coefficients de diffusion                  
                      write (motnld1,'(I1)') J
c                     diffusion 
                      motnl2='DH'
c                     1er indice: numero de la variable d helmholtz
                      motnl3=motnl2//motnl1
c                     2eme indice : direction 
                      motnl4=motnl3//motnld1
                      TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
     #                5+(J-1)*4+1)=motnl4 
c                     print*,motnl4               
                      do K=1,3
                          write (motnla1,'(I1)') K              
c                         direction principales de diffusion 
                          motnlv1='V'
c                         indice numero de la variable
                          motnl2=motnlv1//motnl1
c                         indice direction principale
                          motnl3=motnl2//motnld1
c                         projection sur axe indice de base fixe
                          motnl4=motnl3//motnla1
                          TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
     #                    5+(J-1)*4+1+K)=motnl4
c                         print*,motnl4
                      end do
                    end do
                end do
Cc              on rajoute 2 inutiles pour compenser un decalage constatte dans castem20               
C                 TABFAC(JGOBL0+NB_PARA_HELM+1)='FAC1'
C                 TABFAC(JGOBL0+NB_PARA_HELM+2)='FAC2'              
*  fin modif sellier 
*  cas non prevu
              ELSE
                IRET = 0
                CALL ERREUR(5)
              ENDIF
            ENDIF
          ENDIF
C
          IF (NMAT.LT.2) GOTO 9999
C
       CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
       IF (IADVE.NE.0) THEN
        JG0FAC = JGFAC
C       Cas des Tuyaux 1D (MFR=79)
        IF (MFR .EQ. 79) THEN
          JGFAC   = JG0FAC + 3
          TABFAC(JG0FAC +1)='RHO '
          TABFAC(JG0FAC +2)='C   '
          TABFAC(JG0FAC +3)='VITE'
 
C       Cas Massif (MFR=1)
        ELSEIF (MFR.EQ.1) THEN
C         Bidimensionnel PLAN DPGE
          IF (IFOUR.EQ.-3) THEN
            JGFAC   = JG0FAC +5
            TABFAC(JG0FAC +1)='RHO '
            TABFAC(JG0FAC +2)='C   '
            TABFAC(JG0FAC +3)='VITX'
            TABFAC(JG0FAC +4)='VITY'
            TABFAC(JG0FAC +5)='VITZ'
 
C         Bidimensionnel PLAN (CP/DP)
          ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
            JGFAC   = JG0FAC +4
            TABFAC(JG0FAC +1)='RHO '
            TABFAC(JG0FAC +2)='C   '
            TABFAC(JG0FAC +3)='VITX'
            TABFAC(JG0FAC +4)='VITY'
 
C         Axisymetrie
          ELSEIF (IFOUR .EQ. 0) THEN
            JGFAC   = JG0FAC +4
            TABFAC(JG0FAC +1)='RHO '
            TABFAC(JG0FAC +2)='C   '
            TABFAC(JG0FAC +3)='VITR'
            TABFAC(JG0FAC +4)='VITZ'
 
C         Fourier
          ELSEIF (IFOUR.EQ.1) THEN
            CALL ERREUR(21)
            RETURN
C           JGFAC   = 5C
C           TABFAC(1)='RHO '
C           TABFAC(2)='C   '
C           TABFAC(3)='VITR'
C           TABFAC(4)='VITZ'
C           TABFAC(5)='VITT'
 
C         Tridimensionnel
          ELSEIF (IFOUR.EQ.2) THEN
            JGFAC   = JG0FAC +5
            TABFAC(JG0FAC +1)='RHO '
            TABFAC(JG0FAC +2)='C   '
            TABFAC(JG0FAC +3)='VITX'
            TABFAC(JG0FAC +4)='VITY'
            TABFAC(JG0FAC +5)='VITZ'
 
C         Unidimensionnel (1D)
          ELSE
          ENDIF
        ENDIF
        GOTO 9999
       ENDIF
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE ORTHOTROPE
C       =============================================================
        ELSEIF (IMATEE.EQ.2) THEN
C
*         JOINT UNIDIMENSIONNEL JOI1
          IF (MFR.EQ.75) THEN
            IF(IFOUR.EQ.2)THEN
              JGOBL=12
              TABOBL(1)='V1X '
              TABOBL(2)='V1Y '
              TABOBL(3)='V1Z '
              TABOBL(4)='V2X '
              TABOBL(5)='V2Y '
              TABOBL(6)='V2Z '
              TABOBL(7)='KN  '
              TABOBL(8)='KS1 '
              TABOBL(9)='KS2 '
              TABOBL(10)='QN  '
              TABOBL(11)='QS1 '
              TABOBL(12)='QS2 '
*
              JGFAC=10
              TABFAC(1)='MASS'
              TABFAC(2)='JX  '
              TABFAC(3)='JY  '
              TABFAC(4)='JZ  '
              TABFAC(5)='ALPN'
              TABFAC(6)='ALP1'
              TABFAC(7)='ALP2'
              TABFAC(8)='ALQN'
              TABFAC(9)='ALQ1'
              TABFAC(10)='ALQ2'
            ELSEIF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1)THEN
              JGOBL=5
              TABOBL(1)='V1X '
              TABOBL(2)='V1Y '
              TABOBL(3)='KN  '
              TABOBL(4)='KS  '
              TABOBL(5)='QS  '
*
              JGFAC=6
              TABFAC(1)='MASS'
              TABFAC(2)='JZ'
              TABFAC(4)='ALPN'
              TABFAC(5)='ALPS'
              TABFAC(6)='ALQS'
            ENDIF
*
*         COQUES MINCES
          ELSEIF (MFR.EQ.3) THEN
            JGOBL=6
            TABOBL(1)='YG1 '
            TABOBL(2)='YG2 '
            TABOBL(3)='NU12'
            TABOBL(4)='G12 '
            TABOBL(5)='V1X '
            TABOBL(6)='V1Y '
            IF(IFOUR.EQ.-2) THEN
              JGFAC=4
            ELSE
              JGFAC=3
            ENDIF
            TABFAC(1)='ALP1'
            TABFAC(2)='ALP2'
            TABFAC(3)='RHO '
            IF(IFOUR.EQ.-2) TABFAC(4)='DIM3'
 
*         COQUES AVEC CISAILLEMENT TRANSVERSE
          ELSEIF (MFR.EQ.9.OR.MFR.EQ.5) THEN
            JGOBL=8
            TABOBL(1)='YG1 '
            TABOBL(2)='YG2 '
            TABOBL(3)='NU12'
            TABOBL(4)='G12 '
            TABOBL(5)='G23 '
            TABOBL(6)='G13 '
            TABOBL(7)='V1X '
            TABOBL(8)='V1Y '
*
            JGFAC=3
            TABFAC(1)='RHO '
            TABFAC(2)='ALP1'
            TABFAC(3)='ALP2'
 
*         ELEMENTS MASSIFS
          ELSEIF (MFR.EQ.1.OR.MFR.EQ.31) THEN
            IF(IDIM.EQ.3)THEN
              JGOBL=15
              TABOBL(1)='YG1 '
              TABOBL(2)='YG2 '
              TABOBL(3)='YG3 '
              TABOBL(4)='NU12'
              TABOBL(5)='NU23'
              TABOBL(6)='NU13'
              TABOBL(7)='G12 '
              TABOBL(8)='G23 '
              TABOBL(9)='G13 '
              TABOBL(10)='V1X '
              TABOBL(11)='V1Y '
              TABOBL(12)='V1Z '
              TABOBL(13)='V2X '
              TABOBL(14)='V2Y '
              TABOBL(15)='V2Z '
*
              JGFAC=4
              TABFAC(1)='RHO '
              TABFAC(2)='ALP1'
              TABFAC(3)='ALP2'
              TABFAC(4)='ALP3'
            ELSEIF (IDIM.EQ.2) THEN
              IF(IFOUR.EQ.-2)THEN
                JGOBL=9
                TABOBL(1)='YG1 '
                TABOBL(2)='YG2 '
                TABOBL(3)='NU12'
                TABOBL(4)='G12 '
                TABOBL(5)='V1X '
                TABOBL(6)='V1Y '
                TABOBL(7)='YG3 '
                TABOBL(8)='NU23'
                TABOBL(9)='NU13'
*
                JGFAC=4
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='DIM3'
           ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
                JGOBL=9
                TABOBL(1)='YG1 '
                TABOBL(2)='YG2 '
                TABOBL(3)='YG3 '
                TABOBL(4)='NU12'
                TABOBL(5)='NU23'
                TABOBL(6)='NU13'
                TABOBL(7)='G12 '
                TABOBL(8)='V1X '
                TABOBL(9)='V1Y '
*
                JGFAC=4
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='ALP3'
              ELSEIF (IFOUR.EQ.1) THEN
                JGOBL=15
                TABOBL(1)='YG1 '
                TABOBL(2)='YG2 '
                TABOBL(3)='YG3 '
                TABOBL(4)='NU12'
                TABOBL(5)='NU23'
                TABOBL(6)='NU13'
                TABOBL(7)='G12 '
                TABOBL(8)='G23 '
                TABOBL(9)='G13 '
                TABOBL(10)='V1X '
                TABOBL(11)='V1Y '
                TABOBL(12)='V1Z '
                TABOBL(13)='V2X '
                TABOBL(14)='V2Y '
                TABOBL(15)='V2Z '
*
                JGFAC=4
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='ALP3'
              ENDIF
            ELSEIF (IDIM.EQ.1) THEN
C= Dans le cas UNID SPHErique, on doit avoir YG2=YG3 et NU12=NU13
C= et dans le cas thermomecanique ALP2=ALP3.
              JGOBL=6
              TABOBL(1)='YG1 '
              TABOBL(2)='YG2 '
              TABOBL(3)='YG3 '
              TABOBL(4)='NU12'
              TABOBL(5)='NU23'
              TABOBL(6)='NU13'
              IF (IFOUR.EQ.6) THEN
                JGFAC=2
              ELSEIF (IFOUR.EQ.5 .OR. IFOUR.EQ.10) THEN
                JGFAC=3
                TABFAC(3)='ALP3'
              ELSEIF (IFOUR.EQ.4.OR.IFOUR.EQ.8.OR.
     &                         IFOUR.EQ.13) THEN
                JGFAC=3
                TABFAC(3)='ALP2'
              ELSE
                JGFAC=4
                TABFAC(3)='ALP2'
                TABFAC(4)='ALP3'
              ENDIF
              TABFAC(1)='RHO '
              TABFAC(2)='ALP1'
            ENDIF
 
*         ELEMENTS JOINTS
          ELSEIF (MFR.EQ.35) THEN
            IF (IFOUR.EQ.2) THEN
              JGOBL=5
              TABOBL(1)='KS1 '
              TABOBL(2)='KS2 '
              TABOBL(3)='KN  '
              TABOBL(4)='V1X '
              TABOBL(5)='V1Y '
*
              JGFAC=2
              TABFAC(1)='RHO '
              TABFAC(2)='ALPN'
            ENDIF
          ENDIF
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE ANISOTROPE
C       =============================================================
        ELSEIF (IMATEE.EQ.3)THEN
C
C         JOINT UNIDIMESIONNEL JOI1
          IF(MFR.EQ.75)THEN
            IF (IFOUR.EQ.2) THEN
              JGOBL=27
              TABOBL(1)='V1X '
              TABOBL(2)='V1Y '
              TABOBL(3)='V1Z '
              TABOBL(4)='V2X '
              TABOBL(5)='V2Y '
              TABOBL(6)='V2Z '
              TABOBL(7)='D11 '
              TABOBL(8)='D22 '
              TABOBL(9)='D33 '
              TABOBL(10)='D44 '
              TABOBL(11)='D55 '
              TABOBL(12)='D66 '
              TABOBL(13)='D21 '
              TABOBL(14)='D31 '
              TABOBL(15)='D32 '
              TABOBL(16)='D41 '
              TABOBL(17)='D42 '
              TABOBL(18)='D43 '
              TABOBL(19)='D51 '
              TABOBL(20)='D52 '
              TABOBL(21)='D53 '
              TABOBL(22)='D54 '
              TABOBL(23)='D61 '
              TABOBL(24)='D62 '
              TABOBL(25)='D63 '
              TABOBL(26)='D64 '
              TABOBL(27)='D65 '
*
              JGFAC=10
* MASS: masse totale de l'élément joint
              TABFAC(1)='MASS'
              TABFAC(2)='JX  '
              TABFAC(3)='JY  '
              TABFAC(4)='JZ  '
              TABFAC(5)='ALP1'
              TABFAC(6)='ALP2'
              TABFAC(7)='ALP3'
              TABFAC(8)='ALQ1'
              TABFAC(9)='ALQ2'
              TABFAC(10)='ALQ3'
            ELSEIF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
              JGOBL=8
              TABOBL(1)='V1X '
              TABOBL(2)='V1Y '
              TABOBL(3)='D11 '
              TABOBL(4)='D22 '
              TABOBL(5)='D33 '
              TABOBL(6)='D21 '
              TABOBL(7)='D31 '
              TABOBL(8)='D32 '
*
              JGFAC=5
* MASS: masse totale de l'élément joint
              TABFAC(1)='MASS'
              TABFAC(2)='JZ  '
              TABFAC(3)='ALP1'
              TABFAC(4)='ALP2'
              TABFAC(5)='ALQ3'
            ENDIF
          ENDIF
*
*                   ELEMENTS MASSIFS
*
          IF(MFR.EQ.1.OR.MFR.EQ.31)THEN
            IF(IDIM.EQ.3)THEN
              JGOBL=27
              TABOBL(1)='D11 '
              TABOBL(2)='D21 '
              TABOBL(3)='D22 '
              TABOBL(4)='D31 '
              TABOBL(5)='D32 '
              TABOBL(6)='D33 '
              TABOBL(7)='D41 '
              TABOBL(8)='D42 '
              TABOBL(9)='D43 '
              TABOBL(10)='D44 '
              TABOBL(11)='D51 '
              TABOBL(12)='D52 '
              TABOBL(13)='D53 '
              TABOBL(14)='D54 '
              TABOBL(15)='D55 '
              TABOBL(16)='D61 '
              TABOBL(17)='D62 '
              TABOBL(18)='D63 '
              TABOBL(19)='D64 '
              TABOBL(20)='D65 '
              TABOBL(21)='D66 '
              TABOBL(22)='V1X '
              TABOBL(23)='V1Y '
              TABOBL(24)='V1Z '
              TABOBL(25)='V2X '
              TABOBL(26)='V2Y '
              TABOBL(27)='V2Z '
*
              JGFAC=7
              TABFAC(1)='RHO '
              TABFAC(2)='ALP1'
              TABFAC(3)='ALP2'
              TABFAC(4)='ALP3'
              TABFAC(5)='AL12'
              TABFAC(6)='AL13'
              TABFAC(7)='AL23'
            ELSEIF (IDIM.EQ.2) THEN
              IF (IFOUR.EQ.-2) THEN
                JGOBL=12
                TABOBL(1)='D11 '
                TABOBL(2)='D21 '
                TABOBL(3)='D22 '
                TABOBL(4)='D41 '
                TABOBL(5)='D42 '
                TABOBL(6)='D44 '
                TABOBL(7)='V1X '
                TABOBL(8)='V1Y '
                TABOBL(9)='D31 '
                TABOBL(10)='D32 '
                TABOBL(11)='D33 '
                TABOBL(12)='D43 '
*
                JGFAC=5
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='AL12'
                TABFAC(5)='DIM3'
              ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
                JGOBL=12
                TABOBL(1)='D11 '
                TABOBL(2)='D21 '
                TABOBL(3)='D22 '
                TABOBL(4)='D31 '
                TABOBL(5)='D32 '
                TABOBL(6)='D33 '
                TABOBL(7)='D41 '
                TABOBL(8)='D42 '
                TABOBL(9)='D43 '
                TABOBL(10)='D44 '
                TABOBL(11)='V1X '
                TABOBL(12)='V1Y '
*
                JGFAC=5
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='AL12'
                TABFAC(5)='ALP3'
              ELSEIF (IFOUR.EQ.1) THEN
                JGOBL=15
                TABOBL(1)='D11 '
                TABOBL(2)='D21 '
                TABOBL(3)='D22 '
                TABOBL(4)='D31 '
                TABOBL(5)='D32 '
                TABOBL(6)='D33 '
                TABOBL(7)='D41 '
                TABOBL(8)='D42 '
                TABOBL(9)='D43 '
                TABOBL(10)='D44 '
                TABOBL(11)='D55 '
                TABOBL(12)='D65 '
                TABOBL(13)='D66 '
                TABOBL(14)='V1X '
                TABOBL(15)='V1Y '
*
                JGFAC=5
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='AL12'
                TABFAC(5)='ALP3'
              ENDIF
            ENDIF
          ENDIF
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE POREUX ISOTROPE
C       =============================================================
        ELSEIF (IMATEE.EQ.7) THEN
          IF (MFR.EQ.33) THEN
            JGOBL=12
            TABOBL(1) ='YOUN'
            TABOBL(2) ='NU  '
            TABOBL(3) ='RHO '
            TABOBL(4) ='ALPH'
            TABOBL(5) ='COB '
            TABOBL(6) ='MOB '
            TABOBL(7) ='ALPM'
            TABOBL(8) ='PERM'
            TABOBL(9) ='VISC'
            TABOBL(10)='KF  '
            TABOBL(11)='RHOF'
            TABOBL(12)='ALPF'
            IF (IFOUR.EQ.-2) THEN
              JGFAC=1
              TABFAC(1)='DIM3'
            ENDIF
          ENDIF
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE HOMOGENEISE
C       =============================================================
        ELSEIF (IMATEE.EQ.5) THEN
          IF (MFR.EQ.37) THEN
            JGOBL=16
            TABOBL( 1)='B11 '
            TABOBL( 2)='B22 '
            TABOBL( 3)='B12 '
            TABOBL( 4)='ROF '
            TABOBL( 5)='ROS '
            TABOBL( 6)='YOUN'
            TABOBL( 7)='CSON'
            TABOBL( 8)='RORF'
            TABOBL( 9)='CREF'
            TABOBL(10)='LCAR'
            TABOBL(11)='E111'
            TABOBL(12)='E112'
            TABOBL(13)='E121'
            TABOBL(14)='E122'
            TABOBL(15)='E221'
            TABOBL(16)='E222'
          ENDIF
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE UNIDIRECTIONNEL
C       =============================================================
        ELSEIF (IMATEE.EQ.4) THEN
          IF ((MFR.EQ.1.OR.MFR.EQ.31).AND.IDIM.EQ.3) THEN
            JGOBL=7
            TABOBL(1)='YOUN'
            TABOBL(2)='V1X '
            TABOBL(3)='V1Y '
            TABOBL(4)='V1Z '
            TABOBL(5)='V2X '
            TABOBL(6)='V2Y '
            TABOBL(7)='V2Z '
          ELSE
            JGOBL=3
            TABOBL(1)='YOUN'
            TABOBL(2)='V1X '
            TABOBL(3)='V1Y '
          ENDIF
 
          IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND.
     &                    IFOUR.EQ.-2) THEN
            JGFAC=3
            TABFAC(3)='DIM3'
          ELSE
            JGFAC=2
          ENDIF
          TABFAC(1)='RHO '
          TABFAC(2)='ALPH'
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE SECTION DU MODELE A FIBRE
C       =============================================================
        ELSEIF (IMATEE.EQ.11) THEN
          JGOBL=2
          TABOBL(1)='MODS'
          TABOBL(2)='MATS'
 
          JGFAC=1
          TABFAC(1)='MANO'
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE ARMATURE POUR CABLE PRECONTRAINT
C       =============================================================
        ELSEIF (IMATEE.EQ.8) THEN
          JGOBL=1
          TABOBL(1)='YOUN'
 
          JGFAC=8
          TABFAC(1)='FF'
          TABFAC(2)='PHIF'
          TABFAC(3)='GANC'
          TABFAC(4)='RMU0'
          TABFAC(5)='FPRG'
          TABFAC(6)='RH10'
          TABFAC(7)='ALPH'
          TABFAC(8)='RHO '
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE MODAL
C       =============================================================
        ELSEIF (IMATEE.EQ.9) THEN
          JGOBL=3
          TABOBL(1)='FREQ'
          TABOBL(2)='MASS'
          TABOBL(3)='DEFO'
 
          JGFAC = 9
          TABFAC(1) = 'AMOR'
          TABFAC(2) = 'CGRA'
          TABFAC(3) = 'MADE'
          TABFAC(4) = 'RICR'
          TABFAC(5) = 'MAIB'
          TABFAC(6) = 'MACR'
          TABFAC(7) = 'AMCR'
          TABFAC(8) = 'ALP0'
          TABFAC(9) = 'ECRO'
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE STATIQUE
C       =============================================================
        ELSEIF (IMATEE.EQ.10) THEN
          JGOBL=3
          TABOBL(1)='DEFO'
          TABOBL(2)='RIDE'
          TABOBL(3)='MADE'
 
          JGFAC = 8
          TABFAC(1) = 'AMOR'
          TABFAC(2) = 'RICR'
          TABFAC(3) = 'MAIA'
          TABFAC(4) = 'MAIB'
          TABFAC(5) = 'MACR'
          TABFAC(6) = 'AMCR'
          TABFAC(7) = 'BET0'
          TABFAC(8) = 'ECRO'
C       =============================================================
C       MATERIAU MECANIQUE ELASTIQUE ZONE_COHESIVE
C       =============================================================
        ELSEIF (IMATEE.EQ.12) THEN
          JGOBL=2
          TABOBL(1)='KS'
          TABOBL(2)='KN'
        ENDIF
C
C       LECTURE DES AUTRES MOTS-CLES DE MATMOD (DE INMAT A NMAT)
        IF (NMAT.GE.INMAT) THEN
          CALL MODNLI(MOMODL,NMOD)
          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
          INMAT=INMAT+1
          IRET = 1
C         ===========================================================
C         MATERIAU MECANIQUE ELASTIQUE ... PLASTIQUE (ISOTROPE)
C         ===========================================================
          IF (IPLAC.EQ.1) THEN
            CALL MODPLA(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU MECANIQUE ELASTIQUE ... FLUAGE (NORTON)
C         ===========================================================
          ELSEIF (IPLAC.EQ.2) THEN
            CALL MODFLU(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU MECANIQUE ELASTIQUE ... VISCOPLASTIQUE (ONERA)
C         ===========================================================
          ELSEIF (IPLAC.EQ.3) THEN
            CALL MODVIS(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=2
C
C           Cas particulier des modeles GATT_MONERIE et UO2
C           'RHO ' et 'ALPH' sont obligatoires
            IF (IPLAC.EQ.18.OR.IPLAC.EQ.19) THEN
              JGOBL0 = JGOBL
              JGOBL=JGOBL0+2
              TABOBL(JGOBL0+1)='RHO '
              TABOBL(JGOBL0+2)='ALPH'
 
              JGOBL0=JGFAC
              JGFAC=JGOBL0-2
              IF (JGFAC.GT.0) THEN
                DO 200 I=1,JGFAC
                  TABFAC(I)=TABFAC(I+2)
 200            CONTINUE
              ENDIF
            ENDIF
            CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU MECANIQUE ELASTIQUE ... ENDOMMAGEMENT (MAZARS)
C         ===========================================================
          ELSEIF (IPLAC.EQ.4) THEN
            CALL MODEND(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU MECANIQUE ELASTIQUE ... PLASTIQUE_ENDOM (PSURY)
C         ===========================================================
          ELSEIF (IPLAC.EQ.5) THEN
            CALL MODPLE(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDPLEN(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU MECANIQUE ELASTIQUE ... NON LINEAIRE (EQUIPLAS)
C         ===========================================================
          ELSEIF (IPLAC.EQ.6) THEN
            CALL MODENL(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDELNL(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
          ENDIF
          IF (IRET.EQ.0) GOTO 9999
        ENDIF
C
C       Parametre CRIP
        CALL PLACE(MATMOD,NMAT,ICRIP,'CRIP')
        IF (ICRIP.NE.0) THEN
          JGOBL = JGOBL + 1
          TABOBL(JGOBL) = 'LIMP'
        ENDIF
C
C       Parametre FUSION
        CALL PLACE(MATMOD,NMAT,IFUSI,'FUSION')
        IF (IFUSI.NE.0) THEN
          JGOBL = JGOBL + 1
          TABOBL(JGOBL) = 'TFUS'
        ENDIF
C     =================================================================
C     FORMULATION POREUX
C     =================================================================
      ELSE IF (MFR2.EQ.33) THEN
C
        CALL MODLIN(MOMODL,NMOD)
        CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
        IF (IRET.EQ.0) GOTO 9999
C
        IF (NMAT.LT.2) GOTO 14
C
        CALL MODELA(MOMODL,NMOD)
        CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
C
        INMAT=3
        IF (IPLAC.EQ.0) INMAT=2
C
 14     CONTINUE
C
C       MATERIAU POREUX ELASTIQUE ISOTROPE (PAR DEFAUT)
C       ===============================================
        IF (IMATEE.EQ.1) THEN
          IF(MELE.GE.79.AND.MELE.LE.83)THEN
            JGOBL=4
            TABOBL(1) ='YOUN'
            TABOBL(2) ='NU  '
            TABOBL(3) ='COB '
            TABOBL(4) ='MOB '
*
            JGFAC=8
            IF(IFOUR.EQ.-2) THEN
              JGFAC=9
              TABFAC(9)='DIM3'
            ENDIF
            TABFAC(1)='RHOF'
            TABFAC(2)='ALPF'
            TABFAC(3)='ALPM'
            TABFAC(4)='PERM'
            TABFAC(5)='VISC'
            TABFAC(6)='KF  '
            TABFAC(7)='RHO '
            TABFAC(8)='ALPH'
          ELSEIF(MELE.GE.108.AND.MELE.LE.110)THEN
            JGOBL=4
            TABOBL(1)='KS  '
            TABOBL(2)='KN  '
            TABOBL(3)='COB '
            TABOBL(4)='MOB '
*
            JGFAC=4
            TABFAC(1)='PERT'
            TABFAC(2)='PERH'
            TABFAC(3)='PERB'
            TABFAC(4)='VISC'
          ELSEIF(MELE.GE.173.AND.MELE.LE.177)THEN
            JGOBL=10
            TABOBL(1) ='YOUN'
            TABOBL(2) ='NU  '
            TABOBL(3) ='COP1'
            TABOBL(4) ='COP2'
            TABOBL(5) ='CPP1'
            TABOBL(6) ='CPP2'
            TABOBL(7) ='KK11'
            TABOBL(8) ='KK12'
            TABOBL(9) ='KK21'
            TABOBL(10)='KK22'
*
            JGFAC=8
            IF(IFOUR.EQ.-2) THEN
              JGFAC=9
              TABFAC(9)='DIM3'
            ENDIF
            TABFAC(1)='RHOF'
            TABFAC(2)='ALPF'
            TABFAC(3)='ALPM'
            TABFAC(4)='PK11'
            TABFAC(5)='PK12'
            TABFAC(6)='PK21'
            TABFAC(7)='PK22'
            TABFAC(8)='ALPH'
          ELSEIF(MELE.GE.178.AND.MELE.LE.182)THEN
            JGOBL=17
            TABOBL(1)='YOUN'
            TABOBL(2)='NU  '
            TABOBL(3)='COP1'
            TABOBL(4)='COP2'
            TABOBL(5)='COP3'
            TABOBL(6)='CPP1'
            TABOBL(7)='CPP2'
            TABOBL(8)='CPP3'
            TABOBL(9)='KK11'
            TABOBL(10)='KK12'
            TABOBL(11)='KK13'
            TABOBL(12)='KK21'
            TABOBL(13)='KK22'
            TABOBL(14)='KK23'
            TABOBL(15)='KK31'
            TABOBL(16)='KK32'
            TABOBL(17)='KK33'
*
            JGFAC=14
            IF(IFOUR.EQ.-2) THEN
              JGFAC=15
              TABFAC(15)='DIM3'
            ENDIF
            TABFAC(1)='RHOF'
            TABFAC(2)='ALPF'
            TABFAC(3)='ALPM'
            TABFAC(4)='PK11'
            TABFAC(5)='PK12'
            TABFAC(6)='PK13'
            TABFAC(7)='PK21'
            TABFAC(8)='PK22'
            TABFAC(9)='PK23'
            TABFAC(10)='PK31'
            TABFAC(11)='PK32'
            TABFAC(12)='PK33'
            TABFAC(13)='RHO '
            TABFAC(14)='ALPH'
          ELSEIF(MELE.GE.185.AND.MELE.LE.187)THEN
            JGOBL=10
            TABOBL(1)='KS  '
            TABOBL(2)='KN  '
            TABOBL(3)='COP1'
            TABOBL(4)='COP2'
            TABOBL(5)='CPP1'
            TABOBL(6)='CPP2'
            TABOBL(7)='KK11'
            TABOBL(8)='KK12'
            TABOBL(9)='KK21'
            TABOBL(10)='KK22'
*
            JGFAC=17
            IF(IFOUR.EQ.-2) THEN
              JGFAC=18
              TABFAC(18)='DIM3'
            ENDIF
            TABFAC(1)='RHOF'
            TABFAC(2)='ALPF'
            TABFAC(3)='ALPM'
            TABFAC(4)='PT11'
            TABFAC(5)='PH11'
            TABFAC(6)='PB11'
            TABFAC(7)='PT12'
            TABFAC(8)='PH12'
            TABFAC(9)='PB12'
            TABFAC(10)='PT21'
            TABFAC(11)='PH21'
            TABFAC(12)='PB21'
            TABFAC(13)='PT22'
            TABFAC(14)='PH22'
            TABFAC(15)='PB22'
            TABFAC(16)='RHO '
            TABFAC(17)='ALPH'
          ELSEIF(MELE.GE.188.AND.MELE.LE.190)THEN
            JGOBL=17
            TABOBL(1)='KS  '
            TABOBL(2)='KN  '
            TABOBL(3)='COP1'
            TABOBL(4)='COP2'
            TABOBL(5)='COP3'
            TABOBL(6)='CPP1'
            TABOBL(7)='CPP2'
            TABOBL(8)='CPP3'
            TABOBL(9)='KK11'
            TABOBL(10)='KK12'
            TABOBL(11)='KK13'
            TABOBL(12)='KK21'
            TABOBL(13)='KK22'
            TABOBL(14)='KK23'
            TABOBL(15)='KK31'
            TABOBL(16)='KK32'
            TABOBL(17)='KK33'
*
            JGFAC=32
            IF(IFOUR.EQ.-2) THEN
              JGFAC=33
              TABFAC(33)='DIM3'
            ENDIF
            TABFAC(1)='RHOF'
            TABFAC(2)='ALPF'
            TABFAC(3)='ALPM'
            TABFAC(4)='PT11'
            TABFAC(5)='PH11'
            TABFAC(6)='PB11'
            TABFAC(7)='PT12'
            TABFAC(8)='PH12'
            TABFAC(9)='PB12'
            TABFAC(10)='PT13'
            TABFAC(11)='PH13'
            TABFAC(12)='PB13'
            TABFAC(13)='PT21'
            TABFAC(14)='PH21'
            TABFAC(15)='PB21'
            TABFAC(16)='PT22'
            TABFAC(17)='PH22'
            TABFAC(18)='PB22'
            TABFAC(19)='PT23'
            TABFAC(20)='PH23'
            TABFAC(21)='PB23'
            TABFAC(22)='PT31'
            TABFAC(23)='PH31'
            TABFAC(24)='PB31'
            TABFAC(25)='PT32'
            TABFAC(26)='PH32'
            TABFAC(27)='PB32'
            TABFAC(28)='PT33'
            TABFAC(29)='PH33'
            TABFAC(30)='PB33'
            TABFAC(31)='RHO '
            TABFAC(32)='ALPH'
          ENDIF
C
C       MATERIAU POREUX ELASTIQUE ORTHOTROPE
C       ====================================
        ELSEIF (IMATEE.EQ.2) THEN
          IF (MFR.EQ.33) THEN
            IF(IDIM.EQ.3)THEN
*                      ELEMENTS 3D
              JGOBL=19
              TABOBL(1)='YG1 '
              TABOBL(2)='YG2 '
              TABOBL(3)='YG3 '
              TABOBL(4)='NU12'
              TABOBL(5)='NU23'
              TABOBL(6)='NU13'
              TABOBL(7)='G12 '
              TABOBL(8)='G23 '
              TABOBL(9)='G13 '
              TABOBL(10)='V1X '
              TABOBL(11)='V1Y '
              TABOBL(12)='V1Z '
              TABOBL(13)='V2X '
              TABOBL(14)='V2Y '
              TABOBL(15)='V2Z '
              TABOBL(16)='COB1'
              TABOBL(17)='COB2'
              TABOBL(18)='COB3'
              TABOBL(19)='MOB '
*
              JGFAC=12
              TABFAC(1)='RHO '
              TABFAC(2)='ALP1'
              TABFAC(3)='ALP2'
              TABFAC(4)='ALP3'
              TABFAC(5)='ALPM'
              TABFAC(6)='PER1'
              TABFAC(7)='PER2'
              TABFAC(8)='PER3'
              TABFAC(9)='VISC'
              TABFAC(10)='KF  '
              TABFAC(11)='RHOF'
              TABFAC(12)='ALPF'
            ELSEIF (IDIM.EQ.2) THEN
              IF(IFOUR.EQ.-2)THEN
                JGOBL=12
                TABOBL(1)='YG1 '
                TABOBL(2)='YG2 '
                TABOBL(3)='NU12'
                TABOBL(4)='G12 '
                TABOBL(5)='V1X '
                TABOBL(6)='V1Y '
                TABOBL(7)='COB1'
                TABOBL(8)='COB2'
                TABOBL(9)='MOB '
                TABOBL(10)='YG3 '
                TABOBL(11)='NU23'
                TABOBL(12)='NU13'
*
                JGFAC=11
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='ALPM'
                TABFAC(5)='PER1'
                TABFAC(6)='PER2'
                TABFAC(7)='VISC'
                TABFAC(8)='KF  '
                TABFAC(9)='RHOF'
                TABFAC(10)='ALPF'
                TABFAC(11)='DIM3'
           ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
                JGOBL=13
                TABOBL(1)='YG1 '
                TABOBL(2)='YG2 '
                TABOBL(3)='YG3 '
                TABOBL(4)='NU12'
                TABOBL(5)='NU23'
                TABOBL(6)='NU13'
                TABOBL(7)='G12 '
                TABOBL(8)='V1X '
                TABOBL(9)='V1Y '
                TABOBL(10)='COB1'
                TABOBL(11)='COB2'
                TABOBL(12)='COB3'
                TABOBL(13)='MOB '
*
                JGFAC=11
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='ALP3'
                TABFAC(5)='ALPM'
                TABFAC(6)='PER1'
                TABFAC(7)='PER2'
                TABFAC(8)='VISC'
                TABFAC(9)='KF  '
                TABFAC(10)='RHOF'
                TABFAC(11)='ALPF'
              ELSEIF (IFOUR.EQ.1) THEN
                JGOBL=15
                TABOBL(1)='YG1 '
                TABOBL(2)='YG2 '
                TABOBL(3)='YG3 '
                TABOBL(4)='NU12'
                TABOBL(5)='NU23'
                TABOBL(6)='NU13'
                TABOBL(7)='G12 '
                TABOBL(8)='G23 '
                TABOBL(9)='G13 '
                TABOBL(10)='V1X '
                TABOBL(11)='V1Y '
                TABOBL(12)='COB1'
                TABOBL(13)='COB2'
                TABOBL(14)='COB3'
                TABOBL(15)='MOB '
*
                JGFAC=12
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='ALP3'
                TABFAC(5)='ALPM'
                TABFAC(6)='PER1'
                TABFAC(7)='PER2'
                TABFAC(8)='PER3'
                TABFAC(9)='VISC'
                TABFAC(10)='KF  '
                TABFAC(11)='RHOF'
                TABFAC(12)='ALPF'
              ENDIF
            ENDIF
          ELSE
            IRET = 0
            GOTO 9999
          ENDIF
C
C       MATERIAU POREUX ELASTIQUE ANISOTROPE
C       ====================================
        ELSEIF (IMATEE.EQ.3)THEN
          IF(MFR.EQ.33)THEN
            IF(IDIM.EQ.3)THEN
              JGOBL=34
              TABOBL(1)='D11 '
              TABOBL(2)='D21 '
              TABOBL(3)='D22 '
              TABOBL(4)='D31 '
              TABOBL(5)='D32 '
              TABOBL(6)='D33 '
              TABOBL(7)='D41 '
              TABOBL(8)='D42 '
              TABOBL(9)='D43 '
              TABOBL(10)='D44 '
              TABOBL(11)='D51 '
              TABOBL(12)='D52 '
              TABOBL(13)='D53 '
              TABOBL(14)='D54 '
              TABOBL(15)='D55 '
              TABOBL(16)='D61 '
              TABOBL(17)='D62 '
              TABOBL(18)='D63 '
              TABOBL(19)='D64 '
              TABOBL(20)='D65 '
              TABOBL(21)='D66 '
              TABOBL(22)='V1X '
              TABOBL(23)='V1Y '
              TABOBL(24)='V1Z '
              TABOBL(25)='V2X '
              TABOBL(26)='V2Y '
              TABOBL(27)='V2Z '
              TABOBL(28)='COB1'
              TABOBL(29)='COB2'
              TABOBL(30)='COB3'
              TABOBL(31)='CO12'
              TABOBL(32)='CO13'
              TABOBL(33)='CO23'
              TABOBL(34)='MOB '
*
              JGFAC=18
              TABFAC(1)='RHO '
              TABFAC(2)='ALP1'
              TABFAC(3)='ALP2'
              TABFAC(4)='ALP3'
              TABFAC(5)='AL12'
              TABFAC(6)='AL13'
              TABFAC(7)='AL23'
              TABFAC(8)='ALPM'
              TABFAC(9 )='PER1'
              TABFAC(10)='PER2'
              TABFAC(11)='PER3'
              TABFAC(12)='PE12'
              TABFAC(13)='PE13'
              TABFAC(14)='PE23'
              TABFAC(15)='VISC'
              TABFAC(16)='KF  '
              TABFAC(17)='RHOF'
              TABFAC(18)='ALPF'
            ELSEIF (IDIM.EQ.2) THEN
              IF (IFOUR.EQ.-2) THEN
                JGOBL=16
                TABOBL(1)='D11 '
                TABOBL(2)='D21 '
                TABOBL(3)='D22 '
                TABOBL(4)='D41 '
                TABOBL(5)='D42 '
                TABOBL(6)='D44 '
                TABOBL(7)='V1X '
                TABOBL(8)='V1Y '
                TABOBL(9 )='COB1'
                TABOBL(10)='COB2'
                TABOBL(11)='CO12'
                TABOBL(12)='MOB '
                TABOBL(13)='D31 '
                TABOBL(14)='D32 '
                TABOBL(15)='D33 '
                TABOBL(16)='D43 '
*
                JGFAC=13
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='AL12'
                TABFAC(5)='ALPM'
                TABFAC(6)='PER1'
                TABFAC(7)='PER2'
                TABFAC(8)='PE12'
                TABFAC(9)='VISC'
                TABFAC(10)='KF  '
                TABFAC(11)='RHOF'
                TABFAC(12)='ALPF'
                TABFAC(13)='DIM3'
*
              ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
                JGOBL=17
                TABOBL(1)='D11 '
                TABOBL(2)='D21 '
                TABOBL(3)='D22 '
                TABOBL(4)='D31 '
                TABOBL(5)='D32 '
                TABOBL(6)='D33 '
                TABOBL(7)='D41 '
                TABOBL(8)='D42 '
                TABOBL(9)='D43 '
                TABOBL(10)='D44 '
                TABOBL(11)='V1X '
                TABOBL(12)='V1Y '
                TABOBL(13)='COB1'
                TABOBL(14)='COB2'
                TABOBL(15)='CO12'
                TABOBL(16)='COB3'
                TABOBL(17)='MOB '
*
                JGFAC=13
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='AL12'
                TABFAC(5)='ALP3'
                TABFAC(6)='ALPM'
                TABFAC(7)='PER1'
                TABFAC(8)='PER2'
                TABFAC(9)='PE12'
                TABFAC(10)='VISC'
                TABFAC(11)='KF  '
                TABFAC(12)='RHOF'
                TABFAC(13)='ALPF'
              ELSEIF (IFOUR.EQ.1) THEN
                JGOBL=20
                TABOBL(1)='D11 '
                TABOBL(2)='D21 '
                TABOBL(3)='D22 '
                TABOBL(4)='D31 '
                TABOBL(5)='D32 '
                TABOBL(6)='D33 '
                TABOBL(7)='D41 '
                TABOBL(8)='D42 '
                TABOBL(9)='D43 '
                TABOBL(10)='D44 '
                TABOBL(11)='D55 '
                TABOBL(12)='D65 '
                TABOBL(13)='D66 '
                TABOBL(14)='V1X '
                TABOBL(15)='V1Y '
                TABOBL(16)='COB1'
                TABOBL(17)='COB2'
                TABOBL(18)='CO12'
                TABOBL(19)='COB3'
                TABOBL(20)='MOB '
*
                JGFAC=14
                TABFAC(1)='RHO '
                TABFAC(2)='ALP1'
                TABFAC(3)='ALP2'
                TABFAC(4)='AL12'
                TABFAC(5)='ALP3'
                TABFAC(6)='ALPM'
                TABFAC(7)='PER1'
                TABFAC(8)='PER2'
                TABFAC(9)='PE12'
                TABFAC(10)='PER3'
                TABFAC(11)='VISC'
                TABFAC(12)='KF  '
                TABFAC(13)='RHOF'
                TABFAC(14)='ALPF'
              ENDIF
            ENDIF
          ELSE
            IRET = 0
            GOTO 9999
          ENDIF
C
C       MATERIAU POREUX ELASTIQUE UNIDIRECTIONNEL
C       =========================================
        ELSEIF (IMATEE.EQ.4) THEN
          IF(MFR.EQ.33)THEN
            IF (IDIM.EQ.3) THEN
              JGOBL=9
              TABOBL(1)='YOUN'
              TABOBL(2)='V1X '
              TABOBL(3)='V1Y '
              TABOBL(4)='V1Z '
              TABOBL(5)='V2X '
              TABOBL(6)='V2Y '
              TABOBL(7)='V2Z '
              TABOBL(8)='COB '
              TABOBL(9)='MOB '
            ELSE
              JGOBL=5
              TABOBL(1)='YOUN'
              TABOBL(2)='V1X '
              TABOBL(3)='V1Y '
              TABOBL(4)='COB '
              TABOBL(5)='MOB '
            ENDIF
*
            JGFAC=8
            IF(IFOUR.EQ.-2) THEN
              JGFAC=9
              TABFAC(9)='DIM3'
            ENDIF
            TABFAC(1)='RHO '
            TABFAC(2)='ALPH'
            TABFAC(3)='ALPM'
            TABFAC(4)='PERM'
            TABFAC(5)='VISC'
            TABFAC(6)='KF  '
            TABFAC(7)='RHOF'
            TABFAC(8)='ALPF'
          ELSE
            IRET = 0
            GOTO 9999
          ENDIF
C       =============================================================
        ENDIF
C
        IF (NMAT.LT.2) GOTO 9999
C
C       LECTURE DES AUTRES MOTS-CLES DE MATMOD (DE INMAT A NMAT)
        IF (NMAT.GE.INMAT) THEN
          CALL MODNLI(MOMODL,NMOD)
          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
          INMAT=INMAT+1
          IRET = 1
C         ===========================================================
C         MATERIAU POREUX ELASTIQUE ... PLASTIQUE (ISOTROPE)
C         ===========================================================
          IF (IPLAC.EQ.1) THEN
            CALL MODPLA(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU POREUX ELASTIQUE ... FLUAGE (NORTON)
C         ===========================================================
          ELSEIF (IPLAC.EQ.2) THEN
            CALL MODFLU(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU POREUX ELASTIQUE ... VISCOPLASTIQUE (ONERA)
C         ===========================================================
          ELSEIF (IPLAC.EQ.3) THEN
            CALL MODVIS(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=2
            CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
C         MATERIAU POREUX ELASTIQUE ... ENDOMMAGEMENT (MAZARS)
C         ===========================================================
          ELSEIF (IPLAC.EQ.4) THEN
            CALL MODEND(MOMODL,NMOD)
            CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
            IF (IPLAC.EQ.0) IPLAC=1
            CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
C         ===========================================================
          ENDIF
          IF (IRET.EQ.0) GOTO 9999
C
C         En cas de creation de materiaux combinant plusieurs materiaux
C         deja existant Ex ELASTIQUE ISOTROPE PLASTIQUE PARFAIT FLUAGE N
C         GOTO 40
C
        ENDIF
C     =================================================================
C     FORMULATION CONTACT
C     =================================================================
      ELSE IF (MFR2.EQ.4) THEN
C
        IF (INATUU.EQ.2) THEN
          JGOBL=2
          TABOBL(1)='FF  '
          TABOBL(2)='PHIF'
 
          JGFAC=0
        ELSEIF (INATUU.EQ.1) THEN
          JGOBL=1
          TABOBL(1)='MU  '
 
          JGFAC=3
          TABFAC(1)='COHE'
          TABFAC(2)='ADHE'
          TABFAC(3)='JEU'
        ELSE
          JGOBL=0
 
          JGFAC=2
          TABFAC(1)='JEU'
          TABFAC(2)='ADHE'
        ENDIF
C     =================================================================
C     FORMULATION CONTRAINTE
C     =================================================================
      ELSE IF (MFR2.EQ.6) THEN
C
        IF (CMATEE.EQ.'ROTATION') THEN
          JGOBL=1
          TABOBL(1)='ANGL'
        ELSEIF (CMATEE.EQ.'DEPLACEMENT') THEN
          JGOBL=1
          TABOBL(1)='AMPL'
        ENDIF
C     =================================================================
C     FORMULATION MAGNETODYNAMIQUE
C     =================================================================
      ELSE IF (MFR2.EQ.69) THEN
C
        IF(MFR.EQ.3) THEN
          IF(IMATEE.EQ.1) THEN
            JGOBL=3
            TABOBL(1)='ETA'
            TABOBL(2)='PERM'
            TABOBL(3)='EPAI'
          ELSEIF (IMATEE.EQ.3) THEN
            JGOBL=4
            TABOBL(1)='ETA1'
            TABOBL(2)='ETA2'
            TABOBL(3)='PERM'
            TABOBL(4)='EPAI'
          ENDIF
        ENDIF
C     =================================================================
C     FORMULATION FISSURE
C     =================================================================
      ELSE IF (MFR2.EQ.40) THEN
C
C       si POISEU_BLASIUS ou POISEU_COLEBROOK ou par defaut
        JGOBL=1
        TABOBL(1)='RUGO'
        CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT1')
        CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT2')
        IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
          JGOBL=7
          TABOBL(2)='REC'
          TABOBL(3)='FK'
          TABOBL(4)='FA'
          TABOBL(5)='FB'
          TABOBL(6)='FC'
          TABOBL(7)='FD'
        ELSE
          CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT3')
          CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT4')
          IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
            JGOBL=2
            TABOBL(2)='FK'
          ENDIF
        ENDIF
C     =================================================================
C     FORMULATION NAVIER_STOKES
C     =================================================================
      ELSE IF (MFR2.EQ.52) THEN
C
        IF (IMATEE.EQ.4) THEN
*NLIN          
          JGOBL = 2
          TABOBL(1) ='RHO '
          TABOBL(2) ='MU  '
 
          JGFAC = 1
          TABFAC(1) ='FREQ'
        ENDIF      
C     =================================================================
C     FORMULATION MELANGE
C     =================================================================
      ELSE IF (MFR2.EQ.38) THEN
C
        IF (IMATEE.NE.3.AND.IMATEE.NE.4) THEN
          CALL IDMETA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IMATEE,IRET)
          IF (IRET.EQ.0) GOTO 9999
        ENDIF
C     =================================================================
C     FORMULATION LIAISON
C     =================================================================
      ELSE IF (MFR2.EQ.24) THEN
C
        IF (IMATEE.EQ.0) THEN
          IRET = 0
          GOTO 9999
*
        ELSEIF (IMATEE.EQ.1) THEN
*PO_PL_FL
         JGOBL = 7
         TABOBL(1) ='NORM'
         TABOBL(2) ='INER'
         TABOBL(3) ='CONV'
         TABOBL(4) ='VISC'
         TABOBL(5) ='PELO'
         TABOBL(6) ='PRAP'
         TABOBL(7) ='JFLU'
 
        ELSEIF (IMATEE.EQ.2) THEN
*PO_PL_FR
         JGOBL = 7
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='JEU'
         TABOBL(4) ='GLIS'
         TABOBL(5) ='ADHE'
         TABOBL(6) ='RTAN'
         TABOBL(7) ='ATAN'
 
         JGFAC = 2
         TABFAC(1) = 'AMOR'
         TABFAC(2) = 'LOIC'
 
        ELSEIF (IMATEE.EQ.3) THEN
*PO_PL
         JGOBL= 3
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='JEU'
 
         JGFAC = 4
         TABFAC(1) = 'LOIC'
         TABFAC(2) = 'PERM'
         TABFAC(3) ='SPLA'
         TABFAC(4) ='AMOR'
 
       ELSEIF (IMATEE.EQ.4) THEN
*PO_PO_FR
         JGOBL = 8
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='JEU'
         TABOBL(4) ='POIB'
         TABOBL(5) ='ADHE'
         TABOBL(6) ='RTAN'
         TABOBL(7) ='ATAN'
         TABOBL(8) ='GLIS'
 
         JGFAC =  3
         TABFAC(1) = 'AMOR'
         TABFAC(2) = 'LOIC'
         TABFAC(3) = 'MODE'
 
        ELSEIF (IMATEE.EQ.5) THEN
*PO_PO_DP
         JGOBL = 6
         TABOBL(1) ='NORM'
         TABOBL(2) ='ECRO'
         TABOBL(3) ='JEU'
         TABOBL(4) ='POIB'
         TABOBL(5) ='PERM'
         TABOBL(6) ='LOIC'
 
         JGFAC =  1
         TABFAC(1) = 'AMOR'
 
        ELSEIF (IMATEE.EQ.6) THEN
*PO_PO_RP
         JGOBL = 6
         TABOBL(1) ='AXRO'
         TABOBL(2) ='ECRO'
         TABOBL(3) ='JEU'
         TABOBL(4) ='POIB'
         TABOBL(5) ='PERM'
         TABOBL(6) ='LOIC'
 
         JGFAC =  2
         TABFAC(1) = 'AMOR'
         TABFAC(2) = 'ELAS'
 
        ELSEIF (IMATEE.EQ.7) THEN
*PO_PO
         JGOBL = 5
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='JEU'
         TABOBL(4) ='POIB'
         TABOBL(5) ='PERM'
 
         JGFAC =  2
         TABFAC(1) = 'AMOR'
         TABFAC(2) = 'LOIC'
 
        ELSEIF (IMATEE.EQ.8) THEN
*PO_CE_MO
         JGOBL = 8
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='PCER'
         TABOBL(4) ='RAYO'
         TABOBL(5) ='GLIS'
         TABOBL(6) ='ADHE'
         TABOBL(7) ='RTAN'
         TABOBL(8) ='ATAN'
 
         JGFAC=  2
         TABFAC(1) = 'CINT'
         TABFAC(2) = 'AMOR'
 
        ELSEIF (IMATEE.EQ.9) THEN
*PO_CE_FR
         JGOBL = 8
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='EXCE'
         TABOBL(4) ='RAYO'
         TABOBL(5) ='GLIS'
         TABOBL(6) ='ADHE'
         TABOBL(7) ='RTAN'
         TABOBL(8) ='ATAN'
 
         JGFAC =  2
         TABFAC(1) = 'CINT'
         TABFAC(2) = 'AMOR'
 
        ELSEIF (IMATEE.EQ.10) THEN
*PO_CE
         JGOBL = 4
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='EXCE'
         TABOBL(4) ='RAYO'
 
         JGFAC =  1
         TABFAC(1) = 'AMOR'
 
        ELSEIF (IMATEE.EQ.11) THEN
*CE_PL_FR
         JGOBL = 8
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='JEU'
         TABOBL(4) ='RAYS'
         TABOBL(5) ='GLIS'
         TABOBL(6) ='ADHE'
         TABOBL(7) ='RTAN'
         TABOBL(8) ='ATAN'
 
         JGFAC = 1
         TABFAC(1) = 'AMOR'
 
        ELSEIF (IMATEE.EQ.12) THEN
*CE_CE_FR
         JGOBL = 9
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='EXCE'
         TABOBL(4) ='RAYS'
         TABOBL(5) ='GLIS'
         TABOBL(6) ='ADHE'
         TABOBL(7) ='RTAN'
         TABOBL(8) ='ATAN'
         TABOBL(9) ='RAYB'
 
         JGFAC =  2
         TABFAC(1) = 'AMOR'
         TABFAC(2) = 'CINT'
 
        ELSEIF (IMATEE.EQ.13.OR.IMATEE.EQ.14) THEN
*PR_PR_IN ou PR_PR_EX
         JGOBL = 5
         TABOBL(1) ='NORM'
         TABOBL(2) ='RAID'
         TABOBL(3) ='PFIX'
         TABOBL(4) ='PMOB'
         TABOBL(5) ='ERAI'
 
        ELSEIF (IMATEE.EQ.15) THEN
*LI_LI_FR
         JGOBL = 8
         TABOBL(1) ='NORM'
         TABOBL(2) ='LIMA'
         TABOBL(3) ='LIES'
         TABOBL(4) ='RAID'
         TABOBL(5) ='GLIS'
         TABOBL(6) ='ADHE'
         TABOBL(7) ='RTAN'
         TABOBL(8) ='ATAN'
 
         JGFAC =  4
         TABFAC(1) = 'JEU'
         TABFAC(2) = 'AMOR'
         TABFAC(3) = 'RECH'
         TABFAC(4) = 'SYME'
 
        ELSEIF (IMATEE.EQ.16) THEN
*LI_CE_FR
         JGOBL = 8
         TABOBL(1) ='NORM'
         TABOBL(2) ='LIMA'
         TABOBL(3) ='LIES'
         TABOBL(4) ='RAID'
         TABOBL(5) ='GLIS'
         TABOBL(6) ='ADHE'
         TABOBL(7) ='RTAN'
         TABOBL(8) ='ATAN'
 
         JGFAC =  5
         TABFAC(1) = 'AMOR'
         TABFAC(2) = 'RECH'
         TABFAC(3) = 'RAYO'
         TABFAC(4) = 'ACTN'
         TABFAC(5) = 'INVE'
 
        ELSEIF (IMATEE.EQ.17) THEN
*PA_FL_RO
         JGOBL =10
         TABOBL(1) ='LONG'
         TABOBL(2) ='RAYO'
         TABOBL(3) ='VISC'
         TABOBL(4) ='RHOF'
         TABOBL(5) ='PADM'
         TABOBL(6) ='VROT'
         TABOBL(7) ='EPSI'
         TABOBL(8) ='PHII'
         TABOBL(9) ='AFFI'
         TABOBL(10)='TLOB'
 
         JGFAC =  1
         TABFAC(1) = 'AMOR'
 
        ELSEIF (IMATEE.EQ.23) THEN
*NEWMARK MODAL
         JGOBL =3
         TABOBL(1) ='JEU'
         TABOBL(2) ='MASS'
         TABOBL(3) ='FREQ'
 
         JGFAC =  3
         TABFAC(1) = 'EXCE'
         TABFAC(2) = 'FROT'
         TABFAC(3) = 'MOFR'
        ENDIF
 
C       'SORT' facultatif dans tous les cas
        JGOBL0 = JGFAC
        JGFAC= JGOBL0+1
        TABFAC(JGOBL0+1) = 'SORT'
C     =================================================================
C     FORMULATION ELECTROSTATIQUE
C     =================================================================
      ELSE IF (MFR2.EQ.71) THEN
C
C       Permittivite isotrope
        IF (IMATEE.EQ.1) THEN
          JGOBL = 1
          TABOBL(1)='PEL     '
C
C       Permittivite orthotrope
        ELSEIF (IMATEE.EQ.2) THEN
          IF (IDIM.EQ.2) THEN
            IF (IFOMOD.NE.1) THEN
              JGOBL = 4
              TABOBL(1) = 'PE1     '
              TABOBL(2) = 'PE2     '
              TABOBL(3) = 'V1X     '
              TABOBL(4) = 'V1Y     '
            ELSE
              JGOBL = 5
              TABOBL(1) = 'PE1     '
              TABOBL(2) = 'PE2     '
              TABOBL(3) = 'PE3     '
              TABOBL(4) = 'V1X     '
              TABOBL(5) = 'V1Y     '
            ENDIF
          ELSEIF (IDIM.EQ.3) THEN
            JGOBL = 9
            TABOBL(1) = 'PE1     '
            TABOBL(2) = 'PE2     '
            TABOBL(3) = 'PE3     '
            TABOBL(4) = 'V1X     '
            TABOBL(5) = 'V1Y     '
            TABOBL(6) = 'V1Z     '
            TABOBL(7) = 'V2X     '
            TABOBL(8) = 'V2Y     '
            TABOBL(9) = 'V2Z     '
          ENDIF
C
C       Permittivite anisotrope
        ELSEIF (IMATEE.EQ.3) THEN
          IF (IDIM.EQ.2) THEN
            IF (IFOMOD.NE.1) THEN
              JGOBL = 5
              TABOBL(1) = 'PE11    '
              TABOBL(2) = 'PE22    '
              TABOBL(3) = 'PE21    '
              TABOBL(4) = 'V1X     '
              TABOBL(5) = 'V1Y     '
            ELSE
              JGOBL = 6
              TABOBL(1) = 'PE11    '
              TABOBL(2) = 'PE22    '
              TABOBL(3) = 'PE21    '
              TABOBL(4) = 'PE33    '
              TABOBL(5) = 'V1X     '
              TABOBL(6) = 'V1Y     '
            ENDIF
          ELSEIF (IDIM.EQ.3) THEN
            JGOBL = 12
            TABOBL( 1) = 'PE11    '
            TABOBL( 2) = 'PE22    '
            TABOBL( 3) = 'PE33    '
            TABOBL( 4) = 'PE21    '
            TABOBL( 5) = 'PE31    '
            TABOBL( 6) = 'PE32    '
            TABOBL( 7) = 'V1X     '
            TABOBL( 8) = 'V1Y     '
            TABOBL( 9) = 'V1Z     '
            TABOBL(10) = 'V2X     '
            TABOBL(11) = 'V2Y     '
            TABOBL(12) = 'V2Z     '
          ENDIF
        ENDIF
C     =================================================================
C     FORMULATION DIFFUSION
C     =================================================================
      ELSE IF (MFR2.EQ.73) THEN
C
        JGFAC    = 1
        TABFAC(1)='CDIF'
C
C       Diffusion isotrope
        IF (IMATEE.EQ.1) THEN
          JGOBL     = 1
C*8       TABOBL(1) = 'KDIFF   '
          TABOBL(1) = 'KD      '
C
C       Diffusion orthotrope
        ELSEIF (IMATEE.EQ.2) THEN
          IF (IDIM.EQ.2) THEN
            IF (IFOMOD.NE.1) THEN
              JGOBL     = 4
              TABOBL(1) = 'KD1     '
              TABOBL(2) = 'KD2     '
              TABOBL(3) = 'V1X     '
              TABOBL(4) = 'V1Y     '
            ELSE
              JGOBL     = 5
              TABOBL(1) = 'KD1     '
              TABOBL(2) = 'KD2     '
              TABOBL(3) = 'KD3     '
              TABOBL(4) = 'V1X     '
              TABOBL(5) = 'V1Y     '
            ENDIF
          ELSEIF (IDIM.EQ.3) THEN
            JGOBL     = 9
            TABOBL(1) = 'KD1     '
            TABOBL(2) = 'KD2     '
            TABOBL(3) = 'KD3     '
            TABOBL(4) = 'V1X     '
            TABOBL(5) = 'V1Y     '
            TABOBL(6) = 'V1Z     '
            TABOBL(7) = 'V2X     '
            TABOBL(8) = 'V2Y     '
            TABOBL(9) = 'V2Z     '
          ENDIF
C
C       Diffusion anisotrope
        ELSEIF (IMATEE.EQ.3) THEN
          IF (IDIM.EQ.2) THEN
            IF (IFOMOD.NE.1) THEN
              JGOBL     = 5
              TABOBL(1) = 'KD11    '
              TABOBL(2) = 'KD22    '
              TABOBL(3) = 'KD21    '
              TABOBL(4) = 'V1X     '
              TABOBL(5) = 'V1Y     '
            ELSE
              JGOBL     = 6
              TABOBL(1) = 'KD11    '
              TABOBL(2) = 'KD22    '
              TABOBL(3) = 'KD21    '
              TABOBL(4) = 'KD33    '
              TABOBL(5) = 'V1X     '
              TABOBL(6) = 'V1Y     '
            ENDIF
          ELSEIF (IDIM.EQ.3) THEN
            JGOBL    = 12
            TABOBL( 1) = 'KD11    '
            TABOBL( 2) = 'KD22    '
            TABOBL( 3) = 'KD33    '
            TABOBL( 4) = 'KD21    '
            TABOBL( 5) = 'KD31    '
            TABOBL( 6) = 'KD32    '
            TABOBL( 7) = 'V1X     '
            TABOBL( 8) = 'V1Y     '
            TABOBL( 9) = 'V1Z     '
            TABOBL(10) = 'V2X     '
            TABOBL(11) = 'V2Y     '
            TABOBL(12) = 'V2Z     '
          ENDIF
        ENDIF
C
        JGOBL0 = JGOBL
        JGFAC0 = JGFAC
C
*       Loi FICK :
*       ------------
        IF (INATUU.EQ.0) THEN
*           Rien a ajouter
 
*       Loi SORET :
*       -------------
        ELSEIF (INATUU.EQ.1) THEN
C         -- Diffusion isotrope
          IF (IMATEE.EQ.1) THEN
            JGOBL            = JGOBL+1
C*8         TABOBL(JGOBL0+1) = 'KSORET  '
            TABOBL(JGOBL0+1) = 'KS      '
 
C         -- Diffusion orthotrope
          ELSEIF (IMATEE.EQ.2) THEN
C           ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
            IF (IFOMOD.LE.0) THEN
              JGOBL            = JGOBL + 2
C*8           TABOBL(JGOBL0+1) = 'KSORET1 '
C*8           TABOBL(JGOBL0+2) = 'KSORET2 '
              TABOBL(JGOBL0+1) = 'KS1     '
              TABOBL(JGOBL0+2) = 'KS2     '
 
C             ---- Elements massifs bidimensionnels FOURIER
C             ---- Elements massifs TRIDimensionnels
            ELSE
              JGOBL            = JGOBL + 3
C*8           TABOBL(JGOBL0+1) = 'KSORET1 '
C*8           TABOBL(JGOBL0+2) = 'KSORET2 '
C*8           TABOBL(JGOBL0+3) = 'KSORET3 '
              TABOBL(JGOBL0+1) = 'KS1     '
              TABOBL(JGOBL0+2) = 'KS2     '
              TABOBL(JGOBL0+3) = 'KS3     '
            ENDIF
 
C         -- Diffusion anisotrope
          ELSEIF (IMATEE.EQ.3) THEN
C             ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
            IF (IFOMOD.LE.0) THEN
              JGOBL            = JGOBL + 3
C*8           TABOBL(JGOBL0+1) = 'KSORET11'
C*8           TABOBL(JGOBL0+2) = 'KSORET22'
C*8           TABOBL(JGOBL0+3) = 'KSORET21'
              TABOBL(JGOBL0+1) = 'KS11    '
              TABOBL(JGOBL0+2) = 'KS22    '
              TABOBL(JGOBL0+3) = 'KS21    '
 
C             ---- Elements massifs bidimensionnels FOURIER
            ELSEIF (IFOMOD.EQ.1) THEN
              JGOBL            = JGOBL + 4
C*8           TABOBL(JGOBL0+1) = 'KSORET11'
C*8           TABOBL(JGOBL0+2) = 'KSORET22'
C*8           TABOBL(JGOBL0+3) = 'KSORET21'
C*8           TABOBL(JGOBL0+4) = 'KSORET33'
              TABOBL(JGOBL0+1) = 'KS11    '
              TABOBL(JGOBL0+2) = 'KS22    '
              TABOBL(JGOBL0+3) = 'KS21    '
              TABOBL(JGOBL0+4) = 'KS33    '
 
C             ---- Elements massifs TRIDimensionnels
            ELSEIF (IFOMOD.EQ.2) THEN
              JGOBL            = JGOBL + 6
C*8           TABOBL(JGOBL0+1) = 'KSORET11'
C*8           TABOBL(JGOBL0+2) = 'KSORET22'
C*8           TABOBL(JGOBL0+3) = 'KSORET33'
C*8           TABOBL(JGOBL0+4) = 'KSORET21'
C*8           TABOBL(JGOBL0+5) = 'KSORET31'
C*8           TABOBL(JGOBL0+6) = 'KSORET32'
              TABOBL(JGOBL0+1) = 'KS11    '
              TABOBL(JGOBL0+2) = 'KS22    '
              TABOBL(JGOBL0+3) = 'KS33    '
              TABOBL(JGOBL0+4) = 'KS21    '
              TABOBL(JGOBL0+5) = 'KS31    '
              TABOBL(JGOBL0+6) = 'KS32    '
            ENDIF
 
          ELSE
C           -- Erreur ne devant pas survenir
            CALL ERREUR(261)
            RETURN
          ENDIF
 
*       Loi UTILISATEUR :
*       -------------------
*       ELSEIF (INATUU.EQ.2) THEN
*       Par defaut pas de composantes autres que celles de l'utilisateur
 
        ELSEIF (INATUU.EQ.5) THEN
*       diffusion avec dissipation affine
          JGOBL         = JGOBL+2
          TABOBL(JGOBL0+1)='DAF0'
          TABOBL(JGOBL0+2)='DAF1'
 
        ELSEIF (INATUU.EQ.6) THEN
*       diffusion avec dissipation visqueuse
          JGOBL         = JGOBL+2
          TABOBL(JGOBL0+1)='DVQ0'
          TABOBL(JGOBL0+2)='DVQ1'
 
*       Loi inconnue : ERREUR
*       ----------------
        ELSE
          CALL ERREUR(261)
          RETURN
        ENDIF
 
*       Composantes ADVECTION :
*       -------------------
        CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
        IF (IADVE.NE.0) THEN
 
C         Cas des Tuyaux 1D (MFR=79)
C         ===========================
          IF (MFR .EQ. 79) THEN
            JGOBL             = JGOBL + 1
            TABOBL(JGOBL0 + 1)='VITE'
 
            JGFAC             = JGFAC + 1
            TABFAC(JGFAC0 + 1)='DL'
 
C         Cas Diffusion Massif
C         ============================
          ELSEIF (MFR .EQ. 1) THEN
            IF (IFOUR.EQ.-3) THEN
              JGOBL             = JGOBL + 3
              TABOBL(JGOBL0 + 1)='VITX'
              TABOBL(JGOBL0 + 2)='VITY'
              TABOBL(JGOBL0 + 3)='VITZ'
 
              JGFAC             = JGFAC + 3
              TABFAC(JGFAC0 + 1)='DX'
              TABFAC(JGFAC0 + 2)='DY'
              TABFAC(JGFAC0 + 3)='DZ'
            ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
              JGOBL             = JGOBL + 2
              TABOBL(JGOBL0 + 1)='VITX'
              TABOBL(JGOBL0 + 2)='VITY'
 
              JGFAC             = JGFAC + 2
              TABFAC(JGFAC0 + 1)='DX'
              TABFAC(JGFAC0 + 2)='DY'
            ELSEIF (IFOUR .EQ. 0) THEN
              JGOBL             = JGOBL + 2
              TABOBL(JGOBL0 + 1)='VITR'
              TABOBL(JGOBL0 + 2)='VITZ'
 
              JGFAC             = JGFAC + 2
              TABFAC(JGFAC0 + 1)='DR'
              TABFAC(JGFAC0 + 2)='DZ'
            ELSEIF (IFOUR .EQ. 1) THEN
              CALL ERREUR(21)
              RETURN
C             JGOBL             = JGOBL + 3
C             TABOBL(JGOBL0 + 1)='VITR'
C             TABOBL(JGOBL0 + 2)='VITZ'
C             TABOBL(JGOBL0 + 3)='VITT'
C
C             JGFAC             = JGFAC + 3
C             TABFAC(JGFAC0 + 1)='DR'
C             TABFAC(JGFAC0 + 2)='DZ'
C             TABFAC(JGFAC0 + 3)='DT'
            ELSEIF (IFOUR .EQ. 2) THEN
              JGOBL             = JGOBL + 3
              TABOBL(JGOBL0 + 1)='VITX'
              TABOBL(JGOBL0 + 2)='VITY'
              TABOBL(JGOBL0 + 3)='VITZ'
 
              JGFAC             = JGFAC + 3
              TABFAC(JGFAC0 + 1)='DX'
              TABFAC(JGFAC0 + 2)='DY'
              TABFAC(JGFAC0 + 3)='DZ'
            ELSEIF (IFOUR.GE.3 .AND. IFOUR.LE.15) THEN
              IF (IFOUR.LE.6) THEN
                JGOBL             = JGOBL + 1
                TABOBL(JGOBL0 + 1)='VITX'
 
                JGFAC             = JGFAC + 1
                TABFAC(JGFAC0 + 1)='DX'
 
              ELSEIF (IFOUR.EQ.7 .OR. IFOUR.EQ.8) THEN
                JGOBL             = JGOBL + 2
                TABOBL(JGOBL0 + 1)='VITX'
                TABOBL(JGOBL0 + 2)='VITY'
 
                JGFAC             = JGFAC + 2
                TABFAC(JGFAC0 + 1)='DX'
                TABFAC(JGFAC0 + 2)='DY'
 
              ELSEIF (IFOUR.EQ.9 .OR. IFOUR.EQ.10) THEN
                JGOBL             = JGOBL + 2
                TABOBL(JGOBL0 + 1)='VITX'
                TABOBL(JGOBL0 + 2)='VITZ'
 
                JGFAC             = JGFAC + 2
                TABFAC(JGFAC0 + 1)='DX'
                TABFAC(JGFAC0 + 2)='DZ'
 
              ELSEIF (IFOUR.EQ.11) THEN
                JGOBL             = JGOBL + 3
                TABOBL(JGOBL0 + 1)='VITX'
                TABOBL(JGOBL0 + 2)='VITY'
                TABOBL(JGOBL0 + 3)='VITZ'
 
                JGFAC             = JGFAC + 3
                TABFAC(JGFAC0 + 1)='DX'
                TABFAC(JGFAC0 + 2)='DY'
                TABFAC(JGFAC0 + 3)='DZ'
 
              ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN
                JGOBL             = JGOBL + 1
                TABOBL(JGOBL0 + 1)='VITR'
 
                JGFAC             = JGFAC + 1
                TABFAC(JGFAC0 + 1)='DR'
 
              ELSEIF (IFOUR.EQ.14) THEN
                JGOBL             = JGOBL + 2
                TABOBL(JGOBL0 + 1)='VITR'
                TABOBL(JGOBL0 + 2)='VITZ'
 
                JGFAC             = JGFAC + 2
                TABFAC(JGFAC0 + 1)='DR'
                TABFAC(JGFAC0 + 2)='DZ'
              ENDIF
            ENDIF
          ENDIF
        ENDIF
C     =================================================================
C     FORMULATION CHARGEMENT
C     =================================================================
      ELSE IF (MFR2.EQ.72) THEN
C
        IF (IMATEE.EQ.1) THEN
C         RAJOUTER DIFFERENTIATION EN FONCIONS DE LA FORMULATION MASSIF -> P
C         COQUE PINF PSUP
          JGOBL=1
          TABOBL(1)='PR  '
        ELSE
          IRET = 0
          CALL ERREUR(5)
        ENDIF
C     =================================================================
C     FORMULATION CHANGEMENT_PHASE
C     =================================================================
      ELSE IF (MFR2.EQ.30) THEN
C
        IF (IMODEL.CMATEE.EQ.'CHPH_PAR') THEN
          JGOBL=2
          TABOBL(1)='PRIM'
          TABOBL(2)='DUAL'
 
          JGFAC=0
        ELSEIF (IMODEL.CMATEE.EQ.'CHPH_SOL') THEN
          JGOBL=1
          TABOBL(1)='SOLU'
 
          JGFAC=0
        ENDIF
C     =================================================================
      ENDIF
*--------------------------------------------------------------------
*     CREATION DU NOMID A PARTIR DE TABOBL ET TABFAC
*--------------------------------------------------------------------
 9999 CONTINUE
 
      IF((MFR2.EQ.2) .OR. (MFR2.EQ.33)) THEN
C       CB215821 : Ajout de TREF et TALP à la fin des compostantes facultatives
C                  pour les FORMULATION MECANIQUE, POREUX
C                  Pour ne pas casser l'ordre dans COMP (com2/coml6/comval/comara/...)
        JGFAC=JGFAC+2
        TABFAC(JGFAC-1) ='TREF'
        TABFAC(JGFAC)   ='TALP'
      ENDIF
 
C     Erreur si JGOBL ou JGFAC sont superieurs a ITA
C     (Passage en FORTRAN 77 car la compilation depasse la memoire sur
C     Windows-32bits)
      IF ((JGOBL .GT. ITA) .OR. (JGFAC .GT. ITA)) THEN
        IRET = 0
        CALL ERREUR(5)
        RETURN
      ENDIF
C
C     SI PROBLEME (IRET = 0), ON SORT AVEC IPNOMC A 0
C
      IF (IRET.NE.0) THEN
        NBROBL = JGOBL
        NBRFAC = JGFAC
        SEGINI,NOMID
        IPNOMC=NOMID
C
        DO 100 IO = 1,NBROBL
          NOMID.LESOBL(IO) = TABOBL(IO)
 100    CONTINUE
C
        DO 110 IO=1,NBRFAC
          NOMID.LESFAC(IO) = TABFAC(IO)
110     CONTINUE
C
        SEGACT,NOMID*NOMOD
      ENDIF
      END