Télécharger la plaquette Cast3M


fr en

Télécharger idmatr.eso

Retour à la liste

Numérotation des lignes : 

idmatr

  1. C IDMATR    SOURCE    JK148537  25/02/17    21:15:05     12156          
  2.  
  3.       SUBROUTINE IDMATR(MFR,IPMODL,IPNOMC,NBROBL,NBRFAC)
  4.  
  5. *--------------------------------------------------------------------*
  6. *       Noms de composantes de materiaux                             *
  7. *--------------------------------------------------------------------*
  8. *                                                                    *
  9. *  ENTREES:                                                          *
  10. *     MFR      Numero de formulation                                 *
  11. *     IPMODL   objet modele elementaire ( segment actif )            *
  12. *                                                                    *
  13. *  SORTIES:                                                          *
  14. *     IPNOMC   pointeur sur les listes de noms de composantes        *
  15. *              obligatoires et facultatives                          *
  16. *     NBROBL   nombre de composantes obligatoires                    *
  17. *     NBRFAC   nombre de composantes facultatives                    *
  18. *                                                                    *
  19. *--------------------------------------------------------------------*
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8(A-H,O-Z)
  23.  
  24. -INC PPARAM
  25. -INC CCOPTIO
  26. C==DEB= FORMULATION HHO == Include specifique ==========================
  27. -INC CCHHOPA
  28. C==FIN= FORMULATION HHO ================================================
  29.  
  30. -INC SMLMOTS
  31. -INC SMMODEL
  32.       POINTEUR NOMID1.NOMID
  33.  
  34. ******Modif SELLIER ****************************************************
  35. *      augmentation nombre de caracteristiques      
  36.       PARAMETER (ITA=500)
  37.       CHARACTER*16 MOMODL(500)
  38.       integer I,J,K,L
  39. c     nombre de variables et de parametres par variables
  40. c      include './nombre_helmholtz.h' 
  41. -INC HNBRHEL   
  42.       CHARACTER*1 motnl1,motnlv1,motnld1,motnla1
  43.       CHARACTER*2 motnl2
  44.       CHARACTER*3 motnl3 
  45.       CHARACTER*4 motnl4    
  46. ******fin modif SELLIER ************************************************
  47.  
  48.       CHARACTER*(LOCOMP) TABOBL(ITA),TABFAC(ITA)
  49.  
  50.       IMODEL=IPMODL
  51. * Le segment existe-t-il deja?
  52.       IF (lnomid(6).NE.0) THEN
  53.         nomid = lnomid(6)
  54.         SEGACT,nomid
  55.         nbrobl = lesobl(/2)
  56.         nbrfac = lesfac(/2)
  57.         IF (nbrobl+nbrfac.EQ.0) GOTO 765
  58.         SEGINI,nomid1=nomid
  59.         ipnomc=nomid1
  60.         RETURN
  61.       ENDIF
  62.  765  CONTINUE
  63.  
  64.       JGOBL=0
  65.       JGFAC=0
  66.       ipnomc=0
  67.       IRET = 1
  68. *
  69.       NMAT=MATMOD(/2)
  70.       NFOR=FORMOD(/2)
  71.       MELE=NEFMOD
  72.  
  73. C== FORMULATION HHO == MECANIQUE + idem qu'en MASSIF ==> MFR = 1 =======
  74. C== Les composantes specifiques seront ajoutees dans HHOIDC via INOMID.
  75.       IF (MFR.EQ.HHO_MFR_ELEMENT) MFR = 1
  76. C== FORMULATION HHO ====================================================
  77.  
  78. *--------------------------------------------------------------------
  79. *     CAS DE LA FORMULATION THERMOHYDRIQUE
  80. *--------------------------------------------------------------------
  81.       CALL PLACE(FORMOD,NFOR,ithehy,'THERMOHYDRIQUE')
  82.       IF (ithehy.ne.0) then
  83.         if ( matmod(1).eq.'SCHREFLER') then
  84.            JGOBL = 20
  85.            TABOBL(1)='KGG'
  86.            TABOBL(2)='KGC'
  87.            TABOBL(3)='KGT'
  88.            TABOBL(4)='KCG'
  89.            TABOBL(5)='KCC'
  90.            TABOBL(6)='KCT'
  91.            TABOBL(7)='KTG'
  92.            TABOBL(8)='KTC'
  93.            TABOBL(9)='KTT'
  94.            TABOBL(10)='CGG'
  95.            TABOBL(11)='CGC'
  96.            TABOBL(12)='CGT'
  97.            TABOBL(13)='CCG'
  98.            TABOBL(14)='CCC'
  99.            TABOBL(15)='CCT'
  100.            TABOBL(16)='CTG'
  101.            TABOBL(17)='CTC'
  102.            TABOBL(18)='CTT'
  103.            TABOBL(19)='KTGG'
  104.            TABOBL(20)='KTCG'
  105.         else
  106.           IRET = 0
  107.           call erreur(5)
  108.         endif
  109.         GO TO 9999
  110.       ENDIF
  111.  
  112. *--------------------------------------------------------------------
  113. *     CAS DE LA FORMULATION LIQUIDE ET ELEMENT DE RACCORD LITU
  114. *--------------------------------------------------------------------
  115.       CALL PLACE(FORMOD,NFOR,ILIQU,'LIQUIDE')
  116.       IF (ILIQU.NE.0) THEN
  117.         IF (MFR.NE.41.AND.MFR.NE.43) THEN
  118.            JGOBL = 6
  119.            TABOBL(1)='RHO '
  120.            TABOBL(2)='CSON'
  121.            TABOBL(3)='RORF'
  122.            TABOBL(4)='CREF'
  123.            TABOBL(5)='LCAR'
  124.            TABOBL(6)='G   '
  125.  
  126.         ELSEIF (MFR.EQ.41) THEN
  127.            JGOBL = 5
  128.            TABOBL(1)='RHO '
  129.            TABOBL(2)='CSON'
  130.            TABOBL(3)='RORF'
  131.            TABOBL(4)='CREF'
  132.            TABOBL(5)='LCAR'
  133. c*        ELSEIF (MFR.EQ.43) THEN
  134.         ELSE
  135.            JGOBL = 3
  136.            TABOBL(1)='RHO '
  137.            TABOBL(2)='LCAR'
  138.            TABOBL(3)='RORF'
  139.         ENDIF
  140.         GOTO 9999
  141.       ENDIF
  142.  
  143. *--------------------------------------------------------------------
  144. *     CAS DE LA FORMULATION THERMIQUE
  145. *     CAS DE LA FORMULATION DARCY
  146. *--------------------------------------------------------------------
  147.       CALL PLACE(FORMOD,NFOR,ITHER ,'THERMIQUE')
  148.       CALL PLACE(FORMOD,NFOR,IDARCY,'DARCY'    )
  149.       IF (ITHER.NE.0 .OR. IDARCY.NE.0) THEN
  150.  
  151.         IF (ITHER.NE.0) THEN
  152.  
  153. C         Cas particuliers de THERMIQUE CONVECTION :
  154.           CALL PLACE(MATMOD,NMAT,ICONV,'CONVECTION')
  155.           IF (ICONV.NE.0) THEN
  156. C            IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  157. CC             Coques thermiques
  158. C              JGOBL    =2
  159. C              TABOBL(1)='HINF'
  160. C              TABOBL(2)='HSUP'
  161. C              JGFAC    =2
  162. C              TABFAC(1)='TCINF'
  163. C              TABFAC(2)='TCSUP'
  164. C            else
  165. C             Elements massifs
  166.               JGOBL = 1
  167.               TABOBL(1)='H'
  168.               JGFAC=1
  169.               TABFAC(1)='TC'
  170. C            endif
  171.             GOTO 9999
  172.           ENDIF
  173.  
  174. C         Cas particuliers de THERMIQUE RAYONNEMENT :
  175.           CALL PLACE(MATMOD,NMAT,IRAYE,'RAYONNEMENT')
  176.           IF (IRAYE.NE.0) THEN
  177. C           Cas particuliers de THERMIQUE RAYONNEMENT :
  178.             CALL PLACE(MATMOD,NMAT,ICAVE,'CAVITE')
  179.             CALL PLACE(MATMOD,NMAT,IFACA,'FAC_A_FAC')
  180.             CALL PLACE(MATMOD,NMAT,IINFI,'INFINI')
  181.  
  182.             IF (ICAVE.NE.0) THEN
  183. C             RAYONNEMENT en CAVITE :
  184.               IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  185.                 JGOBL=2
  186.                 TABOBL(1)='EINF'
  187.                 TABOBL(2)='ESUP'
  188.               else
  189.                 JGOBL=1
  190.                 TABOBL(1)='EMIS'
  191.               endif
  192.  
  193.               JGFAC=3
  194.               TABFAC(1)='CABS'
  195.               TABFAC(2)='TABS'
  196.               TABFAC(3)='H'
  197.  
  198.             ELSEIF (IFACA.NE.0) THEN
  199. C             RAYONNEMENT FACE_A_FACE :
  200.               JGOBL=1
  201.               TABOBL(1)='EMIS'
  202.               JGFAC = 1
  203.               TABFAC(1)='H'
  204.  
  205.             ELSEIF (IINFI.NE.0) THEN
  206. C             RAYONNEMENT a l'INFINI :
  207.               IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  208.                 JGOBL=3
  209.                 TABOBL(1)='EINF'
  210.                 TABOBL(2)='ESUP'
  211.                 TABOBL(3)='T_IN'
  212.               ELSE
  213.                 JGOBL=2
  214.                 TABOBL(1)='EMIS'
  215.                 TABOBL(2)='T_IN'
  216.               ENDIF
  217.  
  218.               JGFAC = 2
  219.               TABFAC(1)='E_IN'
  220.               TABFAC(2)='H'
  221.  
  222.             ELSE
  223.               IRET = 0
  224.               CALL ERREUR(5)
  225.             ENDIF
  226.             GOTO 9999
  227.           ENDIF
  228.  
  229. C         Cas particuliers de THERMIQUE SOURCE :
  230.           ISRCE = 0
  231.           CALL PLACE(MATMOD,NMAT,ISRCE,'SOURCE')
  232.           IF (ISRCE.NE.0) THEN
  233.             IF (INATUU.EQ.1.AND.IMATEE.EQ.1) THEN
  234.               JGOBL = 1
  235.               TABOBL(1)='QVOL'
  236.               JGFAC = 2
  237.               TABFAC(1)='QINF'
  238.               TABFAC(2)='QSUP'
  239.               GOTO 9999
  240.             ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.1) THEN
  241.               JGOBL = 3
  242.               TABOBL(1)='QTOT'
  243.               TABOBL(2)='ORIG'
  244.               TABOBL(3)='RGAU'
  245.               GOTO 9999
  246.             ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.2) THEN
  247.               JGOBL = 5
  248.               TABOBL(1)='QTOT'
  249.               TABOBL(2)='ORIG'
  250.               TABOBL(3)='RGAU'
  251.               TABOBL(4)='DIRE'
  252.               TABOBL(5)='ZGAU'
  253.               GOTO 9999
  254.             ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.3) THEN
  255.               JGOBL = 7
  256.               TABOBL(1)='QTOT'
  257.               TABOBL(2)='ORIG'
  258.               TABOBL(3)='RGAU'
  259.               TABOBL(4)='DIRE'
  260.               TABOBL(5)='ZGAU'
  261.               TABOBL(6)='DIRL'
  262.               TABOBL(7)='LGAU'
  263.               GOTO 9999
  264.             ELSE
  265. c     write(6,*) 'INATUU, IMATEE =',INATUU, IMATEE
  266.               WRITE(6,*) ' Dans IDMATR : numero IMATEE non prevu'
  267.               CALL ERREUR(21)
  268.               RETURN
  269.             ENDIF
  270.           ENDIF
  271.  
  272.         ENDIF
  273.  
  274. C* Cas THERMIQUE et DARCY
  275.         CALL PLACE(MATMOD,NMAT,IORTH,'ORTHOTROPE')
  276.         CALL PLACE(MATMOD,NMAT,IANIS,'ANISOTROPE')
  277.         IF (IORTH.EQ.0.AND.IANIS.EQ.0)THEN
  278.           JGOBL = 1
  279.           IF(MFR .EQ. 75)THEN
  280. C           Cas des JOI1 (MFR=75) ==> Ressorts THERMIQUES
  281. C           ====================
  282.             TABOBL(1)='KT'
  283.           ELSE
  284.             TABOBL(1)='K   '
  285.           ENDIF
  286.  
  287.         ELSEIF (IORTH.NE.0) THEN
  288. *         ELEMENTS COQUES
  289.           IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  290.             JGOBL=5
  291.             TABOBL(1)='K1  '
  292.             TABOBL(2)='K2  '
  293.             TABOBL(3)='K3  '
  294.             TABOBL(4)='V1X '
  295.             TABOBL(5)='V1Y '
  296.  
  297.           ELSEIF (MFR.EQ.1.OR.MFR.EQ.45) THEN
  298. *              ELEMENTS MASSIFS et HYBRIDES
  299.             IF (IDIM.EQ.2) THEN
  300.               IF (IFOMOD.EQ.1) THEN
  301. *                   ELEMENT MASSIF DE FOURIER
  302.                 JGOBL=5
  303.                 TABOBL(1)='K1  '
  304.                 TABOBL(2)='K2  '
  305.                 TABOBL(3)='V1X '
  306.                 TABOBL(4)='V1Y '
  307.                 TABOBL(5)='K3  '
  308.  
  309.               ELSE
  310.                 JGOBL=4
  311.                 TABOBL(1)='K1  '
  312.                 TABOBL(2)='K2  '
  313.                 TABOBL(3)='V1X '
  314.                 TABOBL(4)='V1Y '
  315.               ENDIF
  316.  
  317.             ELSEIF (IDIM.EQ.3) THEN
  318.               JGOBL=9
  319.               TABOBL(1)='K1  '
  320.               TABOBL(2)='K2  '
  321.               TABOBL(3)='K3  '
  322.               TABOBL(4)='V1X '
  323.               TABOBL(5)='V1Y '
  324.               TABOBL(6)='V1Z '
  325.               TABOBL(7)='V2X '
  326.               TABOBL(8)='V2Y '
  327.               TABOBL(9)='V2Z '
  328.             ENDIF
  329.           ENDIF
  330.  
  331.         ELSEIF (IANIS.NE.0) THEN
  332. *           ELEMENTS MASSIFS
  333.           IF (MFR.EQ.1.OR.MFR.EQ.45) THEN
  334.             IF (IDIM.EQ.2) THEN
  335.               IF (IFOMOD.EQ.1) THEN
  336. *               ELEMENT MASSIF DE FOURIER
  337.                 JGOBL=6
  338.                 TABOBL(1)='K11 '
  339.                 TABOBL(2)='K22 '
  340.                 TABOBL(3)='K21 '
  341.                 TABOBL(4)='V1X '
  342.                 TABOBL(5)='V1Y '
  343.                 TABOBL(6)='K33 '
  344.  
  345.               ELSE
  346.                 JGOBL=5
  347.                 TABOBL(1)='K11 '
  348.                 TABOBL(2)='K22 '
  349.                 TABOBL(3)='K21 '
  350.                 TABOBL(4)='V1X '
  351.                 TABOBL(5)='V1Y '
  352.               ENDIF
  353.  
  354.             ELSEIF (IDIM.EQ.3) THEN
  355.               JGOBL=12
  356.               TABOBL(1)='K11 '
  357.               TABOBL(2)='K22 '
  358.               TABOBL(3)='K33 '
  359.               TABOBL(4)='K21 '
  360.               TABOBL(5)='K31 '
  361.               TABOBL(6)='K32 '
  362.               TABOBL(7)='V1X '
  363.               TABOBL(8)='V1Y '
  364.               TABOBL(9)='V1Z '
  365.               TABOBL(10)='V2X '
  366.               TABOBL(11)='V2Y '
  367.               TABOBL(12)='V2Z '
  368.             ENDIF
  369.           ENDIF
  370.         ENDIF
  371. C*
  372. C* Cas THERMIQUE CONDUCTION, THERMIQUE PHASE ou THERMIQUE ADVECTION :
  373.         IF (ITHER.NE.0) THEN
  374.  
  375.           CALL PLACE(MATMOD,NMAT,IPHA,'PHASE')
  376.           IF (IPHA.NE.0) THEN
  377.             JGOBL0 = JGOBL
  378.             JGOBL= JGOBL0+4
  379.             TABOBL(JGOBL0+1)='RHO '
  380.             TABOBL(JGOBL0+2)='C   '
  381.             TABOBL(JGOBL0+3)='QLAT'
  382.             TABOBL(JGOBL0+4)='TPHA'
  383.  
  384.             JGFAC=1
  385.             TABFAC(1)='H   '
  386.             GOTO 9999
  387.           ENDIF
  388.  
  389.           CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
  390.           if (iadve .ne. 0) then
  391. C           Cas des Tuyaux 1D (MFR=79)
  392. C           ===========================
  393.             IF (MFR .EQ. 79) THEN
  394.               JGFAC   = 3
  395.               TABFAC(1)='RHO '
  396.               TABFAC(2)='C   '
  397.               TABFAC(3)='VITE'
  398.  
  399. C           Cas Massif (MFR=1)
  400. C           ===========================
  401.             ELSEIF (MFR .EQ. 1) THEN
  402. C             Bidimensionnel PLAN DPGE
  403. C             ===========================
  404.               IF (IFOUR.EQ.-3) THEN
  405.                 JGFAC   = 5
  406.                 TABFAC(1)='RHO '
  407.                 TABFAC(2)='C   '
  408.                 TABFAC(3)='VITX'
  409.                 TABFAC(4)='VITY'
  410.                 TABFAC(5)='VITZ'
  411.  
  412. C             Bidimensionnel PLAN (CP/DP)
  413. C             ===========================
  414.               ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
  415.                 JGFAC   = 4
  416.                 TABFAC(1)='RHO '
  417.                 TABFAC(2)='C   '
  418.                 TABFAC(3)='VITX'
  419.                 TABFAC(4)='VITY'
  420.  
  421. C             Axisymetrie
  422. C             ===========================
  423.               ELSEIF (IFOUR .EQ. 0) THEN
  424.                 JGFAC   = 4
  425.                 TABFAC(1)='RHO '
  426.                 TABFAC(2)='C   '
  427.                 TABFAC(3)='VITR'
  428.                 TABFAC(4)='VITZ'
  429.  
  430. C             Fourier
  431. C             ===========================
  432.               ELSEIF (IFOUR .EQ. 1) THEN
  433.                 CALL ERREUR(21)
  434.                 RETURN
  435. C                JGFAC   = 5C
  436. C                TABFAC(1)='RHO '
  437. C                TABFAC(2)='C   '
  438. C                TABFAC(3)='VITR'
  439. C                TABFAC(4)='VITZ'
  440. C                TABFAC(5)='VITT'
  441.  
  442. C             Tridimensionnel
  443. C             ===========================
  444.               ELSEIF (IFOUR .EQ. 2) THEN
  445.                 JGFAC   = 5
  446.                 TABFAC(1)='RHO '
  447.                 TABFAC(2)='C   '
  448.                 TABFAC(3)='VITX'
  449.                 TABFAC(4)='VITY'
  450.                 TABFAC(5)='VITZ'
  451.  
  452. C             Unidimensionnel (1D)
  453. C             ===========================
  454.               ELSEIF (IFOUR.GE.3 .AND. IFOUR.LE.15) THEN
  455.                 IF (IFOUR.LE.6) THEN
  456.                   JGFAC   = 3
  457.                   TABFAC(1)='RHO '
  458.                   TABFAC(2)='C   '
  459.                   TABFAC(3)='VITX'
  460.  
  461.                 ELSEIF (IFOUR.EQ.7 .OR. IFOUR.EQ.8) THEN
  462.                   JGFAC   = 4
  463. C                 Verifier l'utilite des composantes au dela de 4
  464.                   TABFAC(1)='RHO '
  465.                   TABFAC(2)='C   '
  466.                   TABFAC(3)='VITX'
  467.                   TABFAC(4)='VITY'
  468.  
  469.                 ELSEIF (IFOUR.EQ.9 .OR. IFOUR.EQ.10) THEN
  470.                   JGFAC   = 4
  471. C                 Verifier l'utilite des composantes au dela de 4
  472.                   TABFAC(1)='RHO '
  473.                   TABFAC(2)='C   '
  474.                   TABFAC(3)='VITX'
  475.                   TABFAC(4)='VITZ'
  476.  
  477.                 ELSEIF (IFOUR.EQ.11) THEN
  478.                   JGFAC   = 5
  479. C                 Verifier l'utilite des composantes au dela de 4
  480.                   TABFAC(1)='RHO '
  481.                   TABFAC(2)='C   '
  482.                   TABFAC(3)='VITX'
  483.                   TABFAC(4)='VITY'
  484.                   TABFAC(5)='VITZ'
  485.  
  486.                 ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN
  487.                   JGFAC   = 3
  488.                   TABFAC(1)='RHO '
  489.                   TABFAC(2)='C   '
  490.                   TABFAC(3)='VITR'
  491.  
  492.                 ELSEIF (IFOUR.EQ.14) THEN
  493.                   JGFAC   = 4
  494. C                 Verifier l'utilite des composantes au dela de 4
  495.                   TABFAC(1)='RHO '
  496.                   TABFAC(2)='C   '
  497.                   TABFAC(3)='VITR'
  498.                   TABFAC(4)='VITZ'
  499.                 ENDIF
  500.               ENDIF
  501.             ENDIF
  502.             GOTO 9999
  503.           endif
  504.  
  505. C         composantes facultatives THERMIQUE CONDUCTION seule
  506.           IF(MFR .EQ. 75)THEN
  507. C           Cas des JOI1 (MFR=75) ==> Ressorts THERMIQUES
  508. C           ====================
  509.             JGFAC    = 4
  510.             TABFAC(1)='M'
  511.             TABFAC(2)='C'
  512.             TABFAC(3)='TINI'
  513.           ELSE
  514.             JGFAC    = 4
  515.             TABFAC(1)='RHO '
  516.             TABFAC(2)='C   '
  517.             TABFAC(3)='H   '
  518.             TABFAC(4)='TINI'
  519.           ENDIF
  520.         ENDIF
  521.         GOTO 9999
  522.       ENDIF
  523.  
  524. *--------------------------------------------------------------------
  525. *     CAS DE LA FORMULATION MECANIQUE
  526. *--------------------------------------------------------------------
  527.       CALL PLACE(FORMOD,NFOR,IMECA,'MECANIQUE')
  528.       IF (IMECA.NE.0) THEN
  529. *
  530.        CALL MODLIN(MOMODL,NMOD)
  531.        CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
  532.        IF (IRET.EQ.0) GOTO 9999
  533.        IF (NMAT.GE.2) THEN
  534.  
  535.          CALL MODELA(MOMODL,NMOD)
  536.          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
  537.          if ((IPLAC.EQ.9.OR.IPLAC.EQ.10).AND.NMAT.GT.2) then
  538.             do jm = 1,matmod(/2)
  539.               IF (matmod(jm).eq.'IMPEDANCE') THEN
  540.                 INMAT = 3
  541.                 goto 19
  542.               ENDIF
  543.             enddo
  544.          endif
  545.  
  546.          IF (IPLAC.NE.0) THEN
  547.             INMAT=3
  548.             IF (IPLAC.EQ.1) THEN
  549.                INMAT=3
  550.                GOTO 10
  551.             ENDIF
  552.             IF (IPLAC.EQ.2) THEN
  553. *
  554. *              MATERIAU ELASTIQUE ORTHOTROPE
  555. *
  556.                IF (MFR.EQ.75) THEN
  557. *               JOINT UNIDIMENSIONNEL JOI1
  558. *
  559.                 IF(IFOUR.EQ.2)THEN
  560.                   JGOBL=12
  561.                   TABOBL(1)='V1X '
  562.                   TABOBL(2)='V1Y '
  563.                   TABOBL(3)='V1Z '
  564.                   TABOBL(4)='V2X '
  565.                   TABOBL(5)='V2Y '
  566.                   TABOBL(6)='V2Z '
  567.                   TABOBL(7)='KN  '
  568.                   TABOBL(8)='KS1 '
  569.                   TABOBL(9)='KS2 '
  570.                   TABOBL(10)='QN  '
  571.                   TABOBL(11)='QS1 '
  572.                   TABOBL(12)='QS2 '
  573. *
  574.                   JGFAC=10
  575.                   TABFAC(1)='MASS'
  576.                   TABFAC(2)='JX  '
  577.                   TABFAC(3)='JY  '
  578.                   TABFAC(4)='JZ  '
  579.                   TABFAC(5)='ALPN'
  580.                   TABFAC(6)='ALP1'
  581.                   TABFAC(7)='ALP2'
  582.                   TABFAC(8)='ALQN'
  583.                   TABFAC(9)='ALQ1'
  584.                   TABFAC(10)='ALQ2'
  585. *
  586.                ELSEIF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1)THEN
  587.                   JGOBL=5
  588.                   TABOBL(1)='V1X '
  589.                   TABOBL(2)='V1Y '
  590.                   TABOBL(3)='KN  '
  591.                   TABOBL(4)='KS  '
  592.                   TABOBL(5)='QS  '
  593. *
  594.                   JGFAC=6
  595.                   TABFAC(1)='MASS'
  596.                   TABFAC(2)='JZ'
  597.                   TABFAC(4)='ALPN'
  598.                   TABFAC(5)='ALPS'
  599.                   TABFAC(6)='ALQS'
  600.                ENDIF
  601. *
  602.               ELSEIF (MFR.EQ.3) THEN
  603. *                 COQUES MINCES
  604. *
  605.                   JGOBL=6
  606.                   TABOBL(1)='YG1 '
  607.                   TABOBL(2)='YG2 '
  608.                   TABOBL(3)='NU12'
  609.                   TABOBL(4)='G12 '
  610.                   TABOBL(5)='V1X '
  611.                   TABOBL(6)='V1Y '
  612. *
  613.                   IF(IFOUR.EQ.-2) THEN
  614.                     JGFAC=4
  615.                   ELSE
  616.                     JGFAC=3
  617.                   ENDIF
  618.                   TABFAC(1)='ALP1'
  619.                   TABFAC(2)='ALP2'
  620.                   TABFAC(3)='RHO '
  621.                   IF(IFOUR.EQ.-2) TABFAC(4)='DIM3'
  622.  
  623.                ELSEIF (MFR.EQ.9.OR.MFR.EQ.5) THEN
  624. *                 COQUES    AVEC CISAILLEMENT TRANSVERSE
  625. *
  626.                   JGOBL=8
  627.                   TABOBL(1)='YG1 '
  628.                   TABOBL(2)='YG2 '
  629.                   TABOBL(3)='NU12'
  630.                   TABOBL(4)='G12 '
  631.                   TABOBL(5)='G23 '
  632.                   TABOBL(6)='G13 '
  633.                   TABOBL(7)='V1X '
  634.                   TABOBL(8)='V1Y '
  635. *
  636.                   JGFAC=3
  637.                   TABFAC(1)='RHO '
  638.                   TABFAC(2)='ALP1'
  639.                   TABFAC(3)='ALP2'
  640.  
  641.                ELSEIF (MFR.EQ.1.OR.MFR.EQ.31) THEN
  642. *                 ELEMENTS MASSIFS
  643. *
  644.                   IF(IDIM.EQ.3)THEN
  645. *                    ELEMENTS 3D
  646.                      JGOBL=15
  647.                      TABOBL(1)='YG1 '
  648.                      TABOBL(2)='YG2 '
  649.                      TABOBL(3)='YG3 '
  650.                      TABOBL(4)='NU12'
  651.                      TABOBL(5)='NU23'
  652.                      TABOBL(6)='NU13'
  653.                      TABOBL(7)='G12 '
  654.                      TABOBL(8)='G23 '
  655.                      TABOBL(9)='G13 '
  656.                      TABOBL(10)='V1X '
  657.                      TABOBL(11)='V1Y '
  658.                      TABOBL(12)='V1Z '
  659.                      TABOBL(13)='V2X '
  660.                      TABOBL(14)='V2Y '
  661.                      TABOBL(15)='V2Z '
  662. *
  663.                      JGFAC=4
  664.                      TABFAC(1)='RHO '
  665.                      TABFAC(2)='ALP1'
  666.                      TABFAC(3)='ALP2'
  667.                      TABFAC(4)='ALP3'
  668.  
  669.                   ELSEIF (IDIM.EQ.2) THEN
  670.                      IF(IFOUR.EQ.-2)THEN
  671. *                       CONTRAINTE PLANE
  672.                         JGOBL=9
  673.                         TABOBL(1)='YG1 '
  674.                         TABOBL(2)='YG2 '
  675.                         TABOBL(3)='NU12'
  676.                         TABOBL(4)='G12 '
  677.                         TABOBL(5)='V1X '
  678.                         TABOBL(6)='V1Y '
  679.                         TABOBL(7)='YG3 '
  680.                         TABOBL(8)='NU23'
  681.                         TABOBL(9)='NU13'
  682. *
  683.                         JGFAC=4
  684.                         TABFAC(1)='RHO '
  685.                         TABFAC(2)='ALP1'
  686.                         TABFAC(3)='ALP2'
  687.                         TABFAC(4)='DIM3'
  688.  
  689.          ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  690. *                       DEFORMATION PLANE ,AXISYMETRIE
  691.                         JGOBL=9
  692.                         TABOBL(1)='YG1 '
  693.                         TABOBL(2)='YG2 '
  694.                         TABOBL(3)='YG3 '
  695.                         TABOBL(4)='NU12'
  696.                         TABOBL(5)='NU23'
  697.                         TABOBL(6)='NU13'
  698.                         TABOBL(7)='G12 '
  699.                         TABOBL(8)='V1X '
  700.                         TABOBL(9)='V1Y '
  701. *
  702.                         JGFAC=4
  703.                         TABFAC(1)='RHO '
  704.                         TABFAC(2)='ALP1'
  705.                         TABFAC(3)='ALP2'
  706.                         TABFAC(4)='ALP3'
  707.  
  708.                      ELSEIF (IFOUR.EQ.1) THEN
  709. *                       AXISYMETRIE DE FOURIER
  710.                         JGOBL=15
  711.                         TABOBL(1)='YG1 '
  712.                         TABOBL(2)='YG2 '
  713.                         TABOBL(3)='YG3 '
  714.                         TABOBL(4)='NU12'
  715.                         TABOBL(5)='NU23'
  716.                         TABOBL(6)='NU13'
  717.                         TABOBL(7)='G12 '
  718.                         TABOBL(8)='G23 '
  719.                         TABOBL(9)='G13 '
  720.                         TABOBL(10)='V1X '
  721.                         TABOBL(11)='V1Y '
  722.                         TABOBL(12)='V1Z '
  723.                         TABOBL(13)='V2X '
  724.                         TABOBL(14)='V2Y '
  725.                         TABOBL(15)='V2Z '
  726. *
  727.                         JGFAC=4
  728.                         TABFAC(1)='RHO '
  729.                         TABFAC(2)='ALP1'
  730.                         TABFAC(3)='ALP2'
  731.                         TABFAC(4)='ALP3'
  732.                      ENDIF
  733.                   ELSEIF (IDIM.EQ.1) THEN
  734. C= Dans le cas UNID SPHErique, on doit avoir YG2=YG3 et NU12=NU13
  735. C= et dans le cas thermomecanique ALP2=ALP3.
  736.                     JGOBL=6
  737.                     TABOBL(1)='YG1 '
  738.                     TABOBL(2)='YG2 '
  739.                     TABOBL(3)='YG3 '
  740.                     TABOBL(4)='NU12'
  741.                     TABOBL(5)='NU23'
  742.                     TABOBL(6)='NU13'
  743.  
  744.                     IF (IFOUR.EQ.6) THEN
  745.                       JGFAC=2
  746.  
  747.                     ELSEIF (IFOUR.EQ.5 .OR. IFOUR.EQ.10) THEN
  748.                       JGFAC=3
  749.                       TABFAC(3)='ALP3'
  750.                     ELSEIF (IFOUR.EQ.4.OR.IFOUR.EQ.8.OR.
  751.      &                       IFOUR.EQ.13) THEN
  752.                       JGFAC=3
  753.                       TABFAC(3)='ALP2'
  754.                     ELSE
  755.                       JGFAC=4
  756.                       TABFAC(3)='ALP2'
  757.                       TABFAC(4)='ALP3'
  758.                     ENDIF
  759.                     TABFAC(1)='RHO '
  760.                     TABFAC(2)='ALP1'
  761.                   ENDIF
  762.  
  763.                ELSEIF (MFR.EQ.35) THEN
  764. *                 ELEMENTS JOINTS
  765.                   IF (IFOUR.EQ.2) THEN
  766.                       JGOBL=5
  767.                       TABOBL(1)='KS1 '
  768.                       TABOBL(2)='KS2 '
  769.                       TABOBL(3)='KN  '
  770.                       TABOBL(4)='V1X '
  771.                       TABOBL(5)='V1Y '
  772. *
  773.                       JGFAC=2
  774.                       TABFAC(1)='RHO '
  775.                       TABFAC(2)='ALPN'
  776.                   ENDIF
  777.                ENDIF
  778.  
  779.             ELSEIF (IPLAC.EQ.3)THEN
  780. *              MATERIAU ANISOTROPE ELASTIQUE
  781. *
  782.                IF(MFR.EQ.75)THEN
  783. *                 JOINT UNIDIMESIONNEL JOI1
  784. *
  785.                   IF (IFOUR.EQ.2) THEN
  786.                      JGOBL=27
  787.                      TABOBL(1)='V1X '
  788.                      TABOBL(2)='V1Y '
  789.                      TABOBL(3)='V1Z '
  790.                      TABOBL(4)='V2X '
  791.                      TABOBL(5)='V2Y '
  792.                      TABOBL(6)='V2Z '
  793.                      TABOBL(7)='D11 '
  794.                      TABOBL(8)='D22 '
  795.                      TABOBL(9)='D33 '
  796.                      TABOBL(10)='D44 '
  797.                      TABOBL(11)='D55 '
  798.                      TABOBL(12)='D66 '
  799.                      TABOBL(13)='D21 '
  800.                      TABOBL(14)='D31 '
  801.                      TABOBL(15)='D32 '
  802.                      TABOBL(16)='D41 '
  803.                      TABOBL(17)='D42 '
  804.                      TABOBL(18)='D43 '
  805.                      TABOBL(19)='D51 '
  806.                      TABOBL(20)='D52 '
  807.                      TABOBL(21)='D53 '
  808.                      TABOBL(22)='D54 '
  809.                      TABOBL(23)='D61 '
  810.                      TABOBL(24)='D62 '
  811.                      TABOBL(25)='D63 '
  812.                      TABOBL(26)='D64 '
  813.                      TABOBL(27)='D65 '
  814. *
  815.                      JGFAC=10
  816. * MASS: masse totale de l'élément joint
  817.                      TABFAC(1)='MASS'
  818.                      TABFAC(2)='JX  '
  819.                      TABFAC(3)='JY  '
  820.                      TABFAC(4)='JZ  '
  821.                      TABFAC(5)='ALP1'
  822.                      TABFAC(6)='ALP2'
  823.                      TABFAC(7)='ALP3'
  824.                      TABFAC(8)='ALQ1'
  825.                      TABFAC(9)='ALQ2'
  826.                      TABFAC(10)='ALQ3'
  827. *
  828.                   ELSEIF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  829.                      JGOBL=8
  830.                      TABOBL(1)='V1X '
  831.                      TABOBL(2)='V1Y '
  832.                      TABOBL(3)='D11 '
  833.                      TABOBL(4)='D22 '
  834.                      TABOBL(5)='D33 '
  835.                      TABOBL(6)='D21 '
  836.                      TABOBL(7)='D31 '
  837.                      TABOBL(8)='D32 '
  838. *
  839.                      JGFAC=5
  840. * MASS: masse totale de l'élément joint
  841.                      TABFAC(1)='MASS'
  842.                      TABFAC(2)='JZ  '
  843.                      TABFAC(3)='ALP1'
  844.                      TABFAC(4)='ALP2'
  845.                      TABFAC(5)='ALQ3'
  846.                   ENDIF
  847.                ENDIF
  848. *
  849. *                 ELEMENTS MASSIFS
  850. *
  851.                IF(MFR.EQ.1.OR.MFR.EQ.31)THEN
  852.                   IF(IDIM.EQ.3)THEN
  853. *                    ELEMENTS 3D
  854.                      JGOBL=27
  855.                      TABOBL(1)='D11 '
  856.                      TABOBL(2)='D21 '
  857.                      TABOBL(3)='D22 '
  858.                      TABOBL(4)='D31 '
  859.                      TABOBL(5)='D32 '
  860.                      TABOBL(6)='D33 '
  861.                      TABOBL(7)='D41 '
  862.                      TABOBL(8)='D42 '
  863.                      TABOBL(9)='D43 '
  864.                      TABOBL(10)='D44 '
  865.                      TABOBL(11)='D51 '
  866.                      TABOBL(12)='D52 '
  867.                      TABOBL(13)='D53 '
  868.                      TABOBL(14)='D54 '
  869.                      TABOBL(15)='D55 '
  870.                      TABOBL(16)='D61 '
  871.                      TABOBL(17)='D62 '
  872.                      TABOBL(18)='D63 '
  873.                      TABOBL(19)='D64 '
  874.                      TABOBL(20)='D65 '
  875.                      TABOBL(21)='D66 '
  876.                      TABOBL(22)='V1X '
  877.                      TABOBL(23)='V1Y '
  878.                      TABOBL(24)='V1Z '
  879.                      TABOBL(25)='V2X '
  880.                      TABOBL(26)='V2Y '
  881.                      TABOBL(27)='V2Z '
  882. *
  883.                      JGFAC=7
  884.                      TABFAC(1)='RHO '
  885.                      TABFAC(2)='ALP1'
  886.                      TABFAC(3)='ALP2'
  887.                      TABFAC(4)='ALP3'
  888.                      TABFAC(5)='AL12'
  889.                      TABFAC(6)='AL13'
  890.                      TABFAC(7)='AL23'
  891.  
  892.                   ELSEIF (IDIM.EQ.2) THEN
  893.                      IF (IFOUR.EQ.-2) THEN
  894. *                       CONTRAINTE PLANE
  895.                         JGOBL=12
  896.                         TABOBL(1)='D11 '
  897.                         TABOBL(2)='D21 '
  898.                         TABOBL(3)='D22 '
  899.                         TABOBL(4)='D41 '
  900.                         TABOBL(5)='D42 '
  901.                         TABOBL(6)='D44 '
  902.                         TABOBL(7)='V1X '
  903.                         TABOBL(8)='V1Y '
  904.                         TABOBL(9)='D31 '
  905.                         TABOBL(10)='D32 '
  906.                         TABOBL(11)='D33 '
  907.                         TABOBL(12)='D43 '
  908. *
  909.                         JGFAC=5
  910.                         TABFAC(1)='RHO '
  911.                         TABFAC(2)='ALP1'
  912.                         TABFAC(3)='ALP2'
  913.                         TABFAC(4)='AL12'
  914.                         TABFAC(5)='DIM3'
  915.  
  916.          ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  917. *                    DEFORMATION PLANE ,AXISYMETRIE
  918.                         JGOBL=12
  919.                         TABOBL(1)='D11 '
  920.                         TABOBL(2)='D21 '
  921.                         TABOBL(3)='D22 '
  922.                         TABOBL(4)='D31 '
  923.                         TABOBL(5)='D32 '
  924.                         TABOBL(6)='D33 '
  925.                         TABOBL(7)='D41 '
  926.                         TABOBL(8)='D42 '
  927.                         TABOBL(9)='D43 '
  928.                         TABOBL(10)='D44 '
  929.                         TABOBL(11)='V1X '
  930.                         TABOBL(12)='V1Y '
  931. *
  932.                         JGFAC=5
  933.                         TABFAC(1)='RHO '
  934.                         TABFAC(2)='ALP1'
  935.                         TABFAC(3)='ALP2'
  936.                         TABFAC(4)='AL12'
  937.                         TABFAC(5)='ALP3'
  938.  
  939.                      ELSEIF (IFOUR.EQ.1) THEN
  940. *                       AXISYMETRIE DE FOURIER
  941.                         JGOBL=15
  942.                         TABOBL(1)='D11 '
  943.                         TABOBL(2)='D21 '
  944.                         TABOBL(3)='D22 '
  945.                         TABOBL(4)='D31 '
  946.                         TABOBL(5)='D32 '
  947.                         TABOBL(6)='D33 '
  948.                         TABOBL(7)='D41 '
  949.                         TABOBL(8)='D42 '
  950.                         TABOBL(9)='D43 '
  951.                         TABOBL(10)='D44 '
  952.                         TABOBL(11)='D55 '
  953.                         TABOBL(12)='D65 '
  954.                         TABOBL(13)='D66 '
  955.                         TABOBL(14)='V1X '
  956.                         TABOBL(15)='V1Y '
  957. *
  958.                         JGFAC=5
  959.                         TABFAC(1)='RHO '
  960.                         TABFAC(2)='ALP1'
  961.                         TABFAC(3)='ALP2'
  962.                         TABFAC(4)='AL12'
  963.                         TABFAC(5)='ALP3'
  964.                      ENDIF
  965.                   ENDIF
  966.                ENDIF
  967.  
  968.             ELSEIF (IPLAC.EQ.4) THEN
  969.                IF (MFR.EQ.33) THEN
  970. *                 MILIEU POREUX ISOTROPE
  971. *
  972.                   JGOBL=12
  973.                   TABOBL(1) ='YOUN'
  974.                   TABOBL(2) ='NU  '
  975.                   TABOBL(3) ='RHO '
  976.                   TABOBL(4) ='ALPH'
  977.                   TABOBL(5) ='COB '
  978.                   TABOBL(6) ='MOB '
  979.                   TABOBL(7) ='ALPM'
  980.                   TABOBL(8) ='PERM'
  981.                   TABOBL(9) ='VISC'
  982.                   TABOBL(10)='KF  '
  983.                   TABOBL(11)='RHOF'
  984.                   TABOBL(12)='ALPF'
  985. *
  986.                   IF (IFOUR.EQ.-2) THEN
  987.                     JGFAC=1
  988.                     TABFAC(1)='DIM3'
  989.                   ENDIF
  990.                ENDIF
  991.  
  992.             ELSEIF (IPLAC.EQ.5) THEN
  993.                IF (MFR.EQ.37) THEN
  994. *                 MILIEU HOMOGENEISE
  995. *
  996.                   JGOBL=16
  997.                   TABOBL( 1)='B11 '
  998.                   TABOBL( 2)='B22 '
  999.                   TABOBL( 3)='B12 '
  1000.                   TABOBL( 4)='ROF '
  1001.                   TABOBL( 5)='ROS '
  1002.                   TABOBL( 6)='YOUN'
  1003.                   TABOBL( 7)='CSON'
  1004.                   TABOBL( 8)='RORF'
  1005.                   TABOBL( 9)='CREF'
  1006.                   TABOBL(10)='LCAR'
  1007.                   TABOBL(11)='E111'
  1008.                   TABOBL(12)='E112'
  1009.                   TABOBL(13)='E121'
  1010.                   TABOBL(14)='E122'
  1011.                   TABOBL(15)='E221'
  1012.                   TABOBL(16)='E222'
  1013.                ENDIF
  1014.             ELSEIF (IPLAC.EQ.6) THEN
  1015. C
  1016. C         MATERIAU ELASTIQUE UNIDIRECTIONNEL
  1017. C
  1018.              IF ((MFR.EQ.1.OR.MFR.EQ.31).AND.IDIM.EQ.3) THEN
  1019.               JGOBL=7
  1020.               TABOBL(1)='YOUN'
  1021.               TABOBL(2)='V1X '
  1022.               TABOBL(3)='V1Y '
  1023.               TABOBL(4)='V1Z '
  1024.               TABOBL(5)='V2X '
  1025.               TABOBL(6)='V2Y '
  1026.               TABOBL(7)='V2Z '
  1027.  
  1028.              ELSE
  1029.               JGOBL=3
  1030.               TABOBL(1)='YOUN'
  1031.               TABOBL(2)='V1X '
  1032.               TABOBL(3)='V1Y '
  1033.              ENDIF
  1034.  
  1035.              IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND.
  1036.      .                  IFOUR.EQ.-2) THEN
  1037.                JGFAC=3
  1038.                TABFAC(3)='DIM3'
  1039.              ELSE
  1040.                JGFAC=2
  1041.              ENDIF
  1042.              TABFAC(1)='RHO '
  1043.              TABFAC(2)='ALPH'
  1044.  
  1045.             ELSEIF (IPLAC.EQ.7) THEN
  1046. C      MODELE ET MATERIAU DE LA SECTION DU MODELE A FIBRE
  1047. C
  1048.              JGOBL=2
  1049.              TABOBL(1)='MODS'
  1050.              TABOBL(2)='MATS'
  1051.  
  1052.              JGFAC=1
  1053.              TABFAC(1)='MANO'
  1054. C ---- IPLAC 8: Ajoute par Jiang, 22/08/1995 et fleuret 28/05/96
  1055.  
  1056.             ELSEIF (IPLAC.EQ.8) THEN
  1057. C            MODELE ET MATERIAU POUR LE CABLE PRECONTRAINT
  1058. C
  1059.              JGOBL=1
  1060.              TABOBL(1)='YOUN'
  1061.  
  1062.              JGFAC=8
  1063.              TABFAC(1)='FF'
  1064.              TABFAC(2)='PHIF'
  1065.              TABFAC(3)='GANC'
  1066.              TABFAC(4)='RMU0'
  1067.              TABFAC(5)='FPRG'
  1068.              TABFAC(6)='RH10'
  1069.              TABFAC(7)='ALPH'
  1070.              TABFAC(8)='RHO '
  1071.  
  1072.             ELSEIF (IPLAC.EQ.9) THEN
  1073. C            MODAL
  1074. C
  1075.              JGOBL=3
  1076.              TABOBL(1)='FREQ'
  1077.              TABOBL(2)='MASS'
  1078.              TABOBL(3)='DEFO'
  1079.  
  1080.              JGFAC = 9
  1081.              TABFAC(1) = 'AMOR'
  1082.              TABFAC(2) = 'CGRA'
  1083.              TABFAC(3) = 'MADE'
  1084.              TABFAC(4) = 'RICR'
  1085.              TABFAC(5) = 'MAIB'
  1086.              TABFAC(6) = 'MACR'
  1087.              TABFAC(7) = 'AMCR'
  1088.              TABFAC(8) = 'ALP0'
  1089.              TABFAC(9) = 'ECRO'
  1090.  
  1091.             ELSEIF (IPLAC.EQ.10) THEN
  1092. C            STATIQUE
  1093. C
  1094.              JGOBL=3
  1095.              TABOBL(1)='DEFO'
  1096.              TABOBL(2)='RIDE'
  1097.              TABOBL(3)='MADE'
  1098.  
  1099.              JGFAC = 8
  1100.              TABFAC(1) = 'AMOR'
  1101.              TABFAC(2) = 'RICR'
  1102.              TABFAC(3) = 'MAIA'
  1103.              TABFAC(4) = 'MAIB'
  1104.              TABFAC(5) = 'MACR'
  1105.              TABFAC(6) = 'AMCR'
  1106.              TABFAC(7) = 'BET0'
  1107.              TABFAC(8) = 'ECRO'
  1108.  
  1109.             ELSEIF (IPLAC.EQ.11) THEN
  1110. C            ZONE_COHESIVE
  1111. C
  1112.              JGOBL=2
  1113.              TABOBL(1)='KS'
  1114.              TABOBL(2)='KN'
  1115.             ENDIF
  1116.  
  1117.             INMAT=3
  1118.             GOTO 20
  1119.  
  1120.          ELSE
  1121.            INMAT=2
  1122.            GOTO 10
  1123.          ENDIF
  1124.       ENDIF
  1125.       INMAT=0
  1126.  10   CONTINUE
  1127. *
  1128. *     MATERIAU ELASTIQUE ISOTROPE
  1129. *
  1130.         IF (MFR.EQ.35) THEN
  1131. *         VALABLE EN 2D COMME EN 3D
  1132.           JGOBL=2
  1133.           TABOBL(1)='KS  '
  1134.           TABOBL(2)='KN  '
  1135.  
  1136.           JGFAC=2
  1137.           TABFAC(1)='RHO '
  1138.           TABFAC(2)='ALPN'
  1139. *
  1140.         ELSEIF (MFR.EQ.78) THEN
  1141. *         VALABLE EN 2D COMME EN 3D
  1142.           JGOBL=2
  1143.           TABOBL(1)='KS  '
  1144.           TABOBL(2)='KN  '
  1145. *
  1146. *  JOINT CISAILLEMENT (2D)
  1147. *
  1148.         ELSEIF (MFR.EQ.53) THEN
  1149. *         VALABLE EN 2D "COMME EN 3D"
  1150.  
  1151.           JGOBL=1
  1152.           TABOBL(1)='KS  '
  1153. *
  1154.           JGFAC=2
  1155.           TABFAC(1)='RHO '
  1156.           TABFAC(2)='ALPN'
  1157. *
  1158. * TOUS LES AUTRES CAS
  1159. *
  1160.         ELSE
  1161.           JGOBL=2
  1162.           TABOBL(1)='YOUN'
  1163.           TABOBL(2)='NU  '
  1164. *
  1165.           IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.63).AND.
  1166.      .     IFOUR.EQ.-2) THEN
  1167.             JGFAC=4
  1168.             TABFAC(3)='DIM3'
  1169.             TABFAC(4)='VISQ'
  1170.           ELSE
  1171.             JGFAC=3
  1172.             TABFAC(3)='VISQ'
  1173.           ENDIF
  1174.           TABFAC(1)='RHO '
  1175.           TABFAC(2)='ALPH'
  1176.  
  1177. * AM  CAS FORMULATION NON-LOCALE
  1178.  
  1179.           INLOC = 0
  1180.           MN3 = INFMOD(/1)
  1181.           IF (MN3.GE.13) INLOC=-1*INFMOD(13)
  1182.           IF (INLOC.GT.0) THEN
  1183.             JGOBL0 = JGFAC
  1184. *  moyenne
  1185.             IF(INLOC.EQ.1) THEN
  1186.               JGFAC=JGOBL0+1
  1187.               TABFAC(JGOBL0+1)='LCAR'
  1188. *  stress-based
  1189.             ELSEIF(INLOC.EQ.2) THEN
  1190.               JGFAC=JGOBL0+2
  1191.               TABFAC(JGOBL0+1)='LCAR'
  1192.               TABFAC(JGOBL0+2)='SBFT'
  1193.  
  1194. *  helmholtz
  1195.             ELSE IF(INLOC.EQ.3) THEN
  1196. *  sellier 03/07/20
  1197. *              JGFAC=JGOBL0+6
  1198. *              TABFAC(JGOBL0+1)='LCAR'
  1199. *              TABFAC(JGOBL0+2)='LCF1'
  1200. *              TABFAC(JGOBL0+3)='LCF2'
  1201. *              TABFAC(JGOBL0+4)='LCF3'
  1202. *              TABFAC(JGOBL0+5)='LCF4'
  1203. *              TABFAC(JGOBL0+6)='LCF5'       
  1204.               JGFAC=JGOBL0+NB_PARA_HELM
  1205. c             variables de Helmholtz a diffusions ORThotropes
  1206. c             modif liste para heilmholtz sellier 10//11/22
  1207. c             print*,'Dans IDMATR ', NB_HELM, 'formulationss Helmholtz'
  1208.               do I=1,NB_HELM
  1209. c                 numero de la variable de Helmholtz                 
  1210.                   write (motnl1,'(I1)') I
  1211. c                 Capacité (infinie si BH=1)           
  1212.                   motnl3='CAP' 
  1213.                   motnl4=motnl3//motnl1 
  1214.                   TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+1)=motnl4
  1215. c                 print*,motnl4
  1216. c                 Blocage de Dirichlet si 1 sinon 0           
  1217.                   motnl3='BLO' 
  1218.                   motnl4=motnl3//motnl1          
  1219.                   TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+2)=motnl4
  1220. c                 print*,motnl4
  1221. c                 Valeur imposee a la zone de blocage de Dirichlet             
  1222.                   motnl3='DEP' 
  1223.                   motnl4=motnl3//motnl1          
  1224.                   TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+3)=motnl4
  1225. c                 print*,motnl4
  1226. c                 Valeur initiale de la variable de Helmholtz
  1227.                   motnl3='INI' 
  1228.                   motnl4=motnl3//motnl1          
  1229.                   TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+4)=motnl4
  1230. c                 print*,motnl4
  1231. c                 indicateur logique de lineratite pour ne pas passer par istep=3 si 1
  1232.                   motnl3='LIN' 
  1233.                   motnl4=motnl3//motnl1          
  1234.                   TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+5)=motnl4
  1235. c                 print*,motnl4
  1236.                   do J=1,3
  1237. c                   coefficients de diffusion                  
  1238.                     write (motnld1,'(I1)') J
  1239. c                   diffusion 
  1240.                     motnl2='DH'
  1241. c                   1er indice: numero de la variable d helmholtz
  1242.                     motnl3=motnl2//motnl1
  1243. c                   2eme indice : direction 
  1244.                     motnl4=motnl3//motnld1
  1245.                     TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
  1246.      #              5+(J-1)*4+1)=motnl4 
  1247. c                   print*,motnl4               
  1248.                     do K=1,3
  1249.                         write (motnla1,'(I1)') K              
  1250. c                       direction principales de diffusion 
  1251.                         motnlv1='V'
  1252. c                       indice numero de la variable
  1253.                         motnl2=motnlv1//motnl1
  1254. c                       indice direction principale
  1255.                         motnl3=motnl2//motnld1
  1256. c                       projection sur axe indice de base fixe
  1257.                         motnl4=motnl3//motnla1
  1258.                         TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
  1259.      #                  5+(J-1)*4+1+K)=motnl4
  1260. c                       print*,motnl4
  1261.                     end do
  1262.                   end do
  1263.               end do
  1264. Cc            on rajoute 2 inutiles pour compenser un decalage constatte dans castem20               
  1265. C               TABFAC(JGOBL0+NB_PARA_HELM+1)='FAC1'
  1266. C               TABFAC(JGOBL0+NB_PARA_HELM+2)='FAC2'              
  1267. *  fin modif sellier 
  1268. *  cas non prevu
  1269.             ELSE
  1270.               IRET = 0
  1271.               CALL ERREUR(5)
  1272.             ENDIF
  1273.           ENDIF
  1274. *
  1275.       ENDIF
  1276. *
  1277.       IF (INMAT.EQ.0) GOTO 9999
  1278. *
  1279.  19   CONTINUE
  1280.  20   CONTINUE
  1281.       DO jm = 1,matmod(/2)
  1282.         IF (matmod(jm).eq.'IMPEDANCE') THEN
  1283.           imate = imatee
  1284.           INMAT = INMAT + 1
  1285.           JGOBL = 0
  1286.           JGFAC = 0
  1287.  
  1288.           IF(CMATEE.EQ.'IMPELAST') THEN
  1289. *IMPE_ELAS
  1290.             JGOBL0 = JGOBL
  1291.             JGOBL= JGOBL0+1
  1292.             TABOBL(JGOBL0 + 1) ='RAID'
  1293.             JGOBL0 = JGFAC
  1294.  
  1295.             JGFAC = JGOBL0+4
  1296.             TABFAC(JGOBL0+1) = 'AMOR'
  1297.             TABFAC(JGOBL0+2) = 'ZNU'
  1298.             TABFAC(JGOBL0+3) = 'MASS'
  1299.             TABFAC(JGOBL0+4) = 'ALPH'
  1300.             GOTO 22
  1301.           ELSEIF (CMATEE.EQ.'IMPVOIGT'.or.CMATEE.EQ.'IMPREUSS') THEN
  1302. *IMPE_VOIGT ou IMPE_REUSS
  1303.             JGOBL0 = JGOBL
  1304.             JGOBL= JGOBL0+2
  1305.             TABOBL(JGOBL0 + 1) ='RAID'
  1306.             TABOBL(JGOBL0 + 2) ='VISC'
  1307.  
  1308.             JGOBL0 = JGFAC
  1309.             JGFAC = JGOBL0+2
  1310.             TABFAC(JGOBL0+1) = 'MASS'
  1311.             TABFAC(JGOBL0+2) = 'AMOR'
  1312.             GOTO 22
  1313.           ELSEIF (CMATEE.EQ.'IMPCOMPL') then
  1314. *IMPE_COMPLEXE
  1315.            JGOBL0 = JGOBL
  1316.            JGOBL= JGOBL0+1
  1317.            TABOBL(JGOBL0 + 1) ='RAID'
  1318.  
  1319.            JGOBL0 = JGFAC
  1320.            JGFAC = JGFAC+4
  1321.            TABFAC(JGOBL0+1) = 'MOCO'
  1322.            TABFAC(JGOBL0+2) = 'VISC'
  1323.            TABFAC(JGOBL0+3) = 'MASS'
  1324.            TABFAC(JGOBL0+4) = 'AMOR'
  1325.            GOTO 22
  1326.           ELSE
  1327.           ENDIF
  1328.  22       CONTINUE
  1329. *       if (mele.eq.45) then
  1330.            JGOBL0 = JGFAC
  1331.            JGFAC = JGOBL0+3
  1332.  
  1333.           IF(CMATEE.EQ.'IMPELAST') THEN
  1334.            TABFAC(JGOBL0+1) = TABFAC(JGOBL0-2)
  1335.            TABFAC(JGOBL0+2) = TABFAC(JGOBL0-1)
  1336.            TABFAC(JGOBL0+3) = TABFAC(JGOBL0)
  1337.            TABFAC(JGOBL0-2) = 'CPLE'
  1338.            TABFAC(JGOBL0-1) = 'INER'
  1339.            TABFAC(JGOBL0)   = 'AROT'
  1340.           ELSE
  1341.            TABFAC(JGOBL0+1) = 'CPLE'
  1342.            TABFAC(JGOBL0+2) = 'INER'
  1343.            TABFAC(JGOBL0+3) = 'AROT'
  1344.           ENDIF
  1345. *        endif
  1346.          IF (NMAT.GE.INMAT) inmat = inmat+1
  1347. *            GOTO 9999
  1348.         ENDIF
  1349.       ENDDO
  1350. *
  1351.       IF (NMAT.GE.INMAT) THEN
  1352.          CALL MODNLI(MOMODL,NMOD)
  1353.          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1354.          IF (IPLAC.EQ.1) THEN
  1355.             INMAT=INMAT+1
  1356.             CALL MODPLA(MOMODL,NMOD)
  1357.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1358.             IF (IPLAC.EQ.0) THEN
  1359. *              ISOTROPE
  1360.                IPLAC=1
  1361. *           ELSE
  1362. *              INMAT=INMAT+1
  1363.             ENDIF
  1364.             CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1365.          ELSEIF (IPLAC.EQ.2) THEN
  1366.             INMAT=INMAT+1
  1367.             CALL MODFLU(MOMODL,NMOD)
  1368.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1369.             IF (IPLAC.EQ.0) THEN
  1370. *              NORTON
  1371.                IPLAC=1
  1372. *           ELSE
  1373. *              INMAT=INMAT+1
  1374.             ENDIF
  1375.             CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1376.          ELSEIF (IPLAC.EQ.3) THEN
  1377.             INMAT=INMAT+1
  1378.             CALL MODVIS(MOMODL,NMOD)
  1379.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1380.             IF (IPLAC.EQ.0) THEN
  1381. *              ONERA
  1382.                IPLAC=2
  1383. *           ELSE
  1384. *              INMAT=INMAT+1
  1385.             ENDIF
  1386. C
  1387. C Cas particulier des modeles GATT_MONERIE et UO2
  1388. C   'RHO ' et 'ALPH' sont obligatoires
  1389.             IF (IPLAC.EQ.18.OR.IPLAC.EQ.19) THEN
  1390.               JGOBL0 = JGOBL
  1391.               JGOBL=JGOBL0+2
  1392.               TABOBL(JGOBL0+1)='RHO '
  1393.               TABOBL(JGOBL0+2)='ALPH'
  1394.  
  1395.               JGOBL0=JGFAC
  1396.               JGFAC=JGOBL0-2
  1397.               IF (JGFAC.GT.0) THEN
  1398.                 DO 200 I=1,JGFAC
  1399.                   TABFAC(I)=TABFAC(I+2)
  1400.  200            CONTINUE
  1401.               ENDIF
  1402.  
  1403. C ===
  1404. C   Modeles SYMONDS & COWPER SYCO1 et SYCO2
  1405. C ===
  1406.            ELSEIF (IPLAC.EQ.28) THEN
  1407.              JGOBL0=JGOBL
  1408.              JGOBL=JGOBL0+1
  1409.              TABOBL(JGOBL0+1)='ECRO'
  1410.  
  1411.           ELSEIF (IPLAC.EQ.29) THEN
  1412.              JGOBL0=JGOBL
  1413.              JGOBL=JGOBL0+1
  1414.              TABOBL(JGOBL0+1)='ECRO'
  1415.           ENDIF
  1416.           CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1417. C
  1418.         ELSEIF (IPLAC.EQ.4) THEN
  1419.             INMAT=INMAT+1
  1420.             CALL MODEND(MOMODL,NMOD)
  1421.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1422.             IF (IPLAC.EQ.0) THEN
  1423. *              MAZARS
  1424.                IPLAC=1
  1425. *           ELSE
  1426. *              INMAT=INMAT+1
  1427.             ENDIF
  1428.             CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1429.          ELSEIF (IPLAC.EQ.5) THEN
  1430.             INMAT=INMAT+1
  1431.             CALL MODPLE(MOMODL,NMOD)
  1432.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1433.             IF (IPLAC.EQ.0) THEN
  1434. *              TRIAXIAL P/Y
  1435.                IPLAC=1
  1436. *           ELSE
  1437. *              INMAT=INMAT+1
  1438.             ENDIF
  1439.             CALL IDPLEN(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1440.          ELSEIF (IPLAC.EQ.6) THEN
  1441.             INMAT=INMAT+1
  1442.             CALL MODENL(MOMODL,NMOD)
  1443.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1444.             IF (IPLAC.EQ.0) THEN
  1445. *              ISOTROPE
  1446.                IPLAC=1
  1447. *           ELSE
  1448. *              INMAT=INMAT+1
  1449.             ENDIF
  1450.             CALL IDELNL(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1451.          ELSEIF (IPLAC.EQ.7) THEN
  1452. C        Pas de composantes supplementaires pour une loi 'VISCO_EXTERNE'
  1453.             IRET = 1
  1454.          ELSEIF(IPLAC.EQ.0.AND.IMATEE.EQ.31) THEN
  1455.            GOTO 9999
  1456.          ELSE
  1457.             GOTO 9999
  1458.          ENDIF
  1459. C        IRET code retour de IDPLAS,IDFLUA,IDVISC,IDENDO,IDPLEN,IDELNL
  1460.          IF (IRET.EQ.0) GOTO 9999
  1461.        ENDIF
  1462.  
  1463. *      Parametres CRIP et FUSION :
  1464.        do jma=1,matmod(/2)
  1465.          if(matmod(jma).eq.'CRIP') then
  1466.            JGOBL = JGOBL + 1
  1467.            TABOBL(JGOBL) = 'LIMP'
  1468.  
  1469.          elseif(matmod(jma).eq.'FUSION') then
  1470.            JGOBL = JGOBL + 1
  1471.            TABOBL(JGOBL) = 'TFUS'
  1472.  
  1473.          endif
  1474.        enddo
  1475. *
  1476.         GOTO 9999
  1477.       ENDIF
  1478.  
  1479. *--------------------------------------------------------------------
  1480. *     CAS DE LA FORMULATION POREUX
  1481. *--------------------------------------------------------------------
  1482.       CALL PLACE(FORMOD,NFOR,IPORE,'POREUX')
  1483.       IF (IPORE.NE.0) THEN
  1484. *
  1485.        CALL MODLIN(MOMODL,NMOD)
  1486.        CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
  1487.        IF (IRET.EQ.0) GOTO 9999
  1488. *
  1489. *     D'ABORD :  CAS NON ISOTROPE
  1490. *
  1491.        IF (NMAT.GE.2) THEN
  1492.          CALL MODELA(MOMODL,NMOD)
  1493.          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
  1494.          IF (IPLAC.NE.0) THEN
  1495.             INMAT=3
  1496.             IF (IPLAC.EQ.1) THEN
  1497.                INMAT=3
  1498.                GOTO 30
  1499.             ENDIF
  1500.             IF (IPLAC.EQ.2) THEN
  1501. *
  1502. *              MATERIAU ELASTIQUE ORTHOTROPE
  1503. *
  1504.                IF (MFR.EQ.33) THEN
  1505. *
  1506. *                 ELEMENTS MASSIFS
  1507. *
  1508.                   IF(IDIM.EQ.3)THEN
  1509. *                    ELEMENTS 3D
  1510.                      JGOBL=19
  1511.                      TABOBL(1)='YG1 '
  1512.                      TABOBL(2)='YG2 '
  1513.                      TABOBL(3)='YG3 '
  1514.                      TABOBL(4)='NU12'
  1515.                      TABOBL(5)='NU23'
  1516.                      TABOBL(6)='NU13'
  1517.                      TABOBL(7)='G12 '
  1518.                      TABOBL(8)='G23 '
  1519.                      TABOBL(9)='G13 '
  1520.                      TABOBL(10)='V1X '
  1521.                      TABOBL(11)='V1Y '
  1522.                      TABOBL(12)='V1Z '
  1523.                      TABOBL(13)='V2X '
  1524.                      TABOBL(14)='V2Y '
  1525.                      TABOBL(15)='V2Z '
  1526.                      TABOBL(16)='COB1'
  1527.                      TABOBL(17)='COB2'
  1528.                      TABOBL(18)='COB3'
  1529.                      TABOBL(19)='MOB '
  1530. *
  1531.                      JGFAC=12
  1532.                      TABFAC(1)='RHO '
  1533.                      TABFAC(2)='ALP1'
  1534.                      TABFAC(3)='ALP2'
  1535.                      TABFAC(4)='ALP3'
  1536.                      TABFAC(5)='ALPM'
  1537.                      TABFAC(6)='PER1'
  1538.                      TABFAC(7)='PER2'
  1539.                      TABFAC(8)='PER3'
  1540.                      TABFAC(9)='VISC'
  1541.                      TABFAC(10)='KF  '
  1542.                      TABFAC(11)='RHOF'
  1543.                      TABFAC(12)='ALPF'
  1544.                   ELSEIF (IDIM.EQ.2) THEN
  1545.                      IF(IFOUR.EQ.-2)THEN
  1546. *                       CONTRAINTE PLANE
  1547.                         JGOBL=12
  1548.                         TABOBL(1)='YG1 '
  1549.                         TABOBL(2)='YG2 '
  1550.                         TABOBL(3)='NU12'
  1551.                         TABOBL(4)='G12 '
  1552.                         TABOBL(5)='V1X '
  1553.                         TABOBL(6)='V1Y '
  1554.                         TABOBL(7)='COB1'
  1555.                         TABOBL(8)='COB2'
  1556.                         TABOBL(9)='MOB '
  1557.                         TABOBL(10)='YG3 '
  1558.                         TABOBL(11)='NU23'
  1559.                         TABOBL(12)='NU13'
  1560. *
  1561.                         JGFAC=11
  1562.                         TABFAC(1)='RHO '
  1563.                         TABFAC(2)='ALP1'
  1564.                         TABFAC(3)='ALP2'
  1565.                         TABFAC(4)='ALPM'
  1566.                         TABFAC(5)='PER1'
  1567.                         TABFAC(6)='PER2'
  1568.                         TABFAC(7)='VISC'
  1569.                         TABFAC(8)='KF  '
  1570.                         TABFAC(9)='RHOF'
  1571.                         TABFAC(10)='ALPF'
  1572.                         TABFAC(11)='DIM3'
  1573.          ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  1574. *                       DEFORMATION PLANE ,AXISYMETRIE
  1575.                         JGOBL=13
  1576.                         TABOBL(1)='YG1 '
  1577.                         TABOBL(2)='YG2 '
  1578.                         TABOBL(3)='YG3 '
  1579.                         TABOBL(4)='NU12'
  1580.                         TABOBL(5)='NU23'
  1581.                         TABOBL(6)='NU13'
  1582.                         TABOBL(7)='G12 '
  1583.                         TABOBL(8)='V1X '
  1584.                         TABOBL(9)='V1Y '
  1585.                         TABOBL(10)='COB1'
  1586.                         TABOBL(11)='COB2'
  1587.                         TABOBL(12)='COB3'
  1588.                         TABOBL(13)='MOB '
  1589. *
  1590.                         JGFAC=11
  1591.                         TABFAC(1)='RHO '
  1592.                         TABFAC(2)='ALP1'
  1593.                         TABFAC(3)='ALP2'
  1594.                         TABFAC(4)='ALP3'
  1595.                         TABFAC(5)='ALPM'
  1596.                         TABFAC(6)='PER1'
  1597.                         TABFAC(7)='PER2'
  1598.                         TABFAC(8)='VISC'
  1599.                         TABFAC(9)='KF  '
  1600.                         TABFAC(10)='RHOF'
  1601.                         TABFAC(11)='ALPF'
  1602. *
  1603.                      ELSEIF (IFOUR.EQ.1) THEN
  1604. *                       AXISYMETRIE DE FOURIER
  1605.                         JGOBL=15
  1606.                         TABOBL(1)='YG1 '
  1607.                         TABOBL(2)='YG2 '
  1608.                         TABOBL(3)='YG3 '
  1609.                         TABOBL(4)='NU12'
  1610.                         TABOBL(5)='NU23'
  1611.                         TABOBL(6)='NU13'
  1612.                         TABOBL(7)='G12 '
  1613.                         TABOBL(8)='G23 '
  1614.                         TABOBL(9)='G13 '
  1615.                         TABOBL(10)='V1X '
  1616.                         TABOBL(11)='V1Y '
  1617.                         TABOBL(12)='COB1'
  1618.                         TABOBL(13)='COB2'
  1619.                         TABOBL(14)='COB3'
  1620.                         TABOBL(15)='MOB '
  1621. *
  1622.                         JGFAC=12
  1623.                         TABFAC(1)='RHO '
  1624.                         TABFAC(2)='ALP1'
  1625.                         TABFAC(3)='ALP2'
  1626.                         TABFAC(4)='ALP3'
  1627.                         TABFAC(5)='ALPM'
  1628.                         TABFAC(6)='PER1'
  1629.                         TABFAC(7)='PER2'
  1630.                         TABFAC(8)='PER3'
  1631.                         TABFAC(9)='VISC'
  1632.                         TABFAC(10)='KF  '
  1633.                         TABFAC(11)='RHOF'
  1634.                         TABFAC(12)='ALPF'
  1635.                      ENDIF
  1636.                   ENDIF
  1637.                ELSE
  1638. *
  1639. *                CAS NON PREVU
  1640. *
  1641.                  IRET = 0
  1642.                  GOTO 9999
  1643.                ENDIF
  1644.             ELSEIF (IPLAC.EQ.3)THEN
  1645. *
  1646. *              MATERIAU ANISOTROPE ELASTIQUE
  1647. *
  1648.                IF(MFR.EQ.33)THEN
  1649. *
  1650. *                 ELEMENTS MASSIFS
  1651. *
  1652.                   IF(IDIM.EQ.3)THEN
  1653. *                    ELEMENTS 3D
  1654.                      JGOBL=34
  1655.                      TABOBL(1)='D11 '
  1656.                      TABOBL(2)='D21 '
  1657.                      TABOBL(3)='D22 '
  1658.                      TABOBL(4)='D31 '
  1659.                      TABOBL(5)='D32 '
  1660.                      TABOBL(6)='D33 '
  1661.                      TABOBL(7)='D41 '
  1662.                      TABOBL(8)='D42 '
  1663.                      TABOBL(9)='D43 '
  1664.                      TABOBL(10)='D44 '
  1665.                      TABOBL(11)='D51 '
  1666.                      TABOBL(12)='D52 '
  1667.                      TABOBL(13)='D53 '
  1668.                      TABOBL(14)='D54 '
  1669.                      TABOBL(15)='D55 '
  1670.                      TABOBL(16)='D61 '
  1671.                      TABOBL(17)='D62 '
  1672.                      TABOBL(18)='D63 '
  1673.                      TABOBL(19)='D64 '
  1674.                      TABOBL(20)='D65 '
  1675.                      TABOBL(21)='D66 '
  1676.                      TABOBL(22)='V1X '
  1677.                      TABOBL(23)='V1Y '
  1678.                      TABOBL(24)='V1Z '
  1679.                      TABOBL(25)='V2X '
  1680.                      TABOBL(26)='V2Y '
  1681.                      TABOBL(27)='V2Z '
  1682.                      TABOBL(28)='COB1'
  1683.                      TABOBL(29)='COB2'
  1684.                      TABOBL(30)='COB3'
  1685.                      TABOBL(31)='CO12'
  1686.                      TABOBL(32)='CO13'
  1687.                      TABOBL(33)='CO23'
  1688.                      TABOBL(34)='MOB '
  1689. *
  1690.                      JGFAC=18
  1691.                      TABFAC(1)='RHO '
  1692.                      TABFAC(2)='ALP1'
  1693.                      TABFAC(3)='ALP2'
  1694.                      TABFAC(4)='ALP3'
  1695.                      TABFAC(5)='AL12'
  1696.                      TABFAC(6)='AL13'
  1697.                      TABFAC(7)='AL23'
  1698.                      TABFAC(8)='ALPM'
  1699.                      TABFAC(9 )='PER1'
  1700.                      TABFAC(10)='PER2'
  1701.                      TABFAC(11)='PER3'
  1702.                      TABFAC(12)='PE12'
  1703.                      TABFAC(13)='PE13'
  1704.                      TABFAC(14)='PE23'
  1705.                      TABFAC(15)='VISC'
  1706.                      TABFAC(16)='KF  '
  1707.                      TABFAC(17)='RHOF'
  1708.                      TABFAC(18)='ALPF'
  1709.                   ELSEIF (IDIM.EQ.2) THEN
  1710.                      IF (IFOUR.EQ.-2) THEN
  1711. *                       CONTRAINTE PLANE
  1712.                         JGOBL=16
  1713.                         TABOBL(1)='D11 '
  1714.                         TABOBL(2)='D21 '
  1715.                         TABOBL(3)='D22 '
  1716.                         TABOBL(4)='D41 '
  1717.                         TABOBL(5)='D42 '
  1718.                         TABOBL(6)='D44 '
  1719.                         TABOBL(7)='V1X '
  1720.                         TABOBL(8)='V1Y '
  1721.                         TABOBL(9 )='COB1'
  1722.                         TABOBL(10)='COB2'
  1723.                         TABOBL(11)='CO12'
  1724.                         TABOBL(12)='MOB '
  1725.                         TABOBL(13)='D31 '
  1726.                         TABOBL(14)='D32 '
  1727.                         TABOBL(15)='D33 '
  1728.                         TABOBL(16)='D43 '
  1729. *
  1730.                         JGFAC=13
  1731.                         TABFAC(1)='RHO '
  1732.                         TABFAC(2)='ALP1'
  1733.                         TABFAC(3)='ALP2'
  1734.                         TABFAC(4)='AL12'
  1735.                         TABFAC(5)='ALPM'
  1736.                         TABFAC(6)='PER1'
  1737.                         TABFAC(7)='PER2'
  1738.                         TABFAC(8)='PE12'
  1739.                         TABFAC(9)='VISC'
  1740.                         TABFAC(10)='KF  '
  1741.                         TABFAC(11)='RHOF'
  1742.                         TABFAC(12)='ALPF'
  1743.                         TABFAC(13)='DIM3'
  1744. *
  1745.          ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  1746. *                       DEFORMATION PLANE ,AXISYMETRIE
  1747.                         JGOBL=17
  1748.                         TABOBL(1)='D11 '
  1749.                         TABOBL(2)='D21 '
  1750.                         TABOBL(3)='D22 '
  1751.                         TABOBL(4)='D31 '
  1752.                         TABOBL(5)='D32 '
  1753.                         TABOBL(6)='D33 '
  1754.                         TABOBL(7)='D41 '
  1755.                         TABOBL(8)='D42 '
  1756.                         TABOBL(9)='D43 '
  1757.                         TABOBL(10)='D44 '
  1758.                         TABOBL(11)='V1X '
  1759.                         TABOBL(12)='V1Y '
  1760.                         TABOBL(13)='COB1'
  1761.                         TABOBL(14)='COB2'
  1762.                         TABOBL(15)='CO12'
  1763.                         TABOBL(16)='COB3'
  1764.                         TABOBL(17)='MOB '
  1765. *
  1766.                         JGFAC=13
  1767.                         TABFAC(1)='RHO '
  1768.                         TABFAC(2)='ALP1'
  1769.                         TABFAC(3)='ALP2'
  1770.                         TABFAC(4)='AL12'
  1771.                         TABFAC(5)='ALP3'
  1772.                         TABFAC(6)='ALPM'
  1773.                         TABFAC(7)='PER1'
  1774.                         TABFAC(8)='PER2'
  1775.                         TABFAC(9)='PE12'
  1776.                         TABFAC(10)='VISC'
  1777.                         TABFAC(11)='KF  '
  1778.                         TABFAC(12)='RHOF'
  1779.                         TABFAC(13)='ALPF'
  1780. *
  1781.                      ELSEIF (IFOUR.EQ.1) THEN
  1782. *                       AXISYMETRIE DE FOURIER
  1783.                         JGOBL=20
  1784.                         TABOBL(1)='D11 '
  1785.                         TABOBL(2)='D21 '
  1786.                         TABOBL(3)='D22 '
  1787.                         TABOBL(4)='D31 '
  1788.                         TABOBL(5)='D32 '
  1789.                         TABOBL(6)='D33 '
  1790.                         TABOBL(7)='D41 '
  1791.                         TABOBL(8)='D42 '
  1792.                         TABOBL(9)='D43 '
  1793.                         TABOBL(10)='D44 '
  1794.                         TABOBL(11)='D55 '
  1795.                         TABOBL(12)='D65 '
  1796.                         TABOBL(13)='D66 '
  1797.                         TABOBL(14)='V1X '
  1798.                         TABOBL(15)='V1Y '
  1799.                         TABOBL(16)='COB1'
  1800.                         TABOBL(17)='COB2'
  1801.                         TABOBL(18)='CO12'
  1802.                         TABOBL(19)='COB3'
  1803.                         TABOBL(20)='MOB '
  1804. *
  1805.                         JGFAC=14
  1806.                         TABFAC(1)='RHO '
  1807.                         TABFAC(2)='ALP1'
  1808.                         TABFAC(3)='ALP2'
  1809.                         TABFAC(4)='AL12'
  1810.                         TABFAC(5)='ALP3'
  1811.                         TABFAC(6)='ALPM'
  1812.                         TABFAC(7)='PER1'
  1813.                         TABFAC(8)='PER2'
  1814.                         TABFAC(9)='PE12'
  1815.                         TABFAC(10)='PER3'
  1816.                         TABFAC(11)='VISC'
  1817.                         TABFAC(12)='KF  '
  1818.                         TABFAC(13)='RHOF'
  1819.                         TABFAC(14)='ALPF'
  1820.                      ENDIF
  1821.                   ENDIF
  1822.                ELSE
  1823. *
  1824. *                CAS NON PREVU
  1825. *
  1826.                  IRET = 0
  1827.                  GOTO 9999
  1828.                ENDIF
  1829.  
  1830.             ELSEIF (IPLAC.EQ.6) THEN
  1831. C
  1832. C         MATERIAU ELASTIQUE UNIDIRECTIONNEL
  1833. C
  1834.                IF(MFR.EQ.33)THEN
  1835.                  IF (IDIM.EQ.3) THEN
  1836.                    JGOBL=9
  1837.                    TABOBL(1)='YOUN'
  1838.                    TABOBL(2)='V1X '
  1839.                    TABOBL(3)='V1Y '
  1840.                    TABOBL(4)='V1Z '
  1841.                    TABOBL(5)='V2X '
  1842.                    TABOBL(6)='V2Y '
  1843.                    TABOBL(7)='V2Z '
  1844.                    TABOBL(8)='COB '
  1845.                    TABOBL(9)='MOB '
  1846.                  ELSE
  1847.                    JGOBL=5
  1848.                    TABOBL(1)='YOUN'
  1849.                    TABOBL(2)='V1X '
  1850.                    TABOBL(3)='V1Y '
  1851.                    TABOBL(4)='COB '
  1852.                    TABOBL(5)='MOB '
  1853.                  ENDIF
  1854. *
  1855.                  IF(IFOUR.EQ.-2) THEN
  1856.                     JGFAC=9
  1857.                     TABFAC(9)='DIM3'
  1858.                  ELSE
  1859.                    JGFAC=8
  1860.                  ENDIF
  1861.                  TABFAC(1)='RHO '
  1862.                  TABFAC(2)='ALPH'
  1863.                  TABFAC(3)='ALPM'
  1864.                  TABFAC(4)='PERM'
  1865.                  TABFAC(5)='VISC'
  1866.                  TABFAC(6)='KF  '
  1867.                  TABFAC(7)='RHOF'
  1868.                  TABFAC(8)='ALPF'
  1869.                ELSE
  1870. *
  1871. *                CAS NON PREVU
  1872. *
  1873.                  IRET = 0
  1874.                  GO TO 9999
  1875.                ENDIF
  1876. *
  1877.             ENDIF
  1878.             INMAT=3
  1879.             GOTO 40
  1880.          ELSE
  1881.             INMAT=2
  1882.             GOTO 30
  1883.          ENDIF
  1884.       ENDIF
  1885.       INMAT=0
  1886.  30   CONTINUE
  1887. *
  1888. *     CAS MATERIAU POREUX ELASTIQUE ISOTROPE
  1889. *
  1890.       IF(MELE.GE.79.AND.MELE.LE.83)THEN
  1891. *
  1892.         JGOBL=4
  1893.         TABOBL(1) ='YOUN'
  1894.         TABOBL(2) ='NU  '
  1895.         TABOBL(3) ='COB '
  1896.         TABOBL(4) ='MOB '
  1897. *
  1898.         IF(IFOUR.EQ.-2) THEN
  1899.             JGFAC=9
  1900.             TABFAC(9)='DIM3'
  1901.         ELSE
  1902.             JGFAC=8
  1903.         ENDIF
  1904.         TABFAC(1)='RHOF'
  1905.         TABFAC(2)='ALPF'
  1906.         TABFAC(3)='ALPM'
  1907.         TABFAC(4)='PERM'
  1908.         TABFAC(5)='VISC'
  1909.         TABFAC(6)='KF  '
  1910.         TABFAC(7)='RHO '
  1911.         TABFAC(8)='ALPH'
  1912. *
  1913.       ELSEIF(MELE.GE.108.AND.MELE.LE.110)THEN
  1914. *
  1915. *     CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  1916. *
  1917.         JGOBL=4
  1918.         TABOBL(1)='KS  '
  1919.         TABOBL(2)='KN  '
  1920.         TABOBL(3)='COB '
  1921.         TABOBL(4)='MOB '
  1922. *
  1923.         JGFAC=4
  1924.         TABFAC(1)='PERT'
  1925.         TABFAC(2)='PERH'
  1926.         TABFAC(3)='PERB'
  1927.         TABFAC(4)='VISC'
  1928. *
  1929.       ELSEIF(MELE.GE.173.AND.MELE.LE.177)THEN
  1930. *
  1931.         JGOBL=10
  1932.         TABOBL(1) ='YOUN'
  1933.         TABOBL(2) ='NU  '
  1934.         TABOBL(3) ='COP1'
  1935.         TABOBL(4) ='COP2'
  1936.         TABOBL(5) ='CPP1'
  1937.         TABOBL(6) ='CPP2'
  1938.         TABOBL(7) ='KK11'
  1939.         TABOBL(8) ='KK12'
  1940.         TABOBL(9) ='KK21'
  1941.         TABOBL(10)='KK22'
  1942. *
  1943.         IF(IFOUR.EQ.-2) THEN
  1944.             JGFAC=9
  1945.             TABFAC(9)='DIM3'
  1946.         ELSE
  1947.             JGFAC=8
  1948.         ENDIF
  1949.         TABFAC(1)='RHOF'
  1950.         TABFAC(2)='ALPF'
  1951.         TABFAC(3)='ALPM'
  1952.         TABFAC(4)='PK11'
  1953.         TABFAC(5)='PK12'
  1954.         TABFAC(6)='PK21'
  1955.         TABFAC(7)='PK22'
  1956.         TABFAC(8)='ALPH'
  1957. *
  1958.       ELSEIF(MELE.GE.185.AND.MELE.LE.187)THEN
  1959. *
  1960. *     CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  1961. *
  1962.       JGOBL=10
  1963.       TABOBL(1)='KS  '
  1964.       TABOBL(2)='KN  '
  1965.       TABOBL(3)='COP1'
  1966.       TABOBL(4)='COP2'
  1967.       TABOBL(5)='CPP1'
  1968.       TABOBL(6)='CPP2'
  1969.       TABOBL(7)='KK11'
  1970.       TABOBL(8)='KK12'
  1971.       TABOBL(9)='KK21'
  1972.       TABOBL(10)='KK22'
  1973. *
  1974.       IF(IFOUR.EQ.-2) THEN
  1975.           JGFAC=18
  1976.           TABFAC(18)='DIM3'
  1977.       ELSE
  1978.           JGFAC=17
  1979.       ENDIF
  1980.       TABFAC(1)='RHOF'
  1981.       TABFAC(2)='ALPF'
  1982.       TABFAC(3)='ALPM'
  1983.       TABFAC(4)='PT11'
  1984.       TABFAC(5)='PH11'
  1985.       TABFAC(6)='PB11'
  1986.       TABFAC(7)='PT12'
  1987.       TABFAC(8)='PH12'
  1988.       TABFAC(9)='PB12'
  1989.       TABFAC(10)='PT21'
  1990.       TABFAC(11)='PH21'
  1991.       TABFAC(12)='PB21'
  1992.       TABFAC(13)='PT22'
  1993.       TABFAC(14)='PH22'
  1994.       TABFAC(15)='PB22'
  1995.       TABFAC(16)='RHO '
  1996.       TABFAC(17)='ALPH'
  1997.  
  1998.       ELSEIF(MELE.GE.178.AND.MELE.LE.182)THEN
  1999. *
  2000.       JGOBL=17
  2001.       TABOBL(1)='YOUN'
  2002.       TABOBL(2)='NU  '
  2003.       TABOBL(3)='COP1'
  2004.       TABOBL(4)='COP2'
  2005.       TABOBL(5)='COP3'
  2006.       TABOBL(6)='CPP1'
  2007.       TABOBL(7)='CPP2'
  2008.       TABOBL(8)='CPP3'
  2009.       TABOBL(9)='KK11'
  2010.       TABOBL(10)='KK12'
  2011.       TABOBL(11)='KK13'
  2012.       TABOBL(12)='KK21'
  2013.       TABOBL(13)='KK22'
  2014.       TABOBL(14)='KK23'
  2015.       TABOBL(15)='KK31'
  2016.       TABOBL(16)='KK32'
  2017.       TABOBL(17)='KK33'
  2018. *
  2019.       IF(IFOUR.EQ.-2) THEN
  2020.           JGFAC=15
  2021.           TABFAC(15)='DIM3'
  2022.       ELSE
  2023.           JGFAC=14
  2024.       ENDIF
  2025.       TABFAC(1)='RHOF'
  2026.       TABFAC(2)='ALPF'
  2027.       TABFAC(3)='ALPM'
  2028.       TABFAC(4)='PK11'
  2029.       TABFAC(5)='PK12'
  2030.       TABFAC(6)='PK13'
  2031.       TABFAC(7)='PK21'
  2032.       TABFAC(8)='PK22'
  2033.       TABFAC(9)='PK23'
  2034.       TABFAC(10)='PK31'
  2035.       TABFAC(11)='PK32'
  2036.       TABFAC(12)='PK33'
  2037.       TABFAC(13)='RHO '
  2038.       TABFAC(14)='ALPH'
  2039.  
  2040.       ELSEIF(MELE.GE.188.AND.MELE.LE.190)THEN
  2041. *
  2042. *     CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  2043. *
  2044.       JGOBL=17
  2045.       TABOBL(1)='KS  '
  2046.       TABOBL(2)='KN  '
  2047.       TABOBL(3)='COP1'
  2048.       TABOBL(4)='COP2'
  2049.       TABOBL(5)='COP3'
  2050.       TABOBL(6)='CPP1'
  2051.       TABOBL(7)='CPP2'
  2052.       TABOBL(8)='CPP3'
  2053.       TABOBL(9)='KK11'
  2054.       TABOBL(10)='KK12'
  2055.       TABOBL(11)='KK13'
  2056.       TABOBL(12)='KK21'
  2057.       TABOBL(13)='KK22'
  2058.       TABOBL(14)='KK23'
  2059.       TABOBL(15)='KK31'
  2060.       TABOBL(16)='KK32'
  2061.       TABOBL(17)='KK33'
  2062. *
  2063.       IF(IFOUR.EQ.-2) THEN
  2064.           JGFAC=33
  2065.           TABFAC(33)='DIM3'
  2066.       ELSE
  2067.           JGFAC=32
  2068.       ENDIF
  2069.       TABFAC(1)='RHOF'
  2070.       TABFAC(2)='ALPF'
  2071.       TABFAC(3)='ALPM'
  2072.       TABFAC(4)='PT11'
  2073.       TABFAC(5)='PH11'
  2074.       TABFAC(6)='PB11'
  2075.       TABFAC(7)='PT12'
  2076.       TABFAC(8)='PH12'
  2077.       TABFAC(9)='PB12'
  2078.       TABFAC(10)='PT13'
  2079.       TABFAC(11)='PH13'
  2080.       TABFAC(12)='PB13'
  2081.       TABFAC(13)='PT21'
  2082.       TABFAC(14)='PH21'
  2083.       TABFAC(15)='PB21'
  2084.       TABFAC(16)='PT22'
  2085.       TABFAC(17)='PH22'
  2086.       TABFAC(18)='PB22'
  2087.       TABFAC(19)='PT23'
  2088.       TABFAC(20)='PH23'
  2089.       TABFAC(21)='PB23'
  2090.       TABFAC(22)='PT31'
  2091.       TABFAC(23)='PH31'
  2092.       TABFAC(24)='PB31'
  2093.       TABFAC(25)='PT32'
  2094.       TABFAC(26)='PH32'
  2095.       TABFAC(27)='PB32'
  2096.       TABFAC(28)='PT33'
  2097.       TABFAC(29)='PH33'
  2098.       TABFAC(30)='PB33'
  2099.       TABFAC(31)='RHO '
  2100.       TABFAC(32)='ALPH'
  2101.  
  2102.       ENDIF
  2103. *
  2104.       IF (INMAT.EQ.0) THEN
  2105. **      IRET = 0
  2106.         GOTO 9999
  2107.       ENDIF
  2108. *
  2109.  40   CONTINUE
  2110.       IF (NMAT.GE.INMAT) THEN
  2111.          CALL MODNLI(MOMODL,NMOD)
  2112.          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2113.          IF (IPLAC.EQ.1) THEN
  2114.             INMAT=INMAT+1
  2115.             CALL MODPLA(MOMODL,NMOD)
  2116.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2117.             IF (IPLAC.EQ.0) THEN
  2118. *              ISOTROPE
  2119.                IPLAC=1
  2120. *           ELSE
  2121. *              INMAT=INMAT+1
  2122.             ENDIF
  2123.             CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2124.             IF (IRET.EQ.0) GOTO 9999
  2125.          ELSEIF (IPLAC.EQ.2) THEN
  2126.             INMAT=INMAT+1
  2127.             CALL MODFLU(MOMODL,NMOD)
  2128.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2129.             IF (IPLAC.EQ.0) THEN
  2130. *              NORTON
  2131.                IPLAC=1
  2132. *           ELSE
  2133. *              INMAT=INMAT+1
  2134.             ENDIF
  2135.             CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2136.             IF (IRET.EQ.0) GOTO 9999
  2137.          ELSEIF (IPLAC.EQ.3) THEN
  2138.             INMAT=INMAT+1
  2139.             CALL MODVIS(MOMODL,NMOD)
  2140.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2141.             IF (IPLAC.EQ.0) THEN
  2142. *              ONERA
  2143.                IPLAC=2
  2144. *           ELSE
  2145. *              INMAT=INMAT+1
  2146.             ENDIF
  2147.             CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2148.             IF (IRET.EQ.0) GOTO 9999
  2149.          ELSEIF (IPLAC.EQ.4) THEN
  2150.             INMAT=INMAT+1
  2151.             CALL MODEND(MOMODL,NMOD)
  2152.             CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2153.             IF (IPLAC.EQ.0) THEN
  2154. *              MAZARS
  2155.                IPLAC=1
  2156. *           ELSE
  2157. *              INMAT=INMAT+1
  2158.             ENDIF
  2159.             CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2160.             IF (IRET.EQ.0) GOTO 9999
  2161.          ELSE
  2162. *            GOTO 9999
  2163.          ENDIF
  2164. *
  2165. *         En cas de creation de materiaux combinant plusieurs materiaux
  2166. *         deja existant Ex ELASTIQUE ISOTROPE PLASTIQUE PARFAIT FLUAGE N
  2167. *         GOTO 40
  2168. *
  2169. *          GOTO 9999
  2170.       ENDIF
  2171.       GOTO 9999
  2172.       ENDIF
  2173.  
  2174. *--------------------------------------------------------------------
  2175. *     CAS DE LA FORMULATION CONTACT
  2176. *--------------------------------------------------------------------
  2177.       CALL PLACE(FORMOD,NFOR,ICONT,'CONTACT')
  2178.       IF (ICONT.NE.0) THEN
  2179.         CALL MODFRO(MOMODL,NMOD)
  2180.         iplla =0
  2181.         do iou = 1,NMAT
  2182.           CALL PLACE(MOMODL(1),5,IPLAC,MATMOD(iou))
  2183.           iplla = max(iplac,iplla)
  2184.         enddo
  2185.         IF (iplla.EQ.0) THEN
  2186.           IRET = 0
  2187.           GOTO 9999
  2188.         ENDIF
  2189.  
  2190.         IF ((IPLLA.EQ.1) .OR. (IPLLA.EQ.2) .OR. (IPLLA.EQ.3)) THEN
  2191. *         CONTACT SIMPLE
  2192. *         --------------
  2193.           INMAT=1
  2194.           JGOBL=0
  2195.  
  2196.           JGFAC=2
  2197.           TABFAC(1)='JEU'
  2198.           TABFAC(2)='ADHE'
  2199.  
  2200.         ELSEIF (IPLLA .EQ. 4) THEN
  2201. *         FROTTEMENT DE COULOMB
  2202. *         ---------------------
  2203.           INMAT=1
  2204.           JGOBL=1
  2205.           TABOBL(1)='MU  '
  2206.  
  2207.           JGFAC=3
  2208.           TABFAC(1)='COHE'
  2209.           TABFAC(2)='ADHE'
  2210.           TABFAC(3)='JEU'
  2211.  
  2212.         ELSEIF (IPLLA .EQ. 5) THEN
  2213. *         FROTTEMENTS DE CABLES
  2214. *         ---------------------
  2215.           INMAT=1
  2216.           JGOBL=2
  2217.           TABOBL(1)='FF  '
  2218.           TABOBL(2)='PHIF'
  2219.  
  2220.           JGFAC=0
  2221.  
  2222.         ELSE
  2223.           IRET=0
  2224.           CALL ERREUR (261)
  2225.           RETURN
  2226.         ENDIF
  2227.  
  2228.         GOTO 9999
  2229.       ENDIF
  2230.  
  2231. *--------------------------------------------------------------------
  2232. *     CAS DE LA FORMULATION CONTRAINTE   
  2233. *--------------------------------------------------------------------
  2234.       CALL PLACE(FORMOD,NFOR,ICNTR,'CONTRAINTE')
  2235.       IF (ICNTR.NE.0) THEN
  2236.         CALL MODCLI(MOMODL,NMOD)
  2237.         iplla =0
  2238.         do iou = 1,NMAT
  2239.           CALL PLACE(MOMODL(1),NMOD,IPLAC,MATMOD(iou))
  2240.           iplla = max(iplac,iplla)
  2241.         enddo
  2242.         IF (iplla.EQ.0) THEN
  2243.           IRET = 0
  2244.           GOTO 9999
  2245.         ENDIF
  2246.  
  2247.  
  2248.         IF (IPLLA .EQ. 1) THEN
  2249. *         ROTATION                      
  2250. *         ---------------------
  2251.           INMAT=1
  2252.           JGOBL=1
  2253.           TABOBL(1)='ANGL'
  2254.         ELSEIF (IPLLA .EQ. 2) THEN
  2255. *         DEPLACEMENT                     
  2256. *         ---------------------
  2257.           INMAT=1
  2258.           JGOBL=1
  2259.           TABOBL(1)='AMPL'
  2260.         ENDIF
  2261.       ENDIF
  2262. *--------------------------------------------------------------------
  2263. *     CAS DE LA FORMULATION MAGNETODYNAMIQUE
  2264. *--------------------------------------------------------------------
  2265.       CALL PLACE(FORMOD,NFOR,ICONV,'MAGNETODYNAMIQUE')
  2266.       IF (ICONV.NE.0) THEN
  2267. *
  2268. *       FORMULATION EN COQUES
  2269.         IF(MFR.EQ.3) THEN
  2270.           CALL PLACE(MATMOD,NMAT,ISOT,'ISOTROPE')
  2271.           IF(ISOT.NE.0) THEN
  2272.             JGOBL=3
  2273.  
  2274.             TABOBL(1)='ETA'
  2275.             TABOBL(2)='PERM'
  2276.             TABOBL(3)='EPAI'
  2277.           ELSE
  2278.             CALL PLACE(MATMOD,NMAT,IORTH,'ORTHOTROPE')
  2279.             IF(IORTH.NE.0) THEN
  2280.               JGOBL=4
  2281.  
  2282.               TABOBL(1)='ETA1'
  2283.               TABOBL(2)='ETA2'
  2284.               TABOBL(3)='PERM'
  2285.               TABOBL(4)='EPAI'
  2286.             ENDIF
  2287.           ENDIF
  2288.         ENDIF
  2289. *
  2290.         GOTO 9999
  2291.       ENDIF
  2292.  
  2293. *--------------------------------------------------------------------
  2294. *     CAS DE LA FORMULATION FISSURE
  2295. *--------------------------------------------------------------------
  2296.       CALL PLACE(FORMOD,NFOR,ICONV,'FISSURE')
  2297.       IF (ICONV.NE.0) THEN
  2298. *
  2299. * si POISEU_BLASIUS ou POISEU_COLEBROOK ou par defaut
  2300.             JGOBL=1
  2301.  
  2302.             TABOBL(1)='RUGO'
  2303. * si FROTTEMENT1 ou FROTTEMENT2
  2304.             CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT1')
  2305.             CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT2')
  2306.             IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
  2307.               JGOBL=7
  2308.  
  2309.               TABOBL(2)='REC'
  2310.               TABOBL(3)='FK'
  2311.               TABOBL(4)='FA'
  2312.               TABOBL(5)='FB'
  2313.               TABOBL(6)='FC'
  2314.               TABOBL(7)='FD'
  2315.             ELSE
  2316. * si FROTTEMENT3 ou FROTTEMENT4
  2317.               CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT3')
  2318.               CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT4')
  2319.               IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
  2320.                 JGOBL=2
  2321.  
  2322.                 TABOBL(2)='FK'
  2323.               ENDIF
  2324.             ENDIF
  2325.         GOTO 9999
  2326.       ENDIF
  2327.  
  2328. *--------------------------------------------------------------------
  2329. *     CAS DE LA FORMULATION NAVIER_STOKES
  2330. *--------------------------------------------------------------------
  2331.       CALL PLACE(FORMOD,NFOR,INAST,'NAVIER_STOKES')
  2332.       IF (INAST.GT.0) THEN
  2333.         iplac = imatee
  2334.         if (iplac.eq.0) then
  2335.           iret = 0
  2336.           goto 9999
  2337.         endif
  2338. *
  2339.         JGOBL = 0
  2340.         JGFAC = 0
  2341.  
  2342.         if (iplac.eq.4) then
  2343. *NLIN        
  2344.          JGOBL0 = JGOBL
  2345.          JGOBL = JGOBL0+2
  2346.          TABOBL(JGOBL0 + 1) ='RHO '
  2347.          TABOBL(JGOBL0 + 2) ='MU  '
  2348.          JGOBL0 = JGFAC
  2349.          JGFAC = JGOBL0+1
  2350.          TABFAC(JGOBL0 + 1) ='FREQ'
  2351.         endif      
  2352.         GOTO 9999
  2353.       ENDIF
  2354.  
  2355. *--------------------------------------------------------------------
  2356. *     CAS DE LA FORMULATION MELANGE
  2357. *--------------------------------------------------------------------
  2358.       CALL PLACE(FORMOD,NFOR,ICONV,'MELANGE')
  2359.       IF (ICONV.NE.0) THEN
  2360.          CALL MODMEL(MOMODL,NMOD)
  2361.          CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(1))
  2362.          if (iplac.eq.0) then
  2363.            goto 9999
  2364.          endif
  2365.          INMAT=1
  2366. *
  2367. * a priori elements massifs ou coques : pas de verif
  2368. *
  2369.         IF (IPLAC.NE.3.AND.IPLAC.NE.4) THEN
  2370.           JGOBL = 0
  2371.           JGFAC = 0
  2372.           CALL IDMETA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2373.           IF (IRET.EQ.0) GOTO 9999
  2374.         ELSE
  2375.         ENDIF
  2376.         GOTO 9999
  2377.       ENDIF
  2378.  
  2379. *--------------------------------------------------------------------
  2380. *     CAS DE LA FORMULATION LIAISON
  2381. *--------------------------------------------------------------------
  2382.       CALL PLACE(FORMOD,NFOR,ICONV,'LIAISON')
  2383.       IF (ICONV.NE.0) THEN
  2384.         iplac = imatee
  2385.         if (iplac.eq.0) then
  2386.           iret = 0
  2387.           goto 9999
  2388.         endif
  2389. *
  2390.         JGOBL = 0
  2391.         JGFAC = 0
  2392.  
  2393. * 'SORT' facultatif dans tous les cas
  2394.         IF (iplac.EQ.1) THEN
  2395. *PO_PL_FL
  2396.          JGOBL0 = JGOBL
  2397.          JGOBL = JGOBL0+7
  2398.          TABOBL(JGOBL0 + 1) ='NORM'
  2399.          TABOBL(JGOBL0 + 2) ='INER'
  2400.          TABOBL(JGOBL0 + 3) ='CONV'
  2401.          TABOBL(JGOBL0 + 4) ='VISC'
  2402.          TABOBL(JGOBL0 + 5) ='PELO'
  2403.          TABOBL(JGOBL0 + 6) ='PRAP'
  2404.          TABOBL(JGOBL0 + 7) ='JFLU'
  2405.  
  2406.         ELSEIF (iplac.eq.2) then
  2407. *PO_PL_FR
  2408.          JGOBL0 = JGOBL
  2409.          JGOBL = JGOBL0+7
  2410.          TABOBL(JGOBL0 + 1) ='NORM'
  2411.          TABOBL(JGOBL0 + 2) ='RAID'
  2412.          TABOBL(JGOBL0 + 3) ='JEU'
  2413.          TABOBL(JGOBL0 + 4) ='GLIS'
  2414.          TABOBL(JGOBL0 + 5) ='ADHE'
  2415.          TABOBL(JGOBL0 + 6) ='RTAN'
  2416.          TABOBL(JGOBL0 + 7) ='ATAN'
  2417.  
  2418.          JGOBL0 = JGFAC
  2419.          JGFAC = JGOBL0 +2
  2420.          TABFAC(JGOBL0+1) = 'AMOR'
  2421.          TABFAC(JGOBL0+2) = 'LOIC'
  2422.  
  2423.         ELSEIF (iplac.eq.3) then
  2424. *PO_PL
  2425.          JGOBL0 = JGOBL
  2426.          JGOBL= JGOBL0+3
  2427.          TABOBL(JGOBL0 + 1) ='NORM'
  2428.          TABOBL(JGOBL0 + 2) ='RAID'
  2429.          TABOBL(JGOBL0 + 3) ='JEU'
  2430.  
  2431.          JGOBL0 = JGFAC
  2432.          JGFAC = JGOBL0+4
  2433.          TABFAC(JGOBL0+1) = 'LOIC'
  2434.          TABFAC(JGOBL0+2) = 'PERM'
  2435.          TABFAC(JGOBL0 + 3) ='SPLA'
  2436.          TABFAC(JGOBL0 + 4) ='AMOR'
  2437.  
  2438.        ELSEIF (iplac.eq.4) then
  2439. *PO_PO_FR
  2440.          JGOBL0 = JGOBL
  2441.          JGOBL= JGOBL0 + 8
  2442.          TABOBL(JGOBL0 + 1) ='NORM'
  2443.          TABOBL(JGOBL0 + 2) ='RAID'
  2444.          TABOBL(JGOBL0 + 3) ='JEU'
  2445.          TABOBL(JGOBL0 + 4) ='POIB'
  2446.          TABOBL(JGOBL0 + 5) ='ADHE'
  2447.          TABOBL(JGOBL0 + 6) ='RTAN'
  2448.          TABOBL(JGOBL0 + 7) ='ATAN'
  2449.          TABOBL(JGOBL0 + 8) ='GLIS'
  2450.  
  2451.          JGOBL0 = JGFAC
  2452.          JGFAC = JGOBL0 + 3
  2453.          TABFAC(JGOBL0+1) = 'AMOR'
  2454.          TABFAC(JGOBL0+2) = 'LOIC'
  2455.          TABFAC(JGOBL0+3) = 'MODE'
  2456.  
  2457.         ELSEIF (iplac.eq.5) then
  2458. *PO_PO_DP
  2459.          JGOBL0 = JGOBL
  2460.          JGOBL= JGOBL0 + 6
  2461.          TABOBL(JGOBL0 + 1) ='NORM'
  2462.          TABOBL(JGOBL0 + 2) ='ECRO'
  2463.          TABOBL(JGOBL0 + 3) ='JEU'
  2464.          TABOBL(JGOBL0 + 4) ='POIB'
  2465.          TABOBL(JGOBL0 + 5) ='PERM'
  2466.          TABOBL(JGOBL0 + 6) ='LOIC'
  2467.  
  2468.          JGOBL0 = JGFAC
  2469.          JGFAC = JGOBL0 + 1
  2470.          TABFAC(JGOBL0+1) = 'AMOR'
  2471.  
  2472.         ELSEIF (iplac.eq.6) then
  2473. *PO_PO_RP
  2474.          JGOBL0 = JGOBL
  2475.          JGOBL= JGOBL0 + 6
  2476.          TABOBL(JGOBL0 + 1) ='AXRO'
  2477.          TABOBL(JGOBL0 + 2) ='ECRO'
  2478.          TABOBL(JGOBL0 + 3) ='JEU'
  2479.          TABOBL(JGOBL0 + 4) ='POIB'
  2480.          TABOBL(JGOBL0 + 5) ='PERM'
  2481.          TABOBL(JGOBL0 + 6) ='LOIC'
  2482.  
  2483.          JGOBL0 = JGFAC
  2484.          JGFAC = JGOBL0 + 2
  2485.          TABFAC(JGOBL0+1) = 'AMOR'
  2486.          TABFAC(JGOBL0+2) = 'ELAS'
  2487.  
  2488.         ELSEIF (iplac.eq.7) then
  2489. *PO_PO
  2490.          JGOBL0 = JGOBL
  2491.          JGOBL= JGOBL0 + 5
  2492.          TABOBL(JGOBL0 + 1) ='NORM'
  2493.          TABOBL(JGOBL0 + 2) ='RAID'
  2494.          TABOBL(JGOBL0 + 3) ='JEU'
  2495.          TABOBL(JGOBL0 + 4) ='POIB'
  2496.          TABOBL(JGOBL0 + 5) ='PERM'
  2497.  
  2498.          JGOBL0 = JGFAC
  2499.          JGFAC = JGOBL0 + 2
  2500.          TABFAC(JGOBL0+1) = 'AMOR'
  2501.          TABFAC(JGOBL0+2) = 'LOIC'
  2502.  
  2503.         ELSEIF (iplac.eq.8) then
  2504. *PO_CE_MO
  2505.          JGOBL0 = JGOBL
  2506.          JGOBL= JGOBL0 + 8
  2507.          TABOBL(JGOBL0 + 1) ='NORM'
  2508.          TABOBL(JGOBL0 + 2) ='RAID'
  2509.          TABOBL(JGOBL0 + 3) ='PCER'
  2510.          TABOBL(JGOBL0 + 4) ='RAYO'
  2511.          TABOBL(JGOBL0 + 5) ='GLIS'
  2512.          TABOBL(JGOBL0 + 6) ='ADHE'
  2513.          TABOBL(JGOBL0 + 7) ='RTAN'
  2514.          TABOBL(JGOBL0 + 8) ='ATAN'
  2515.  
  2516.          JGOBL0 = JGFAC
  2517.          JGFAC= JGOBL0 + 2
  2518.          TABFAC(JGOBL0+1) = 'CINT'
  2519.          TABFAC(JGOBL0+2) = 'AMOR'
  2520.  
  2521.         ELSEIF (iplac.eq.9) then
  2522. *PO_CE_FR
  2523.          JGOBL0 = JGOBL
  2524.          JGOBL= JGOBL0 + 8
  2525.          TABOBL(JGOBL0 + 1) ='NORM'
  2526.          TABOBL(JGOBL0 + 2) ='RAID'
  2527.          TABOBL(JGOBL0 + 3) ='EXCE'
  2528.          TABOBL(JGOBL0 + 4) ='RAYO'
  2529.          TABOBL(JGOBL0 + 5) ='GLIS'
  2530.          TABOBL(JGOBL0 + 6) ='ADHE'
  2531.          TABOBL(JGOBL0 + 7) ='RTAN'
  2532.          TABOBL(JGOBL0 + 8) ='ATAN'
  2533.  
  2534.          JGOBL0 = JGFAC
  2535.          JGFAC = JGOBL0 + 2
  2536.          TABFAC(JGOBL0+1) = 'CINT'
  2537.          TABFAC(JGOBL0+2) = 'AMOR'
  2538.  
  2539.         ELSEIF (iplac.eq.10) then
  2540. *PO_CE
  2541.          JGOBL0 = JGOBL
  2542.          JGOBL= JGOBL0 + 4
  2543.          TABOBL(JGOBL0 + 1) ='NORM'
  2544.          TABOBL(JGOBL0 + 2) ='RAID'
  2545.          TABOBL(JGOBL0 + 3) ='EXCE'
  2546.          TABOBL(JGOBL0 + 4) ='RAYO'
  2547.  
  2548.          JGOBL0 = JGFAC
  2549.          JGFAC = JGOBL0 + 1
  2550.          TABFAC(JGOBL0+1) = 'AMOR'
  2551.  
  2552.         ELSEIF (iplac.eq.11) then
  2553. *CE_PL_FR
  2554.          JGOBL0 = JGOBL
  2555.          JGOBL= JGOBL0 + 8
  2556.          TABOBL(JGOBL0 + 1) ='NORM'
  2557.          TABOBL(JGOBL0 + 2) ='RAID'
  2558.          TABOBL(JGOBL0 + 3) ='JEU'
  2559.          TABOBL(JGOBL0 + 4) ='RAYS'
  2560.          TABOBL(JGOBL0 + 5) ='GLIS'
  2561.          TABOBL(JGOBL0 + 6) ='ADHE'
  2562.          TABOBL(JGOBL0 + 7) ='RTAN'
  2563.          TABOBL(JGOBL0 + 8) ='ATAN'
  2564.  
  2565.          JGOBL0 = JGFAC
  2566.          JGFAC = JGOBL0+1
  2567.          TABFAC(JGOBL0+1) = 'AMOR'
  2568.  
  2569.         ELSEIF (iplac.eq.12) then
  2570. *CE_CE_FR
  2571.          JGOBL0 = JGOBL
  2572.          JGOBL= JGOBL0 + 9
  2573.          TABOBL(JGOBL0 + 1) ='NORM'
  2574.          TABOBL(JGOBL0 + 2) ='RAID'
  2575.          TABOBL(JGOBL0 + 3) ='EXCE'
  2576.          TABOBL(JGOBL0 + 4) ='RAYS'
  2577.          TABOBL(JGOBL0 + 5) ='GLIS'
  2578.          TABOBL(JGOBL0 + 6) ='ADHE'
  2579.          TABOBL(JGOBL0 + 7) ='RTAN'
  2580.          TABOBL(JGOBL0 + 8) ='ATAN'
  2581.          TABOBL(JGOBL0 + 9) ='RAYB'
  2582.  
  2583.          JGOBL0 = JGFAC
  2584.          JGFAC = JGOBL0 + 2
  2585.          TABFAC(JGOBL0+1) = 'AMOR'
  2586.          TABFAC(JGOBL0+2) = 'CINT'
  2587.  
  2588.         ELSEIF (iplac.eq.13.or.iplac.eq.14) then
  2589. *PR_PR_IN ou PR_PR_EX
  2590.          JGOBL0 = JGOBL
  2591.          JGOBL= JGOBL0 + 5
  2592.          TABOBL(JGOBL0 + 1) ='NORM'
  2593.          TABOBL(JGOBL0 + 2) ='RAID'
  2594.          TABOBL(JGOBL0 + 3) ='PFIX'
  2595.          TABOBL(JGOBL0 + 4) ='PMOB'
  2596.          TABOBL(JGOBL0 + 5) ='ERAI'
  2597.  
  2598.         ELSEIF (iplac.eq.15) then
  2599. *LI_LI_FR
  2600.          JGOBL0 = JGOBL
  2601.          JGOBL= JGOBL0 + 8
  2602.          TABOBL(JGOBL0 + 1) ='NORM'
  2603.          TABOBL(JGOBL0 + 2) ='LIMA'
  2604.          TABOBL(JGOBL0 + 3) ='LIES'
  2605.          TABOBL(JGOBL0 + 4) ='RAID'
  2606.          TABOBL(JGOBL0 + 5) ='GLIS'
  2607.          TABOBL(JGOBL0 + 6) ='ADHE'
  2608.          TABOBL(JGOBL0 + 7) ='RTAN'
  2609.          TABOBL(JGOBL0 + 8) ='ATAN'
  2610.  
  2611.          JGOBL0 = JGFAC
  2612.          JGFAC = JGOBL0 + 4
  2613.          TABFAC(JGOBL0+1) = 'JEU'
  2614.          TABFAC(JGOBL0+2) = 'AMOR'
  2615.          TABFAC(JGOBL0+3) = 'RECH'
  2616.          TABFAC(JGOBL0+4) = 'SYME'
  2617.  
  2618.         ELSEIF (iplac.eq.16) then
  2619. *LI_CE_FR
  2620.          JGOBL0 = JGOBL
  2621.          JGOBL= JGOBL0 + 8
  2622.          TABOBL(JGOBL0 + 1) ='NORM'
  2623.          TABOBL(JGOBL0 + 2) ='LIMA'
  2624.          TABOBL(JGOBL0 + 3) ='LIES'
  2625.          TABOBL(JGOBL0 + 4) ='RAID'
  2626.          TABOBL(JGOBL0 + 5) ='GLIS'
  2627.          TABOBL(JGOBL0 + 6) ='ADHE'
  2628.          TABOBL(JGOBL0 + 7) ='RTAN'
  2629.          TABOBL(JGOBL0 + 8) ='ATAN'
  2630.  
  2631.          JGOBL0 = JGFAC
  2632.          JGFAC = JGOBL0 + 5
  2633.          TABFAC(JGOBL0+1) = 'AMOR'
  2634.          TABFAC(JGOBL0+2) = 'RECH'
  2635.          TABFAC(JGOBL0+3) = 'RAYO'
  2636.          TABFAC(JGOBL0+4) = 'ACTN'
  2637.          TABFAC(JGOBL0+5) = 'INVE'
  2638.  
  2639.         ELSEIF (iplac.eq.17) then
  2640. *PA_FL_RO
  2641.          JGOBL0 = JGOBL
  2642.          JGOBL= JGOBL0+10
  2643.          TABOBL(JGOBL0 + 1) ='LONG'
  2644.          TABOBL(JGOBL0 + 2) ='RAYO'
  2645.          TABOBL(JGOBL0 + 3) ='VISC'
  2646.          TABOBL(JGOBL0 + 4) ='RHOF'
  2647.          TABOBL(JGOBL0 + 5) ='PADM'
  2648.          TABOBL(JGOBL0 + 6) ='VROT'
  2649.          TABOBL(JGOBL0 + 7) ='EPSI'
  2650.          TABOBL(JGOBL0 + 8) ='PHII'
  2651.          TABOBL(JGOBL0 + 9) ='AFFI'
  2652.          TABOBL(JGOBL0 + 10)='TLOB'
  2653.  
  2654.          JGOBL0 = JGFAC
  2655.          JGFAC = JGOBL0 + 1
  2656.          TABFAC(JGOBL0+1) = 'AMOR'
  2657.  
  2658.         ELSEIF (iplac.eq.23) then
  2659. *NEWMARK MODAL
  2660.          JGOBL0 = JGOBL
  2661.          JGOBL= JGOBL0+3
  2662.          TABOBL(JGOBL0 + 1) ='JEU'
  2663.          TABOBL(JGOBL0 + 2) ='MASS'
  2664.          TABOBL(JGOBL0 + 3) ='FREQ'
  2665.  
  2666.          JGOBL0 = JGFAC
  2667.          JGFAC = JGOBL0 + 3
  2668.          TABFAC(JGOBL0+1) = 'EXCE'
  2669.          TABFAC(JGOBL0+2) = 'FROT'
  2670.          TABFAC(JGOBL0+3) = 'MOFR'
  2671.         ENDIF
  2672.  
  2673. * 'SORT' facultatif dans tous les cas
  2674.         JGOBL0 = JGFAC
  2675.         JGFAC= JGOBL0+1
  2676.         TABFAC(JGOBL0+1) = 'SORT'
  2677.         GOTO 9999
  2678.       ENDIF
  2679.  
  2680. *--------------------------------------------------------------------
  2681. *     CAS DE LA FORMULATION ELECTROSTATIQUE
  2682. *--------------------------------------------------------------------
  2683.       CALL PLACE(FORMOD,NFOR,IELEC,'ELECTROSTATIQUE')
  2684.       IF (IELEC.NE.0) THEN
  2685. C -- Permittivite isotrope
  2686.         IF (IMATEE.EQ.1) THEN
  2687. C*      IF (CMATEE.EQ.'ISOTROPE') THEN
  2688. C*      IF (MATMOD(1).EQ.'ISOTROPE        ') THEN
  2689.           JGOBL = 1
  2690.           TABOBL(1)='PEL     '
  2691. C -- Permittivite orthotrope
  2692.         ELSEIF (IMATEE.EQ.2) THEN
  2693. C*      ELSEIF (CMATEE.EQ.'ORTHOTRO') THEN
  2694. C*      ELSEIF (MATMOD(1).EQ.'ORTHOTROPE      ') THEN
  2695. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2696.           IF (IDIM.EQ.2) THEN
  2697.             IF (IFOMOD.NE.1) THEN
  2698.               JGOBL = 4
  2699.               TABOBL(1) = 'PE1     '
  2700.               TABOBL(2) = 'PE2     '
  2701.               TABOBL(3) = 'V1X     '
  2702.               TABOBL(4) = 'V1Y     '
  2703.  
  2704. C ---- Elements massifs bidimensionnels FOURIER
  2705.             ELSE
  2706.               JGOBL = 5
  2707.               TABOBL(1) = 'PE1     '
  2708.               TABOBL(2) = 'PE2     '
  2709.               TABOBL(3) = 'PE3     '
  2710.               TABOBL(4) = 'V1X     '
  2711.               TABOBL(5) = 'V1Y     '
  2712.             ENDIF
  2713. C ---- Elements massifs TRIDimensionnels
  2714.           ELSEIF (IDIM.EQ.3) THEN
  2715.             JGOBL = 9
  2716.             TABOBL(1) = 'PE1     '
  2717.             TABOBL(2) = 'PE2     '
  2718.             TABOBL(3) = 'PE3     '
  2719.             TABOBL(4) = 'V1X     '
  2720.             TABOBL(5) = 'V1Y     '
  2721.             TABOBL(6) = 'V1Z     '
  2722.             TABOBL(7) = 'V2X     '
  2723.             TABOBL(8) = 'V2Y     '
  2724.             TABOBL(9) = 'V2Z     '
  2725.           ENDIF
  2726.  
  2727. C -- Permittivite anisotrope
  2728.         ELSEIF (IMATEE.EQ.3) THEN
  2729. C*      ELSEIF (CMATEE.EQ.'ANISOTRO') THEN
  2730. C*      ELSEIF (MATMOD(1).EQ.'ANISOTROPE      ') THEN
  2731. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2732.           IF (IDIM.EQ.2) THEN
  2733.             IF (IFOMOD.NE.1) THEN
  2734.               JGOBL = 5
  2735.               TABOBL(1) = 'PE11    '
  2736.               TABOBL(2) = 'PE22    '
  2737.               TABOBL(3) = 'PE21    '
  2738.               TABOBL(4) = 'V1X     '
  2739.               TABOBL(5) = 'V1Y     '
  2740.  
  2741. C ---- Elements massifs bidimensionnels FOURIER
  2742.             ELSE
  2743.               JGOBL = 6
  2744.               TABOBL(1) = 'PE11    '
  2745.               TABOBL(2) = 'PE22    '
  2746.               TABOBL(3) = 'PE21    '
  2747.               TABOBL(4) = 'PE33    '
  2748.               TABOBL(5) = 'V1X     '
  2749.               TABOBL(6) = 'V1Y     '
  2750.             ENDIF
  2751.  
  2752. C ---- Elements massifs TRIDimensionnels
  2753.           ELSEIF (IDIM.EQ.3) THEN
  2754.             JGOBL = 12
  2755.             TABOBL( 1) = 'PE11    '
  2756.             TABOBL( 2) = 'PE22    '
  2757.             TABOBL( 3) = 'PE33    '
  2758.             TABOBL( 4) = 'PE21    '
  2759.             TABOBL( 5) = 'PE31    '
  2760.             TABOBL( 6) = 'PE32    '
  2761.             TABOBL( 7) = 'V1X     '
  2762.             TABOBL( 8) = 'V1Y     '
  2763.             TABOBL( 9) = 'V1Z     '
  2764.             TABOBL(10) = 'V2X     '
  2765.             TABOBL(11) = 'V2Y     '
  2766.             TABOBL(12) = 'V2Z     '
  2767.           ENDIF
  2768.  
  2769.         ELSE
  2770.           IRET = 0
  2771.           CALL ERREUR(5)
  2772.         ENDIF
  2773.         GOTO 9999
  2774.       ENDIF
  2775.  
  2776. *-----------------------------------------------------------------------
  2777. *     CAS DE LA FORMULATION 'DIFFUSION'
  2778. *-----------------------------------------------------------------------
  2779.       CALL PLACE(FORMOD,NFOR,IDIFF,'DIFFUSION')
  2780.       IF (IDIFF.NE.0) THEN
  2781.         JGFAC    = 1
  2782.         TABFAC(1)='CDIF'
  2783.  
  2784. C-      Diffusion isotrope
  2785.         IF (IMATEE .EQ. 1) THEN
  2786.           JGOBL     = 1
  2787. C*8       TABOBL(1) = 'KDIFF   '
  2788.           TABOBL(1) = 'KD      '
  2789.  
  2790. C-      Diffusion orthotrope
  2791.         ELSEIF (IMATEE .EQ. 2) THEN
  2792. C         --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2793.           IF (IDIM.EQ.2) THEN
  2794.             IF (IFOMOD.NE.1) THEN
  2795.               JGOBL     = 4
  2796.               TABOBL(1) = 'KD1     '
  2797.               TABOBL(2) = 'KD2     '
  2798.               TABOBL(3) = 'V1X     '
  2799.               TABOBL(4) = 'V1Y     '
  2800.  
  2801. C           --- Elements massifs bidimensionnels FOURIER
  2802.             ELSE
  2803.               JGOBL     = 5
  2804.               TABOBL(1) = 'KD1     '
  2805.               TABOBL(2) = 'KD2     '
  2806.               TABOBL(3) = 'KD3     '
  2807.               TABOBL(4) = 'V1X     '
  2808.               TABOBL(5) = 'V1Y     '
  2809.             ENDIF
  2810.  
  2811. C         --- Elements massifs TRIDimensionnels
  2812.           ELSEIF (IDIM.EQ.3) THEN
  2813.             JGOBL     = 9
  2814.             TABOBL(1) = 'KD1     '
  2815.             TABOBL(2) = 'KD2     '
  2816.             TABOBL(3) = 'KD3     '
  2817.             TABOBL(4) = 'V1X     '
  2818.             TABOBL(5) = 'V1Y     '
  2819.             TABOBL(6) = 'V1Z     '
  2820.             TABOBL(7) = 'V2X     '
  2821.             TABOBL(8) = 'V2Y     '
  2822.             TABOBL(9) = 'V2Z     '
  2823.           ENDIF
  2824.  
  2825. C       - Diffusion anisotrope
  2826.         ELSEIF (IMATEE.EQ.3) THEN
  2827. C         --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2828.           IF (IDIM.EQ.2) THEN
  2829.             IF (IFOMOD.NE.1) THEN
  2830.               JGOBL     = 5
  2831.               TABOBL(1) = 'KD11    '
  2832.               TABOBL(2) = 'KD22    '
  2833.               TABOBL(3) = 'KD21    '
  2834.               TABOBL(4) = 'V1X     '
  2835.               TABOBL(5) = 'V1Y     '
  2836.  
  2837. C           --- Elements massifs bidimensionnels FOURIER
  2838.             ELSE
  2839.               JGOBL     = 6
  2840.               TABOBL(1) = 'KD11    '
  2841.               TABOBL(2) = 'KD22    '
  2842.               TABOBL(3) = 'KD21    '
  2843.               TABOBL(4) = 'KD33    '
  2844.               TABOBL(5) = 'V1X     '
  2845.               TABOBL(6) = 'V1Y     '
  2846.             ENDIF
  2847.  
  2848. C         --- Elements massifs TRIDimensionnels
  2849.           ELSEIF (IDIM.EQ.3) THEN
  2850.             JGOBL    = 12
  2851.             TABOBL( 1) = 'KD11    '
  2852.             TABOBL( 2) = 'KD22    '
  2853.             TABOBL( 3) = 'KD33    '
  2854.             TABOBL( 4) = 'KD21    '
  2855.             TABOBL( 5) = 'KD31    '
  2856.             TABOBL( 6) = 'KD32    '
  2857.             TABOBL( 7) = 'V1X     '
  2858.             TABOBL( 8) = 'V1Y     '
  2859.             TABOBL( 9) = 'V1Z     '
  2860.             TABOBL(10) = 'V2X     '
  2861.             TABOBL(11) = 'V2Y     '
  2862.             TABOBL(12) = 'V2Z     '
  2863.           ENDIF
  2864.         ENDIF
  2865.  
  2866.         JGOBL0 = JGOBL
  2867.         JGFAC0 = JGFAC
  2868.  
  2869. *        Loi FICK :
  2870. *       ------------
  2871.         IF (INATUU.EQ.0) THEN
  2872. *         Rien a ajouter
  2873.  
  2874. *        Loi SORET :
  2875. *       -------------
  2876.         ELSEIF (INATUU.EQ.1) THEN
  2877. C         -- Diffusion isotrope
  2878.           IF (IMATEE.EQ.1) THEN
  2879.             JGOBL            = JGOBL+1
  2880. C*8         TABOBL(JGOBL0+1) = 'KSORET  '
  2881.             TABOBL(JGOBL0+1) = 'KS      '
  2882.  
  2883. C         -- Diffusion orthotrope
  2884.           ELSEIF (IMATEE.EQ.2) THEN
  2885. C           ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2886.             IF (IFOMOD.LE.0) THEN
  2887.               JGOBL            = JGOBL + 2
  2888. C*8           TABOBL(JGOBL0+1) = 'KSORET1 '
  2889. C*8           TABOBL(JGOBL0+2) = 'KSORET2 '
  2890.               TABOBL(JGOBL0+1) = 'KS1     '
  2891.               TABOBL(JGOBL0+2) = 'KS2     '
  2892.  
  2893. C           ---- Elements massifs bidimensionnels FOURIER
  2894. C           ---- Elements massifs TRIDimensionnels
  2895.             ELSE
  2896.               JGOBL            = JGOBL + 3
  2897. C*8           TABOBL(JGOBL0+1) = 'KSORET1 '
  2898. C*8           TABOBL(JGOBL0+2) = 'KSORET2 '
  2899. C*8           TABOBL(JGOBL0+3) = 'KSORET3 '
  2900.               TABOBL(JGOBL0+1) = 'KS1     '
  2901.               TABOBL(JGOBL0+2) = 'KS2     '
  2902.               TABOBL(JGOBL0+3) = 'KS3     '
  2903.             ENDIF
  2904.  
  2905. C         -- Diffusion anisotrope
  2906.           ELSEIF (IMATEE.EQ.3) THEN
  2907. C           ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2908.             IF (IFOMOD.LE.0) THEN
  2909.               JGOBL            = JGOBL + 3
  2910. C*8           TABOBL(JGOBL0+1) = 'KSORET11'
  2911. C*8           TABOBL(JGOBL0+2) = 'KSORET22'
  2912. C*8           TABOBL(JGOBL0+3) = 'KSORET21'
  2913.               TABOBL(JGOBL0+1) = 'KS11    '
  2914.               TABOBL(JGOBL0+2) = 'KS22    '
  2915.               TABOBL(JGOBL0+3) = 'KS21    '
  2916.  
  2917. C           ---- Elements massifs bidimensionnels FOURIER
  2918.             ELSEIF (IFOMOD.EQ.1) THEN
  2919.               JGOBL            = JGOBL + 4
  2920. C*8           TABOBL(JGOBL0+1) = 'KSORET11'
  2921. C*8           TABOBL(JGOBL0+2) = 'KSORET22'
  2922. C*8           TABOBL(JGOBL0+3) = 'KSORET21'
  2923. C*8           TABOBL(JGOBL0+4) = 'KSORET33'
  2924.               TABOBL(JGOBL0+1) = 'KS11    '
  2925.               TABOBL(JGOBL0+2) = 'KS22    '
  2926.               TABOBL(JGOBL0+3) = 'KS21    '
  2927.               TABOBL(JGOBL0+4) = 'KS33    '
  2928.  
  2929. C           ---- Elements massifs TRIDimensionnels
  2930.             ELSEIF (IFOMOD.EQ.2) THEN
  2931.               JGOBL            = JGOBL + 6
  2932. C*8           TABOBL(JGOBL0+1) = 'KSORET11'
  2933. C*8           TABOBL(JGOBL0+2) = 'KSORET22'
  2934. C*8           TABOBL(JGOBL0+3) = 'KSORET33'
  2935. C*8           TABOBL(JGOBL0+4) = 'KSORET21'
  2936. C*8           TABOBL(JGOBL0+5) = 'KSORET31'
  2937. C*8           TABOBL(JGOBL0+6) = 'KSORET32'
  2938.               TABOBL(JGOBL0+1) = 'KS11    '
  2939.               TABOBL(JGOBL0+2) = 'KS22    '
  2940.               TABOBL(JGOBL0+3) = 'KS33    '
  2941.               TABOBL(JGOBL0+4) = 'KS21    '
  2942.               TABOBL(JGOBL0+5) = 'KS31    '
  2943.               TABOBL(JGOBL0+6) = 'KS32    '
  2944.             ENDIF
  2945.  
  2946.           ELSE
  2947. C         -- Erreur ne devant pas survenir
  2948.             CALL ERREUR(261)
  2949.             RETURN
  2950.           ENDIF
  2951.  
  2952. *        Loi UTILISATEUR :
  2953. *       -------------------
  2954. *       ELSEIF (INATUU.EQ.2) THEN
  2955. *       Par defaut pas de composantes autres que celles de l'utilisateur
  2956.  
  2957.         ELSEIF (INATUU.EQ.5) THEN
  2958. *       diffusion avec dissipation affine
  2959.           JGOBL         = JGOBL+2
  2960.           TABOBL(JGOBL0+1)='DAF0'
  2961.           TABOBL(JGOBL0+2)='DAF1'
  2962.  
  2963.         ELSEIF (INATUU.EQ.6) THEN
  2964. *       diffusion avec dissipation visqueuse
  2965.           JGOBL         = JGOBL+2
  2966.           TABOBL(JGOBL0+1)='DVQ0'
  2967.           TABOBL(JGOBL0+2)='DVQ1'
  2968.  
  2969. *        Loi inconnue : ERREUR
  2970. *       ----------------
  2971.         ELSE
  2972.           CALL ERREUR(261)
  2973.           RETURN
  2974.         ENDIF
  2975.  
  2976. *        Composantes ADVECTION :
  2977. *       -------------------
  2978.         CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
  2979.         if (iadve .ne. 0) then
  2980.  
  2981. C         Cas des Tuyaux 1D (MFR=79)
  2982. C         ===========================
  2983.           IF (MFR .EQ. 79) THEN
  2984.             JGOBL             = JGOBL + 1
  2985.             TABOBL(JGOBL0 + 1)='VITE'
  2986.  
  2987.             JGFAC             = JGFAC + 1
  2988.             TABFAC(JGFAC0 + 1)='DL'
  2989.  
  2990. C         Cas Diffusion Massif (MFR=73)
  2991. C         ============================
  2992.           ELSEIF (MFR .EQ. 73) THEN
  2993. C           Bidimensionnel PLAN DPGE
  2994. C           ===========================
  2995.             IF (IFOUR.EQ.-3) THEN
  2996.               JGOBL             = JGOBL + 3
  2997.               TABOBL(JGOBL0 + 1)='VITX'
  2998.               TABOBL(JGOBL0 + 2)='VITY'
  2999.               TABOBL(JGOBL0 + 3)='VITZ'
  3000.  
  3001.               JGFAC             = JGFAC + 3
  3002.               TABFAC(JGFAC0 + 1)='DX'
  3003.               TABFAC(JGFAC0 + 2)='DY'
  3004.               TABFAC(JGFAC0 + 3)='DZ'
  3005.  
  3006. C           Bidimensionnel PLAN (CP/DP)
  3007. C           ===========================
  3008.             ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
  3009.               JGOBL             = JGOBL + 2
  3010.               TABOBL(JGOBL0 + 1)='VITX'
  3011.               TABOBL(JGOBL0 + 2)='VITY'
  3012.  
  3013.               JGFAC             = JGFAC + 2
  3014.               TABFAC(JGFAC0 + 1)='DX'
  3015.               TABFAC(JGFAC0 + 2)='DY'
  3016.  
  3017. C           Axisymetrie
  3018. C           ===========================
  3019.             ELSEIF (IFOUR .EQ. 0) THEN
  3020.               JGOBL             = JGOBL + 2
  3021.               TABOBL(JGOBL0 + 1)='VITR'
  3022.               TABOBL(JGOBL0 + 2)='VITZ'
  3023.  
  3024.               JGFAC             = JGFAC + 2
  3025.               TABFAC(JGFAC0 + 1)='DR'
  3026.               TABFAC(JGFAC0 + 2)='DZ'
  3027.  
  3028. C           Fourier
  3029. C           ===========================
  3030.             ELSEIF (IFOUR .EQ. 1) THEN
  3031.               CALL ERREUR(21)
  3032.               RETURN
  3033. C             JGOBL             = JGOBL + 3
  3034. C             TABOBL(JGOBL0 + 1)='VITR'
  3035. C             TABOBL(JGOBL0 + 2)='VITZ'
  3036. C             TABOBL(JGOBL0 + 3)='VITT'
  3037. C
  3038. C             JGFAC             = JGFAC + 3
  3039. C             TABFAC(JGFAC0 + 1)='DR'
  3040. C             TABFAC(JGFAC0 + 2)='DZ'
  3041. C             TABFAC(JGFAC0 + 3)='DT'
  3042.  
  3043. C           Tridimensionnel
  3044. C           ===========================
  3045.             ELSEIF (IFOUR .EQ. 2) THEN
  3046.               JGOBL             = JGOBL + 3
  3047.               TABOBL(JGOBL0 + 1)='VITX'
  3048.               TABOBL(JGOBL0 + 2)='VITY'
  3049.               TABOBL(JGOBL0 + 3)='VITZ'
  3050.  
  3051.               JGFAC             = JGFAC + 3
  3052.               TABFAC(JGFAC0 + 1)='DX'
  3053.               TABFAC(JGFAC0 + 2)='DY'
  3054.               TABFAC(JGFAC0 + 3)='DZ'
  3055.  
  3056. C           Unidimensionnel (1D)
  3057. C           ===========================
  3058.             ELSEIF (IFOUR.GE.3 .AND. IFOUR.LE.15) THEN
  3059.               IF (IFOUR.LE.6) THEN
  3060.                 JGOBL             = JGOBL + 1
  3061.                 TABOBL(JGOBL0 + 1)='VITX'
  3062.  
  3063.               JGFAC             = JGFAC + 1
  3064.               TABFAC(JGFAC0 + 1)='DX'
  3065.  
  3066.               ELSEIF (IFOUR.EQ.7 .OR. IFOUR.EQ.8) THEN
  3067.                 JGOBL             = JGOBL + 2
  3068.                 TABOBL(JGOBL0 + 1)='VITX'
  3069.                 TABOBL(JGOBL0 + 2)='VITY'
  3070.  
  3071.               JGFAC             = JGFAC + 2
  3072.               TABFAC(JGFAC0 + 1)='DX'
  3073.               TABFAC(JGFAC0 + 2)='DY'
  3074.  
  3075.               ELSEIF (IFOUR.EQ.9 .OR. IFOUR.EQ.10) THEN
  3076.                 JGOBL             = JGOBL + 2
  3077.                 TABOBL(JGOBL0 + 1)='VITX'
  3078.                 TABOBL(JGOBL0 + 2)='VITZ'
  3079.  
  3080.               JGFAC             = JGFAC + 2
  3081.               TABFAC(JGFAC0 + 1)='DX'
  3082.               TABFAC(JGFAC0 + 2)='DZ'
  3083.  
  3084.               ELSEIF (IFOUR.EQ.11) THEN
  3085.                 JGOBL             = JGOBL + 3
  3086.                 TABOBL(JGOBL0 + 1)='VITX'
  3087.                 TABOBL(JGOBL0 + 2)='VITY'
  3088.                 TABOBL(JGOBL0 + 3)='VITZ'
  3089.  
  3090.                 JGFAC             = JGFAC + 3
  3091.                 TABFAC(JGFAC0 + 1)='DX'
  3092.                 TABFAC(JGFAC0 + 2)='DY'
  3093.                 TABFAC(JGFAC0 + 3)='DZ'
  3094.  
  3095.               ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN
  3096.                 JGOBL             = JGOBL + 1
  3097.                 TABOBL(JGOBL0 + 1)='VITR'
  3098.  
  3099.                 JGFAC             = JGFAC + 1
  3100.                 TABFAC(JGFAC0 + 1)='DR'
  3101.  
  3102.               ELSEIF (IFOUR.EQ.14) THEN
  3103.                 JGOBL             = JGOBL + 2
  3104.                 TABOBL(JGOBL0 + 1)='VITR'
  3105.                 TABOBL(JGOBL0 + 2)='VITZ'
  3106.  
  3107.                 JGFAC             = JGFAC + 2
  3108.                 TABFAC(JGFAC0 + 1)='DR'
  3109.                 TABFAC(JGFAC0 + 2)='DZ'
  3110.               ENDIF
  3111.             ENDIF
  3112.           ENDIF
  3113.         ENDIF
  3114.  
  3115. *       WRITE(6,*) 'TABOBL',JGOBL,(TABOBL(i),i=1,JGOBL)
  3116. *       WRITE(6,*) 'TABFAC',JGFAC,(TABFAC(i),i=1,JGFAC)
  3117.         GOTO 9999
  3118.       ENDIF
  3119.  
  3120. *-----------------------------------------------------------------------
  3121. *     CAS DE LA FORMULATION 'CHARGEMENT'
  3122. *-----------------------------------------------------------------------
  3123.       CALL PLACE(FORMOD,NFOR,ICHAR,'CHARGEMENT')
  3124.       IF (ICHAR.NE.0) THEN
  3125.         IPLAC = IMATEE
  3126.         IF (IPLAC.EQ.1) THEN
  3127. C         RAJOUTER DIFFERENTIATION EN FONCIONS DE LA FORMULATION MASSIF -> P
  3128. C         COQUE PINF PSUP
  3129.           JGOBL=1
  3130.           TABOBL(1)='PR  '
  3131.         ELSE
  3132.           IRET = 0
  3133.           CALL ERREUR(5)
  3134.         ENDIF
  3135.         GOTO 9999
  3136.       ENDIF
  3137.  
  3138. *--------------------------------------------------------------------
  3139.  9999 CONTINUE
  3140.  
  3141.       IF(IMECA.GT.0 .OR. IPORE.GT.0)THEN
  3142. C       CB215821 : Ajout de TREF et TALP à la fin des compostantes facultatives
  3143. C                  pour les FORMULATION MECANIQUE, POREUX
  3144. C                  Pour ne pas casser l'ordre dans COMP (com2/coml6/comval/comara/...)
  3145.         JGFAC=JGFAC+2
  3146.         TABFAC(JGFAC-1) ='TREF'
  3147.         TABFAC(JGFAC)   ='TALP'
  3148.       ENDIF
  3149.  
  3150. C     Erreur si JGOBL ou JGFAC sont superieurs a ITA
  3151. C    (Passage en FORTRAN 77 car la compilation depasse la memoire sur
  3152. C     Windows-32bits)
  3153.       IF ((JGOBL .GT. ITA) .OR. (JGFAC .GT. ITA)) THEN
  3154.         IRET = 0
  3155.         CALL ERREUR(5)
  3156.         RETURN
  3157.       ENDIF
  3158.  
  3159.       NBROBL = JGOBL
  3160.       NBRFAC = JGFAC
  3161.       NOMID  = 0
  3162.  
  3163. *     SI PROBLEME (IRET = 0), ON SORT AVEC IPNOMC A 0
  3164. *
  3165.       IF (IRET.NE.0) THEN
  3166.         SEGINI,NOMID
  3167.         DO 100 IO = 1,NBROBL
  3168.           NOMID.LESOBL(IO) = TABOBL(IO)
  3169.  100    CONTINUE
  3170.  
  3171.         DO 110 IO=1,NBRFAC
  3172.           NOMID.LESFAC(IO) = TABFAC(IO)
  3173. 110     CONTINUE
  3174.         SEGACT,NOMID*NOMOD
  3175.       ENDIF
  3176.  
  3177.       IPNOMC = NOMID
  3178.  
  3179.       END
  3180.  
  3181.  
  3182.  
  3183.  

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