Sources Esope | Cast3M

prquoi
C PRQUOI    SOURCE    OF166741  26/02/23    21:15:26     12480          
C=======================================================================
C=                  Sous-programme PRQUOI                              =
C=                  ---------------------                              =
C=                                                                     =
C= Ce sous programme fournit les caracteristiques des elements finis.  =
C=                                                                     =
C= Entrees :                                                           =
C= ---------                                                           =
C=   IMODEL    : pointeur sur un IMODEL, zone elementaire d'un MMODEL  =
C=                                                                     =
C= Sorties :                                                           =
C= ---------                                                           =
C=  L'element INFELE du segment IMODEL a ete renseigne et contient :   =
C=   INFELE( 1) : numero de l'element fini                             =
C=   INFELE( 2) : nb de points d'integration en contrainte multicouche =
C=   INFELE( 3) : nb de points d'integration pour MASSE                =
C=   INFELE( 4) : nb de points d'integration pour SIGMA BSIGMA KSIGMA  =
C=   INFELE( 5) : nb de caracteristiques                               =
C=   INFELE( 6) : nb de points d'integration pour RIGI/COND            =
C=   INFELE( 7) : longueur d'un tableau de travail pour l'element      =
C=   INFELE( 8) : nombre de fonctions de forme                         =
C=   INFELE( 9) : nb de d.d.l. dans la matrice de RIGIDITE             =
C=   INFELE(10) : taille de la matrice de Hooke                        =
C=   INFELE(11) : pointeur sur le segment d'integration                =
C=   INFELE(12) : inutilise                                            =
C=   INFELE(13) : numero de la formulation de l'element fini (NUMMFR)  =
C=   INFELE(14) : numero de l'element geometrique associe    (NUMGEO)  =
C=   INFELE(15) : nb maximal de d.d.l. par noeud                       =
C=   INFELE(16) : nb de composantes de contraintes ou de deformations  =
C=======================================================================
      SUBROUTINE PRQUOI(IMODEL)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
-INC PPARAM
-INC CCOPTIO
-INC CCGEOME
-INC CCHAMP
C= DEB === Formulation HHO ===== Include specifique ====================
-INC CCHHOPA
C= FIN === Formulation HHO =============================================
-INC SMMODEL
-INC SMELEME
-INC SMLMOTS
C
      logical dcmate,dcmat2,bmodi
C
      IF (IMODEL.INFMOD(/1).LT.1) THEN
        WRITE(IOIMP,*)
        WRITE(IOIMP,*) 'PRQUOI : INFMOD(/1) incorrect',INFMOD(/1),' &lt; 1'
        CALL ERREUR(5)
      ENDIF
C
C (0) petite verrue pour la formulation CONTACT
      IF (IMODEL.FORMOD(1) .EQ. 'CONTACT         ') RETURN
 
C     Recuperation des donnees de INFELE en entree
      MELE = IMODEL.NEFMOD
      IELE = NUMGEO(MELE)
      MFR  = NUMMFR(MELE)
C
C (1) petite verrue pour les MULT
      IF (MELE .EQ. 22) RETURN
C (2) petite verrue pour les SURE
C     Travail fait par ailleurs :
      IF (MELE .EQ. 259) RETURN
C==== Formulation HHO ===== Appel specifique a HHOELQ ==================
      IF (MELE .EQ. HHO_NUM_ELEMENT) THEN
C test equivalent      IF (MFR .EQ. HHO_MFR_ELEMENT) THEN
        CALL HHOELQ(IMODEL)
        RETURN
      END IF
C==== Formulation HHO ==================================================
 
      dcmate = .false.
      dcmat2 = .false.
      bmodi=.false.
      do im = 1, imodel.matmod(/2)
        if (imodel.matmod(im).eq.'IMPEDANCE') then
          dcmate = .true.
          meleme = imodel.imamod
          ityp1 = meleme.itypel
          mfr = nummfr(mele)
          if (ityp1.eq.1) then
            bmodi=.true.
            mele = 45
            melep=imodel.nefmod
            imodel.nefmod=mele
          endif
          iele = numgeo(mele)
        endif
      enddo
C
C     Nb de points d'integration dans l'epaisseur
      NPINT = IMODEL.INFMOD(1)
C
C     =================================================================
C     INFELE(1)  : Numero de l'element fini
C     =================================================================
      INFELE(1) = MELE
C     =================================================================
C     INFELE(8)  : Nombre de fonctions de forme de l'element
C     =================================================================
      INFELE(8) = NBNNE(IELE)
C     =================================================================
C     INFELE(13) : Numero de la formulation de l'element fini
C     =================================================================
      INFELE(13) = MFR
      IF (FORMOD(1).EQ.'LIAISON'.AND.MFR.EQ.27)    INFELE(13) = 24
      DO IM = 1,MATMOD(/2)
       IF (MFR.EQ.27.AND.MATMOD(IM).EQ.'MODAL')    INFELE(13) = 26
       IF (MFR.EQ.27.AND.MATMOD(IM).EQ.'STATIQUE') INFELE(13) = 28
      ENDDO
      IF (FORMOD(1).EQ.'CHARGEMENT') THEN
        IF (MFR.EQ.1) INFELE(13) = 72
        IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.7.OR.MFR.EQ.9.OR.
     &      MFR.EQ.13) INFELE(13) = 74
      ENDIF
      MFR = INFELE(13)
      IF (MELE.GE.99.AND.MELE.LE.103) INFELE(13)=1
 
C     ELECTROSTATIQUE : MFR=71=INFELE(13) et NUMMFR(mele)=1 (MASSIF)
      IF (FORMOD(1).EQ.'ELECTROSTATIQUE') INFELE(13)=71
 
C     DIFFUSION : MFR=73=INFELE(13) et NUMMFR(mele)=1 (MASSIF)
      IF (FORMOD(1).EQ.'DIFFUSION') THEN
        IF (MFR.EQ.1) INFELE(13)=73
      ENDIF
C     =================================================================
C     INFELE(14) : Numero de l'element geometrique
C     =================================================================
      INFELE(14) = NUMGEO(MELE)
 
C     Redirection suivant l'element fini
C     ----------------------------------------------------------
      IF(MELE.GE.1.AND.MELE.LE.100) THEN
C            CABL SEG2 SEG3 TRI3 TRI4 TRI6 TRI7 QUA4 QUA5 QUA8
      GOTO (  300,   2,   3,   4, 300,   6,   7,   8, 300,  10
C            QUA9 RAC2 RAC3 CUB8 CU20 PRI6 PR15 LIA3 LIA4 LIA6
     1      ,  11,  12,   3,  14,  15,  16,  17,  18,  19,   6
C            LIA8 MULT TET4 TE10 PYR5 PY13 COQ3  DKT POUT LISP
     2      ,  10, 300,  23,  24,  25,  26,  27,  28,  29,  30
C            FAC3 FAC4 FAC6 FAC8 LTR3 LQU4 LCU8 LPR6 LTE4 LPY5
     3      ,  31,  32,  33,  34,  35,  36,  37,  38,  39,  40
C            COQ8 TUYA TUFI COQ2 POI1 BARR RACO LSU2 COQ4 LISM
     4      ,  41,  42,  43,  44,  45,  46,  47,  48,  49,  50
C            COF3 RES2 LSU3 LSU4 LICO COQ6 CVS2 CVS3 CVT3 CVT6
     5      ,  51, 300,  53,  54,  55,  56,  57,  58,  59,  60
C            CVQ4 CVQ8 THP5 TH13 THP6 TH15 THC8 TH20 ICT3 ICQ4
     6      ,  61,  62,  63,  64,  65,  66,  67,  68,  69,  70
C            ICT6 ICQ8 ICC8 ICT4 ICP6 IC20 IC10 IC15 TRIP QUAP
     7      ,  71,  72,  73,  74,  75,  76,  77,  78,  79,  80
C            CUBP TETP PRIP TIMO JOI2 JOI3 JOT3 JOI4 JOI6 JOI8
     8      ,  81,  82,  83,  84,  85,  86,  87,  88,  89,  90
C            LISC TRIH  DST LIC4 CERC TUYO LSE2 LITU HYT3 HYQ4
     9      ,  91,  92,  93,  94,  95,  96,  97,  98,   4,   8)
C
     &      ,MELE
      ELSEIF(MELE.GE.101.AND.MELE.LE.200) THEN
C            HYT4 HYP6 HYC8 TRIS QUAS POIS FOR3 JOP3 JOP6 JOP8
      GOTO (   23,  16,  14, 104, 105, 106, 300, 108, 109, 110
C            POL3 POL4 POL5 POL6 POL7 POL8 POL9 PO10 PO11 PO12
     1      , 111, 111, 111, 111, 111, 111, 111, 111, 111, 111
C            PO13 PO14 BAR3 BAEX LIA2 QUAH CUBH ROT3 SEF2 TRF3
     2      , 111, 111, 123, 124, 125, 126, 127, 128,   2,   4
C            QUF4 CUF8 PRF6 TEF4 PYF5 MSE3 MTR6 MQU9 MC27 MP18
     3      ,   8,  14,  16,  23,  25,   3,   6,  11, 275, 300
C            MT10 MP14 SEF3 TRF7 QUF9 CF27 PF21 TF15 PF19 SEG6
     4      ,  24, 300,   3,   7,  11, 275, 276, 277, 278, 300
C            TR21 QU36 C216 P126 TE56 PY91 TRH6 BSE2 BTR4 BQU5
     5      , 300, 300, 300, 300, 300, 300, 157, 300, 300, 300
C            BCU9 BPR7 BTE5 BPY6 FRO4 SEGS POJS JCT3 JCI4 JGI2
     6      , 300, 300, 300, 300, 300, 166, 167, 168, 169, 170
C            JGT3 JGI4 TRIQ QUAQ CUBQ TETQ PRIQ TRIR QUAR CUBR
     7      , 171, 172, 173, 174, 175, 176, 177, 178, 179, 180
C            TETR PRIR Q4RI Q8RI JOQ3 JOQ6 JOQ8 JOR3 JOR6 JOR8
     8      , 181, 182, 183, 184, 185, 186, 187, 188, 189, 190
C            T1D2 T1D3 M1D2 M1D3 LC03 LC07 LC09 LC27 LC21 LC15
     9      , 191, 192, 191, 192, 195, 196, 197, 300, 300, 300)
C
     &      ,MELE-100
      ELSEIF(MELE.GE.201.AND.MELE.LE.300) THEN
C            LC19 LS03 LS07 LS09 LS27 LS21 LS15 LS19 BS03 BS07
      GOTO (  300, 195, 196, 197, 300, 300, 300, 300, 195, 196
C            BS09 BS27 BS21 BS15 BS19 MC03 MC07 MC09 MC27 MC21
     1      , 197, 300, 300, 300, 300, 216, 217, 218, 300, 300
C            MC15 MC19 M103 M107 M109 M127 M121 M115 M119 MS03
     2      , 300, 300, 223, 224, 225, 300, 300, 300, 300, 223
C            MS07 MS09 MS27 MS21 MS15 MS19 QC03 QC07 QC09 QC27
     3      , 224, 225, 300, 300, 300, 300, 237, 238, 239, 300
C            QC21 QC15 QC19 Q103 Q107 Q109 Q127 Q121 Q115 Q119
     4      , 300, 300, 300, 237, 238, 239, 300, 300, 300, 300
C            QS03 QS07 QS09 QS27 QS21 QS15 QS19 CIFL SURE SHB8
     5      , 237, 238, 239, 300, 300, 300, 300, 300, 300, 260
C            CAF2 CAF3 XQ4R XC8R JOI1 ZCO2 ZCO3 ZCO4 TUY2 TUY3
     6      , 300, 300, 263, 264, 265, 266, 267, 268, 269, 270
C            COS2 COA2 ICY5 IC13 CU27 PR21 TE15 PY19 C20R P15R
     7      , 271, 272,  25,  26, 275, 276, 277, 278, 279, 280
C            HHO  .... .... .... .... .... .... .... .... ....
     8      , 999, 999, 999, 999, 999, 999, 999, 999, 999, 999
C            .... .... .... .... .... .... .... .... .... ....
     9      , 999, 999, 999, 999, 999, 999, 999, 999, 999, 999)
C
     &      ,MELE-200
      ENDIF
C
C     =================================================================
C     ERREUR
 999  CONTINUE
      CALL ERREUR(5)
C
C     =================================================================
C     SEG2
C     =================================================================
 2    CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(6)=2
      INFELE(7)=0
      INFELE(9)=4
      INFELE(10)=1
      INFELE(15)=2
      INFELE(16)=1
      IF (IFOUR.EQ.1) THEN
        INFELE(9)=6
        INFELE(15)=3
      ENDIF
      GOTO 300
C     =================================================================
C     SEG3
C     =================================================================
 3    CONTINUE
      INFELE(2)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(7)=0
      INFELE(9)=6
      INFELE(10)=1
      INFELE(15)=2
      INFELE(16)=1
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
        IF (IFOUR.EQ.1) THEN
          INFELE(9)=9
          INFELE(15)=3
        ENDIF
      ENDIF
      IF (MELE.EQ.99) THEN
        INFELE(3)=4
        INFELE(6)=4
      ENDIF
      GOTO 300
C     =================================================================
C     TRI3
C     =================================================================
 4    CONTINUE
      INFELE(2)=1
      INFELE(3)=4
      INFELE(4)=1
      INFELE(6)=1
      INFELE(9)=6
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1)  INFELE(3)=7
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=9
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     TRI6
C     =================================================================
 6    CONTINUE
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=7
        INFELE(4)=7
        INFELE(6)=7
      ELSE
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
      ENDIF
      INFELE(3)=7
      INFELE(9)=12
      IF (IFOUR.EQ.-3) INFELE(9)=15
      IF (IFOUR.EQ.1)  INFELE(9)=18
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     TRI7
C     =================================================================
 7    CONTINUE
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=7
        INFELE(4)=7
        INFELE(6)=7
      ELSE
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
      ENDIF
      INFELE(3)=7
      INFELE(9)=14
      IF (IFOUR.EQ.-3) INFELE(9)=17
      IF (IFOUR.EQ.1)  INFELE(9)=21
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     QUA4
C     =================================================================
 8    CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=4
C**   INFELE(4)=5    CORRECTION  MILL  LE 8/8/90
      INFELE(6)=4
C**   IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9
      INFELE(9)=8
      IF (IFOUR.EQ.-3) INFELE(9)=11
      IF (IFOUR.EQ.1)  INFELE(9)=12
      IF (MELE.EQ.100) THEN
        INFELE(3)=4
        INFELE(6)=4
      ENDIF
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     QUA8
C     =================================================================
 10   CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16
      INFELE(4)=9
      INFELE(6)=9
C     IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(9)=16
      IF (IFOUR.EQ.-3) INFELE(9)=19
      IF (IFOUR.EQ.1)  INFELE(9)=24
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     QUA9
C     =================================================================
 11   CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16
      INFELE(4)=9
      INFELE(6)=9
