Télécharger la plaquette Cast3M


fr en

Télécharger bdata.eso

Retour à la liste

Numérotation des lignes : 

bdata

  1. C BDATA     SOURCE    JK148537  23/11/20    21:15:03     11790          
  2.       BLOCK DATA
  3.  
  4. -INC PPARAM
  5. -INC CCOPTIO
  6. -INC CCPRECO
  7. -INC CCGEOME
  8. -INC CCREEL
  9. -INC CCNOYAU
  10. -INC CCHAMP
  11. -INC CCVQUA4
  12. -INC CCASSIS
  13. -INC CCFXDR
  14. -INC CCTRACE
  15. -INC CCPERF
  16.  
  17.       PARAMETER (UNTIER= 0.3333333333333333)
  18.       PARAMETER (UNSIXI= 0.1666666666666667)
  19.       PARAMETER (UNSIXM=-0.1666666666666667)
  20.       PARAMETER (UNDOUZ= 0.0833333333333333)
  21.  
  22.  
  23.       COMMON /CTDRAN/ r(32),xa,xm,xg,xval,xgerme,iflag,ibuff
  24.       REAL*8 xa,xm,xg,xval,xgerme,r
  25.  
  26. C  -----------------------
  27. C  Valeurs de PPARAM.INC :
  28. C  -----------------------
  29. C     TAG pour GREP : LOCHAI,LOCOMP
  30.       DATA          FMCHAI/'(A512)'/
  31.       DATA          FMCOMP/'(A8)'/
  32.  
  33. C  -----------------------
  34. C  Valeurs de CCPERF.INC :
  35. C  -----------------------
  36.       DATA ITPSPR,ITPSBO,IPTEMP /0,0,0/
  37.  
  38. C  ------------------------
  39. C  Valeurs de CCPRECO.INC :
  40. C  ------------------------
  41.  
  42. C  Preconditionnement numerotation
  43.       DATA PRENUM / npreco*0/
  44.  
  45. C  Preconditionnement COMP
  46.       DATA PREORI / 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
  47.      .              0,0,0,0,0,0,0,0,0,0 /
  48.       DATA PREPRE / 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
  49.      .              0,0,0,0,0,0,0,0,0,0 /
  50.       DATA PRECLE / ' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',
  51.      .              ' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',
  52.      .              ' ',' ',' ',' ',' ',' ',' ',' ',' ',' ' /
  53. C  Preconditionnement crechp
  54.       PARAMETER (nbemet= nbemel * 65       )
  55.       DATA PREMEL /nbemet    *0/
  56. *     DATA PREMEL / 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
  57. *    >              0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
  58. *    .              0,0,0,0,0,0,0,0,0,0 /
  59.  
  60. C  Preconditionnement part6
  61.       DATA PARMOD / NMOPAR*0 /
  62.       DATA PESCLA / NMOPAR*0 /
  63.  
  64.  
  65. C  Preconditionnement reduaf (Echelle globale MMODEL et MCHELM)
  66.       PARAMETER(ITAIL1=NBASMA+1, ITAIL2=ITAIL1*NPREDU)
  67.       DATA NBPRRE / ITAIL1*0 /
  68.       DATA PRECMO / ITAIL2*0 /
  69.       DATA PRECM1 / ITAIL2*0 /
  70.       DATA PRECM2 / ITAIL2*0 /
  71.       DATA PRECM3 / ITAIL2*0 /
  72.       DATA PRECM4 / ITAIL2*0 /
  73.       DATA PRECM5 / ITAIL2*0 /
  74.  
  75.  
  76. C  Preconditionnement extension de MMODEL (modete.eso pour l'instant)
  77.       DATA NBMOMO / ITAIL1*0 /
  78.       DATA PMOMO1 / ITAIL2*0 /
  79.       DATA PMOMO2 / ITAIL2*0 /
  80.  
  81.  
  82. C  Preconditionnement reduaf (Intersections des MELEME)
  83.       PARAMETER(ITAIL3=ITAIL1*NTRIPL)
  84.       DATA NINTSA / ITAIL1*0 /
  85.       DATA PMAMOD / ITAIL3*0 /
  86.       DATA PMAMOH / ITAIL3*0 /
  87.       DATA PMACHA / ITAIL3*0 /
  88.       DATA PMACHH / ITAIL3*0 /
  89.       DATA PMLENT / ITAIL3*0 /
  90.  
  91.  
  92. C  Preconditionnement actobj (pour faire les SEGACT par paquets)
  93.       DATA PACTOB / ITAIL1*0 /
  94.  
  95. C  Preconditionnement chame1 (pour se rappeler les chpts)
  96.       PARAMETER (ITAIL4=NPRCHA*(NBASMA+1))
  97.       DATA iprma / ITAIL4*0 /
  98.  
  99.  
  100. C  Nombres d'assistants et parametres assistants
  101.       DATA NBESC / 0 /
  102.       DATA LODESL,LOTRMA,LODEFE / .FALSE.,.FALSE.,.FALSE. /
  103.       DATA INILU,NOMLU / 0,0 /
  104.       DATA IEGAL / 1 /
  105.  
  106.  
  107. C  ------------------------
  108. C  Valeurs de CCOPTIO.INC :
  109. C  ------------------------
  110.  
  111.       DATA MOTERR/' '/
  112.       DATA TITREE/'CAST3M FECIT
  113.      .          ' /
  114.       DATA IERPER,IERMAX,IERR / 2,0,0 /
  115.       DATA IOLEC,IOIMP,IOCAR,IOPER,MEPSIL,IOSAU,IORES / 5,6,4,7,0,8,2 /
  116.       DATA IOACQ / 9 /
  117.       DATA IOSCR / 0 /
  118.       DATA IOTER / 5 /
  119.       DATA LUPARA /0/
  120.       DATA MISAUP /0/
  121.       DATA ZLOPRO /.FALSE./
  122. C  Pour IBM : a surcharger dans DEPMAC)
  123.       DATA IZROSF / 4 /
  124.       DATA IOSPI / 0 /
  125.       DATA IIMPI / 0 /
  126.       DATA IDIM / 0 /
  127.       DATA IFOMOD,IFOUR,NIFOUR,NSDPGE / -1,-1, 0, 0 /
  128.       DATA DIOCAD / 30. /
  129.  
  130.       DATA NOMSAU,NOMRES / '        ','        ' /
  131.  
  132. C ... Lorsque IONIVE est modifié il ne faut pas oublier de modifier la
  133. C     valeur maximale (en dur) dans OPTION  ...
  134.       DATA IONIVE / 24 /
  135.  
  136.       DATA NGMAXY/50000/
  137.       DATA ISOTYP/1/
  138.       DATA IOMBRE/1/
  139.       DATA NORVAL/0/
  140.       DATA NORVAD/0/
  141.       DATA NORINC/-1/
  142.       DATA NORIND/0/
  143.       DATA NUCROU/0/
  144.       DATA LANGUE/'FRAN' /
  145.       DATA IOEIL/0/
  146.       DATA Ilgni/0/
  147. C  Pour l'impression Postscript en vertical
  148.       DATA ZHORIZ /.TRUE./
  149.       DATA ZINIPS /.TRUE./
  150.  
  151. C  Dimension d'une file en nombre de flottant simple precision
  152.       DATA DIMFIC / 4.8E8 /
  153.  
  154. C  Epaisseur du trait pour le tracé
  155.       DATA IEPTR / 1 /
  156.  
  157. C  pour la surveillance et l'attente dans GEMAT
  158.       DATA MSURVE/0/
  159.       DATA MATTEN/0/
  160.  
  161. C  ------------------------
  162. C  Valeurs de CCGEOME.INC :
  163. C  -------------------------
  164.  
  165.       PARAMETER (NDIR1=426,NDIR1P=NDIR1+1,NDIR2=310)
  166.       PARAMETER (NDIR2P=NDIR1+NDIR2+1,NDIR3=764)
  167.       PARAMETER (NFAC1=371,NFAC1P=NFAC1+1,NFAC2=339)
  168.       PARAMETER (NFAC2P=NFAC1+NFAC2+1,NFAC3=290)
  169.       INTEGER KSEGM1(NDIR1),KSEGM2(NDIR2),KSEGM3(NDIR3)
  170.       INTEGER LFAC1(NFAC1),LFAC2(NFAC2),LFAC3(NFAC3)
  171.       EQUIVALENCE (KSEGM(1),KSEGM1(1)),(KSEGM(NDIR1P),KSEGM2(1)),
  172.      .            (KSEGM(NDIR2P),KSEGM3(1))
  173.       EQUIVALENCE (LFAC(1),LFAC1(1)),(LFAC(NFAC1P),LFAC2(1)),
  174.      .            (LFAC(NFAC2P),LFAC3(1))
  175.       INTEGER LISNO1(378),LISNO2(378),LISNO3(340)
  176.       INTEGER LIELM1(356),LIELM2(368)
  177.       INTEGER LPOS11(2,105),LPOS12(2,109)
  178.       EQUIVALENCE (LPOS1(1,1),LPOS11(1,1)),(LPOS1(1,106),LPOS12(1,1))
  179.       EQUIVALENCE (LISNOE(1),LISNO1(1)),(LISNOE(379),LISNO2(1)),
  180.      .            (LISNOE(757),LISNO3(1))
  181.       EQUIVALENCE (LIELM(1),LIELM1(1)),(LIELM(357),LIELM2(1))
  182.  
  183.       PARAMETER (NTABM1=248,NTABM4=217)
  184.       INTEGER ITABM1(NTABM1),ITABM2(NTABM1)
  185.       INTEGER ITABM3(NTABM1),ITABM4(NTABM4)
  186.       EQUIVALENCE (ITABM(0,0),ITABM1(1)),(ITABM(0,8),ITABM2(1))
  187.       EQUIVALENCE (ITABM(0,16),ITABM3(1)),(ITABM(0,24),ITABM4(1))
  188.  
  189. C  Couleur par defaut, indice du noir
  190. C  Tableau des couleurs (indices de 0 à 100)
  191.       DATA IDCOUL,IDNOIR,NBCOUL / 0,8,31 /
  192.       DATA NCOUL / 'DEFA','BLEU','ROUG','ROSE','VERT','TURQ','JAUN',
  193.      .             'BLAN','NOIR','VIOL','ORAN','AZUR','OCEA','CYAN',
  194.      .             'OLIV','GRIS','POUR','BRUN','BRIQ','CORA','BEIG',
  195.      .             'OR'  ,'MARI','BOUT','LIME','LAVA','BRON','KAKI',
  196.      .             'PEAU','CARA','INDI',70*'    ' /
  197.  
  198. C  Tableau de mélange des couleurs deux à deux [0:15,0:15]
  199. *      DATA ITABM / 0,5,3,2,6,1,0,0,
  200. *     .             5,1,3,3,5,5,4,5,
  201. *     .             3,3,2,3,6,3,6,3,
  202. *     .             2,3,3,3,2,3,6,2,
  203. *     .             6,5,6,2,4,5,6,6,
  204. *     .             1,5,3,4,5,5,4,1,
  205. *     .             0,4,6,6,6,4,6,7,
  206. *     .             0,5,3,2,6,1,7,7 /
  207. c       DATA ITABM / 0 ,1 ,2 ,3 ,4 ,5 ,6 ,7 ,0 ,9 ,10,11,12,13,14,15,
  208. c      .             1 ,1 ,11,11,14,5 ,4 ,7 ,1 ,1 ,6 ,11,12,13,14,1 ,
  209. c      .             2 ,11,2 ,3 ,6 ,11,10,7 ,2 ,11,2 ,11,11,11,14,2 ,
  210. c      .             3 ,11,3 ,3 ,2 ,3 ,6 ,7 ,3 ,11,2 ,9 ,13,2 ,15,3 ,
  211. c      .             4 ,14,6 ,2 ,4 ,5 ,6 ,7 ,4 ,1 ,1 ,7 ,13,7 ,2 ,4 ,
  212. c      .             5 ,5 ,3 ,4 ,5 ,5 ,4 ,7 ,5 ,4 ,12,1 ,13,5 ,15,5 ,
  213. c      .             6 ,4 ,6 ,6 ,6 ,4 ,6 ,7 ,6 ,4 ,10,7 ,10,10,6 ,6 ,
  214. c      .             7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,7 ,
  215. c      .             0 ,1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10,11,12,13,14,15,
  216. c      .             9 ,1 ,0 ,11,1 ,4 ,4 ,7 ,9 ,9 ,14,1 ,1 ,9 ,15,9 ,
  217. c      .             10,10,0 ,2 ,1 ,12,10,7 ,10,14,10,7 ,5 ,14,14,10,
  218. c      .             11,11,0 ,9 ,7 ,1 ,7 ,7 ,11,1 ,7 ,11,12,1 ,15,11,
  219. c      .             12,12,0 ,13,13,13,10,7 ,12,1 ,5 ,12,12,1 ,5 ,12,
  220. c      .             13,13,0 ,2 ,7 ,5 ,10,7 ,13,9 ,14,1 ,1 ,13,5 ,13,
  221. c      .             14,14,0 ,15,2 ,15,6 ,7 ,14,15,14,15,5 ,5 ,14,14,
  222. c      .             15,1 ,2 ,3 ,4 ,5 ,6 ,7 ,15,9 ,10,11,12,13,14,15 /
  223.       DATA ITABM1 /  0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,11,12,13,14,
  224.      &               15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,
  225.      .                1, 1, 30,9, 12,11,12,13,22,1, 16,1, 11,11,12,
  226.      &               11,9, 30,30,9, 13,12,1, 22,12,11,30,5, 9, 30,1,
  227.      .                2, 30,2, 16,17,29,10,19,18,16,2, 16,17,29,29,
  228.      &               19,18,18,2, 2, 19,10,18,17,26,19,29,19,19,18,18,
  229.      .                3, 9, 16,3, 12,9, 19,28,30,3, 19,9, 16,28,16,
  230.      &               28,3, 16,16,16,28,19,9, 30,28,28,16,28,3, 16,9,
  231.      .                4, 12,17,12,4, 12,24,24,23,12,14,12,4, 12,4,
  232.      &               12,12,23,17,14,14,14,23,4, 14,12,14,14,14,14,12,
  233.      .                5, 11,29,9, 12,5, 14,13,12,11,14,5, 5, 5, 12,
  234.      &               13,13,12,12,15,15,14,12,12,12,13,14,15,15,12,11,
  235.      .                6, 12,10,19,24,14,6, 27,17,29,21,14,14,27,24,
  236.      &               27,19,26,26,26,27,6, 17,14,6, 27,21,6, 27,26,29,
  237.      .                7, 13,19,28,24,13,27,7, 12,28,27,13,15,25,27,
  238.      &               25,28,29,28,28,25,27,13,12,27,0, 27,20,25,20,15/
  239.       DATA ITABM2 /  8, 22,18,30,23,12,17,12,8, 30,17,22,23,12,17,
  240.      &               12,17,8, 17,17,17,17,8, 8, 17,12,17,17,17,17,22,
  241.      .                9, 1, 16,3, 12,11,29,28,30,9, 16,9, 30,9, 16,
  242.      &               16,9, 16,16,16,28,29,30,30,29,15,16,16,16,16,9,
  243.      .                10,16,2, 19,14,14,21,27,17,16,10,12,14,27,26,
  244.      &               27,19,29,29,26,27,10,17,17,26,27,10,26,19,26,17,
  245.      .                11,1, 16,9, 12,5, 14,13,22,9, 12,11,5, 5, 12,
  246.      &               13,9, 12,30,12,13,14,1, 12,12,13,12,13,13,12,11,
  247.      .                12,11,17,16,4, 5, 14,15,23,30,14,5, 12,5, 12,
  248.      &               5, 29,12,17,29,15,14,12,12,14,13,14,14,15,14,12,
  249.      .                13,11,29,28,12,5, 27,25,12,9, 27,5, 5, 13,15,
  250.      &               13,15,12,29,15,15,27,11,12,27,15,27,15,15,15,11,
  251.      .                14,12,29,16,4, 12,24,27,17,16,26,12,12,15,14,
  252.      &               27,29,29,29,29,27,26,12,12,14,27,14,14,27,14,17,
  253.      .                15,11,19,28,12,13,27,25,12,16,27,13,5, 13,27,
  254.      &               15,28,29,29,28,15,27,12,12,27,15,27,20,15,27,16/
  255.       DATA ITABM3 / 16,9, 18,3, 12,13,19,28,17,9, 19,9, 29,15,29,
  256.      &               28,16,18,16,16,28,19,30,17,29,28,29,19,16,16,9,
  257.      .                17,30,18,16,23,12,26,29,8, 16,29,12,12,12,29,
  258.      &               29,18,17,17,29,29,26,30,17,14,29,29,29,29,17,17,
  259.      .                18,30,2, 16,17,12,26,28,17,16,29,30,17,29,29,
  260.      &               29,16,17,18,29,29,26,30,17,29,29,29,29,19,18,18,
  261.      .                19,9, 2, 16,14,15,26,28,17,16,26,12,29,15,29,
  262.      &               28,16,29,29,19,27,26,17,17,26,28,29,19,19,19,16,
  263.      .                20,13,19,28,14,15,27,25,17,28,27,13,15,15,27,
  264.      &               15,28,29,29,27,20,27,12,14,27,15,27,20,20,27,16,
  265.      .                21,12,10,19,14,14,6, 27,17,29,10,14,14,27,26,
  266.      &               27,19,26,26,26,27,21,17,14,6, 27,10,21,27,26,29,
  267.      .                22,1, 18,9, 23,12,17,13,8, 30,17,1, 12,11,12,
  268.      &               12,30,30,30,17,12,17,22,8, 12,5, 17,12,16,17,22,
  269.      .                23,22,17,30,4, 12,14,12,8, 30,17,12,12,12,12,
  270.      &               12,17,17,17,17,14,14,8, 23,14,12,17,14,12,17,22/
  271.       DATA ITABM4 / 24,12,26,28,14,12,6, 27,17,29,26,12,14,27,14,
  272.      &               27,29,14,29,26,27,6, 12,14,24,27,14,24,27,14,12,
  273.      .                25,11,19,28,12,13,27,0, 12,15,27,13,13,15,27,
  274.      &               15,28,29,29,28,15,27,5, 12,27,25,27,20,15,28,13,
  275.      .                26,30,29,16,14,14,21,27,17,16,10,12,14,27,14,
  276.      &               27,29,29,29,29,27,10,17,17,14,27,26,26,19,26,17,
  277.      .                27,5, 19,28,14,15,6, 20,17,16,26,13,14,15,14,
  278.      &               20,19,29,29,19,20,21,12,14,24,20,26,27,20,27,29,
  279.      .                28,9, 19,3, 14,15,27,25,17,16,19,13,15,15,27,
  280.      &               15,16,29,19,19,20,27,16,12,27,15,19,20,28,19,16,
  281.      .                29,30,18,16,14,12,26,20,17,16,26,12,14,15,14,
  282.      &               27,16,17,18,19,27,26,17,17,14,28,26,27,19,29,16,
  283.      .                30,1, 18,9, 12,11,29,15,22,9, 17,11,12,11,17,
  284.      &               16,9, 17,18,16,16,29,22,22,12,13,17,29,16,16,30/
  285.  
  286.  
  287.       DATA DENSIT / 0. /
  288.       DATA ILCOUR / 0 /
  289.  
  290. C  Noms des elements support geometrique
  291. C  ATTENTION : nombre LIMITE a 100 dans CCGEOME.INC)
  292.       DATA NOMBR / 49 /
  293.       DATA NOMS( 1),NOMS( 2),NOMS( 3),NOMS( 4),NOMS( 5),NOMS( 6),
  294.      .     NOMS( 7),NOMS( 8),NOMS( 9),NOMS(10),NOMS(11),NOMS(12),
  295.      .     NOMS(13),NOMS(14),NOMS(15),NOMS(16),NOMS(17),NOMS(18),
  296.      .     NOMS(19),NOMS(20),NOMS(21),NOMS(22),NOMS(23),NOMS(24),
  297.      .     NOMS(25),NOMS(26),NOMS(27),NOMS(28),NOMS(29),NOMS(30),
  298.      .     NOMS(31),NOMS(32),NOMS(33),NOMS(34),NOMS(35),NOMS(36),
  299.      .     NOMS(37),NOMS(38),NOMS(39),NOMS(40),NOMS(41),NOMS(42),
  300.      .     NOMS(43),NOMS(44),NOMS(45),NOMS(46),NOMS(47),NOMS(48),
  301.      .     NOMS(49)
  302.      . /   'POI1','SEG2','SEG3','TRI3','TRI4','TRI6',
  303.      .     'TRI7','QUA4','QUA5','QUA8','QUA9','RAC2',
  304.      .     'RAC3','CUB8','CU20','PRI6','PR15','LIA3',
  305.      .     'LIA4','LIA6','LIA8','MULT','TET4','TE10',
  306.      .     'PYR5','PY13','ATTA','SUPE','RAP3','LIP6',
  307.      .     'LIP8','POLY','CU27','PR21','TE15','PY19',
  308.      .     'SEG4','QU16','TR12','PR18','SEG6','TR21',
  309.      .     'QU36','C216','P126','TE56','PY91','SURE',
  310.      .     'DOUB' /
  311. C  Dimension des elements geometriques
  312. C  -1 si pas de sens
  313.       DATA LDLR( 1),LDLR( 2),LDLR( 3),LDLR( 4),LDLR( 5),LDLR( 6),
  314.      .     LDLR( 7),LDLR( 8),LDLR( 9),LDLR(10),LDLR(11),LDLR(12),
  315.      .     LDLR(13),LDLR(14),LDLR(15),LDLR(16),LDLR(17),LDLR(18),
  316.      .     LDLR(19),LDLR(20),LDLR(21),LDLR(22),LDLR(23),LDLR(24),
  317.      .     LDLR(25),LDLR(26),LDLR(27),LDLR(28),LDLR(29),LDLR(30),
  318.      .     LDLR(31),LDLR(32),LDLR(33),LDLR(34),LDLR(35),LDLR(36),
  319.      .     LDLR(37),LDLR(38),LDLR(39),LDLR(40),LDLR(41),LDLR(42),
  320.      .     LDLR(43),LDLR(44),LDLR(45),LDLR(46),LDLR(47),LDLR(48)
  321.      . /        0,     1,     1,     2,     2,     2,
  322.      .          2,     2,     2,     2,     2,     1,
  323.      .          1,     3,     3,     3,     3,     2,
  324.      .          2,     2,     2,    -1,     3,     3,
  325.      .          3,     3,    -1,    -1,     2,     3,
  326.      .          3,     2,     3,     3,     3,     3,
  327.      .          1,     2,     2,     3,     1,     2,
  328.      .          2,     3,     3,     3,     3,    -1 /
  329.  
  330. C  Nombre de noeuds par element geometrique
  331.       DATA NBNNE( 1),NBNNE( 2),NBNNE( 3),NBNNE( 4),NBNNE( 5),NBNNE( 6),
  332.      .     NBNNE( 7),NBNNE( 8),NBNNE( 9),NBNNE(10),NBNNE(11),NBNNE(12),
  333.      .     NBNNE(13),NBNNE(14),NBNNE(15),NBNNE(16),NBNNE(17),NBNNE(18),
  334.      .     NBNNE(19),NBNNE(20),NBNNE(21),NBNNE(22),NBNNE(23),NBNNE(24),
  335.      .     NBNNE(25),NBNNE(26),NBNNE(27),NBNNE(28),NBNNE(29),NBNNE(30),
  336.      .     NBNNE(31),NBNNE(32),NBNNE(33),NBNNE(34),NBNNE(35),NBNNE(36),
  337.      .     NBNNE(37),NBNNE(38),NBNNE(39),NBNNE(40),NBNNE(41),NBNNE(42),
  338.      .     NBNNE(43),NBNNE(44),NBNNE(45),NBNNE(46),NBNNE(47),NBNNE(48)
  339.      .   /   1,  2,  3,  3,  4,  6,  7,  4,  5,  8,  9,  4,  6,  8, 20,
  340.      .       6, 15,  6,  8, 12, 16,  0,  4, 10,  5, 13,  0,  0,  8, 15,
  341.      .      20,  0, 27, 21, 15, 19,  4, 16, 12, 18,  6, 21, 36,216,126,
  342.      .      56, 91, 0 /
  343. C  Degre de l'interpolation geometrique (nombre de noeuds par segment)
  344.       DATA KDEGRE( 1),KDEGRE( 2),KDEGRE( 3),KDEGRE( 4),KDEGRE( 5),
  345.      .     KDEGRE( 6),KDEGRE( 7),KDEGRE( 8),KDEGRE( 9),KDEGRE(10),
  346.      .     KDEGRE(11),KDEGRE(12),KDEGRE(13),KDEGRE(14),KDEGRE(15),
  347.      .     KDEGRE(16),KDEGRE(17),KDEGRE(18),KDEGRE(19),KDEGRE(20),
  348.      .     KDEGRE(21),KDEGRE(22),KDEGRE(23),KDEGRE(24),KDEGRE(25),
  349.      .     KDEGRE(26),KDEGRE(27),KDEGRE(28),KDEGRE(29),KDEGRE(30),
  350.      .     KDEGRE(31),KDEGRE(32),KDEGRE(33),KDEGRE(34),KDEGRE(35),
  351.      .     KDEGRE(36),KDEGRE(37),KDEGRE(38),KDEGRE(39),KDEGRE(40),
  352.      .     KDEGRE(41),KDEGRE(42),KDEGRE(43),KDEGRE(44),KDEGRE(45),
  353.      .     KDEGRE(46),KDEGRE(47),KDEGRE(48)
  354. *          1                              11                        20
  355.      .   /  1, 2, 3, 2, 2, 3, 3, 2, 2, 3, 3, 2, 3, 2, 3, 2, 3, 2, 2, 3,
  356. *          21                             31                        40
  357.      .      3, 0, 2, 3, 2, 3, 0, 0, 3, 3, 3, 2, 3, 3, 3, 3, 4, 4, 4, 2,
  358.      .      6, 6, 6, 6, 6, 6, 6, 0 /
  359. C  Faces des elements (= numero element support des faces dans NOMS)
  360.       DATA KSURF( 1),KSURF( 2),KSURF( 3),KSURF( 4),KSURF( 5),KSURF( 6),
  361.      .     KSURF( 7),KSURF( 8),KSURF( 9),KSURF(10),KSURF(11),KSURF(12),
  362.      .     KSURF(13),KSURF(14),KSURF(15),KSURF(16),KSURF(17),KSURF(18),
  363.      .     KSURF(19),KSURF(20),KSURF(21),KSURF(22),KSURF(23),KSURF(24),
  364.      .     KSURF(25),KSURF(26),KSURF(27),KSURF(28),KSURF(29),KSURF(30),
  365.      .     KSURF(31),KSURF(32),KSURF(33),KSURF(34),KSURF(35),KSURF(36),
  366.      .     KSURF(37),KSURF(38),KSURF(39),KSURF(40),KSURF(41),KSURF(42),
  367.      .     KSURF(43),KSURF(44),KSURF(45),KSURF(46),KSURF(47),KSURF(48)
  368. *          1                             11                         20
  369.      .   /  0, 0, 0, 4, 5, 6, 7, 8, 9,10,11, 8, 0, 8,10, 8,10, 8, 6, 8,
  370. *          21                            31                         40
  371.      .     10, 0, 4, 6, 8,10, 0, 0,29, 6,10,32,11,11, 7,11,37,38,39,40,
  372.      .      0,42,43,43,43,42,43, 0 /
  373. C  Nombre de "sommets" de l'element geometrique
  374.       DATA NBSOM( 1),NBSOM( 2),NBSOM( 3),NBSOM( 4),NBSOM( 5),NBSOM( 6),
  375.      .     NBSOM( 7),NBSOM( 8),NBSOM( 9),NBSOM(10),NBSOM(11),NBSOM(12),
  376.      .     NBSOM(13),NBSOM(14),NBSOM(15),NBSOM(16),NBSOM(17),NBSOM(18),
  377.      .     NBSOM(19),NBSOM(20),NBSOM(21),NBSOM(22),NBSOM(23),NBSOM(24),
  378.      .     NBSOM(25),NBSOM(26),NBSOM(27),NBSOM(28),NBSOM(29),NBSOM(30),
  379.      .     NBSOM(31),NBSOM(32),NBSOM(33),NBSOM(34),NBSOM(35),NBSOM(36),
  380.      .     NBSOM(37),NBSOM(38),NBSOM(39),NBSOM(40),NBSOM(41),NBSOM(42),
  381.      .     NBSOM(43),NBSOM(44),NBSOM(45),NBSOM(46),NBSOM(47),NBSOM(48)
  382. *          1                             11                         20
  383.      .   /  1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 2, 8, 8, 6, 6, 6, 8, 3,
  384. *          21                            31                         40
  385.      .      4, 0, 4, 4, 5, 5, 0, 0, 4, 6, 8, 0, 8, 6, 4, 5, 2, 4, 3, 6,
  386.      .      2, 3, 4, 8, 6, 4, 5, 0 /
  387. C  "Position" dans le tableau IBSOM pour avoir numeros locaux sommets
  388.       DATA NSPOS( 1),NSPOS( 2),NSPOS( 3),NSPOS( 4),NSPOS( 5),NSPOS( 6),
  389.      .     NSPOS( 7),NSPOS( 8),NSPOS( 9),NSPOS(10),NSPOS(11),NSPOS(12),
  390.      .     NSPOS(13),NSPOS(14),NSPOS(15),NSPOS(16),NSPOS(17),NSPOS(18),
  391.      .     NSPOS(19),NSPOS(20),NSPOS(21),NSPOS(22),NSPOS(23),NSPOS(24),
  392.      .     NSPOS(25),NSPOS(26),NSPOS(27),NSPOS(28),NSPOS(29),NSPOS(30),
  393.      .     NSPOS(31),NSPOS(32),NSPOS(33),NSPOS(34),NSPOS(35),NSPOS(36),
  394.      .     NSPOS(37),NSPOS(38),NSPOS(39),NSPOS(40),NSPOS(41),NSPOS(42),
  395.      .     NSPOS(43),NSPOS(44),NSPOS(45),NSPOS(46),NSPOS(47),NSPOS(48)
  396. *            1                                      11              15
  397.      .   /   1,  2,  4,  6,  9, 12, 15, 18, 22, 26, 30, 18,  4, 34, 42,
  398. *           16                  21                                  30
  399.      .      50, 56, 50, 34, 12, 26,  0, 62, 66, 70, 75,  0,  0, 80, 84,
  400. *           31                                      41              45
  401.      .      90, 98, 42, 56, 66, 75,126,126,126, 56,130,132,135,139,147,
  402.      .     153,157,  0 /
  403. C  Tableau des sommets par element (numero local des sommets)
  404.       DATA IBSOM
  405. *         1  2    4    6      9      12     15     18       22
  406.      .  / 1, 1,2, 1,3, 1,2,3, 1,2,3, 1,3,5, 1,3,5, 1,2,3,4, 1,2,3,4,
  407. *         26       30       34               42
  408.      .    1,3,5,7, 1,3,5,7, 1,2,3,4,5,6,7,8, 1,3,5,7,13,15,17,19,
  409. *         50           56              62       66        70
  410.      .    1,2,3,4,5,6, 1,3,5,10,12,14 ,1,2,3,4, 1,3,5,10, 1,2,3,4,5,
  411. *         75          80       84            90
  412.      .    1,3,5,7,13, 1,3,4,6, 1,3,5,7,9,11, 1,3,5,7,9,11,13,15,
  413. *         98                                112
  414.      .    1,2,3,4,5,6,7,8,9,10,11,12,13,14, 1,3,5,7,13,15,17,19,
  415. *         120             126       130  132     135
  416.      .    1,3,5,10,12,14, 1,4,7,10, 1,6, 1,6,11, 1,6,11,16,
  417. *         139                        147                 153
  418.      .    1,6,11,16,181,186,191,196, 1,6,11,106,111,116, 1,6,11,56,
  419. *         157
  420.      .    1,6,11,16,91, 139*0 /
  421. C  Nombre de segments par element (entre 2 noeuds successifs)
  422.       DATA LPL( 1),LPL( 2),LPL( 3),LPL( 4),LPL( 5),LPL( 6),LPL( 7),
  423.      .     LPL( 8),LPL( 9),LPL(10),LPL(11),LPL(12),LPL(13),LPL(14),
  424.      .     LPL(15),LPL(16),LPL(17),LPL(18),LPL(19),LPL(20),LPL(21),
  425.      .     LPL(22),LPL(23),LPL(24),LPL(25),LPL(26),LPL(27),LPL(28),
  426.      .     LPL(29),LPL(30),LPL(31),LPL(32),LPL(33),LPL(34),LPL(35),
  427.      .     LPL(36),LPL(37),LPL(38),LPL(39),LPL(40),LPL(41),LPL(42),
  428.      .     LPL(43),LPL(44),LPL(45),LPL(46),LPL(47),LPL(48)
  429. *           1                            11                         20
  430.      .   /  0, 1, 2, 3, 3, 6, 6, 4, 4, 8, 8, 4, 2,12,24, 9,18, 9,12, 6,
  431. *           21                           31                         40
  432.      .      8, 0, 6,12, 8,16, 0, 0, 8,18,24, 0,24,18,12,16,12, 9, 3,18,
  433.      .      5,15,20,60,45,30,40, 0 /
  434. C  "Position" dans tableau KSEGM de la liste des segments d'un element
  435.       DATA LPT( 1),LPT( 2),LPT( 3),LPT( 4),LPT( 5),LPT( 6),LPT( 7),
  436.      .     LPT( 8),LPT( 9),LPT(10),LPT(11),LPT(12),LPT(13),LPT(14),
  437.      .     LPT(15),LPT(16),LPT(17),LPT(18),LPT(19),LPT(20),LPT(21),
  438.      .     LPT(22),LPT(23),LPT(24),LPT(25),LPT(26),LPT(27),LPT(28),
  439.      .     LPT(29),LPT(30),LPT(31),LPT(32),LPT(33),LPT(34),LPT(35),
  440.      .     LPT(36),LPT(37),LPT(38),LPT(39),LPT(40),LPT(41),LPT(42),
  441.      .     LPT(43),LPT(44),LPT(45),LPT(46),LPT(47),LPT(48)
  442. *            1                                      11              15
  443.      .   /   0,  1,  3,  7, 13, 19, 31, 43, 51, 59, 75, 43,  3, 91,115,
  444. *           16                  21                                  30
  445.      .     163,181,163, 91, 19, 59,  0,217,229,253,269,  0,  0,301,317,
  446. *           31                                      41              45
  447.      .     353,401,115,181,229,269,401,401,401,181,427,437,467,467,587,
  448.      .     677,737,  0 /
  449.       DATA KSEGM1 /
  450.      .  1,2, 1,2,2,3, 1,2,2,3,3,1, 1,2,2,3,3,1, 1,2,2,3,3,4,4,5,5,6,6,1,
  451.      .  1,2,2,3,3,4,4,5,5,6,6,1, 1,2,2,3,3,4,4,1, 1,2,2,3,3,4,4,1,
  452.      .  1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,1,  1,2,2,3,3,4,4,5,5,6,6,7,7,8,
  453.      .  8,1,  1,2,2,3,3,4,4,1,5,6,6,7,7,8,8,5,1,5,2,6,3,7,4,8,
  454.      .  1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,1,
  455.      .  1,9,9,13,3,10,10,15,5,11,11,17,7,12,12,19,13,14,14,15,15,
  456.      .  16,16,17,17,18,18,19,19,20,20,13,
  457.      .  1,2,2,3,3,1,4,5,5,6,6,4,1,4,2,5,3,6, 1,2,2,3,3,4,4,5,5,6,6,1,
  458.      .  10,11,11,12,12,13,13,14,14,15,15,10,
  459.      .  1,7,7,10,3,8,8,12,5,9,9,14, 1,2,2,3,3,1,1,4,2,4,3,4,
  460.      .  1,2,2,3,3,4,4,5,5,6,6,1,1,7,7,10,3,8,8,10,5,9,9,10,
  461.      .  1,2,2,3,3,4,4,1,1,5,2,5,3,5,4,5,
  462.      .  1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,1,1,9,9,13,3,10,10,13,5,11,11,
  463.      .  13,7,12,12,13, 1,2,2,3,3,8,8,4,4,5,5,6,6,7,7,1,
  464.      .  1,2,2,3,3,4,4,5,5,6,6,1, 1,13,13,7, 3,14,14,9, 5,15,15,11,
  465.      .  7,8,8,9,9,10,10,11,11,12,12,7, 1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,1,
  466.      .  1,17,17,9, 3,18,18,11, 5,19,19,13,
  467.      .  7,20,20,15,9,10,10,11,11,12,12,13,13,14,14,15,15,16,16,9,
  468.      .  1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13,14 /
  469.       DATA KSEGM2 /
  470.      . 1,2,2,3,3,4,4,5,5,6, 1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,
  471.      . 11,11,12,12,13,13,14,14,15,15,1, 1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,
  472.      . 9,9,10,10,11,11,12,12,13,13,14,14,15,15,16,16,17,17,18,18,19,19,
  473.      . 20,20,1, 181,182,182,183,183,184,184,185,185,186,186,187,187,
  474.      . 188,188,189,189,190,190,191,191,192,192,193,193,194,194,195,195,
  475.      . 196,196,197,197,198,198,199,199,200,200,181, 1,37,37,73,73,109,
  476.      . 109,145,145,181, 6,42,42,78,78,114,114,150,150,186, 11,47,47,83,
  477.      . 83,119,119,155,155,191, 16,52,52,88,88,124,124,160,160,196,
  478.      . 1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13,14,14,
  479.      . 15,15,1, 106,107,107,108,108,109,109,110,110,111,111,112,112,113,
  480.      . 113,114,114,115,115,116,116,117,117,118,118,119,119,120,120,106 ,
  481.      . 1,22,22,43,43,64,64,85,85,106, 6,27,27,48,48,69,69,90,90,111,
  482.      . 11,32,32,53,53,74,74,95,95,116,
  483.      . 1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13,14,14,
  484.      . 15,15,1, 1,22,22,37,37,47,47,53,53,56, 6,26,26,40,40,49,49,54,54,
  485.      . 56, 11,30,30,43,43,51,51,55,55,56 /
  486.       DATA KSEGM3 /
  487.      . 1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13,14,14,
  488.      . 15,15,16,16,17,17,18,18,19,19,20,20,1, 1,37,37,62,62,78,78,87,87,
  489.      . 91, 6,41,41,65,65,80,80,88,88,91, 11,45,45,68,68,82,82,89,89,91,
  490.      . 16,49,49,71,71,84,84,90,90,91, 684*0 /
  491. C  KDFAC : 1) nombre de points de la face
  492. C          2) adresse de description des triangles dans KFAC
  493. C          3) nombre de triangles dans la face
  494.       DATA KDFAC /3,1,1, 4,4,2, 6,10,4, 8,22,6, 6,40,4, 0,0,0, 7,52,6,
  495.      .     9,70,8,  21,94,25,  36,169,50, 30*0/
  496. C  FACES : T3 Q4 T6 Q8 ?? POLY T7 Q9 T21 Q36
  497. C  KFAC   DESCRIPTION DES TRIANGLES DANS LA FACE
  498. *                 1      4            10
  499.       DATA KFAC / 1,2,3, 1,2,3,1,3,4, 1,2,6, 2,3,4, 4,5,6, 4,6,2,
  500. *          22
  501.      .     1,2,8, 2,3,4, 4,5,6, 7,8,6, 2,4,8, 4,6,8,
  502. *          40
  503.      .     1,2,5, 1,5,6, 2,3,4, 2,4,5,
  504. *          52
  505.      .     1,2,7, 2,3,7, 3,4,7, 4,5,7, 5,6,7, 6,1,7,
  506. *          70
  507.      .     1,2,9, 2,3,9, 3,4,9, 4,5,9, 5,6,9, 6,7,9, 7,8,9, 8,1,9,
  508. *          94
  509.      .      1, 2,15,  2, 3,16,  3, 4,17,  4, 5,18,  5, 6, 7,
  510.      .      2,16,15,  3,17,16,  4,18,17,  5, 7,18, 15,16,14,
  511.      .     16,17,19, 17,18,20, 18, 7, 8, 16, 9,14, 17,20,19,
  512.      .     18, 8,20, 14,19,13, 19,20,21, 20, 8, 9, 19,21,13,
  513.      .     20, 9,21, 13,21,12, 21, 9,10, 21,10,12, 12,10,11,
  514. *          169
  515.      .      1,21,20,  1, 2,21,  2, 3,21,  3,22,21,  3,23,22,
  516.      .      3, 4,23,  4, 5,23,  5,24,23,  5, 7,24,  5, 6, 7,
  517.      .     20,21,19, 21,25,19, 21,26,25, 21,22,26, 22,23,26,
  518.      .     23,27,26, 23,28,27, 23,24,28, 24, 7,28,  7, 8,28,
  519.      .     19,29,18, 19,25,29, 25,26,29, 26,30,29, 26,31,30,
  520.      .     26,27,31, 27,28,31, 28,32,31, 28, 9,32, 28, 8, 9,
  521.      .     18,29,17, 29,33,17, 29,34,33, 29,30,34, 30,31,34,
  522.      .     31,35,34, 31,36,35, 31,32,36, 32, 9,36,  9,10,36,
  523.      .     17,15,16, 17,33,15, 33,34,15, 34,14,15, 34,13,14,
  524.      .     34,35,13, 35,36,13, 36,12,13, 36,11,12, 36,10,11, 682*0 /
  525. C
  526. C  LTEL : 1) nombre de faces dans l'element
  527. C         2) adresse dans LDEL de la description des faces
  528. C               POI1 SEG2 SEG3 TRI3 TRI4 TRI6 TRI7  QUA4 QUA5 QUA8 QUA9
  529.       DATA LTEL/0,0, 0,0, 0,0, 1,1, 1,1, 1,2, 1,93, 1,3, 1,3, 1,4, 1,94,
  530. C      RAC2 RAC3 CUB8 CU20  PRI6  PR15  LIA3  LIA4 LIA6 LIA8 MULT TET4
  531.      . 1,3, 0,0, 6,5, 6,11, 5,17, 5,22, 5,17, 6,5, 1,2, 1,4, 0,0, 4,27,
  532. C      TE10  PYR5  PY13  ATTA SUPE RAP3  LIP6  LIP8  POLY
  533.      . 4,31, 5,35, 5,40, 0,0, 0,0, 1,45, 2,46, 2,48, 1,50,
  534. C      CU27  PR21  TE15  PY19
  535.      . 6,73, 5,79, 4,84, 5,88,
  536. C      SEG4 QU16 TR12 PR18 SEG6 TR21  QU36  C216  P126  TE56  PY91
  537.      . 0,0, 0,0, 0,0, 0,0, 0,0, 1,51, 1,52, 6,53, 5,59, 4,64, 5,68,
  538. C      SURE
  539.      . 0,0, 0,0, 0,0, 100*0 /
  540.  
  541. C  LDEL : 1) type de la face
  542. C         2) adresse de la face dans LFAC
  543. c                 T3/4 T6/7 Q4/5 Q8/9 <(5)----------CUB8------------>
  544. c       DATA LDEL / 1,1, 3,1, 2,1, 4,1, 2,1, 2,5, 2,9, 2,13, 2,17, 2,21,
  545.       DATA LDEL / 1,1, 3,1, 2,1, 4,1, 2,1,2,684,2,688,2,692, 2,17, 2,21,
  546. c     <(11)-------------CU20-----------> <(17)---------PRI6-------->
  547.      . 4,1, 4,25, 4,33, 4,41, 4,49, 4,57, 1,1,1,696,2,699,2,703,2,707,
  548. c     <(22)-----------PR15-------> <(27)--------TET4------->
  549.      . 3,1, 3,77, 4,83, 4,91, 4,99, 1,1, 1,107, 1,110, 1,113,
  550. c     <(31)--------TE10-------> <(35)------------PYR5---------->
  551.      . 3,1, 3,116, 3,122, 3,128, 2,1, 1,134, 1,137, 1,140, 1,143,
  552. c     <(40)-----------PY13----------->  RAP3  <(46)-LIP6> <(48)-LIP8>
  553.      . 4,1, 3,146, 3,152, 3,158, 3,164, 4,170, 3,1, 3,178, 4,1, 4,184,
  554. c     <POLY>  TR21   QU36
  555.      . 6,192, 9,192, 10,192,
  556. C     <-----------------------C216------------------->
  557.      . 10,192, 10,228, 10,264, 10,300, 10,336, 10,372,
  558. c     <--------------P126----------------->  <-----------TE56--------->
  559.      . 7,192, 7,408, 10,429, 10,465, 10,501, 7,192, 7,537, 7,558, 7,579,
  560. c     <(68)----------PY91--------------->  <(73)------CU27----------
  561.      . 10,192, 9,600, 9,621, 9,642, 9,663, 8,720,  8,729,  8,738, 8,747,
  562. c      ---------->  <(79)-----------PR21------------>  <(84)----TE15-
  563.      . 8,756,8,765,  7,774, 7,781, 8,788, 8,797, 8,806, 7,815,  7,822,
  564. c      ----------->  <(88)------------PY19----------->  (93) T7 (94) Q9
  565.      . 7,829, 7,836, 8,843, 7,852, 7,859, 7,866, 7,873, 7,711 , 8,711,
  566.      .   12*0/
  567.  
  568. C  LFAC : description de la face
  569. c BP2016 : on oriente les faces vers l'interieur des elements
  570. c          (sauf pour les elements C216, P126 TE56 et PY91)
  571. c rem : chiffres 1 a 8 = FACE1 du CUB8 CU20 PRI6 PR15 TET4 TE10 PYR5 ...
  572. c 1:24        CUB8 <FACE1>           <FACE3>  <FACE4>  <FACE5>  <FACE6>
  573.       DATA LFAC1 / 1,2,3,4, 5,6,7,8, 1,2,6,5, 2,3,7,6, 4,3,7,8, 1,4,8,5,
  574. c 25:40
  575. c CU20 <     FACE2           >  <      FACE3      >
  576.      . 13,20,19,18,17,16,15,14, 1,9,13,14,15,10,3,2,
  577. c 41:64
  578. c CU20 <     FACE4        >  <      FACE5       > <      FACE6       >
  579.      . 3,10,15,16,17,11,5,4, 7,6,5,11,17,18,19,12, 1,8,7,12,19,20,13,9,
  580. c 65:82
  581. c PRI6 <FACE3>  <FACE4>  <FACE5>   PR15 <     FAC2      >
  582.      . 1,2,5,4, 2,3,6,5, 3,1,4,6,       10,15,14,13,12,11,
  583. c 83:106
  584. c PR15 <     FACE3      >  <      FACE4     >  <      FACE5     >
  585.      . 1,7,10,11,12,8,3,2, 3,8,12,13,14,9,5,4, 5,9,14,15,10,7,1,6,
  586. c 107:133
  587. c TET4<FAC2><FAC3><FAC4>  TE10 <  FAC2 >  <   FAC3   > <    FAC4   >
  588.      . 1,4,2,2,4,3,3,4,1,   1,7,10,8,3,2, 3,8,10,9,5,4, 5,9,10,7,1,6,
  589. c 134:157
  590. c PYR5 <FAC2><FAC3><FAC4><FAC5>  PY13 <    FAC2   >  <    FAC3    >
  591.      . 1,5,2,2,5,3,3,5,4,4,5,1,       1,9,13,10,3,2, 3,10,13,11,5,4,
  592. c 158:169
  593. c PY13  <    FAC4   >  <    FAC5   >
  594.      . 5,11,13,12,7,6, 7,12,13,9,1,8,
  595. c 170:191 <RAP3 FACE1>  < LIP6 FACE2 >  <   LIP8     FACE2   >
  596.      . 1,2,3,8,4,5,6,7, 7,8,9,10,11,12, 9,10,11,12,13,14,15,16,
  597. c 192:215 <  liste 1 a 36  pour les FACES 1  des  POLY TR21 QU36
  598.      . 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,
  599. c 216:233  C216 P126 TE56  PY91 ...      >
  600.      . 25,26,27,28,29,30,31,32,33,34,35,36, 181,182,183,184,185,186,
  601. c 234:249
  602.      . 187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,
  603. c 250:263
  604.      . 203,204,205,206,207,208,209,210,211,212,213,214,215,216,
  605. c 264:283
  606.      . 1,2,3,4,5,6,42,78,114,150,186,185,184,183,182,181,145,109,73,37,
  607. c 284:299
  608.      . 38,39,40,41,74,75,76,77,110,111,112,113,146,147,148,149,
  609. c 300:318
  610.      . 6,7,8,9,10,11,47,83,119,155,191,190,189,188,187,186,150,114,78,
  611. c 319:335
  612.      . 42,43,44,45,46,79,80,81,82,115,116,117,118,151,152,153,154,
  613. c 336:353
  614.      . 11,12,13,14,15,16,52,88,124,160,196,195,194,193,192,191,155,119,
  615. c 354:371
  616.      . 83,47,48,49,50,51,84,85,86,87,120,121,122,123,156,157,158,159 /
  617. C  Tableau LFAC2 :
  618. c 372:386
  619.       DATA LFAC2 / 16,17,18,19,20,1,37,73,109,145,181,200,199,198,197,
  620. c 387:404
  621.      . 196,160,124,88,52,53,54,55,56,89,90,91,92,125,126,127,128,161,
  622. c 405:407
  623.      . 162,163,164,
  624. c 48:423
  625.      . 106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,
  626. c 424:442
  627.      . 122,123,124,125,126, 1,2,3,4,5,6,27,48,69,90,111,110,109,108,
  628. c 443:462
  629.      . 107,106,85,64,43,22,23,24,25,26,44,45,46,47,65,66,67,68,86,87,
  630. c 463:482
  631.      . 88,89, 6,7,8,9,10,11,32,53,74,95,116,115,114,113,112,111,90,69,
  632. c 483:500
  633.      . 48,27,28,29,30,31,49,50,51,52,70,71,72,73,91,92,93,94,
  634. c 501:519
  635.      . 11,12,13,14,15,1,22,43,64,85,106,120,119,118,117,116,95,74,53,
  636. c 520:536
  637.      . 32,33,34,35,36,54,55,56,57,75,76,77,78,96,97,98,99,
  638. c 537:557
  639.      . 1,2,3,4,5,6,26,40,49,54,56,53,47,37,22,23,24,25,38,39,48,
  640. c 558:578
  641.      . 6,7,8,9,10,11,30,43,51,55,56,54,49,40,26,27,28,29,41,42,50,
  642. c 579:599
  643.      . 11,12,13,14,15,1,22,37,47,53,56,55,51,43,30,31,32,33,44,45,52,
  644. c 600:620
  645.      . 1,2,3,4,5,6,41,65,80,88,91,87,78,62,37,38,39,40,63,64,79,
  646. c 621:641
  647.      . 6,7,8,9,10,11,45,68,82,89,91,88,80,65,41,42,43,44,66,67,81,
  648. c 642:662
  649.      . 11,12,13,14,15,16,49,71,84,90,91,89,82,68,45,46,47,48,69,70,83,
  650. c 663:683
  651.      . 16,17,18,19,20,1,37,62,78,87,91,90,84,71,49,50,51,52,72,73,85,
  652. c     684       688      692
  653.      . 5,8,7,6, 1,5,6,2, 2,6,7,3,
  654. c     696     699      703      707
  655.      . 4,6,5, 1,4,5,2, 2,5,6,3, 3,6,4,1 /
  656. C  Tableau LFAC3 :
  657. c 711:728
  658. C                 T7/Q9               CU27  < FACE1     >
  659.       DATA LFAC3 /1,2,3,4,5,6,7,8,9,  1,2,3,4,5,6,7,8,25,
  660. c 729:746
  661. c CU27 <     FACE2               >   <      FACE3          >
  662.      . 13,20,19,18,17,16,15,14,26,   1,9,13,14,15,10,3,2,21,
  663. c 747:764
  664. c CU27 <     FACE4            >   <      FACE5           >
  665.      . 3,10,15,16,17,11,5,4,22,   7,6,5,11,17,18,19,12,23,
  666. c 765:787
  667. c CU27 <      FACE6          >  PR21 <FACE1   >  <      FACE2        >
  668.      . 1,8,7,12,19,20,13,9,24,  1,2,3,4,5,6,19,  10,15,14,13,12,11,20,
  669. c 788:805
  670. c PR21 <     FACE3          >  <      FACE4         >
  671.      . 1,7,10,11,12,8,3,2,16,  3,8,12,13,14,9,5,4,17,
  672. c 806:821
  673. c PR21 <      FACE5         >  TE15  <  FAC1 >
  674.      . 5,9,14,15,10,7,1,6,18,  1,2,3,4,5,6,11,
  675. c 822:842
  676. c TE15 <  FAC2        > <   FAC3       >  <    FAC4      >
  677.      . 1,7,10,8,3,2,12, 3,8,10,9,5,4,13,  5,9,10,7,1,6,14,
  678. c 843:865
  679. c PY19 <FAC1             > <    FAC2       > <    FAC3        >
  680.      . 1,2,3,4,5,6,7,8,14, 1,9,13,10,3,2,15, 3,10,13,11,5,4,16,
  681. c 866:879
  682. c PY19 <    FAC4        > <    FAC5       >
  683.      . 5,11,13,12,7,6,17, 7,12,13,9,1,8,18,
  684. c 880:1000
  685.      . 121*0 /
  686. cbp2016C  KSIF :  un point interne a un element doit voir les faces de cet
  687. cbp2016C  element sous un angle solide positif
  688. cbp2016      DATA KSIF / -1, 1, 1, 1,-1,-1,-1, 1, 1, 1, 1 /
  689. C
  690. C Tableau des noeuds (de l'element initial) auxquels il faudra
  691. C affecter les coefficients des fonctions de forme (XCOEEF) pour
  692. C obtenir les coordonnes des nouveaux noeuds
  693. C Nb valeurs = 21x52 + 4 = 1096
  694.       DATA LISNO1 /
  695. c       <-----TRI3-----><-(7)-----------------------TRI6-----------(21)
  696.      .  1, 2, 2, 3, 1, 3, 2, 1, 3, 2, 3, 1, 4, 3, 5, 4, 5, 3, 6, 5, 1,
  697. c      -------------------------------------------------------------->
  698.      .  6, 1, 5, 2, 6, 4, 3, 5, 1, 2, 4, 6, 1, 5, 3, 4, 6, 2, 1, 3, 5,
  699. c      <-----------QUA4-------------------><----------------QUA8---(63)
  700.      .  1, 2, 2, 3, 3, 4, 1, 4, 1, 2, 3, 4, 2, 1, 3, 2, 3, 1, 4, 3, 5,
  701. c      ------------------------------------------------------------(84)
  702.      .  4, 5, 3, 6, 5, 7, 6, 7, 5, 8, 7, 1, 8, 1, 7, 2, 4, 8, 6, 1, 3,
  703. c      -----------------------------------------------------------(105)
  704.      .  5, 7, 2, 4, 6, 8, 1, 3, 5, 7, 8, 2, 6, 4, 1, 3, 5, 7, 4, 2, 6,
  705. c      --------------------------------------><-------------------(126)
  706.      .  8, 1, 3, 5, 7, 6, 4, 8, 2, 1, 3, 5, 7, 1, 2, 1, 2, 3, 4, 1, 4,
  707.      .  1, 5, 1, 2, 5, 6, 1, 2, 3, 4, 5, 6, 7, 8, 1, 4, 5, 8, 2, 3, 2,
  708.      .  6, 2, 3, 6, 7, 3, 4, 3, 7, 3, 4, 7, 8, 4, 8, 5, 6, 5, 6, 7, 8,
  709.      .  5, 8, 6, 7, 7, 8, 2, 1, 3, 2, 4, 8, 6, 1, 3, 5, 7, 2, 4, 6, 8,
  710.      .  1, 3, 5, 7, 8, 6, 2, 4, 1, 3, 5, 7, 8, 1, 7, 9, 1,13, 2, 9,10,
  711.      . 14, 1, 3,13,15, 2, 4, 6, 8, 9,10,11,12,14,16,18,20,13,15,17,19,
  712.      .  1, 3, 5, 7, 8, 9,12,20, 1, 7,13,19, 9, 2,14,10, 1, 3,13,15, 2,
  713.      .  9,10,14, 1, 3,13,15, 2, 9,10,14, 4, 8,16,20, 6,11,12,18, 5, 7,
  714.      . 17,19, 1, 3,13,15, 2, 4, 6, 8, 9,10,11,12,14,16,18,20, 1, 3, 5,
  715.      .  7,13,15,17,19, 8, 9,12,20, 2, 6,14,18, 4,10,11,16, 3, 5,15,17,
  716.      .  1, 7,13,19, 8, 9,12,20, 1, 7,13,19, 9, 8,20,12, 1, 7,13,19, 2,
  717.      .  3, 1, 4, 3, 5, 4, 2, 6, 8, 1, 3, 5, 7,10, 3,15, 4,10,11,16, 3,
  718.      .  5,15,17,10, 2,14, 9, 1, 3,13,15,10, 4,16,11, 5, 3,15,17, 4,10 /
  719. C
  720.       DATA LISNO2 /
  721.      . 11,16, 3, 5,15,17, 4,10,11,16, 2, 6,14,18, 8, 9,12,20, 1, 7,13,
  722.      . 19, 3, 5,15,17, 4, 5, 3, 6, 5, 7, 6, 4, 8, 2, 1, 3, 5, 7,11, 5,
  723.      . 17, 6,11,12,18, 5, 7,17,19,11, 4,16,10, 3, 5,15,17,11, 6,18,12,
  724.      .  5, 7,17,19, 6,11,12,18, 5, 7,17,19, 6,11,12,18, 4, 8,16,20, 2,
  725.      .  9,10,14, 1, 3,13,15, 5, 7,17,19, 6, 7, 5, 8, 7, 1,12, 7,19,12,
  726.      .  6,18,11, 5, 7,17,19,12, 8,20, 9, 1, 7,13,19, 9,13, 1,14, 9,10,
  727.      .  2, 1, 3,13,15,14,16,18,20, 9,10,11,12, 2, 4, 6, 8, 1, 3, 5, 7,
  728.      . 13,15,17,19,20, 9,12, 8, 1, 7,13,19,14,13,15,14,16,20,18,13,15,
  729.      . 17,19,14,16,18,20,13,15,17,19,20,14,18,16,13,15,17,19,20,13,19,
  730.      . 10,15, 3,16,10,11, 4, 3, 5,15,17,14,15,13,16,15,17,16,14,18,20,
  731.      . 13,15,17,19,11,17, 5,18,11,12, 6, 5, 7,17,19,16,17,15,18,17,19,
  732.      . 18,16,20,14,13,15,17,19,12,19, 7,18,19,17,20,19,13, 1, 2, 2, 3,
  733.      .  1, 3, 4, 5, 5, 6, 4, 6, 2, 1, 3, 2, 3, 1, 4, 3, 5, 4, 5, 3, 6,
  734.      .  5, 1, 6, 1, 5, 2, 6, 4, 3, 5, 1, 2, 4, 6, 1, 5, 3, 4, 6, 2, 1,
  735.      .  3, 5, 2, 7, 8,11, 1, 3,10,12, 4, 8, 9,13, 3, 5,12,14, 6, 7, 9,
  736.      . 15, 1, 5,10,14,11,10,12,11,12,10,13,12,14,13,14,12,15,14,10,15,
  737.      . 10,14,11,15,13,12,14,10,11,13,15,10,14,12,15,13,11,12,10,14, 1,
  738.      .  2, 2, 3, 1, 3, 1, 4, 2, 4, 3, 4, 2, 1, 3, 2, 6, 4, 3, 5, 1, 6 /
  739. C
  740.       DATA LISNO3 /
  741.      .  1, 5, 7, 1,10, 2, 7, 8, 3,10, 1, 6, 7, 9, 5,10, 1, 2, 3, 1, 4,
  742.      .  3, 5, 2, 4, 6, 1, 5, 3, 2, 8, 7, 1,10, 3, 8, 3,10, 4, 8, 9, 5,
  743.      . 10, 3, 4, 5, 3, 6, 5, 1, 4, 6, 2, 1, 3, 5, 4, 9, 8, 3,10, 5, 9,
  744.      .  5,10, 6, 9, 7, 1,10, 3, 7, 8, 2, 1, 3,10, 8, 9, 4, 3, 5,10, 7,
  745.      .  9, 6, 1, 5,10, 7,10, 1, 8,10, 3, 9,10, 5, 2, 4, 6, 7, 8, 9, 1,
  746.      .  3, 5,10, 1, 2, 2, 3, 3, 4, 1, 4, 1, 2, 3, 4, 1, 5, 2, 5, 3, 5,
  747.      .  4, 5, 2, 1, 3, 2, 4, 8, 6, 1, 3, 5, 7, 2, 4, 6, 8, 1, 3, 5, 7,
  748.      .  8, 6, 2, 4, 1, 3, 5, 7, 8, 1, 7, 9, 1,13, 2, 9,10, 3,13, 1, 2,
  749.      .  8, 9, 4, 6,10,12,11, 5,13, 1, 3, 7, 8, 9,12, 7,13, 1, 2, 3, 1,
  750.      .  4, 3, 5, 4, 2, 6, 8, 1, 3, 5, 7, 2,10, 9, 1,13, 3,10, 3,13, 4,
  751.      . 10,11, 5,13, 3, 2, 4,10, 6, 8, 9,11,12, 7,13, 1, 3, 5, 4, 5, 3,
  752.      .  6, 5, 7, 6, 8, 4, 2, 1, 3, 5, 7, 4, 6,11, 2, 8,10,12, 9, 1,13,
  753.      .  3, 5, 7, 4,11,10, 3,13, 5,11, 5,13, 6,11,12, 7,13, 5, 6, 7, 5,
  754.      .  8, 7, 1, 8,12, 9, 1,13, 7, 6, 8,12, 2, 4, 9,11,10, 3,13, 1, 5,
  755.      .  7, 6,12,11, 5,13, 7,12, 7,13, 9,10, 2, 1, 3,13,10,11, 4, 3, 5,
  756.      . 13,11,12, 6, 5, 7,13, 9,12, 8, 1, 7,13, 9,13, 1,10,13, 3,11,13,
  757.      .  5,12,13, 7 /
  758. C
  759. C LPOS1 : 1) Position dans le tableau XCOEFF
  760. C         2) Nombre d'elements de la liste pointee par 1)
  761. C Nb valeurs = 21x20 + 8 = 428
  762.       DATA LPOS11 /
  763. c      <------TRI3------><-----------------------TRI6-------------(21)
  764.      .  1, 2, 1, 2, 1, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,10, 6,10,
  765. c      -------------------------------------->
  766. c               <------------QUA4---(32)-----><------QUA8---------(42)
  767.      .  6,10, 6, 1, 2, 1, 2, 1, 2, 1, 2, 6, 4, 3, 3, 3, 3, 3, 3, 3, 3,
  768. c      -----------------------(51)------------------------->      (63)
  769.      .  3, 3, 3, 3, 3, 3, 3, 3,16, 8,24, 8,16, 8,16, 8,16, 8, 1, 2, 6,
  770.      .  4, 1, 2, 1, 2, 6, 4,32, 8, 6, 4, 1, 2, 1, 2, 6, 4, 1, 2, 1, 2,
  771.      .  6, 4, 1, 2, 1, 2, 6, 4, 1, 2, 1, 2, 1, 2, 3, 3,16, 8,24, 8,16,
  772.      .  8, 3, 3, 3, 3,16, 8,63,20,16, 8,16, 8,24, 8,63,20,83,20,63,20,
  773.      . 24, 8,16, 8, 3, 3, 3, 3,16, 8, 3, 3,16, 8,16, 8,16, 8,24, 8,63,
  774.      . 20, 3, 3, 3, 3,16, 8, 3, 3,16, 8,16, 8,16, 8,24, 8,63,20, 3, 3,
  775.      .  3, 3, 3, 3,16, 8,16, 8, 3, 3,16, 8,63,20,16, 8, 3, 3,16, 8,24,
  776.      .  8,16, 8, 3, 3, 3, 3,16, 8, 3, 3, 3, 3,16, 8, 3, 3,16, 8, 3, 3 /
  777.       DATA LPOS12 /
  778.      .  3, 3,16, 8, 3, 3, 3, 3, 3, 3, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1,
  779.      .  2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,10, 6,10, 6,10, 6,24, 8,
  780.      . 24, 8,24, 8, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,10, 6,10, 6,10,
  781.      .  6, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 3, 3,10, 6, 3, 3, 3, 3,
  782.      . 10, 6,10, 6, 3, 3, 3, 3,10, 6,10, 6, 3, 3,10, 6, 3, 3, 3, 3,10,
  783.      .  6,10, 6, 3, 3,10, 6,10, 6,10, 6,10, 6, 3, 3, 3, 3, 3, 3,40,10,
  784.      .  1, 2, 1, 2, 1, 2, 1, 2, 6, 4, 1, 2, 1, 2, 1, 2, 1, 2, 3, 3,16,
  785.      .  8,24, 8,16, 8, 3, 3, 3, 3,10, 6,50,13,10, 6, 3, 3, 3, 3,16, 8,
  786.      . 10, 6, 3, 3,10, 6,50,13, 3, 3, 3, 3,16, 8,50,13,10, 6, 3, 3,10,
  787.      .  6, 3, 3, 3, 3,10, 6,50,13,10, 6, 3, 3,10, 6,10, 6,10, 6,10, 6,
  788.      .  3, 3, 3, 3, 3, 3, 3, 3 /
  789. C
  790. C Position de l'element dans LPOS1
  791.       DATA LPOS2 /
  792. c          TRI3     TRI6     QUA4     QUA8     CUB8 CU20 PRI6 PR15
  793.      . 3*0,  1 ,1*0, 4, 1*0,  13, 1*0, 18, 3*0, 31, 50,  111, 117, 5*0,
  794. c      TET4 TE10 PYR5 PY13  etc...
  795.      . 138, 144, 169, 178,  74*0 /
  796. C
  797. C Position de l'element dans LISNOE
  798.       DATA LPOS3 /
  799. c          TRI3     TRI6     QUA4     QUA8     CUB8 CU20 PRI6 PR15
  800.      . 3*0, 1, 1*0,  7, 1*0, 43, 1*0, 55, 3*0, 119, 175, 627, 639, 5*0,
  801. c      TET4 TE10 PYR5 PY13  etc...
  802.      . 735, 747, 865, 885,  74*0 /
  803. C
  804. C Liste des coefficients a affecter aux noeuds (valeurs obtenues par
  805. C le calcul des fonctions de forme, selon le type de noeud a creer)
  806.       DATA XCOEFF /
  807. C Nb valeurs = 7x14 + 4 = 102
  808.      +      0.5,     0.5,    0.75,   0.375,  -0.125,    0.25,    0.25,
  809.      +     0.25,    0.25,     0.5,     0.5,    0.25,  -0.125,  -0.125,
  810.      +        0,    0.75,   0.375,   0.375,    0.25, -0.1875, -0.1875,
  811.      +  -0.1875, -0.1875,     0.5,     0.5,     0.5,     0.5,   -0.25,
  812.      +    -0.25,   -0.25,   -0.25,   0.125,   0.125,   0.125,   0.125,
  813.      +    0.125,   0.125,   0.125,   0.125,    0.25,    0.25,    0.25,
  814.      +     0.25,    0.25,    0.25,  -0.125,  -0.125,  -0.125,  -0.125,
  815.      + UNTIER  , UNTIER , UNTIER , UNSIXI , UNSIXI , UNSIXI , UNSIXI ,
  816.      +  UNDOUZ ,  -0.125,  -0.125, UNSIXM , UNSIXM , UNSIXM ,   0.375,
  817.      +    0.375,   0.375,   0.375,  0.1875,  0.1875,  0.1875,  0.1875,
  818.      +    0.125,   0.125,   0.125,   0.125,-0.15625,-0.15625,-0.15625,
  819.      + -0.15625,-0.28125,-0.28125,-0.28125,-0.28125,    0.25,    0.25,
  820.      +     0.25,    0.25,    0.25,    0.25,    0.25,    0.25,    0.25,
  821.      +     0.25,    0.25,    0.25,   -0.25,   -0.25,   -0.25,   -0.25,
  822.      +    -0.25,   -0.25,   -0.25,   -0.25 /
  823. C
  824. C Nombre de points a creer pour raffiner l'element
  825.       DATA NBRAF /
  826. c           TRI3     TRI6     QUA4     QUA8     CUB8 CU20 PRI6 PR15
  827.      .  3*0, 3  ,1*0, 9  ,1*0, 5  ,1*0, 13 ,3*0, 19 , 61 , 6,   21,
  828. c           TET4 TE10 PYR5 PY13  etc...
  829.      .  5*0, 6  , 25,  9,   37,  74*0 /
  830. C
  831. C LNELM : 1) nombre d'elements a creer
  832. C         2) type d'elements a creer
  833. C Indexe par (ITYPEL-1)*2+1
  834. C Nb valeurs = 108
  835.       DATA LNELM /
  836.      .  12*0,4,4,6*0,4,6,6*0,4,8,6*0,4,10,14*0,8,14,2*0,8,15,2*0,4,16,
  837.      .  2*0,4,17,22*0,8,23,2*0,8,24,2*0,6,25,4,23,6,26,4,24,4*0 /
  838. C
  839. C Liste des noeuds formant les elements a creer
  840. C Rq: pour les tetraedres il y a deux types de decomposition selon
  841. C la longueur des diagonales 7-4 et 9-2
  842. C LIELM1 : TRI3-TRI6-QUA4-QUA8-CU20-TET4 (21x16+20=356 valeurs)
  843. C LIELM2 : TE10-PRI6-PR15-PYR5-PY13 (21x17+11=368 valeurs)
  844.       DATA LIELM1 /
  845.      .  1, 4, 6, 2, 5, 4, 3, 6, 5, 4, 5, 6, 1, 7, 2,13, 6,12, 2, 8, 3,
  846.      .  9, 4,14, 4,10, 5,11, 6,15, 2,14, 4,15, 6,13, 1, 5, 9, 8, 2, 6,
  847.      .  9, 5, 3, 7, 9, 6, 4, 8, 9, 7, 1, 9, 2,17,18,19, 8,16, 2,10, 3,
  848.      . 11, 4,20,18,17, 4,12, 5,13, 6,21,18,20, 6,14, 7,15, 8,19,18,21,
  849.      .  1, 9,10,11,12,13,14,15, 9, 2,16,10,13,17,18,14,10,16, 3,19,14,
  850.      . 18,20,21,11,10,19, 4,15,14,21,22,12,13,14,15, 5,23,24,25,13,17,
  851.      . 18,14,23, 6,26,24,14,18,20,21,24,26, 7,27,15,14,21,22,25,24,27,
  852.      .  8, 1,21, 2,22,23,24, 8,25,26,27,28,29, 9,30,31,32,33,34,35,36,
  853.      .  2,37, 3,38, 4,39,23,22,27,40,41,28,31,42,10,43,44,45,33,32,23,
  854.      . 39, 4,46, 5,47, 6,48,28,41,49,50,33,45,44,51,11,52,53,54, 8,24,
  855.      . 23,48, 6,55, 7,56,29,28,50,57,35,34,33,54,53,58,12,59, 9,30,31,
  856.      . 32,33,34,35,36,60,61,62,63,13,64,14,65,66,67,20,68,31,42,10,43,
  857.      . 44,45,33,32,61,69,70,62,14,71,15,72,16,73,66,65,33,45,44,51,11,
  858.      . 52,53,54,62,70,74,75,66,73,16,76,17,77,18,78,35,34,33,54,53,58,
  859.      . 12,59,63,62,75,79,20,67,66,78,18,80,19,81, 1, 5, 7, 8, 5, 2, 6,
  860.      .  9, 7, 6, 3,10, 8, 9,10, 4, 7, 6, 9, 5, 7, 9, 8, 5, 8, 9, 7,10,
  861.      .  9, 6, 7,10, 5, 9, 7, 8, 9,10, 7, 8,10, 9, 7, 6, 9, 5, 7, 6 /
  862.  
  863.       DATA LIELM2 /
  864.      .  1,11, 2,12, 6,13,14,15,16, 7, 2,17, 3,18, 4,19,20,21,22, 8, 4,
  865.      . 23, 5,24, 6,25,26,27,28, 9, 7,29, 8,30, 9,31,32,33,34,10, 8,22,
  866.      .  4,25, 6,35,30,26,28, 9, 7,29, 8,35, 6,16,31,30,28, 9, 6,35, 8,
  867.      . 29, 7,16,12,20,15, 2, 6,25, 4,22, 8,35,12,19,20, 2, 8,30, 9,28,
  868.      .  6,35,29,31,16, 7, 2,20, 8,35, 6,12,15,29,16, 7, 9,30, 8,35, 6,
  869.      . 28,26,22,25, 4, 8,20, 2,12, 6,35,22,19,25, 4, 1, 7, 9, 4,10,12,
  870.      .  7, 2, 8,10, 5,11, 8, 3, 9,11, 6,12, 7, 8, 9,10,11,12, 1,16, 2,
  871.      . 22, 6,21, 7,25,27,10,28,11,34,15,33, 2,17, 3,18, 4,23,25, 8,26,
  872.      . 11,29,12,30,13,35, 6,24, 4,19, 5,20,27,26, 9,15,36,13,31,14,32,
  873.      .  2,23, 4,24, 6,22,25,26,27,11,35,13,36,15,34, 1, 6,10, 9,11, 6,
  874.      .  2, 7,10,12, 7, 3, 8,10,13, 9,10, 8, 4,14,11,12,13,14, 5,14,13,
  875.      . 12,11,10,11,10,12, 6,12,10,13, 7,13,10,14, 8,14,10,11, 9, 1,14,
  876.      .  2,15,16,17, 8,18,19,20,21,22, 9, 2,23, 3,24, 4,25,16,15,26,27,
  877.      . 28,29,10,16,25, 4,30, 5,31, 6,32,33,34,35,36,11, 8,17,16,32, 6,
  878.      . 37, 7,38,39,40,41,42,12, 9,43,10,44,11,45,12,46,47,48,49,50,13,
  879.      .  9,46,12,45,11,44,10,43,21,40,33,29,16, 9,21,16,29,10,43,20,15,
  880.      . 26, 2,10,29,16,33,11,44,28,25,34, 4,11,33,16,40,12,45,36,32,41,
  881.      . 6,12,40,16,21, 9,46,39,17,22, 8 /
  882. C
  883. C Position de l'element dans LIELM
  884.  
  885. C      DATA XXXX / 3*0,TRI3,1*0,TRI6,1*0,QUA4,1*0,QUA8,3*0,CUB8,CU20,
  886. C     . PRI6,PR15,5*0,TET4,TE10,PYR5,PY13,74*0 /
  887.       DATA LPOS4 / 3*0,1,1*0,13,1*0,37,1*0,53,3*0,85,149,477,501,5*0,
  888.      . 309,357,561,607,74*0 /
  889. C
  890. C Nombre de noeuds a creer internes a l'element
  891.       DATA NBINTE /
  892. c         TRI3  TRI6  QUA4  QUA8  CUB8 CU20 PRI6 PR15
  893.      . 3*0,0,1*0,3,1*0,1,1*0,5,3*0,1,   7,   0,   0, 5*0,0,1,0,4,74*0 /
  894. C
  895. C Type du noeud a creer
  896. C Taille : 31 x 6 + 28 = 214
  897.       DATA JTYPNO /
  898.      .  0,0,0,0,0,0,0,0,0,7,7,7,0,0,0,0,7,0,0,0,0,0,0,0,0,7,7,7,7,7,0,
  899.      .  1,0,0,3,7,6,0,0,4,0,0,5,0,0,2,0,0,0,0,1,1,1,0,0,3,7,6,3,3,7,7,
  900.      .  7,6,6,0,0,1,0,4,3,4,4,7,0,0,1,0,5,4,5,5,7,0,0,0,5,6,0,3,7,6,0,
  901.      .  2,2,2,0,0,4,0,0,2,0,5,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,
  902.      .  1,3,4,5,0,0,0,0,0,0,2,2,2,0,0,0,0,0,0,0,1,0,0,2,4,0,0,1,2,0,3,
  903.      .  0,0,1,3,0,4,2,3,4,0,0,0,7,0,0,0,0,1,0,0,0,0,0,1,1,1,0,0,2,7,5,
  904.      .  0,0,1,2,0,3,7,0,0,1,7,3,0,4,0,0,5,7,4,0,2,3,4,5,0,0,0,0 /
  905.  
  906. C Position de l'element dans JTYPNO
  907.       DATA LPOS5 / 3*0,1,1*0,4,1*0,13,1*0,18,3*0,31,50,111,117,5*0,
  908.      . 138,144,169,178,74*0 /
  909. C
  910. C  Valeurs de CCHAMP.INC :
  911. C  -----------------------
  912. C  NOMTP stocke tous les noms d ELEMENTS FINIS possibles
  913. C  LNOMTP : nombre d'elements finis ACTUELLEMENT utilisables
  914. C attention nomtp fonctionne en relation avec NUMMFR et NUMGEO
  915. C  ATTENTION : NOMTP est LIMITE a 300 dans CCHAMP.INC
  916.       DATA LNOMTP / 280 /
  917.  
  918.       DATA NOMTP( 1) / 'CABL' /
  919.       DATA NOMTP( 2) / 'SEG2' /
  920.       DATA NOMTP( 3) / 'SEG3' /
  921.       DATA NOMTP( 4) / 'TRI3' /
  922.       DATA NOMTP( 5) / 'TRI4' /
  923.       DATA NOMTP( 6) / 'TRI6' /
  924.       DATA NOMTP( 7) / 'TRI7' /
  925.       DATA NOMTP( 8) / 'QUA4' /
  926.       DATA NOMTP( 9) / 'QUA5' /
  927.       DATA NOMTP(10) / 'QUA8' /
  928.       DATA NOMTP(11) / 'QUA9' /
  929.       DATA NOMTP(12) / 'RAC2' /
  930.       DATA NOMTP(13) / 'RAC3' /
  931.       DATA NOMTP(14) / 'CUB8' /
  932.       DATA NOMTP(15) / 'CU20' /
  933.       DATA NOMTP(16) / 'PRI6' /
  934.       DATA NOMTP(17) / 'PR15' /
  935.       DATA NOMTP(18) / 'LIA3' /
  936.       DATA NOMTP(19) / 'LIA4' /
  937.       DATA NOMTP(20) / 'LIA6' /
  938.       DATA NOMTP(21) / 'LIA8' /
  939.       DATA NOMTP(22) / 'MULT' /
  940.       DATA NOMTP(23) / 'TET4' /
  941.       DATA NOMTP(24) / 'TE10' /
  942.       DATA NOMTP(25) / 'PYR5' /
  943.       DATA NOMTP(26) / 'PY13' /
  944.       DATA NOMTP(27) / 'COQ3' /
  945.       DATA NOMTP(28) / 'DKT ' /
  946.       DATA NOMTP(29) / 'POUT' /
  947.       DATA NOMTP(30) / 'LISP' /
  948.       DATA NOMTP(31),NOMTP(32) / 'FAC3','FAC4' /
  949.       DATA NOMTP(33),NOMTP(34) / 'FAC6','FAC8' /
  950.       DATA NOMTP(35) / 'LTR3' /
  951.       DATA NOMTP(36) / 'LQU4' /
  952.       DATA NOMTP(37) / 'LCU8' /
  953.       DATA NOMTP(38) / 'LPR6' /
  954.       DATA NOMTP(39) / 'LTE4' /
  955.       DATA NOMTP(40) / 'LPY5' /
  956.       DATA NOMTP(41) / 'COQ8' /
  957.       DATA NOMTP(42) / 'TUYA' /
  958.       DATA NOMTP(43) / 'TUFI' /
  959.       DATA NOMTP(44) / 'COQ2' /
  960.       DATA NOMTP(45) / 'POI1' /
  961.       DATA NOMTP(46) / 'BARR' /
  962.       DATA NOMTP(47) / 'RACO' /
  963.       DATA NOMTP(48) / 'LSU2' /
  964.       DATA NOMTP(49) / 'COQ4' /
  965.       DATA NOMTP(50) / 'LISM' /
  966.       DATA NOMTP(51) / 'COF3' /
  967.       DATA NOMTP(52) / 'RES2' /
  968.       DATA NOMTP(53) / 'LSU3' /
  969.       DATA NOMTP(54) / 'LSU4' /
  970.       DATA NOMTP(55) / 'LICO' /
  971.       DATA NOMTP(56) / 'COQ6' /
  972.       DATA NOMTP(57),NOMTP(58),NOMTP(59) / 'CVS2','CVS3','CVT3' /
  973.       DATA NOMTP(60),NOMTP(61),NOMTP(62) / 'CVT6','CVQ4','CVQ8' /
  974.       DATA NOMTP(63),NOMTP(64),NOMTP(65) / 'THP5','TH13','THP6' /
  975.       DATA NOMTP(66),NOMTP(67),NOMTP(68) / 'TH15','THC8','TH20' /
  976.       DATA NOMTP(69),NOMTP(70),NOMTP(71) / 'ICT3','ICQ4','ICT6' /
  977.       DATA NOMTP(72),NOMTP(73) / 'ICQ8','ICC8' /
  978.       DATA NOMTP(74),NOMTP(75),NOMTP(76) / 'ICT4','ICP6','IC20' /
  979.       DATA NOMTP(77),NOMTP(78) / 'IC10','IC15' /
  980.       DATA NOMTP(79),NOMTP(80),NOMTP(81) / 'TRIP','QUAP','CUBP' /
  981.       DATA NOMTP(82),NOMTP(83),NOMTP(84) / 'TETP','PRIP','TIMO' /
  982.       DATA NOMTP(85),NOMTP(86),NOMTP(87) / 'JOI2','JOI3','JOT3' /
  983.       DATA NOMTP(88),NOMTP(89),NOMTP(90) / 'JOI4','JOI6','JOI8' /
  984.       DATA NOMTP(91) / 'LISC' /
  985.       DATA NOMTP(92) / 'TRIH' /
  986.       DATA NOMTP(93),NOMTP(94) / 'DST ','LIC4' /
  987.       DATA NOMTP(95),NOMTP(96) / 'CERC','TUYO' /
  988.       DATA NOMTP(97),NOMTP(98) / 'LSE2','LITU' /
  989.       DATA NOMTP(99),NOMTP(100) / 'HYT3','HYQ4' /
  990.       DATA NOMTP(101),NOMTP(102),NOMTP(103) / 'HYT4','HYP6','HYC8' /
  991.       DATA NOMTP(104),NOMTP(105),NOMTP(106) / 'TRIS','QUAS','POIS' /
  992.       DATA NOMTP(107),NOMTP(108),NOMTP(109) / 'FRO3','JOP3','JOP6' /
  993.       DATA NOMTP(110) / 'JOP8' /
  994.       DATA NOMTP(111),NOMTP(112),NOMTP(113) / 'POLY','POLY','POLY' /
  995.       DATA NOMTP(114),NOMTP(115),NOMTP(116) / 'POLY','POLY','POLY' /
  996.       DATA NOMTP(117),NOMTP(118),NOMTP(119) / 'POLY','POLY','POLY' /
  997.       DATA NOMTP(120),NOMTP(121),NOMTP(122) / 'POLY','POLY','POLY' /
  998.       DATA NOMTP(123),NOMTP(124),NOMTP(125) / 'BAR3','BAEX','LIA2' /
  999.       DATA NOMTP(126),NOMTP(127) / 'QUAH','CUBH' /
  1000.       DATA NOMTP(128) / 'ROT3' /
  1001.       DATA NOMTP(129),NOMTP(130),NOMTP(131),NOMTP(132),NOMTP(133),
  1002.      .     NOMTP(134),NOMTP(135)
  1003.      .   / 'SEF2','TRF3','QUF4','CUF8','PRF6','TEF4','PYF5' /
  1004.       DATA NOMTP(136),NOMTP(137),NOMTP(138),NOMTP(139),NOMTP(140),
  1005.      .     NOMTP(141),NOMTP(142)
  1006.      .   / 'MSE3','MTR6','MQU9','MC27','MP18','MT10','MP14' /
  1007.       DATA NOMTP(143),NOMTP(144),NOMTP(145),NOMTP(146),NOMTP(147),
  1008.      .     NOMTP(148),NOMTP(149)
  1009.      .   / 'SEF3','TRF7','QUF9','CF27','PF21','TF15','PF19' /
  1010.       DATA NOMTP(150),NOMTP(151),NOMTP(152) / 'SEG6','TR21','QU36' /
  1011.       DATA NOMTP(153),NOMTP(154),NOMTP(155) / 'C216','P126','TE56' /
  1012.       DATA NOMTP(156),NOMTP(157) / 'PY91','TRH6' /
  1013.       DATA NOMTP(158),NOMTP(159),NOMTP(160),NOMTP(161),NOMTP(162),
  1014.      .     NOMTP(163),NOMTP(164),NOMTP(165)
  1015.      .   / 'BSE2','BTR4','BQU5','BCU9','BPR7','BTE5','BPY6','FRO4' /
  1016.       DATA NOMTP(166),NOMTP(167) / 'SEGS','POJS' /
  1017.       DATA NOMTP(168),NOMTP(169),NOMTP(170) / 'JCT3','JCI4','JGI2' /
  1018.       DATA NOMTP(171),NOMTP(172),NOMTP(173) / 'JGT3','JGI4','TRIQ' /
  1019.       DATA NOMTP(174),NOMTP(175),NOMTP(176) / 'QUAQ','CUBQ','TETQ' /
  1020.       DATA NOMTP(177),NOMTP(178),NOMTP(179) / 'PRIQ','TRIR','QUAR' /
  1021.       DATA NOMTP(180),NOMTP(181),NOMTP(182) / 'CUBR','TETR','PRIR' /
  1022.       DATA NOMTP(183),NOMTP(184),NOMTP(185) / 'Q4RI','Q8RI','JOQ3' /
  1023.       DATA NOMTP(186),NOMTP(187),NOMTP(188) / 'JOQ6','JOQ8','JOR3' /
  1024.       DATA NOMTP(189),NOMTP(190) / 'JOR6','JOR8' /
  1025. C  10/2003 : Elements finis (MODE UNIDimensionnel - IDIM=1)
  1026.       DATA NOMTP(191),NOMTP(192) / 'T1D2','T1D3' /
  1027.       DATA NOMTP(193),NOMTP(194) / 'M1D2','M1D3' /
  1028. C  10/05/2004 : Elements finis mécanique fluides
  1029. C               éléments décrivant la dicrétisation du couple
  1030. C               vitesse/pression utilisée
  1031. C     Vitesse      Pression        Nb de pts
  1032. C     L: LINE      C: CENTRE   03: SEG3  21: PR21
  1033. C     B: LINEB     1: CENTREP1 07: TRI7  15: TE15
  1034. C     M: MACRO     S: MSOMMET  09: QUA9  19: PY19
  1035. C     Q: QUAF                  27: CU27
  1036. C     Noms des éléments en LINE/CENTRE pour les éléments géométriques
  1037. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1038.       DATA NOMTP(195),NOMTP(196),NOMTP(197),NOMTP(198),NOMTP(199),
  1039.      .     NOMTP(200),NOMTP(201)
  1040.      .   / 'LC03','LC07','LC09','LC27','LC21','LC15','LC19' /
  1041. C     Noms des éléments en LINE/MSOMMET pour les éléments géométriques
  1042. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1043.       DATA NOMTP(202),NOMTP(203),NOMTP(204),NOMTP(205),NOMTP(206),
  1044.      .     NOMTP(207),NOMTP(208)
  1045.      .   / 'LS03','LS07','LS09','LS27','LS21','LS15','LS19' /
  1046. C     Noms des éléments en LINEB/MSOMMET pour les éléments géométriques
  1047. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1048.       DATA NOMTP(209),NOMTP(210),NOMTP(211),NOMTP(212),NOMTP(213),
  1049.      .     NOMTP(214),NOMTP(215)
  1050.      .   / 'BS03','BS07','BS09','BS27','BS21','BS15','BS19' /
  1051. C     Noms des éléments en MACRO/CENTRE pour les éléments géométriques
  1052. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1053.       DATA NOMTP(216),NOMTP(217),NOMTP(218),NOMTP(219),NOMTP(220),
  1054.      .     NOMTP(221),NOMTP(222)
  1055.      .   / 'MC03','MC07','MC09','MC27','MC21','MC15','MC19' /
  1056. C     Noms des éléments en MACRO/CENTREP1 pour les éléments géométriques
  1057. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1058.       DATA NOMTP(223),NOMTP(224),NOMTP(225),NOMTP(226),NOMTP(227),
  1059.      .     NOMTP(228),NOMTP(229)
  1060.      .   / 'M103','M107','M109','M127','M121','M115','M119' /
  1061. C     Noms des éléments en MACRO/MSOMMET pour les éléments géométriques
  1062. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1063.       DATA NOMTP(230),NOMTP(231),NOMTP(232),NOMTP(233),NOMTP(234),
  1064.      .     NOMTP(235),NOMTP(236)
  1065.      .   / 'MS03','MS07','MS09','MS27','MS21','MS15','MS19' /
  1066. C     Noms des éléments en QUAF/CENTRE pour les éléments géométriques
  1067. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1068.       DATA NOMTP(237),NOMTP(238),NOMTP(239),NOMTP(240),NOMTP(241),
  1069.      .     NOMTP(242),NOMTP(243)
  1070.      .   / 'QC03','QC07','QC09','QC27','QC21','QC15','QC19' /
  1071. C     Noms des éléments en QUAF/CENTREP1 pour les éléments géométriques
  1072. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1073.       DATA NOMTP(244),NOMTP(245),NOMTP(246),NOMTP(247),NOMTP(248),
  1074.      .     NOMTP(249),NOMTP(250)
  1075.      .   / 'Q103','Q107','Q109','Q127','Q121','Q115','Q119' /
  1076. C     Noms des éléments en QUAF/MSOMMET pour les éléments géométriques
  1077. C     SEG3, TRI7, QUA9, CU27, PR21, TE15, PY19
  1078.       DATA NOMTP(251),NOMTP(252),NOMTP(253),NOMTP(254),NOMTP(255),
  1079.      .     NOMTP(256),NOMTP(257),NOMTP(258)
  1080.      .   / 'QS03','QS07','QS09','QS27','QS21','QS15','QS19','CIFL'/
  1081. C    Ajout du type d'element SURE (SUpport de RElations)
  1082. C    Ajout des elements XFEM (mecanique de la rupture)
  1083.       DATA NOMTP(259),NOMTP(260),NOMTP(261),NOMTP(262),NOMTP(263),
  1084.      .     NOMTP(264)
  1085.      . / 'SURE' ,'SHB8','CAF2','CAF3','XQ4R','XC8R'/
  1086. C    Element joint unidimesionnel
  1087.       DATA NOMTP(265) / 'JOI1' /
  1088. C    Elements de zones cohesives
  1089.       DATA NOMTP(266),NOMTP(267),NOMTP(268) / 'ZCO2','ZCO3','ZCO4' /
  1090. C    elements tuy2 et tuy3 pour advection thermique
  1091.       DATA NOMTP(269),nomtp(270)/'TUY2','TUY3'/
  1092. c    element coaxial COS2 (3D pour liaison acier-beton analytique)
  1093.       DATA NOMTP(271)/'COS2'/
  1094. c    element coaxial COA2 (3D pour liaison acier-beton numerique)
  1095.       DATA NOMTP(272)/'COA2'/
  1096. c    element BBAR PYR5 et PY13
  1097.       DATA NOMTP(273),NOMTP(274)/'ICY5','IC13'/
  1098. c    element QUAF 3D
  1099.       DATA NOMTP(275),NOMTP(276),NOMTP(277),NOMTP(278)
  1100.      .     /'CU27','PR21','TE15','PY19'/
  1101. c    element C20R et P15R
  1102.       DATA NOMTP(279),NOMTP(280)/'C20R','P15R'/
  1103.  
  1104.  
  1105. C  Stockage du mot-cle indiquant comment calculer la matrice de RIGIDITE
  1106. C  ATTENTION : NOMIN est DIMENSIONNE a 5 dans CCHAMP.INC.
  1107. C  => Tableau NOMIN non utilise actuellement !
  1108.       DATA LNOMIN / 2 /
  1109.       DATA NOMIN(1) / 'ANALYTIQ' /
  1110.       DATA NOMIN(2) / 'NUMERIQU' /
  1111.  
  1112. C  Stockage des noms de FORMULATIONS possibles
  1113. C  ATTENTION : NOMFR est DIMENSIONNE a 100 dans CCHAMP.INC.
  1114.       DATA LNOMFR / 31 /
  1115.       DATA NOMFR( 1) / 'MASSIF  ' /
  1116.       DATA NOMFR( 2) / 'COQUE   ' /
  1117.       DATA NOMFR( 3) / 'COQUEPAI' /
  1118.       DATA NOMFR( 4) / 'POUTRE  ' /
  1119.       DATA NOMFR( 5) / 'CISATRAN' /
  1120.       DATA NOMFR( 6) / 'LIQUIDE ' /
  1121.       DATA NOMFR( 7) / 'TUYAU   ' /
  1122.       DATA NOMFR( 8) / 'LINESPRI' /
  1123.       DATA NOMFR( 9) / 'TUFISSUR' /
  1124.       DATA NOMFR(10) / 'RACCMASS' /
  1125.       DATA NOMFR(11) / 'RACCCOQU' /
  1126.       DATA NOMFR(12) / 'SURFLIBR' /
  1127.       DATA NOMFR(13) / 'MEMBRANE' /
  1128.       DATA NOMFR(14) / 'UNIAXIAL' /
  1129.       DATA NOMFR(15) / 'THERMIQU' /
  1130.       DATA NOMFR(16) / 'INCOPRES' /
  1131.       DATA NOMFR(17) / 'POREUX  ' /
  1132.       DATA NOMFR(18) / 'RACCORD ' /
  1133.       DATA NOMFR(19) / 'HOMOGENE' /
  1134.       DATA NOMFR(20) / 'TUYACOUS' /
  1135.       DATA NOMFR(21) / 'DARCY   ' /
  1136.       DATA NOMFR(22) / 'FIBRE   ' /
  1137.       DATA NOMFR(23) / 'RAYONNEM' /
  1138.       DATA NOMFR(24) / 'BAREXENT' /
  1139.       DATA NOMFR(25) / 'LIAISON2' /
  1140.       DATA NOMFR(26) / 'RACCORDC' /
  1141.       DATA NOMFR(27) / 'RACCORDG' /
  1142.       DATA NOMFR(28) / 'MACRO_EL' /
  1143.       DATA NOMFR(29) / 'THERMOHY' /
  1144.       DATA NOMFR(30) / 'ELECTROS' /
  1145.       DATA NOMFR(31) / 'DIFFUSIO' /
  1146.  
  1147. C  Stockage du nom des types de CHAMPS PAR ELEMENTS possibles
  1148. C  ATTENTION : NOMCH est DIMENSIONNE a 100 dans CCHAMP.INC.
  1149. C  => Tableau NOMCH non utilise actuellement !
  1150.       DATA LNOMCH / 21 /
  1151.       DATA NOMCH( 1) / 'NOEUD   ' /
  1152.       DATA NOMCH( 2) / 'GRAVITE ' /
  1153.       DATA NOMCH( 3) / 'RIGIDITE' /
  1154.       DATA NOMCH( 4) / 'MASSE   ' /
  1155.       DATA NOMCH( 5) / 'STRESSES' /
  1156.       DATA NOMCH( 6) / 'DEPLACEM' /
  1157.       DATA NOMCH( 7) / 'FORCES  ' /
  1158.       DATA NOMCH( 8) / 'REACTUAL' /
  1159.       DATA NOMCH( 9) / 'FVOLUMIQ' /
  1160.       DATA NOMCH(10) / 'GRADIENT' /
  1161.       DATA NOMCH(11) / 'CONTRAIN' /
  1162.       DATA NOMCH(12) / 'DEFORMAT' /
  1163.       DATA NOMCH(13) / 'MATERIAU' /
  1164.       DATA NOMCH(14) / 'CARACTER' /
  1165.       DATA NOMCH(15) / 'TEMPERAT' /
  1166.       DATA NOMCH(16) / 'PRINCIPA' /
  1167.       DATA NOMCH(17) / 'MAHOOKE ' /
  1168.       DATA NOMCH(18) / 'HOTANGEN' /
  1169.       DATA NOMCH(19) / 'DILATATI' /
  1170.       DATA NOMCH(20) / 'VARINTER' /
  1171.       DATA NOMCH(21) / 'GRAFLEXI' /
  1172.  
  1173. C  Stockage du nom des MATERIAUx possibles
  1174. C  ATTENTION : NOMAT est DIMENSIONNE a 100 dans CCHAMP.INC.
  1175.       DATA LNOMAT / 10 /
  1176.       DATA NOMAT(1) / 'STANDARD' /
  1177.       DATA NOMAT(2) / 'ORTHOTRO' /
  1178.       DATA NOMAT(3) / 'ANISOTRO' /
  1179.       DATA NOMAT(4) / 'COMPOSIT' /
  1180.       DATA NOMAT(5) / 'LIQUIDE ' /
  1181.       DATA NOMAT(6) / 'NONLIN  ' /
  1182.       DATA NOMAT(7) / 'POREUX  ' /
  1183.       DATA NOMAT(8) / 'THERMIQU' /
  1184.       DATA NOMAT(9) / 'HOMOGENE' /
  1185.       DATA NOMAT(10) / 'SCHREFLE' /
  1186.  
  1187. C  Stockage du nom des composantes de DEPLACEMENT (PRIMALES)
  1188. C  ATTENTION : NOMDD est DIMENSIONNE a 100 dans CCHAMP.INC.
  1189.       DATA LNOMDD / 67 /
  1190.       DATA NOMDD( 1),NOMDD( 2),NOMDD( 3) / 'UX  ','UY  ','UZ  ' /
  1191.       DATA NOMDD( 4),NOMDD( 5),NOMDD( 6) / 'RX  ','RY  ','RZ  ' /
  1192.       DATA NOMDD( 7),NOMDD( 8),NOMDD( 9) / 'UR  ','UT  ','RT  ' /
  1193.       DATA NOMDD(10),NOMDD(11),NOMDD(12) / 'LX  ','ALFA','BETA' /
  1194.       DATA NOMDD(13),NOMDD(14),NOMDD(15) / 'FBET','P   ','PI  ' /
  1195.       DATA NOMDD(16),NOMDD(17),NOMDD(18) / 'T   ','RR  ','TINF' /
  1196.       DATA NOMDD(19),NOMDD(20),NOMDD(21) / 'TSUP','TH  ','FC  ' /
  1197.       DATA NOMDD(22),NOMDD(23),NOMDD(24) / 'IUX ','IUY ','IUZ ' /
  1198.       DATA NOMDD(25),NOMDD(26),NOMDD(27) / 'IRX ','IRY ','IRZ ' /
  1199.       DATA NOMDD(28),NOMDD(29),NOMDD(30) / 'IUR ','IUT ','IRT ' /
  1200.       DATA NOMDD(31),NOMDD(32),NOMDD(33) / 'PQ  ','TP  ','IP  ' /
  1201.       DATA NOMDD(34)                     / 'IPI ' /
  1202.       DATA NOMDD(35),NOMDD(36),NOMDD(37) / 'AX  ','AY  ','AZ  ' /
  1203.       DATA NOMDD(38),NOMDD(39),NOMDD(40) / 'B1X ','B1Y ','B1Z ' /
  1204.       DATA NOMDD(41),NOMDD(42),NOMDD(43) / 'C1X ','C1Y ','C1Z ' /
  1205.       DATA NOMDD(44),NOMDD(45),NOMDD(46) / 'D1X ','D1Y ','D1Z ' /
  1206.       DATA NOMDD(47),NOMDD(48),NOMDD(49) / 'E1X ','E1Y ','E1Z ' /
  1207.       DATA NOMDD(50),NOMDD(51),NOMDD(52) / 'B2X ','B2Y ','B2Z ' /
  1208.       DATA NOMDD(53),NOMDD(54),NOMDD(55) / 'C2X ','C2Y ','C2Z ' /
  1209.       DATA NOMDD(56),NOMDD(57),NOMDD(58) / 'D2X ','D2Y ','D2Z ' /
  1210.       DATA NOMDD(59),NOMDD(60),NOMDD(61) / 'E2X ','E2Y ','E2Z ' /
  1211.       DATA NOMDD(62),NOMDD(63)           / 'PG  ','PC  ' /
  1212.       DATA NOMDD(64)                     / 'VEL ' /
  1213.       DATA NOMDD(65),NOMDD(66),NOMDD(67) / 'LX1 ','LX2 ','LX3 ' /
  1214.  
  1215.  
  1216. C  Stockage du nom des composantes des "FORCES" (DUALES)
  1217. C  ATTENTION : NOMDU est DIMENSIONNE a 100 dans CCHAMP.INC.
  1218.       DATA LNOMDU / 67 /
  1219.       DATA NOMDU( 1),NOMDU( 2),NOMDU( 3) / 'FX  ','FY  ','FZ  ' /
  1220.       DATA NOMDU( 4),NOMDU( 5),NOMDU( 6) / 'MX  ','MY  ','MZ  ' /
  1221.       DATA NOMDU( 7),NOMDU( 8),NOMDU( 9) / 'FR  ','FT  ','MT  ' /
  1222.       DATA NOMDU(10),NOMDU(11),NOMDU(12) / 'FLX ','FALF','FBET' /
  1223.       DATA NOMDU(13),NOMDU(14),NOMDU(15) / 'BETA','FP  ','FPI ' /
  1224.       DATA NOMDU(16),NOMDU(17),NOMDU(18) / 'Q   ','MR  ','QINF' /
  1225.       DATA NOMDU(19),NOMDU(20),NOMDU(21) / 'QSUP','FLUX','ED  ' /
  1226.       DATA NOMDU(22),NOMDU(23),NOMDU(24) / 'IFX ','IFY ','IFZ ' /
  1227.       DATA NOMDU(25),NOMDU(26),NOMDU(27) / 'IMX ','IMY ','IMZ ' /
  1228.       DATA NOMDU(28),NOMDU(29),NOMDU(30) / 'IFR ','IFT ','IMT ' /
  1229.       DATA NOMDU(31),NOMDU(32),NOMDU(33) / 'FPQ ','FTP ','IFP ' /
  1230.       DATA NOMDU(34)                     / 'IFPI' /
  1231.       DATA NOMDU(35),NOMDU(36),NOMDU(37) / 'FAX ','FAY ','FAZ ' /
  1232.       DATA NOMDU(38),NOMDU(39),NOMDU(40) / 'FB1X','FB1Y','FB1Z' /
  1233.       DATA NOMDU(41),NOMDU(42),NOMDU(43) / 'FC1X','FC1Y','FC1Z' /
  1234.       DATA NOMDU(44),NOMDU(45),NOMDU(46) / 'FD1X','FD1Y','FD1Z' /
  1235.       DATA NOMDU(47),NOMDU(48),NOMDU(49) / 'FE1X','FE1Y','FE1Z' /
  1236.       DATA NOMDU(50),NOMDU(51),NOMDU(52) / 'FB2X','FB2Y','FB2Z' /
  1237.       DATA NOMDU(53),NOMDU(54),NOMDU(55) / 'FC2X','FC2Y','FC2Z' /
  1238.       DATA NOMDU(56),NOMDU(57),NOMDU(58) / 'FD2X','FD2Y','FD2Z' /
  1239.       DATA NOMDU(59),NOMDU(60),NOMDU(61) / 'FE2X','FE2Y','FE2Z' /
  1240.       DATA NOMDU(62),NOMDU(63)           / 'QG  ','QC  ' /
  1241.       DATA NOMDU(64)                     / 'QEL ' /
  1242.       DATA NOMDU(65),NOMDU(66),NOMDU(67) / 'FLX1','FLX2','FLX3' /
  1243.  
  1244. C  Stockage du nom des modeles de PLASTICITE
  1245. C  ATTENTION : NOMAC est DIMENSIONNE a 100 dans CCHAMP.INC.
  1246. C  => Nom peu explicite !
  1247.       DATA LNOMAC / 40 /
  1248.       DATA NOMAC( 1),NOMAC( 2),NOMAC( 3) / 'PARF','LISP','DRUP' /
  1249.       DATA NOMAC( 4),NOMAC( 5),NOMAC( 6) / 'CINE','ISOT','CINL' /
  1250.       DATA NOMAC( 7),NOMAC( 8),NOMAC( 9) / 'CHA1','EXPL','BETO' /
  1251.       DATA NOMAC(10),NOMAC(11),NOMAC(12) / 'MUCO','CHA2','CHA3' /
  1252.       DATA NOMAC(13),NOMAC(14),NOMAC(15) / 'CHA4','TUFP','DRUC' /
  1253.       DATA NOMAC(16),NOMAC(17),NOMAC(18) / 'GAUV','GUIO','TFEC' /
  1254.       DATA NOMAC(19),NOMAC(20),NOMAC(21) / 'NORT','BLAC','POLY' /
  1255.       DATA NOMAC(22),NOMAC(23),NOMAC(24) / 'RCC1','RCC2','LEMA' /
  1256.       DATA NOMAC(25),NOMAC(26),NOMAC(27) / 'CHAB','ENDO','LISC' /
  1257.       DATA NOMAC(28),NOMAC(29),NOMAC(30) / 'UBIQ','VDOM','MAZA' /
  1258.       DATA NOMAC(31),NOMAC(32),NOMAC(33) / 'UNIL','STEI','CAMC' /
  1259.       DATA NOMAC(34),NOMAC(35),NOMAC(36) / 'MOHR','JOID','BETA' /
  1260.       DATA NOMAC(37),NOMAC(38),NOMAC(39) / 'ROTC','GURS','BETU' /
  1261.       DATA NOMAC(40)                     / 'ACIE' /
  1262.  
  1263. C  Stockage du noms des composantes de VITESSES
  1264. C  Stockage du noms des composantes de VITESSES
  1265. C  ATTENTION : NOMVI est DIMENSIONNE a 100 dans CCHAMP.INC.
  1266. C  => Tableau NOMVI non utilise actuellement !
  1267.       DATA LNOMVI / 10 /
  1268.       DATA NOMVI( 1),NOMVI( 2),NOMVI( 3) / 'VX  ','VY  ','VZ  ' /
  1269.       DATA NOMVI( 4),NOMVI( 5),NOMVI( 6) / 'VMX ','VMY ','VMZ ' /
  1270.       DATA NOMVI( 7),NOMVI( 8),NOMVI( 9) / 'VR  ','VT  ','VMT ' /
  1271.       DATA NOMVI(10)                     / 'VFLX' /
  1272.  
  1273. C  Stockage du nom des compsantes du tenseur des CONTRAINTES
  1274. C  ATTENTION : NOMST est DIMENSIONNE a 100 dans CCHAMP.INC.
  1275. C  => Tableau NOMST non utilise actuellement !
  1276. C    SM = contraintes de MEMBRANE  ,  SF = contraintes de FLEXION
  1277. C    X  Y  Z  elements MASSIFS
  1278. C    S  T  N  elements de PLAQUE
  1279. C    R  Z  T  elements AXISYMETRIQUES
  1280.       DATA LNOMST / 59 /
  1281.       DATA NOMST( 1) / 'SMXX' /
  1282.       DATA NOMST( 2) / 'SMYY' /
  1283.       DATA NOMST( 3) / 'SMXY' /
  1284.       DATA NOMST( 4) / 'SFXX' /
  1285.       DATA NOMST( 5) / 'SFYY' /
  1286.       DATA NOMST( 6) / 'SFXY' /
  1287.       DATA NOMST( 7) / 'SMXZ' /
  1288.       DATA NOMST( 8) / 'SMYZ' /
  1289.       DATA NOMST( 9) / 'SMZZ' /
  1290.       DATA NOMST(10) / 'SFXZ' /
  1291.       DATA NOMST(11) / 'SFYZ' /
  1292.       DATA NOMST(12) / 'SFZZ' /
  1293.       DATA NOMST(13) / 'SMSS' /
  1294.       DATA NOMST(14) / 'SMTT' /
  1295.       DATA NOMST(15) / 'SMST' /
  1296.       DATA NOMST(16) / 'SFSS' /
  1297.       DATA NOMST(17) / 'SFTT' /
  1298.       DATA NOMST(18) / 'SFST' /
  1299.       DATA NOMST(19) / 'SMSN' /
  1300.       DATA NOMST(20) / 'SMTN' /
  1301.       DATA NOMST(21) / 'SMNN' /
  1302.       DATA NOMST(22) / 'SFSN' /
  1303.       DATA NOMST(23) / 'SFTN' /
  1304.       DATA NOMST(24) / 'SFNN' /
  1305.       DATA NOMST(25) / 'SMRR' /
  1306.       DATA NOMST(26) / 'SMZZ' /
  1307.       DATA NOMST(27) / 'SMRZ' /
  1308.       DATA NOMST(28) / 'SFRR' /
  1309.       DATA NOMST(29) / 'SFZZ' /
  1310.       DATA NOMST(30) / 'SFRZ' /
  1311.       DATA NOMST(31) / 'SMRT' /
  1312.       DATA NOMST(32) / 'SMZT' /
  1313.       DATA NOMST(33) / 'SMTT' /
  1314.       DATA NOMST(34) / 'SFRT' /
  1315.       DATA NOMST(35) / 'SFZT' /
  1316.       DATA NOMST(36) / 'SFTT' /
  1317.       DATA NOMST(37) / 'EFFX' /
  1318.       DATA NOMST(38) / 'EFFY' /
  1319.       DATA NOMST(39) / 'EFFZ' /
  1320.       DATA NOMST(40) / 'MOMX' /
  1321.       DATA NOMST(41) / 'MOMY' /
  1322.       DATA NOMST(42) / 'MOMZ' /
  1323.       DATA NOMST(43) / 'N11 ' /
  1324.       DATA NOMST(44) / 'N22 ' /
  1325.       DATA NOMST(45) / 'N12 ' /
  1326.       DATA NOMST(46) / 'M11 ' /
  1327.       DATA NOMST(47) / 'M22 ' /
  1328.       DATA NOMST(48) / 'M12 ' /
  1329.       DATA NOMST(49) / 'NZZ ' /
  1330.       DATA NOMST(50) / 'MZZ ' /
  1331.       DATA NOMST(51) / 'V1  ' /
  1332.       DATA NOMST(52) / 'V2  ' /
  1333.       DATA NOMST(53) / 'NXZ ' /
  1334.       DATA NOMST(54) / 'NYZ ' /
  1335.       DATA NOMST(55) / 'MXX ' /
  1336.       DATA NOMST(56) / 'DELX' /
  1337.       DATA NOMST(57) / 'DELY' /
  1338.       DATA NOMST(58) / 'DELZ' /
  1339.       DATA NOMST(59) / 'DELR' /
  1340.  
  1341. C  Stockage du nom des composantes du tenseur des DEFORMATIONS
  1342. C  ATTENTION : NOMDF est DIMENSIONNE a 100 dans CCHAMP.INC.
  1343. C  => Tableau NOMDF non utilise actuellement !
  1344. C  => LNOMDF est different de LNOMST ????
  1345. C    Meme convention que ci-dessus SAUF que :
  1346. C    SM devient EP  et  SF devient RT
  1347.       DATA LNOMDF / 48 /
  1348.       DATA NOMDF( 1) / 'EPXX' /
  1349.       DATA NOMDF( 2) / 'EPYY' /
  1350.       DATA NOMDF( 3) / 'GAXY' /
  1351.       DATA NOMDF( 4) / 'RTXX' /
  1352.       DATA NOMDF( 5) / 'RTYY' /
  1353.       DATA NOMDF( 6) / 'RTXY' /
  1354.       DATA NOMDF( 7) / 'GAXZ' /
  1355.       DATA NOMDF( 8) / 'GAYZ' /
  1356.       DATA NOMDF( 9) / 'EPZZ' /
  1357.       DATA NOMDF(10) / 'RTXZ' /
  1358.       DATA NOMDF(11) / 'RTYZ' /
  1359.       DATA NOMDF(12) / 'RTZZ' /
  1360.       DATA NOMDF(13) / 'EPSS' /
  1361.       DATA NOMDF(14) / 'EPTT' /
  1362.       DATA NOMDF(15) / 'GAST' /
  1363.       DATA NOMDF(16) / 'RTSS' /
  1364.       DATA NOMDF(17) / 'RTTT' /
  1365.       DATA NOMDF(18) / 'RTST' /
  1366.       DATA NOMDF(19) / 'GASN' /
  1367.       DATA NOMDF(20) / 'GATN' /
  1368.       DATA NOMDF(21) / 'EPNN' /
  1369.       DATA NOMDF(22) / 'RTSN' /
  1370.       DATA NOMDF(23) / 'RTTN' /
  1371.       DATA NOMDF(24) / 'RTNN' /
  1372.       DATA NOMDF(25) / 'EPRR' /
  1373.       DATA NOMDF(26) / 'EPZZ' /
  1374.       DATA NOMDF(27) / 'GARZ' /
  1375.       DATA NOMDF(28) / 'RTRR' /
  1376.       DATA NOMDF(29) / 'RTZZ' /
  1377.       DATA NOMDF(30) / 'RTRZ' /
  1378.       DATA NOMDF(31) / 'GART' /
  1379.       DATA NOMDF(32) / 'GAZT' /
  1380.       DATA NOMDF(33) / 'EPTT' /
  1381.       DATA NOMDF(34) / 'RTRT' /
  1382.       DATA NOMDF(35) / 'RTZT' /
  1383.       DATA NOMDF(36) / 'RTTT' /
  1384.       DATA NOMDF(37) / 'EPS ' /
  1385.       DATA NOMDF(38) / 'GXY ' /
  1386.       DATA NOMDF(39) / 'GXZ ' /
  1387.       DATA NOMDF(40) / 'CX  ' /
  1388.       DATA NOMDF(41) / 'CY  ' /
  1389.       DATA NOMDF(42) / 'CZ  ' /
  1390.       DATA NOMDF(43) / 'GASN' /
  1391.       DATA NOMDF(44) / 'GATN' /
  1392.       DATA NOMDF(45) / 'EELX' /
  1393.       DATA NOMDF(46) / 'EELY' /
  1394.       DATA NOMDF(47) / 'EELZ' /
  1395.       DATA NOMDF(48) / 'EELR' /
  1396.  
  1397. C Stockage du nom des composantes des champs de MATERIAU
  1398. C  ATTENTION : NOMYO est DIMENSIONNE a 200 dans CCHAMP.INC.
  1399. C  => Tableau NOMYO non utilise actuellement !
  1400.       DATA LNOMYO / 115 /
  1401.       DATA NOMYO(1) / 'YOUN' /
  1402.       DATA NOMYO(2) / 'NU  ' /
  1403.       DATA NOMYO(3) / 'RHO ' /
  1404.       DATA NOMYO(4) / 'ALPH' /
  1405.       DATA NOMYO(5) / 'YG1 ' /
  1406.       DATA NOMYO(6) / 'YG2 ' /
  1407.       DATA NOMYO(7) / 'YG3 ' /
  1408.       DATA NOMYO(8) / 'NU12' /
  1409.       DATA NOMYO(9) / 'NU23' /
  1410.       DATA NOMYO(10) / 'NU13' /
  1411.       DATA NOMYO(11) / 'G12 ' /
  1412.       DATA NOMYO(12) / 'G23 ' /
  1413.       DATA NOMYO(13) / 'G13 ' /
  1414.       DATA NOMYO(14) / 'SIGY' /
  1415.       DATA NOMYO(15) / 'H   ' /
  1416.       DATA NOMYO(16) / 'TRAC' /
  1417.       DATA NOMYO(17) / 'CSON' /
  1418.       DATA NOMYO(18) / 'RORF' /
  1419.       DATA NOMYO(19) / 'CREF' /
  1420.       DATA NOMYO(20) / 'LCAR' /
  1421.       DATA NOMYO(21) / 'ALP1' /
  1422.       DATA NOMYO(22) / 'ALP2' /
  1423.       DATA NOMYO(23) / 'ALP3' /
  1424.       DATA NOMYO(24) / 'G   ' /
  1425.       DATA NOMYO(25) / 'V1X ' /
  1426.       DATA NOMYO(26) / 'V1Y ' /
  1427.       DATA NOMYO(27) / 'V1Z ' /
  1428.       DATA NOMYO(28) / 'V2X ' /
  1429.       DATA NOMYO(29) / 'V2Y ' /
  1430.       DATA NOMYO(30) / 'V2Z ' /
  1431.       DATA NOMYO(31) / 'V3X ' /
  1432.       DATA NOMYO(32) / 'V3Y ' /
  1433.       DATA NOMYO(33) / 'V3Z ' /
  1434.       DATA NOMYO(34) / 'K   ' /
  1435.       DATA NOMYO(35) / 'ALF ' /
  1436.       DATA NOMYO(36) / 'J1C ' /
  1437.       DATA NOMYO(37) / 'T   ' /
  1438.       DATA NOMYO(38) / 'ALFA' /
  1439.       DATA NOMYO(39) / 'BETA' /
  1440.       DATA NOMYO(40) / 'GAMM' /
  1441.       DATA NOMYO(41) / 'DELT' /
  1442.       DATA NOMYO(42) / 'ETA ' /
  1443.       DATA NOMYO(43) / 'MU  ' /
  1444.       DATA NOMYO(44) / 'KL  ' /
  1445.       DATA NOMYO(45) / 'PORO' /
  1446.       DATA NOMYO(46) / 'PERM' /
  1447.       DATA NOMYO(47) / 'VISC' /
  1448.       DATA NOMYO(48) / 'COMP' /
  1449.       DATA NOMYO(49) / 'NORM' /
  1450.       DATA NOMYO(50) / 'A   ' /
  1451.       DATA NOMYO(51) / 'C   ' /
  1452.       DATA NOMYO(52) / 'R   ' /
  1453.       DATA NOMYO(53) / 'R0  ' /
  1454.       DATA NOMYO(54) / 'RM  ' /
  1455.       DATA NOMYO(55) / 'B   ' /
  1456.       DATA NOMYO(56) / 'A1  ' /
  1457.       DATA NOMYO(57) / 'C1  ' /
  1458.       DATA NOMYO(58) / 'A2  ' /
  1459.       DATA NOMYO(59) / 'C2  ' /
  1460.       DATA NOMYO(60) / 'PSI ' /
  1461.       DATA NOMYO(61) / 'OMEG' /
  1462.       DATA NOMYO(62) / 'LT1 ' /
  1463.       DATA NOMYO(63) / 'LT2 ' /
  1464.       DATA NOMYO(64) / 'LT3 ' /
  1465.       DATA NOMYO(65) / 'EPT1' /
  1466.       DATA NOMYO(66) / 'EPT2' /
  1467.       DATA NOMYO(67) / 'EPT3' /
  1468.       DATA NOMYO(68) / 'LCS ' /
  1469.       DATA NOMYO(69) / 'LCT ' /
  1470.       DATA NOMYO(70) / 'LCD ' /
  1471.       DATA NOMYO(71) / 'PCT ' /
  1472.       DATA NOMYO(72) / 'PCD ' /
  1473.       DATA NOMYO(73) / 'EPCS' /
  1474.       DATA NOMYO(74) / 'EPCD' /
  1475.       DATA NOMYO(75) / 'LPH ' /
  1476.       DATA NOMYO(76) / 'PENT' /
  1477.       DATA NOMYO(77) / 'BETA' /
  1478.       DATA NOMYO(78) / 'B11 ' /
  1479.       DATA NOMYO(79) / 'B22 ' /
  1480.       DATA NOMYO(80) / 'B12 ' /
  1481.       DATA NOMYO(81) / 'ROF ' /
  1482.       DATA NOMYO(82) / 'ROS ' /
  1483.       DATA NOMYO(83) / 'SIGF' /
  1484.       DATA NOMYO(84) / 'LTR ' /
  1485.       DATA NOMYO(85) / 'EPTR' /
  1486.       DATA NOMYO(86) / 'ENDO' /
  1487.       DATA NOMYO(87) / 'OUV1' /
  1488.       DATA NOMYO(88) / 'OUV2' /
  1489.       DATA NOMYO(89) / 'OUV3' /
  1490.       DATA NOMYO(90) / 'IFIS' /
  1491.       DATA NOMYO(91) / 'ANGL' /
  1492.       DATA NOMYO(92) / 'LBIC' /
  1493.       DATA NOMYO(93) / 'JDA ' /
  1494.       DATA NOMYO(94) / 'EMIS' /
  1495.       DATA NOMYO(95) / 'EINF' /
  1496.       DATA NOMYO(96) / 'ESUP' /
  1497.       DATA NOMYO(97) / 'ETA1' /
  1498.       DATA NOMYO(98) / 'ETA2' /
  1499.       DATA NOMYO(99) / 'RUGO' /
  1500.       DATA NOMYO(100) / 'REC ' /
  1501.       DATA NOMYO(101) / 'FK  ' /
  1502.       DATA NOMYO(102) / 'FA  ' /
  1503.       DATA NOMYO(103) / 'FB  ' /
  1504.       DATA NOMYO(104) / 'FC  ' /
  1505.       DATA NOMYO(105) / 'FD  ' /
  1506.       DATA NOMYO(106) / 'PEL ' /
  1507.       DATA NOMYO(107) / 'PE1 ' /
  1508.       DATA NOMYO(108) / 'PE2 ' /
  1509.       DATA NOMYO(109) / 'PE3 ' /
  1510.       DATA NOMYO(110) / 'PE11' /
  1511.       DATA NOMYO(111) / 'PE22' /
  1512.       DATA NOMYO(112) / 'PE33' /
  1513.       DATA NOMYO(113) / 'PE21' /
  1514.       DATA NOMYO(114) / 'PE31' /
  1515.       DATA NOMYO(115) / 'PE32' /
  1516.  
  1517. C  Stockage du nom des CARACTERISTIQUES geometriques
  1518. C  ATTENTION : NOMCR est DIMENSIONNE a 100 dans CCHAMP.INC.
  1519. C  => Tableau NOMCR non utilise actuellement !
  1520.       DATA LNOMCR / 29 /
  1521.       DATA NOMCR( 1) / 'EPAI' /
  1522.       DATA NOMCR( 2) / 'TORS' /
  1523.       DATA NOMCR( 3) / 'INRY' /
  1524.       DATA NOMCR( 4) / 'INRZ' /
  1525.       DATA NOMCR( 5) / 'SECT' /
  1526.       DATA NOMCR( 6) / 'SECY' /
  1527.       DATA NOMCR( 7) / 'SECZ' /
  1528.       DATA NOMCR( 8) / 'INRZ' /
  1529.       DATA NOMCR( 9) / 'RAYO' /
  1530.       DATA NOMCR(10) / 'VX  ' /
  1531.       DATA NOMCR(11) / 'VY  ' /
  1532.       DATA NOMCR(12) / 'VZ  ' /
  1533.       DATA NOMCR(13) / 'FISS' /
  1534.       DATA NOMCR(14) / 'RACO' /
  1535.       DATA NOMCR(15) / 'VXF ' /
  1536.       DATA NOMCR(16) / 'VYF ' /
  1537.       DATA NOMCR(17) / 'VZF ' /
  1538.       DATA NOMCR(18) / 'ANGL' /
  1539.       DATA NOMCR(19) / 'CALF' /
  1540.       DATA NOMCR(20) / 'DX  ' /
  1541.       DATA NOMCR(21) / 'DY  ' /
  1542.       DATA NOMCR(22) / 'DZ  ' /
  1543.       DATA NOMCR(23) / 'EXCE' /
  1544.       DATA NOMCR(24) / 'SCEL' /
  1545.       DATA NOMCR(25) / 'SFLU' /
  1546.       DATA NOMCR(26) / 'EPS ' /
  1547.       DATA NOMCR(27) / 'PRES' /
  1548.       DATA NOMCR(28) / 'NOF1' /
  1549.       DATA NOMCR(29) / 'NOF2' /
  1550.  
  1551. C  Stockage du nom des coefficients de la matrice de Hooke
  1552. C  ATTENTION : NOMHO est DIMENSIONNE a 100 dans CCHAMP.INC.
  1553. C  => Tableau NOMHO non utilise actuellement !
  1554.       DATA LNOMHO / 73 /
  1555.       DATA NOMHO(1 ),NOMHO(2 ),NOMHO(3 ),NOMHO(4 ),NOMHO(5 ),NOMHO(6 ),
  1556.      .     NOMHO(7 ),NOMHO(8 ),NOMHO(9 ),NOMHO(10),NOMHO(11),NOMHO(12),
  1557.      .     NOMHO(13),NOMHO(14),NOMHO(15),NOMHO(16),NOMHO(17),NOMHO(18),
  1558.      .     NOMHO(19),NOMHO(20),NOMHO(21),NOMHO(22),NOMHO(23),NOMHO(24),
  1559.      .     NOMHO(25),NOMHO(26),NOMHO(27),NOMHO(28),NOMHO(29),NOMHO(30),
  1560.      .     NOMHO(31),NOMHO(32),NOMHO(33),NOMHO(34),NOMHO(35),NOMHO(36)
  1561.      .   / 'HO11','HO21','HO31','HO41','HO51','HO61',
  1562.      .     'HO12','HO22','HO32','HO42','HO52','HO62',
  1563.      .     'HO13','HO23','HO33','HO43','HO53','HO63',
  1564.      .     'HO14','HO24','HO34','HO44','HO54','HO64',
  1565.      .     'HO15','HO25','HO35','HO45','HO55','HO65',
  1566.      .     'HO16','HO26','HO36','HO46','HO56','HO66' /
  1567.       DATA NOMHO(37),NOMHO(38),NOMHO(39) / 'KXX','KYX','KZX' /
  1568.       DATA NOMHO(40),NOMHO(41),NOMHO(42) / 'KXY','KYY','KZY' /
  1569.       DATA NOMHO(43),NOMHO(44),NOMHO(45) / 'KXZ','KYZ','KZZ' /
  1570.       DATA NOMHO(46),NOMHO(47),NOMHO(48),NOMHO(49),NOMHO(50),NOMHO(51),
  1571.      .     NOMHO(52),NOMHO(53),NOMHO(54),NOMHO(55),NOMHO(56),NOMHO(57),
  1572.      .     NOMHO(58),NOMHO(59),NOMHO(60),NOMHO(61),NOMHO(62),NOMHO(63),
  1573.      .     NOMHO(64),NOMHO(65),NOMHO(66),NOMHO(67),NOMHO(68),NOMHO(69),
  1574.      .     NOMHO(70),NOMHO(71),NOMHO(72),NOMHO(73)
  1575.      .   / 'HO71','HO81','HO72','HO82','HO73','HO83','HO74','HO84',
  1576.      .     'HO75','HO85','HO76','HO86',
  1577.      .     'HO17','HO27','HO37','HO47','HO57','HO67','HO77','HO87',
  1578.      .     'HO18','HO28','HO38','HO48','HO58','HO68','HO78','HO88' /
  1579.  
  1580. C  Stockage des noms des VARIABLES INTERNES
  1581. C  ATTENTION : NOMVRI est DIMENSIONNE a 100 dans CCHAMP.INC.
  1582. C  => Tableau NOMVRI non utilise actuellement !
  1583.       DATA LNOVRI / 1 /
  1584.       DATA NOMVRI( 1) / 'EPSE' /
  1585.  
  1586. C  NNAVI stoke tous les couples vitesse/pression utilisable
  1587. C  ATTENTION : NNAVI est DIMENSIONNE a 20 dans CCHAMP.INC.
  1588.       DATA LNNAVI / 9 /
  1589.       DATA NNAVI(1),NNAVI(2),NNAVI(3),NNAVI(4),NNAVI(5),NNAVI(6),
  1590.      .     NNAVI(7),NNAVI(8),NNAVI(9)
  1591.      .   / 'LICE','LIMS','LBMS','MCCE','MCP1','MCMS',
  1592.      .     'QFCE','QFP1','QFMS' /
  1593. C  ELEMENT U/P de Navier-Stokes COURANT
  1594.       DATA ILNAVI / 0 /
  1595.  
  1596. C  Valeurs de CCVQUA4.INC :
  1597. C  ------------------------
  1598. C  Tabulation de SIGMA
  1599.       DATA VSIGMA / 1.D0,1.D0,1.D0,1.D0 /
  1600. C  Tabulation de LAMBDA1
  1601.       DATA VLAMD1 / -1.D0,1.D0,1.D0,-1.D0 /
  1602. C  Tabulation de LAMBDA2
  1603.       DATA VLAMD2 / -1.D0,-1.D0,1.D0,1.D0 /
  1604. C  Tabulation de GAMA
  1605.       DATA VGAMA / 1.D0,-1.D0,1.D0,-1D0 /
  1606. C           TABULATION DE LAMBDA1 * LAMBDA1T
  1607.       DATA VL1L1T / 1.D0,-1.D0,-1.D0, 1.D0,
  1608.      .             -1.D0, 1.D0, 1.D0,-1.D0,
  1609.      .             -1.D0, 1.D0, 1.D0,-1.D0,
  1610.      .              1.D0,-1.D0,-1.D0, 1.D0 /
  1611. C  Tabulation de LAMBDA2 * LAMBDA2T
  1612.       DATA VL2L2T / 1.D0, 1.D0,-1.D0,-1.D0,
  1613.      .              1.D0, 1.D0,-1.D0,-1.D0,
  1614.      .             -1.D0,-1.D0, 1.D0, 1.D0,
  1615.      .             -1.D0,-1.D0, 1.D0, 1.D0 /
  1616. C  Tabulation de GAMA * GAMAT
  1617.       DATA VGGT / 1.D0,-1.D0, 1.D0,-1.D0,
  1618.      .           -1.D0, 1.D0,-1.D0, 1.D0,
  1619.      .            1.D0,-1.D0, 1.D0,-1.D0,
  1620.      .           -1.D0, 1.D0,-1.D0, 1.D0 /
  1621. C  Tabulation de SIGMA * SIGMAT
  1622.       DATA VSST / 1.D0, 1.D0, 1.D0, 1.D0,
  1623.      .            1.D0, 1.D0, 1.D0, 1.D0,
  1624.      .            1.D0, 1.D0, 1.D0, 1.D0,
  1625.      .            1.D0, 1.D0, 1.D0, 1.D0 /
  1626. C  Tabulation de LAMBDA1 GAMAT + GAMA LAMBDA1T / 2
  1627.       DATA VL1GT2 / -1.D0, 1.D0, 0.D0, 0.D0,
  1628.      .               1.D0,-1.D0, 0.D0, 0.D0,
  1629.      .               0.D0, 0.D0, 1.D0,-1.D0,
  1630.      .               0.D0, 0.D0,-1.D0, 1.D0 /
  1631. C  Tabulation de LAMBDA2 GAMAT + GAMA LAMBDA2T / 2
  1632.       DATA VL2GT2 / -1.D0, 0.D0, 0.D0, 1.D0,
  1633.      .               0.D0, 1.D0,-1.D0, 0.D0,
  1634.      .               0.D0,-1.D0, 1.D0, 0.D0,
  1635.      .               1.D0, 0.D0, 0.D0,-1.D0 /
  1636. C  Valeur de CCFXDR
  1637.       DATA IXDRR,IXDRW /0,0/
  1638.  
  1639. C  Valeurs du COMMON CTDRAN
  1640. C     XA=16807.D0,XM=2147483647.D0,XG=123.D0
  1641.       DATA IBUFF/32/,IFLAG/0/
  1642.       DATA XA/16807.d0/,XM/2147483647.D0/,XG/123.D0/,R/32*0.D0/
  1643. c valeurs de CCTRACE
  1644. c par defaut police courier 14 => IOPOTR=4
  1645.       DATA ICOSC/1/,IOPOLI/1/
  1646.       DATA IOPOTR/5/
  1647. *sg Ancienne version ?
  1648. *      DATA ICOTAB/3,1,5,4,6,2,7,11,9,13,12,14,10,15,8,16/
  1649.       DATA ICOTAB /11,3,9,1,13,5,12,4,14,6,10,2,15,7,8/
  1650.       DATA ISOTAB/9,0,0,0,    0,0,0,0,      0,0,0,0,       0,0,0,
  1651.      &            1,15,0,0,   0,0,0,0,      0,0,0,0,       0,0,0,
  1652.      &            1,8,15,0,   0,0,0,0,      0,0,0,0,       0,0,0,
  1653.      &            1,6,11,15,  0,0,0,0,      0,0,0,0,       0,0,0,
  1654.      &            1,5,9,12,  15,0,0,0,      0,0,0,0,       0,0,0,
  1655.      &            1,4,7,10,  13,15,0,0,     0,0,0,0,       0,0,0,
  1656.      &            1,4,7,9,  11,14,15,0,     0,0,0,0,       0,0,0,
  1657.      &            1,3,5,7,   9,12,14,15,    0,0,0,0,       0,0,0,
  1658.      &            1,3,5,7,   9,11,13,14,   15,0,0,0,       0,0,0,
  1659.      &            1,2,4,6,   8,10,12,13,   14,15,0,0,      0,0,0,
  1660.      &            1,3,4,6,   7,9,10,12,   13,14,15,0,      0,0,0,
  1661.      &            1,2,3,5,   6,7,9,11,   12,13,14,15,      0,0,0,
  1662.      &            1,2,4,   5,6,7,      9,10,11,12,  13,14,15,0,0,
  1663.      &            1,2,3,4,   6,7,8,      9,10,11,12,  13,14,15,0,
  1664.      &            1,2,3,4,   5,6,7,      8,9,10,11,  12,13,14,15/
  1665. csg : pour opti isov lign, on évite blan noir gris i.e. 15,7,8
  1666. c   : formule pour générer la nieme ligne du tableau (sauf la 1ere)
  1667. c     cc(i=1,n) = enti(1.5+ (i-1)*(12-1)/n )
  1668.       DATA ISOTA0/7,0,0,0,    0,0,0,0,      0,0,0,0,
  1669.      &            1,12,0,0,   0,0,0,0,      0,0,0,0,
  1670.      &            1,7,12,0,   0,0,0,0,      0,0,0,0,
  1671.      &            1,5,8,12,   0,0,0,0,      0,0,0,0,
  1672.      &            1,4,7, 9,  12,0,0,0,      0,0,0,0,
  1673.      &            1,3,5, 8,  10,12,0,0,     0,0,0,0,
  1674.      &            1,3,5, 7,   8,10,12,0,     0,0,0,0,
  1675.      &            1,3,4,6,   7, 9,10,12,    0,0,0,0,
  1676.      &            1,2,4,5,   7, 8, 9,11,   12,0,0,0,
  1677.      &            1,2,3,5,   6, 7, 8,10,   11,12,0,0,
  1678.      &            1,2,3,4,   5,7, 8, 9,   10,11,12,0,
  1679.      &            1,2,3,4,   5,6,7,8,     9,10,11,12/
  1680.       END
  1681.  
  1682.  
  1683.  
  1684.  
  1685.  
  1686.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales