Sources Esope | Cast3M

nummfr
C NUMMFR    SOURCE    OF166741  24/05/06    21:15:24     11082          
 
C---------------------------------------------------------------------
C                             N U M M F R
C                             -----------
C
C  Fonction :
C ------------
C  Donner le numero de FORMULATION associee a un element FINI
C
C  Parametres :
C --------------
C  NUMEF   (E)   Numero de l'element FINI (place dans NOMTP)
C                voir bdata.eso et CCHAMP.INC
C  NUMMFR  (S)   Numero de la FORMULATION associee
C
C  Constantes :
C --------------
C  NUMFR : On trouve a la N-ieme position du tableau le numero de la
C          la FORMULATION pour le N-ieme element FINI
C
C-----------------------------------------------------------------------
C  FORMULATION de l element fini    1 = MASSIF
C                                   3 = COQUE
C                                   5 = COQUE EPAISSE
C                                   7 = POUTRE
C                                   9 = CISATRAN
C                                  11 = LIQUIDE
C                                  13 = TUYAU
C                                  15 = LINESPRING
C                                  17 = TUYAU FISSURE
C                                  19 = RACCORD MASSIF
C                                  21 = RACCORD COQUE
C                                  23 = SURFACE LIBRE
C                                  25 = MEMBRANE
C                                  27 = UNIAXIALE
C                                  29 = THERMIQUE
C                                  31 = INCOMPRESSIBLES
C                                  33 = POREUX
C                                  35 = JOINT
C                                  37 = HOMOGENEISE
C                                  39 = TUYO
C                                  41 = TUYAU ACOUSTIQUE PURE
C                                  43 = RACCORD TUYAU FLUIDE
C                                  45 = HYBRIDE
C                                  47 = FIBRE
C                                  49 = BARRE 3D EXCENTREE
C                                  51 = LIASON A 2 NOEUDS
C                                  52 = NAVIER_STOKES
C                                  53 = JOINT CISAILLEMENT 2D
C                                  55 = JOINT GENERALISE
C                                  57 = POREUX TYPE Q
C                                  59 = POREUX TYPE R
C                                  61 = MACRO-ELEMENT
C                                  63 = XFEM meca rupt plan ou 3D
C                                  65 = THERMOHYDRIQUE
C                                  69 = MAGNETODYNAMIQUE
C                                  71 = ELECTROSTATIQUE (base MASSIF=1)
C                                  73 = DIFFUSION       (base MASSIF=1)
C                                  75 = JOINT UNIDIMENSIONNEL
C                                  77 = ZONE COHESIVE
c                                  78 = element coaxial COS2
C                                  79 = TUYAU pour l'advection TUY2 et TUY3
C                                  89 = HHO (base MASSIF=1)
C-----------------------------------------------------------------------
 
      FUNCTION NUMMFR(NUMEF)
 
C  LNUMFR DEVRAIT etre egale a LNOMTP definie dans bdata.eso
C  MAIS CE N'EST PAS LE CAS ..!!
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
c cccccc
      PARAMETER (LNUMFR=281)
c cccccc
      DIMENSION NUMFR(LNUMFR)
 
      DATA     NUMFR( 1),NUMFR( 2),NUMFR( 3),NUMFR( 4),NUMFR( 5)
     .        ,NUMFR( 6),NUMFR( 7),NUMFR( 8),NUMFR( 9),NUMFR(10)
     .        ,NUMFR(11),NUMFR(12),NUMFR(13),NUMFR(14),NUMFR(15)
     .        ,NUMFR(16),NUMFR(17),NUMFR(18),NUMFR(19),NUMFR(20)
     .        ,NUMFR(21),NUMFR(22),NUMFR(23),NUMFR(24),NUMFR(25)
     .        ,NUMFR(26),NUMFR(27),NUMFR(28),NUMFR(29),NUMFR(30)
     .        ,NUMFR(31),NUMFR(32),NUMFR(33),NUMFR(34),NUMFR(35)
     .        ,NUMFR(36),NUMFR(37),NUMFR(38),NUMFR(39),NUMFR(40)
     .        ,NUMFR(41),NUMFR(42),NUMFR(43),NUMFR(44),NUMFR(45)
     .        ,NUMFR(46),NUMFR(47),NUMFR(48),NUMFR(49),NUMFR(50)
     .        ,NUMFR(51),NUMFR(52),NUMFR(53),NUMFR(54),NUMFR(55)