C     IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(9)=18
      IF (IFOUR.EQ.-3) INFELE(9)=21
      IF (IFOUR.EQ.1)  INFELE(9)=27
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     RAC2 (raccord a 4 noeuds en 2 dimensions)
C     =================================================================
 12   CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(5)=2
      INFELE(6)=2
      INFELE(9)=8
      IF (IFOUR.EQ.1) INFELE(9)=10
      INFELE(10)=6
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     CUB8
C     =================================================================
 14   CONTINUE
      INFELE(2)=8
      INFELE(3)=8
      INFELE(4)=8
      INFELE(6)=8
      INFELE(9)=24
      IF (MELE.EQ.103) THEN
        INFELE(3)=8
        INFELE(6)=8
      ENDIF
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     CU20
C     =================================================================
 15   CONTINUE
      INFELE(2)=27
      INFELE(3)=27
      INFELE(4)=27
      INFELE(6)=27
      INFELE(9)=60
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     PRI6
C     =================================================================
 16   CONTINUE
      INFELE(2)=6
      INFELE(3)=8
      INFELE(4)=6
      INFELE(6)=6
      INFELE(9)=18
      IF (MELE.EQ.102) THEN
        INFELE(3)=6
        INFELE(6)=6
      ENDIF
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     PR15
C     =================================================================
 17   CONTINUE
      INFELE(2)=12
      INFELE(3)=21
      INFELE(4)=12
      INFELE(6)=12
      INFELE(9)=45
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     LIA3 (raccord a 6 noeudes en 3 dimensions)
C     =================================================================
 18   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(5)=3
      INFELE(6)=4
      INFELE(9)=15
      INFELE(10)=6
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     LIA4 (raccord a 8 noeuds en 3 dimensions)
C     =================================================================
 19   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(5)=3
      INFELE(6)=4
      INFELE(9)=20
      INFELE(10)=6
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     TET4
C     =================================================================
 23   CONTINUE
      INFELE(2)=1
      INFELE(3)=4
      INFELE(4)=1
      INFELE(6)=1
      INFELE(9)=12
      IF (MELE.EQ.101) THEN
        INFELE(3)=4
        INFELE(6)=4
      ENDIF
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     TE10
C     =================================================================
 24   CONTINUE
      INFELE(2)=5
cbp      INFELE(3)=5
      INFELE(3)=15
      INFELE(4)=5
      INFELE(6)=5
      INFELE(9)=30
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     PYR5
C     =================================================================
 25   CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(9)=15
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     PY13
C     =================================================================
 26   CONTINUE
      INFELE(2)=10
      INFELE(3)=17
      INFELE(4)=10
      INFELE(6)=10
      INFELE(9)=39
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      IF (MELE.EQ.274) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ENDIF
      GOTO 300
C     =================================================================
C     COQ3
C     =================================================================
 27   CONTINUE
      INFELE(2)=3
      INFELE(4)=1
      INFELE(5)=3
      INFELE(6)=1
      INFELE(7)=600
      INFELE(9)=18
      IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=9
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     DKT
C     =================================================================
 28   CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(5)=3
      INFELE(7)=600
      INFELE(9)=18
      INFELE(15)=6
      IF (NPINT.EQ.0) THEN
C  Element DKT COQUE GLOBALE
        INFELE(4)=3
        INFELE(6)=3
        INFELE(10)=6
        INFELE(16)=6
      ELSE
C  Element DKT COQUE INTEGREE dans l'epaisseur
        INFELE(4)=3*NPINT
        INFELE(6)=3*NPINT
        INFELE(10)=6
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     POUTRE
C     =================================================================
 29   CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(6)=2
      INFELE(7)=200
      IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
C Element POUTRE 2D plan
        INFELE(5)=6
        INFELE(9)=6
        INFELE(10)=3
        INFELE(15)=3
        INFELE(16)=3
      ELSE
C Element POUTRE 3D
        INFELE(5)=12
        INFELE(9)=12
        INFELE(10)=6
        INFELE(15)=6
        INFELE(16)=6
      ENDIF
      GOTO 300
C     =================================================================
C     LINE SPRING LISP
C     =================================================================
 30   CONTINUE
      INFELE(2)=3
      INFELE(4)=3
      INFELE(5)=5
      INFELE(6)=3
      INFELE(7)=100
      INFELE(9)=24
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     FAC3 pour les PRESSIONS uniquement
C     =================================================================
 31   CONTINUE
      INFELE(2)=1
      INFELE(4)=1
      INFELE(6)=1
      INFELE(9)=9
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     FAC4 pour les PRESSIONS uniquement
C     =================================================================
 32   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=12
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     FAC6 pour les PRESSIONS uniquement
C     =================================================================
 33   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=18
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     FAC8 pour les PRESSIONS uniquement
C     =================================================================
 34   CONTINUE
      INFELE(2)=9
      INFELE(4)=9
      INFELE(6)=9
      INFELE(9)=24
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     LTR3 : element LIQUIDE triangle a 3 noeuds
C     =================================================================
 35   CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=7
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=6
      INFELE(10)=5
      INFELE(15)=2
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     LQU4 : element LIQUIDE quadrangle a 4 noeuds
C     =================================================================
 36   CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=8
      INFELE(10)=5
      INFELE(15)=2
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     LCU8 : Element LIQUIDE CUBE a 8 noeuds
C     =================================================================
 37   CONTINUE
      INFELE(2)=8
      INFELE(3)=8
      INFELE(4)=8
      INFELE(6)=8
      INFELE(9)=16
      INFELE(10)=5
      INFELE(15)=2
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     LPR6 : Element LIQUIDE PRISME a 6 noeuds
C     =================================================================
 38   CONTINUE
      INFELE(2)=6
      INFELE(3)=6
      INFELE(4)=6
      INFELE(6)=6
      INFELE(9)=12
      INFELE(10)=5
      INFELE(15)=2
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     LTE4 : Element LIQUIDE TETRAEDRE a 4 noeuds
C     =================================================================
 39   CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=8
      INFELE(10)=5
      INFELE(15)=2
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     LPY5 : Element LIQUIDE PYRAMIDE a 5 noeuds
C     =================================================================
 40   CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(9)=10
      INFELE(10)=5
      INFELE(15)=2
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     COQ8 : Element COQUE EPAISSE a 8 noeuds
C     =================================================================
 41   CONTINUE
      INFELE(2)=NPINT
      IF (NPINT.EQ.0) INFELE(2)=8
      INFELE(3)=8
      INFELE(4)=NPINT
      IF (NPINT.EQ.0) INFELE(4)=8
      INFELE(5)=2
      INFELE(6)=8
      INFELE(7)=1700
      INFELE(9)=48
      IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=24
      INFELE(10)=5
      INFELE(15)=6
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     TUYA : Element TUYAU
C     =================================================================
 42   CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(5)=7
      INFELE(6)=2
      INFELE(7)=200
      INFELE(9)=12
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     TUFI : Element TUYAU FISSURE
C     =================================================================
 43   CONTINUE
      INFELE(2)=1
      INFELE(3)=0
      INFELE(4)=1
      INFELE(5)=9
      INFELE(6)=1
      INFELE(7)=200
      INFELE(9)=12
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=8
      GOTO 300
C     =================================================================
C     COQ2
C     =================================================================
 44   CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(5)=3
      INFELE(6)=3
      INFELE(7)=368
      INFELE(9)=8
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(9)=6
      IF (IFOUR.EQ.-3) INFELE(9)=9
      INFELE(10)=6
      IF (IFOUR.LE.0) INFELE(10)=4
      INFELE(15)=4
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(15)=3
      INFELE(16)=6
      IF (IFOUR.LE.0) INFELE(16)=4
      GOTO 300
C     =================================================================
C     POI1 (MODEs en DEFOrmations PLANes GENEralisees)
C     =================================================================
 45   CONTINUE
C  Valeurs a adapter pour certains MODE UNID (1D)
      INFELE(2)=1
      INFELE(3)=1
      INFELE(4)=1
      INFELE(5)=1
      INFELE(6)=1
      INFELE(7)=50
      INFELE(15)=2
      IF (IFOUR.EQ.11) THEN
        INFELE(9)=3
        INFELE(10)=2
        INFELE(16)=2
      ELSE
        INFELE(9)=5
        IF (IFOUR.GE.3.AND.IFOUR.LE.15) INFELE(9)=2
        INFELE(10)=1
        INFELE(16)=1
      ENDIF
      if (mfr.eq.26.or.mfr.eq.28) INFELE(9) = 1
      GOTO 300
C     =================================================================
C     BARRE
C     =================================================================
 46   CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(5)=1
      INFELE(6)=2
      INFELE(7)=200
      INFELE(9)=6
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(9)=4
      IF (IFOUR.EQ.-3) INFELE(9)=7
      IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=2
      INFELE(10)=1
      INFELE(15)=3
      IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(15)=1
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(15)=2
      INFELE(16)=1
      GOTO 300
C     =================================================================
C     RACO : Element de RACCORD LIQUIDE-COQUE a 4 noeuds en 2D
C     =================================================================
 47   CONTINUE
      INFELE(2)=1
      INFELE(4)=1
      INFELE(5)=2
      INFELE(6)=1
      INFELE(9)=10
      IF (IFOUR.EQ.1) INFELE(9)=12
      INFELE(15)=4
      INFELE(10)=6
      GOTO 300
C     =================================================================
C     LSU2 : Element de SURFACE LIBRE a 2 noeuds en 2D
C     =================================================================
 48   CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(6)=2
      INFELE(9)=6
      INFELE(10)=6
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     COQ4
C     =================================================================
 49   CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(5)=3
      INFELE(6)=5
      INFELE(7)=100
      INFELE(9)=24
      IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=12
      INFELE(10)=8
      INFELE(15)=6
      INFELE(16)=8
      GOTO 300
C     =================================================================
C     LINE SPRING LISM
C     =================================================================
 50   CONTINUE
      INFELE(2)=1
      INFELE(4)=1
      INFELE(5)=5
      INFELE(6)=3
      INFELE(7)=100
      INFELE(9)=24
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     COF3
C     =================================================================
 51   CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(5)=3
      INFELE(6)=3
      INFELE(7)=368
      INFELE(9)=12
      INFELE(10)=6
      INFELE(15)=4
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     LSU3 : Element de SURFACE LIBRE a 3 noeuds en 3D
C     =================================================================
 53   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=9
      INFELE(10)=6
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     LSU4 : Element de SURFACE LIBRE a 4 noeuds en 3D
C     =================================================================
 54   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(9)=12
      INFELE(10)=6
      INFELE(15)=3
      GOTO 300
C     =================================================================
C     LICO : Element de RACCORD LIQUIDE-COQUE a 6 noeuds en 3D
C     =================================================================
 55   CONTINUE
      INFELE(2)=3
      INFELE(4)=3
      INFELE(5)=3
      INFELE(6)=3
      INFELE(9)=24
      INFELE(10)=6
      INFELE(15)=6
      GOTO 300
C     =================================================================
C     COQ6 : Element COQUE EPAISSE a 6 noeuds
C     =================================================================
 56   CONTINUE
      INFELE(2)=NPINT
      IF(NPINT.EQ.0) INFELE(2)=6
      INFELE(3)=6
      INFELE(4)=NPINT
      IF (NPINT.EQ.0) INFELE(4)=6
      INFELE(5)=2
      INFELE(6)=6
      INFELE(7)=1000
      INFELE(9)=36
      IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=18
      INFELE(10)=5
      INFELE(15)=6
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     THT3 (thermique)
C     =================================================================
 57   CONTINUE
      INFELE(2)=1
      INFELE(6)=1
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=4
      INFELE(9)=3
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=9
      INFELE(10)=2
      INFELE(15)=1
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     THT6
C     =================================================================
 58   CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(6)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=7
      INFELE(9)=6
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=18
      INFELE(10)=2
      INFELE(15)=1
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     THQ4
C     =================================================================
 59   CONTINUE
      INFELE(2)=4
      INFELE(6)=4
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9
      INFELE(9)=4
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=12
      INFELE(10)=2
      INFELE(15)=1
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     THQ8
C     =================================================================
 60   CONTINUE
      INFELE(2)=9
      INFELE(6)=9
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(9)=8
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=24
      INFELE(10)=2
      INFELE(15)=1
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     THT4
C     =================================================================
 61   CONTINUE
      INFELE(2)=1
      INFELE(6)=1
      INFELE(9)=4
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     TH10
C     =================================================================
 62   CONTINUE
      INFELE(2)=5
      INFELE(6)=5
      INFELE(9)=10
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     THP5
C     =================================================================
 63   CONTINUE
      INFELE(2)=5
      INFELE(6)=5
      INFELE(9)=5
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     TH13
C     =================================================================
 64   CONTINUE
      INFELE(2)=27
      INFELE(6)=27
      INFELE(9)=13
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     THP6
C     =================================================================
 65   CONTINUE
      INFELE(2)=2
      INFELE(6)=2
      INFELE(9)=6
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     TH15
C     =================================================================
 66   CONTINUE
      INFELE(2)=12
      INFELE(6)=12
      INFELE(9)=15
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     THC8
C     =================================================================
 67   CONTINUE
      INFELE(2)=8
      INFELE(6)=8
      INFELE(9)=8
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     TH20
C     =================================================================
 68   CONTINUE
      INFELE(2)=27
      INFELE(6)=27
      INFELE(9)=20
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     ICT3
C     =================================================================
 69   CONTINUE
      INFELE(2)=1
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=7
      INFELE(4)=1
      INFELE(6)=1
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=4
      INFELE(9)=6
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=9
      INFELE(10)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4
      INFELE(15)=2
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3
      INFELE(16)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4
      GOTO 300
C     =================================================================
C     ICQ4
C     =================================================================
 70   CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=4
C**   INFELE(4)=5    CORRECTION  MILL  LE 8/8/90
      INFELE(6)=4
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9
      INFELE(9)=8
      IF (IFOUR.EQ.-3) INFELE(9)=11
      IF (IFOUR.EQ.1) INFELE(9)=12
      INFELE(10)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4
      INFELE(15)=2
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3
      INFELE(16)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4
      GOTO 300
C     =================================================================
C     ICT6
C     =================================================================
 71   CONTINUE
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=7
        INFELE(4)=7
        INFELE(6)=7
      ELSE
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
      ENDIF
      INFELE(3)=7
      INFELE(9)=12
      IF (IFOUR.EQ.-3) INFELE(9)=15
      IF (IFOUR.EQ.1) INFELE(9)=18
      INFELE(10)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4
      INFELE(15)=2
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3
      INFELE(16)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4
      GOTO 300
C     =================================================================
C     ICQ8
C     =================================================================
 72   CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16
      INFELE(4)=9
      INFELE(6)=9
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(9)=16
      IF (IFOUR.EQ.-3) INFELE(9)=19
      IF (IFOUR.EQ.1) INFELE(9)=24
      INFELE(10)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4
      INFELE(15)=2
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3
      INFELE(16)=4
C**   IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4
      GOTO 300
C     =================================================================
C     ICC8
C     =================================================================
 73   CONTINUE
      INFELE(2)=8
      INFELE(3)=8
      INFELE(4)=8
      INFELE(6)=8
      INFELE(9)=24
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     ICT4
C     =================================================================
 74   CONTINUE
      INFELE(2)=1
      INFELE(3)=4
      INFELE(4)=1
      INFELE(6)=1
      INFELE(9)=12
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     ICP6
C     =================================================================
 75   CONTINUE
      INFELE(2)=6
      INFELE(3)=8
      INFELE(4)=6
      INFELE(6)=6
      INFELE(9)=18
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     IC20
C     =================================================================
 76   CONTINUE
      INFELE(2)=27
      INFELE(3)=27
      INFELE(4)=27
      INFELE(6)=27
      INFELE(9)=60
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     IC10
C     =================================================================
 77   CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(9)=30
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     IC15
C     =================================================================
 78   CONTINUE
      INFELE(2)=12
      INFELE(3)=21
      INFELE(4)=12
      INFELE(6)=12
      INFELE(9)=45
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     TRIP
C     =================================================================
 79   CONTINUE
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=7
        INFELE(4)=7
        INFELE(6)=7
      ELSE
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
      ENDIF
      INFELE(3)=7
      INFELE(8)=9
      INFELE(9)=15
      IF (IFOUR.EQ.-3.OR.IFOUR.EQ.1) INFELE(9)=21
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN
        INFELE(10)=6
        INFELE(15)=4
        INFELE(16)=7
      ELSE
        INFELE(10)=4
        INFELE(15)=3
        INFELE(16)=5
      ENDIF
      GOTO 300