C                 CABL SEG2 SEG3 TRI3 TRI4 TRI6 TRI7 QUA4 QUA5 QUA8 QUA9
     .           /  27,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1
C                 RAC2 RAC3 CUB8 CU20 PRI6 PR15 LIA3 LIA4 LIA6 LIA8 MULT
     .           ,  19,   0,   1,   1,   1,   1,  19,  19,   0,   0,   0
C                 TET4 TE10 PYR5 PY13 COQ3  DKT POUT LISP FAC3 FAC4 FAC6
     .           ,   1,   1,   1,   1,   3,   3,   7,  15,   1,   1,   1
C                 FAC8 LTR3 LQU4 LCU8 LPR6 LTE4 LPY5 COQ8 TUYA TUFI COQ2
     .           ,   1,  11,  11,  11,  11,  11,  11,   5,  13,  17,   3
C                 POI1 BARR RACO LSU2 COQ4 LISM COF3 RES2 LSU3 LSU4 LICO
     .           ,  27,  27,  21,  23,   9,  15,   3,   0,  23,  23, 21/
 
      DATA     NUMFR(56),NUMFR(57),NUMFR(58),NUMFR(59),NUMFR(60)
     .        ,NUMFR(61),NUMFR(62),NUMFR(63),NUMFR(64),NUMFR(65)
     .        ,NUMFR(66),NUMFR(67),NUMFR(68),NUMFR(69),NUMFR(70)
     .        ,NUMFR(71),NUMFR(72),NUMFR(73),NUMFR(74),NUMFR(75)
     .        ,NUMFR(76),NUMFR(77),NUMFR(78),NUMFR(79),NUMFR(80)
     .        ,NUMFR(81),NUMFR(82),NUMFR(83),NUMFR(84),NUMFR(85)
     .        ,NUMFR(86),NUMFR(87),NUMFR(88),NUMFR(89),NUMFR(90)
     .        ,NUMFR(91),NUMFR(92),NUMFR(93),NUMFR(94),NUMFR(95)
     .        ,NUMFR(96),NUMFR(97),NUMFR(98),NUMFR(99),NUMFR(100)
     .        ,NUMFR(101),NUMFR(102),NUMFR(103),NUMFR(104),NUMFR(105)
     .        ,NUMFR(106),NUMFR(107),NUMFR(108),NUMFR(109),NUMFR(110)
C                 COQ6 CVS2 CVS3 CVT3 CVT6 CVQ4 CVQ8 THP5 TH13 THP6 TH15
     .           /   5,  29,  29,  29,  29,  29,  29,  29,  29,  29,  29
C                 THC8 TH20 ICT3 ICQ4 ICT6 ICQ8 ICC8 ICT4 ICP6 IC20 IC10
     .           ,  29,  29,  31,  31,  31,  31,  31,  31,  31,  31,  31
C                 IC15 TRIP QUAP CUBP TETP PRIP TIMO JOI2 JOI3 JOT3 JOI4
     .           ,  31,  33,  33,  33,  33,  33,   7,  35,  35,  35,  35
C                 JOI6 JOI8 LISC TRIH  DST LIC4 CERC TUYO LSE2 LITU HYT3
     .           ,  35,  35,  15,  37,   9,  21,  27,  39,  41,  43,  45