C     =================================================================
C     QUAP
C     =================================================================
 80   CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=9
      INFELE(6)=9
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(8)=12
      INFELE(9)=20
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=28
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN
        INFELE(10)=6
        INFELE(15)=4
        INFELE(16)=7
      ELSE
        INFELE(10)=4
        INFELE(15)=3
        INFELE(16)=5
      ENDIF
      GOTO 300
C     =================================================================
C     CUBP
C     =================================================================
 81   CONTINUE
      INFELE(2)=27
      INFELE(3)=27
      INFELE(4)=27
      INFELE(6)=27
      INFELE(8)=28
      INFELE(9)=68
      INFELE(10)=6
      INFELE(15)=4
      INFELE(16)=7
      GOTO 300
C     =================================================================
C     TETP
C     =================================================================
 82   CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(8)=14
      INFELE(9)=34
      INFELE(10)=6
      INFELE(15)=4
      INFELE(16)=7
      GOTO 300
C     =================================================================
C     PRIP
C     =================================================================
 83   CONTINUE
      INFELE(2)=12
      INFELE(3)=12
      INFELE(4)=12
      INFELE(6)=12
      INFELE(8)=21
      INFELE(9)=51
      INFELE(10)=6
      INFELE(15)=4
      INFELE(16)=7
      GOTO 300
C     =================================================================
C     TIMO : Element POUTRE de TIMOSHENCKO
C     =================================================================
 84   CONTINUE
      INFELE(2)=1
      INFELE(4)=1
      INFELE(6)=1
      INFELE(7)=200
      IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
C  Element TIMO 2D plan
        INFELE(5)=6
        INFELE(9)=6
        INFELE(10)=3
        INFELE(15)=3
        INFELE(16)=3
      ELSE
C  Element TIMO 3D
        INFELE(5)=12
        INFELE(9)=12
        INFELE(10)=6
        INFELE(15)=6
        INFELE(16)=6
      ENDIF
      GOTO 300
C     =================================================================
C     JOI2
C     =================================================================
 85   CONTINUE
      INFELE(2)=4
      INFELE(3)=2
      INFELE(4)=2
      INFELE(5)=0
      INFELE(6)=2
      INFELE(7)=300
      INFELE(9)=8
      INFELE(10)=2
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     JOI3
C     =================================================================
 86   CONTINUE
      INFELE(2)=3
      INFELE(3)=2
      INFELE(4)=3
      INFELE(5)=0
      INFELE(6)=3
      INFELE(7)=100
      INFELE(9)=12
      INFELE(10)=2
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     JOT3
C     =================================================================
 87   CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(5)=0
      INFELE(6)=3
      INFELE(7)=100
      INFELE(9)=18
      INFELE(10)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     JOI4
C     =================================================================
 88   CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(5)=0
      INFELE(6)=4
      INFELE(7)=100
      INFELE(9)=24
      INFELE(10)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     JOI6
C     =================================================================
 89   CONTINUE
      INFELE(2)=5
      INFELE(4)=5
      INFELE(6)=4
      INFELE(9)=36
      INFELE(10)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     JOI8
C     =================================================================
 90   CONTINUE
      INFELE(2)=9
      INFELE(4)=9
      INFELE(6)=9
      INFELE(9)=48
      INFELE(10)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     LINE SPRING LISC
C     =================================================================
 91   CONTINUE
      INFELE(2)=3
      INFELE(4)=3
      INFELE(5)=5
      INFELE(6)=5
      INFELE(7)=100
      INFELE(9)=36
      INFELE(10)=2
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     TRIH
C     =================================================================
 92   CONTINUE
      INFELE(2)=4
      INFELE(6)=4
      INFELE(7)= 0
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN
        INFELE(3)=7
        INFELE(5)=5
        IF (IFOUR.EQ.0) INFELE(9)=12
        IF (IFOUR.EQ.1) INFELE(9)=18
        INFELE(15)=6
      ELSE
        INFELE(3)=4
        INFELE(5)=5
        INFELE(9)=12
        INFELE(15)=4
      ENDIF
      INFELE(10)=10
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     DST
C     =================================================================
 93   CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(5)=3
      INFELE(6)=3
      INFELE(7)=600
      INFELE(9)=18
      INFELE(10)=8
      INFELE(15)=6
      INFELE(16)=8
      GOTO 300
C     =================================================================
C     LIC4 : Element de RACCORD LIQUIDE-COQUE a 8 noeuds en 3D
C     =================================================================
 94   CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(5)=3
      INFELE(6)=4
      INFELE(9)=32
      INFELE(10)=6
      INFELE(15)=6
      GOTO 300
C     =================================================================
C     CERCE
C     =================================================================
 95   CONTINUE
      INFELE(2)=1
      INFELE(3)=1
      INFELE(4)=1
      INFELE(5)=1
      INFELE(6)=1
      INFELE(7)=50
      INFELE(9)=2
      IF (IFOUR.EQ.1) INFELE(9)=3
      INFELE(10)=1
      INFELE(15)=2
      IF (IFOUR.EQ.1) INFELE(15)=3
      INFELE(16)=1
      GOTO 300
C     =================================================================
C     TUYO
C     =================================================================
 96   CONTINUE
      NTETA=40
      INFELE(2)=2*NTETA
      INFELE(3)=2*NTETA
      INFELE(4)=2*NTETA
      INFELE(5)=7
      INFELE(6)=2*NTETA
      INFELE(7)=500
      INFELE(8)=2+4
      INFELE(9)=12+2*9
      INFELE(10)=6
      INFELE(15)=6+9
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     LSE2 : TUYAU ACOUSTIQUE PURE
C     =================================================================
 97   CONTINUE
      INFELE(2)=2
      INFELE(3)=2
      INFELE(4)=2
      INFELE(5)=7
      INFELE(6)=2
      INFELE(7)=200
      INFELE(9)=4
      INFELE(10)=2
      INFELE(15)=2
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     LITU : COUPLAGE LIQUIDE TUYAU
C     =================================================================
 98   CONTINUE
      INFELE(2)=2
      INFELE(3)=2
      INFELE(4)=2
      INFELE(5)=6
      INFELE(6)=2
      INFELE(7)=200
      INFELE(9)=10
      INFELE(10)=2
      INFELE(15)=5
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     HYT3 : Element HYBRIDE TRI3
C     =================================================================
 99   CONTINUE
      INFELE(3)=4
      INFELE(6)=4
      GOTO 300
C     =================================================================
C     HYQ4 : Element HYBRIDE QUA4
C     =================================================================
 100  CONTINUE
      INFELE(3)=4
      INFELE(6)=4
      GOTO 300
C     =================================================================
C     HYT4 : Element HYBRIDE TET4
C     =================================================================
 101  CONTINUE
      INFELE(3)=4
      INFELE(6)=4
      GOTO 300
C     =================================================================
C     HYP6 : Element HYBRIDE PRI6
C     =================================================================
 102  CONTINUE
      INFELE(3)=6
      INFELE(6)=6
      GOTO 300
C     =================================================================
C     HYC8 : Element HYBRIDE CUB8
C     =================================================================
 103  CONTINUE
      INFELE(3)=8
      INFELE(6)=8
      GOTO 300
C     =================================================================
C     TRIS (TRI3-SECTION)
C     =================================================================
 104  CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(5)=2
      INFELE(6)=4
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
        INFELE(10)=2
        INFELE(16)=2
      ELSE IF (IFOUR.EQ.2) THEN
        INFELE(10)=3
        INFELE(16)=3
      ENDIF
      GOTO 300
C     =================================================================
C     QUAS (QUA4-SECTION)
C     =================================================================
 105  CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(5)=2
      INFELE(6)=4
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
        INFELE(10)=2
        INFELE(16)=2
      ELSE IF (IFOUR.EQ.2) THEN
        INFELE(10)=3
        INFELE(16)=3
      ENDIF
      GOTO 300
C     =================================================================
C     POIS (POI1-SECTION)
C     =================================================================
 106  CONTINUE
      INFELE(2)=1
      INFELE(3)=1
      INFELE(4)=1
      INFELE(5)=1
      INFELE(6)=1
      INFELE(10)=1
      INFELE(16)=1
      GOTO 300
C     =================================================================
C     JOP3 : JOINTS POREUX
C     =================================================================
 108  CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(7)=100
      INFELE(8)=12
      INFELE(9)=18
      INFELE(10)=2
      INFELE(15)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     JOP6 : JOINTS POREUX
C     =================================================================
 109  CONTINUE
      INFELE(2)=7
      INFELE(3)=7
      INFELE(4)=7
      INFELE(6)=7
      INFELE(7)=100
      INFELE(8)=21
      INFELE(9)=45
      INFELE(10)=3
      INFELE(15)=4
      INFELE(16)=4
      GOTO 300
C     =================================================================
C     JOP8 : JOINTS POREUX
C     =================================================================
 110  CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      INFELE(4)=9
      INFELE(6)=9
      INFELE(7)=100
      INFELE(8)=28
      INFELE(9)=60
      INFELE(10)=3
      INFELE(15)=4
      INFELE(16)=4
      GOTO 300
C     =================================================================
C     POLYGONE a N cotes
C     =================================================================
 111  CONTINUE
      INFELE(8) = MELE - 108
      NBSH=INFELE(8)
      INFELE(2)=NBSH
      INFELE(3)=NBSH
      INFELE(4)=NBSH
      INFELE(6)=NBSH
      INFELE(9)=2*NBSH
      IF (IFOUR.EQ.1)  INFELE(9)=3*NBSH
      IF (IFOUR.EQ.-3) INFELE(9)=2*NBSH+3
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     BAR3 (barre a 3 noeuds)
C     =================================================================
 123  CONTINUE
      NPOINT=3
      INFELE(2)=NPOINT
      INFELE(4)=NPOINT
      INFELE(5)=1
      INFELE(6)=NPOINT
      INFELE(7)=200
      INFELE(9)=9
      IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) INFELE(9)=6
      INFELE(10)=1
      INFELE(15)=3
      IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) INFELE(15)=2
      INFELE(16)=1
      GOTO 300
C     =================================================================
C     BAEX (barre excentree 3D a 2 noeuds)
C     =================================================================
 124  CONTINUE
      INFELE(2)=2
      INFELE(3)=2
      INFELE(4)=2
      INFELE(5)=6
      INFELE(6)=2
      INFELE(7)=700
      INFELE(9)=12
      INFELE(10)=1
      INFELE(15)=6
      INFELE(16)=1
      GOTO 300
C     =================================================================
C     LIA2 : element de liaison a 2 noeuds (6 ddl par noeuds)
C     =================================================================
 125  CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(5)=9
      INFELE(6)=2
      INFELE(7)=200
      INFELE(9)=12
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     QUAH
C     =================================================================
 126  CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(5)=5
      INFELE(6)=4
      INFELE(7)=0
      IF (IFOUR.EQ.0) INFELE(9)=16
      IF (IFOUR.EQ.1) INFELE(9)=24
      INFELE(10)=10
      INFELE(15)=6
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     CUBH
C     =================================================================
 127  CONTINUE
      INFELE(2)=8
      INFELE(3)=27
      INFELE(5)=5
      INFELE(6)=27
      INFELE(7)=0
      INFELE(9)=48
      INFELE(10)=10
      INFELE(15)=6
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     ROT3 (1 pt GAUSS pr MASSE, 7 pts GAUSS pr RIGIDITE)
C     =================================================================
 128  CONTINUE
      INFELE(3)=1
      INFELE(4)=INFELE(3)
      INFELE(6)=7
      INFELE(9)=3
      INFELE(15)=1
      GOTO 300
C     =================================================================
C     TR6H
C     =================================================================
 157  CONTINUE
      INFELE(2)=4
      INFELE(6)=6
      INFELE(7)=0
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN
        INFELE(3)=7
        INFELE(5)=5
        IF (IFOUR.EQ.0) INFELE(9)=12
        IF (IFOUR.EQ.1) INFELE(9)=18
        INFELE(15)=6
      ELSE
        INFELE(3)=6
        INFELE(5)=5
        INFELE(9)=24
        INFELE(15)=4
      ENDIF
      INFELE(10)=10
      INFELE(16)=0
      GOTO 300
C     =================================================================
C     SEGS (SEG2-SECTION 2D)
C     =================================================================
 166  CONTINUE
      INFELE(2)=2
      INFELE(3)=2
      INFELE(4)=2
      INFELE(5)=3
      INFELE(6)=2
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
        INFELE(10)=2
        INFELE(16)=2
      ELSE IF (IFOUR.EQ.2) THEN
        INFELE(10)=3
        INFELE(16)=3
      ENDIF
      GOTO 300
C     =================================================================
C     POJS (POI1-SECTION)
C     =================================================================
 167  CONTINUE
      INFELE(2)=1
      INFELE(3)=1
      INFELE(4)=1
      INFELE(5)=3
      INFELE(6)=1
      IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
        INFELE(10)=2
        INFELE(16)=2
      ELSE IF (IFOUR.EQ.2) THEN
        INFELE(10)=3
        INFELE(16)=3
      ENDIF
      GOTO 300
C     =================================================================
C     JCT3
C     =================================================================
 168  CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(7)=300
      INFELE(9)=INFELE(8)*IDIM
      INFELE(10)=2
      INFELE(15)=IDIM
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     JCI4
C     =================================================================
 169  CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(7)=300
      INFELE(9)=INFELE(8)*IDIM
      INFELE(10)=2
      INFELE(15)=IDIM
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     JGI2
C     =================================================================
 170  CONTINUE
      INFELE(2)=4
      INFELE(3)=2
      INFELE(4)=2
      INFELE(6)=2
      IF (IFOUR.EQ.-3) INFELE(9)=INFELE(9)+3
      INFELE(7)=300
      INFELE(9)=INFELE(8)*IDIM
      INFELE(10)=4
      INFELE(15)=IDIM
      INFELE(16)=4
      GOTO 300
C     =================================================================
C     JGT3
C     =================================================================
 171  CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(7)=300
      INFELE(9)=INFELE(8)*IDIM
      INFELE(10)=6
      INFELE(15)=IDIM
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     JGI4
C     =================================================================
 172  CONTINUE
      INFELE(2)=4
      INFELE(3)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(7)=300
      INFELE(9)=INFELE(8)*IDIM
      INFELE(10)=6
      INFELE(15)=IDIM
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     TRIQ
C     =================================================================
 173  CONTINUE
      INFELE(3)=7
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=7
        INFELE(4)=7
        INFELE(6)=7
      ELSE
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
      ENDIF
      INFELE(8)=9
      INFELE(9)=18
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=24
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=8
        INFELE(15)=5
        INFELE(16)=8
      ELSE
        INFELE(10)=6
        INFELE(15)=4
        INFELE(16)=6
      ENDIF
      GOTO 300
C     =================================================================
C     QUAQ
C     =================================================================
 174  CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=9
      INFELE(6)=9
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(8)=12
      INFELE(9)=24
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=32
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=8
        INFELE(15)=5
        INFELE(16)=8
      ELSE
        INFELE(10)=6
        INFELE(15)=4
        INFELE(16)=6
      ENDIF
      GOTO 300
C     =================================================================
C     CUBQ
C     =================================================================
 175  CONTINUE
      INFELE(2)=27
      INFELE(3)=27
      INFELE(4)=27
      INFELE(6)=27
      INFELE(8)=28
      INFELE(9)=76
      INFELE(10)=8
      INFELE(15)=5
      INFELE(16)=8
      GOTO 300