C                 HYQ4 HYT4 HYP6 HYC8 TRIS QUAS POIS FRO3 JOP3 JOP6 JOP8
     .           ,  45,  45,  45,  45,  47,  47,  27,   0,  33,  33, 33/
 
      DATA     NUMFR(111),NUMFR(112),NUMFR(113),NUMFR(114),NUMFR(115)
     .        ,NUMFR(116),NUMFR(117),NUMFR(118),NUMFR(119),NUMFR(120)
     .        ,NUMFR(121),NUMFR(122),NUMFR(123),NUMFR(124),NUMFR(125)
     .        ,NUMFR(126),NUMFR(127),NUMFR(128),NUMFR(129),NUMFR(130)
     .        ,NUMFR(131),NUMFR(132),NUMFR(133),NUMFR(134),NUMFR(135)
     .        ,NUMFR(136),NUMFR(137),NUMFR(138),NUMFR(139),NUMFR(140)
     .        ,NUMFR(141),NUMFR(142),NUMFR(143),NUMFR(144),NUMFR(145)
     .        ,NUMFR(146),NUMFR(147),NUMFR(148),NUMFR(149),NUMFR(150)
     .        ,NUMFR(151),NUMFR(152),NUMFR(153),NUMFR(154),NUMFR(155)
     .        ,NUMFR(156),NUMFR(157),NUMFR(158),NUMFR(159),NUMFR(160)
     .        ,NUMFR(161),NUMFR(162),NUMFR(163),NUMFR(164),NUMFR(165)
C                 POL3 POL4 POL5 POL6 POL7 POL8 POL9 PO10 PO11 PO12 PO13
     .           /   1,   1,   1,   1,   1,   1,   1,   1,   1,   1,   1
C                 PO14 BAR3 BAEX LIA2 QUAH CUBH ROT3 SEF2 TRF3 QUF4 CUF8
     .           ,   1,  27,  49,  51,  37,  37,   3,  52,  52,  52,  52
C                 PRF6 TEF4 PYF5 MSE3 MTR6 MQU9 MC27 MP18 MT10 MP14 SEF3
     .           ,  52,  52,  52,  52,  52,  52,  52,  52,  52,  52,  52
C                 TRF7 QUF9 CF27 PF21 TF15 PF19 SEG6 TR21 QU36 C216 P126
     .           ,  52,  52,  52,  52,  52,  52,   1,   1,   1,   1,   1
C                 TE56 PY91 TRH6 BSE2 BTR4 BQU5 BCU9 BPR7 BTE5 BPY6 FRO4
     .           ,   1,   1,  37,  52,  52,  52,  52,  52,  52,  52,  1/
 
      DATA     NUMFR(166),NUMFR(167),NUMFR(168),NUMFR(169),NUMFR(170)
     .        ,NUMFR(171),NUMFR(172),NUMFR(173),NUMFR(174),NUMFR(175)
     .        ,NUMFR(176),NUMFR(177),NUMFR(178),NUMFR(179),NUMFR(180)
     .        ,NUMFR(181),NUMFR(182),NUMFR(183),NUMFR(184),NUMFR(185)
     .        ,NUMFR(186),NUMFR(187),NUMFR(188),NUMFR(189),NUMFR(190)
     .        ,NUMFR(191),NUMFR(192),NUMFR(193),NUMFR(194)
C                 SEGS POJS JCT3 JCI4 JGI2 JGT3 JGI4 TRIQ QUAQ CUBQ TETQ
     .           /  47,  47,  53,  53,  55,  55,  55,  57,  57,  57,  57
C                 PRIQ TRIR QUAR CUBR TETR PRIR Q4RI Q8RI JOQ3 JOQ6 JOQ8
     .           ,  57,  59,  59,  59,  59,  59,   1,  1,   57,  57,  57