C     =================================================================
C     TETQ
C     =================================================================
 176  CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(8)=14
      INFELE(9)=38
      INFELE(10)=8
      INFELE(15)=5
      INFELE(16)=8
      GOTO 300
C     =================================================================
C     PRIQ
C     =================================================================
 177  CONTINUE
      INFELE(2)=12
      INFELE(3)=12
      INFELE(4)=12
      INFELE(6)=12
      INFELE(8)=21
      INFELE(9)=57
      INFELE(10)=8
      INFELE(15)=5
      INFELE(16)=8
      GOTO 300
C     =================================================================
C     TRIR
C     =================================================================
 178  CONTINUE
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
        INFELE(2)=7
        INFELE(6)=7
        INFELE(4)=7
      ELSE
        INFELE(2)=4
        INFELE(4)=4
        INFELE(6)=4
      ENDIF
      INFELE(3)=7
      INFELE(8)=9
      INFELE(9)=21
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=27
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=9
        INFELE(15)=6
        INFELE(16)=9
      ELSE
        INFELE(10)=7
        INFELE(15)=5
        INFELE(16)=7
      ENDIF
      GOTO 300
C     =================================================================
C     QUAR
C     =================================================================
 179  CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=9
      INFELE(6)=9
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(8)=12
      INFELE(9)=28
      IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=36
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=9
        INFELE(15)=6
        INFELE(16)=9
      ELSE
        INFELE(10)=7
        INFELE(15)=5
        INFELE(16)=7
      ENDIF
      GOTO 300
C     =================================================================
C     CUBR
C     =================================================================
 180  CONTINUE
      INFELE(2)=27
      INFELE(3)=27
      INFELE(4)=27
      INFELE(6)=27
      INFELE(8)=28
      INFELE(9)=84
      INFELE(10)=9
      INFELE(15)=6
      INFELE(16)=9
      GOTO 300
C     =================================================================
C     TETR
C     =================================================================
 181  CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(8)=14
      INFELE(9)=42
      INFELE(10)=9
      INFELE(15)=6
      INFELE(16)=9
      GOTO 300
C     =================================================================
C     PRIR
C     =================================================================
 182  CONTINUE
      INFELE(2)=12
      INFELE(3)=12
      INFELE(4)=12
      INFELE(6)=12
      INFELE(8)=21
      INFELE(9)=61
      INFELE(10)=9
      INFELE(15)=6
      INFELE(16)=9
      GOTO 300
C     =================================================================
C     Q4RI (QUA4 with 1x1 Gauss points)
C     =================================================================
 183  CONTINUE
      INFELE(2)=4
      INFELE(3)=1
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9
      INFELE(4)=1
C**   INFELE(4)=5    CORRECTION  MILL  LE 8/8/90
      INFELE(6)=1
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9
      INFELE(9)=8
      IF (IFOUR.EQ.-3) INFELE(9)=11
      IF (IFOUR.EQ.1)  INFELE(9)=12
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     Q8RI (QUA8 with 2x2 Gauss points)
C     =================================================================
 184  CONTINUE
      INFELE(2)=9
      INFELE(3)=4
      IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16
      INFELE(4)=4
      INFELE(6)=4
C**   IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16
      INFELE(9)=16
      IF (IFOUR.EQ.-3) INFELE(9)=19
      IF (IFOUR.EQ.1)  INFELE(9)=24
      IF (IFOUR.EQ.1) THEN
        INFELE(10)=6
        INFELE(15)=3
        INFELE(16)=6
      ELSE
        INFELE(10)=4
        INFELE(15)=2
        INFELE(16)=4
      ENDIF
      GOTO 300
C     =================================================================
C     JOQ3
C     =================================================================
 185  CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(7)=150
      INFELE(8)=12
      INFELE(9)=24
      INFELE(10)=4
      INFELE(15)=4
      INFELE(16)=4
      GOTO 300
C     =================================================================
C     JOQ6
C     =================================================================
 186  CONTINUE
      INFELE(2)=7
      INFELE(3)=7
      INFELE(4)=7
      INFELE(6)=7
      INFELE(7)=150
      INFELE(8)=21
      INFELE(9)=54
      INFELE(10)=5
      INFELE(15)=5
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     JOQ8
C     =================================================================
 187  CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      INFELE(4)=9
      INFELE(6)=9
      INFELE(7)=150
      INFELE(8)=28
      INFELE(9)=72
      INFELE(10)=5
      INFELE(15)=5
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     JOR3
C     =================================================================
 188  CONTINUE
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(7)=200
      INFELE(8)=12
      INFELE(9)=30
      INFELE(10)=5
      INFELE(15)=5
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     JOR6
C     =================================================================
 189  CONTINUE
      INFELE(2)=7
      INFELE(3)=7
      INFELE(4)=7
      INFELE(6)=7
      INFELE(7)=200
      INFELE(8)=21
      INFELE(9)=63
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     JOR8
C     =================================================================
 190  CONTINUE
      INFELE(2)=9
      INFELE(3)=9
      INFELE(4)=9
      INFELE(6)=9
      INFELE(7)=200
      INFELE(8)=28
      INFELE(9)=84
      INFELE(10)=6
      INFELE(15)=6
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     T1D2/M1D2 : Element THERMIQUE/MECANIQUE MASSIF 1D a 2 noeuds
C     =================================================================
 191  CONTINUE
C  On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1)
      INFELE(2)=2
      INFELE(3)=2
      INFELE(4)=2
      INFELE(5)=3
      INFELE(6)=2
      INFELE(7)=300
      INFELE(9)=2
      IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) INFELE(9)=3
      IF (IFOUR.EQ.11) INFELE(9)=4
      INFELE(10)=3
      INFELE(15)=1
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     T1D3/M1D3 : Element THERMIQUE/MECANIQUE MASSIF 1D a 3 noeuds
C     =================================================================
 192  CONTINUE
C  On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1)
      INFELE(2)=3
      INFELE(3)=3
      INFELE(4)=3
      INFELE(6)=3
      INFELE(9)=3
      IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) INFELE(9)=4
      IF (IFOUR.EQ.11) INFELE(9)=5
        INFELE(10)=3
        INFELE(15)=1
        INFELE(16)=3
      GOTO 300
C     =================================================================
C     LC03 LS03 BS03
C     =================================================================
 195  CONTINUE
      INFELE(6)=2
      GOTO 300
C     =================================================================
C     LC07 LS07 BS07
C     =================================================================
 196  CONTINUE
      INFELE(6)=7
      GOTO 300
C     =================================================================
C     LC09 LS09 BS09
C     =================================================================
 197  CONTINUE
      INFELE(6)=4
      GOTO 300
C     =================================================================
C     MC03
C     =================================================================
 216  CONTINUE
      INFELE(6)=2
      GOTO 300
C     =================================================================
C     MC07
C     =================================================================
 217  CONTINUE
      INFELE(6)=7
      GOTO 300
C     =================================================================
C     MC09
C     =================================================================
 218  CONTINUE
      INFELE(6)=4
      GOTO 300
C     =================================================================
C     M103 MS03
C     =================================================================
 223  CONTINUE
      INFELE(6)=3
      GOTO 300
C     =================================================================
C     M107 MS07
C     =================================================================
 224  CONTINUE
      INFELE(6)=7
      GOTO 300
C     =================================================================
C     M109 MS09
C     =================================================================
 225  CONTINUE
      INFELE(6)=16
      GOTO 300
C     =================================================================
C     QC03 Q103 QS03
C     =================================================================
 237  CONTINUE
      INFELE(6)=3
      GOTO 300
C     =================================================================
C     QC07 Q107 QS07
C     =================================================================
 238  CONTINUE
      INFELE(6)=7
      GOTO 300
C     =================================================================
C     QC09 Q109 QS09
C     =================================================================
 239  CONTINUE
      INFELE(6)=16
      GOTO 300
C     =================================================================
C     CIFL
C     =================================================================
 258  CONTINUE
      INFELE(2)=1
      INFELE(3)=1
      INFELE(4)=1
      INFELE(6)=1
      INFELE(8)=0
      INFELE(9)=8
      INFELE(10)=5
      INFELE(15)=3
      INFELE(16)=5
      GOTO 300
C     =================================================================
C     SHB8
C     =================================================================
 260  CONTINUE
      INFELE(2)=5
      INFELE(3)=5
      INFELE(4)=5
      INFELE(6)=5
      INFELE(9)=24
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     XQ4R
C     =================================================================
 263  CONTINUE
      INFELE(2)=0
      INFELE(3)=64
      INFELE(4)=64
      INFELE(6)=64
      INFELE(8)=40
      INFELE(10)=4
      INFELE(15)=20
      INFELE(16)=4
      INFELE(9)=4*INFELE(15)
      GOTO 300
C     =================================================================
C     XC8R
C     =================================================================
 264  CONTINUE
      INFELE(2)=0
      INFELE(3)=64
      INFELE(4)=64
      INFELE(6)=64
      INFELE(8)=80
      INFELE(10)=6
      INFELE(15)=30
      INFELE(16)=6
      INFELE(9)=8*INFELE(15)
      GOTO 300
C     =================================================================
C     JOI1 : element de liaison a 2 noeuds (6 ddl par noeuds)
C     =================================================================
 265  CONTINUE
      INFELE(2)=1
      INFELE(3)=1
      INFELE(4)=1
      INFELE(5)=0
      INFELE(6)=1
      INFELE(7)=700
      IF(IFOUR.EQ.2) THEN
        INFELE(9)=12
        INFELE(10)=6
        INFELE(15)=6
        INFELE(16)=6
      ELSE IF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
        INFELE(9)=6
        INFELE(10)=3
        INFELE(15)=3
        INFELE(16)=3
      ENDIF
      GOTO 300
C     =================================================================
C     ZCO2 : zone cohesive a 2 noeuds (2 ddls par noeud)
C     =================================================================
 266  CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(6)=2
      INFELE(7)=100
      INFELE(9)=4
      INFELE(10)=2
      INFELE(15)=2
      INFELE(16)=2
      GOTO 300
C     =================================================================
C     ZCO3 : zone cohesive a 3 noeuds (2 ddls par noeud)
C     =================================================================
 267  CONTINUE
      INFELE(2)=1
      INFELE(4)=1
      INFELE(6)=1
      INFELE(7)=100
      INFELE(9)=9
      INFELE(10)=3
      INFELE(15)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     ZCO4 : zone cohesive a 4 noeuds (2 ddls par noeud)
C     =================================================================
 268  CONTINUE
      INFELE(2)=4
      INFELE(4)=4
      INFELE(6)=4
      INFELE(7)=100
      INFELE(9)=12
      INFELE(10)=3
      INFELE(15)=3
      INFELE(16)=3
      GOTO 300
C     =================================================================
C     TUY2 : tuyau pour modele d'advection en thermique
C     =================================================================
 269  CONTINUE
      INFELE(2)=2
      INFELE(4)=2
      INFELE(5)=5
      INFELE(6)=2
      INFELE(7)=200
      INFELE(9)=2
      INFELE(10)=1
      INFELE(15)=1
      INFELE(16)=1
      GOTO 300
C     =================================================================
C     TUY3 : tuyau pour modele d'advection en thermique
C     =================================================================
 270  CONTINUE
      INFELE(2)=3
      INFELE(4)=3
      INFELE(5)=5
      INFELE(6)=3
      INFELE(7)=200
      INFELE(9)=2
      INFELE(10)=1
      INFELE(15)=1
      INFELE(16)=1
      GOTO 300
C     =================================================================
c     COS2 (3D pour liaison acier-beton)
C     =================================================================
 271  CONTINUE
      INFELE(2)=1
      INFELE(4)=2
      INFELE(5)=2
      INFELE(6)=2
      INFELE(7)=100
      INFELE(9)=4*idim
      INFELE(10)=3
      INFELE(15)=IDIM
      INFELE(16)=3
      GOTO 300
C     =================================================================
c     COA2 (3D pour liaison acier-beton)
C     =================================================================
 272  CONTINUE
      INFELE(2)=1
      INFELE(4)=2
      INFELE(5)=2
      INFELE(6)=2
      INFELE(7)=100
      INFELE(9)=4*idim
      INFELE(10)=3
      INFELE(15)=IDIM
      INFELE(16)=3
      GOTO 300
C     =================================================================
c     CU27, comme le CU20
C     =================================================================
 275  CONTINUE
      INFELE(2)=27
      INFELE(3)=27
      INFELE(4)=27
      INFELE(6)=27
      INFELE(9)=81
      GOTO 300
C     =================================================================
c     PR21
C     =================================================================
 276  CONTINUE
      INFELE(2)=21
cbp      INFELE(3)=12
      INFELE(3)=21
      INFELE(4)=21
      INFELE(6)=21
      INFELE(9)=63
      GOTO 300
C     =================================================================
c     TE15
C     =================================================================
 277  CONTINUE
      INFELE(2)=15
cbp      INFELE(3)=5
      INFELE(3)=15
      INFELE(4)=15
      INFELE(6)=15
      INFELE(9)=45
      GOTO 300
C     =================================================================
c     PY19
C     =================================================================
 278  CONTINUE
      INFELE(2)=17
      INFELE(3)=17
      INFELE(4)=17
      INFELE(6)=17
      INFELE(9)=57
      GOTO 300
C     =================================================================
C     C20R
C     =================================================================
 279  CONTINUE
      INFELE(2)=27
      INFELE(3)=8
      INFELE(4)=8
      INFELE(6)=8
      INFELE(9)=60
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
C     P15R
C     =================================================================
 280  CONTINUE
      INFELE(2)=12
      INFELE(3)=8
      INFELE(4)=6
      INFELE(6)=6
      INFELE(9)=45
      INFELE(10)=6
      INFELE(15)=3
      INFELE(16)=6
      GOTO 300
C     =================================================================
 300  CONTINUE
C
      IF (INFELE(3).EQ.0) INFELE(3)=INFELE(6)
C
C     Cas particuliers des autres formulations massives
      IF (FORMOD(1).EQ.'ELECTROSTATIQUE') THEN
        INFELE(9) = NBSH
        INFELE(10)=IDIM
        INFELE(15)=1
        INFELE(16)=IDIM
      ENDIF
 
* nombre composante contraintes/defo, materiau IMPEDANCE
      if (dcmate) then
         MFR3=INFELE(13)
         call idprim(IMODEL,MFR3,IPNOMC,NBROBL,NBRFAC)
        if(tymode(/2).gt.0) then
         if (tymode(1).eq.'LISTMOTS') then
          dcmat2 = .true.
          mlmot5 = ivamod(1)
          segact mlmot5
          nbrobl = mlmot5.mots(/2)*2
         endif
        endif
        if (dcmat2.and.mele.eq.2) then
          INFELE(9) = nbrobl
          INFELE(10) = nbrobl/2
          INFELE(15) = nbrobl/2
         INFELE(16) = nbrobl/2
        else
         INFELE(15) = nbrobl
          if (nbrfac.gt.0) INFELE(15) = INFELE(15) + nbrfac
         INFELE(16) = INFELE(15)
         if (mele.eq.45) INFELE(9) = nbrobl
         if (mele.eq.2) INFELE(9) = nbrobl*2
        endif
      endif
CJK148537 en attendant les fonctions de forme pr18, pr14, py19
      if (MELE.EQ.140.OR.MELE.EQ.142.OR.MELE.EQ.149) GOTO 950
C
      INTTYP = -5
      CALL ELQUOI(IMODEL,INTTYP,IPTR)
      INFELE(11) = IPTR
C
      if (bmodi) imodel.nefmod=melep
C
 950  CONTINUE
      END