C                 JOR3 JOR6 JOR8 T1D2 T1D3 M1D2 M1D3
     .           ,  59,  59,  59,   1,   1,   1,   1 /
 
      DATA     NUMFR(195)
     .        ,NUMFR(196),NUMFR(197),NUMFR(198),NUMFR(199),NUMFR(200)
     .        ,NUMFR(201),NUMFR(202),NUMFR(203),NUMFR(204),NUMFR(205)
     .        ,NUMFR(206),NUMFR(207),NUMFR(208),NUMFR(209),NUMFR(210)
     .        ,NUMFR(211),NUMFR(212),NUMFR(213),NUMFR(214),NUMFR(215)
     .        ,NUMFR(216),NUMFR(217),NUMFR(218),NUMFR(219),NUMFR(220)
     .        ,NUMFR(221),NUMFR(222),NUMFR(223),NUMFR(224),NUMFR(225)
     .        ,NUMFR(226),NUMFR(227),NUMFR(228),NUMFR(229),NUMFR(230)
     .        ,NUMFR(231),NUMFR(232),NUMFR(233),NUMFR(234),NUMFR(235)
     .        ,NUMFR(236)
C                 LC03 LC07 LC09 LC27 LC21 LC15 LC19
     .           /  52,  52,  52,  52,  52,  52,  52
C                 LS03 LS07 LS09 LS27 LS21 LS15 LS19
     .           ,  52,  52,  52,  52,  52,  52,  52
C                 BS03 BS07 BS09 BS27 BS21 BS15 BS19
     .           ,  52,  52,  52,  52,  52,  52,  52
C                 MC03 MC07 MC09 MC27 MC21 MC15 MC19
     .           ,  52,  52,  52,  52,  52,  52,  52
C                 M103 M107 M109 M127 M121 M115 M119
     .           ,  52,  52,  52,  52,  52,  52,  52
C                 MS03 MS07 MS09 MS27 MS21 MS15 MS19
     .           ,  52,  52,  52,  52,  52,  52,  52 /
 
      DATA     NUMFR(237),NUMFR(238),NUMFR(239),NUMFR(240),NUMFR(241)
     .        ,NUMFR(242),NUMFR(243),NUMFR(244),NUMFR(245),NUMFR(246)
     .        ,NUMFR(247),NUMFR(248),NUMFR(249),NUMFR(250),NUMFR(251)
     .        ,NUMFR(252),NUMFR(253),NUMFR(254),NUMFR(255),NUMFR(256)
     .        ,NUMFR(257),NUMFR(258),NUMFR(259),NUMFR(260),NUMFR(261)
     .        ,NUMFR(262),NUMFR(263),NUMFR(264),NUMFR(265),NUMFR(266)
c cccccc
     .        ,NUMFR(267),NUMFR(268),NUMFR(269),NUMFR(270),numfr(271)
     .        ,NUMFR(272),NUMFR(273),NUMFR(274),NUMFR(275),NUMFR(276)
     .        ,NUMFR(277),NUMFR(278),NUMFR(279),NUMFR(280)
c cccccc
C                 QC03 QC07 QC09 QC27 QC21 QC15 QC19
     .           /  52,  52,  52,  52,  52,  52,  52
C                 Q103 Q107 Q109 Q127 Q121 Q115 Q119
     .           ,  52,  52,  52,  52,  52,  52,  52
C                 QS03 QS07 QS09 QS27 QS21 QS15 QS19
     .           ,  52,  52,  52,  52,  52,  52,  52
C                MACRO SURE SHB8 CAF2 CAF3 XQ4R XC8R
     .           ,  61,   1,   1,   0,   0,  63,  63
c cccccc
C                JOI1  ZCO2 ZCO3 ZCO4 TUY2 TUY3 COS2
     .           ,  75,  77,  77,  77, 79,  79,   78
C                COA2  ICY5 IC13 CU27 PR21 TE15  PY19
     .           ,  78,  31,  31,   1,   1,   1,    1
C                C20R  P15R
     .           ,   1,   1/
c cccccc
 
C FORMULATION HHO (voir INCLUDE CCHHOPA) :
C NUMFR(HHO_NUM_ELEMENT) = HHO_MFR_ELEMENT
      DATA NUMFR(281) / 89 /
 
      IF ((NUMEF.GT.LNUMFR).OR.(NUMEF.LE.0)) THEN
            NUMMFR=0
      ELSE
         NUMMFR=NUMFR(NUMEF)
      ENDIF
 
C      RETURN
      END