Télécharger idgrad.eso

Retour à la liste

idgrad
C IDGRAD    SOURCE    MB234859  25/08/04    21:15:17     12339          
 
C=======================================================================
C=         RECHERCHE DES NOMS DE COMPOSANTES DE GRADIENTS              =
C=         ----------------------------------------------              =
C=                                                                     =
C=  Entrees :                                                          =
C=    IPMODE    Pointeur sur un MMODEL.KMODEL                          =
C=    IFOU      valeur de IFOUR de CCOPTIO                             =
C=                                                                     =
C=  Sorties :                                                          =
C=    IPNOMC    pointeur de type NOMID sur les listes de noms de       =
C=              composantes OBLigatoires et FACultatives               =
C=    NBROBL    nombre de composantes OBLigatoires                     =
C=    NBRFAC    nombre de composantes FACultatives                     =
C=                                                                     =
C=  Remarque : Voir INOMID ou CCOPTIO pour signification IFOUR         =
C=======================================================================
C
      SUBROUTINE IDGRAD(IPMODE,IFOU,IPNOMC,NBROBL,NBRFAC)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
-INC PPARAM
-INC CCOPTIO
C==DEB= FORMULATION HHO == Include specifique ==========================
-INC CCHHOPA
C==FIN= FORMULATION HHO ================================================
-INC SMLMOTS
-INC SMMODEL
C
      CHARACTER*(LOCOMP) CCOMP
      EXTERNAL LONG
C
      NBROBL=0
      NBRFAC=0
C
      IMODEL=IPMODE
      NOMID =IMODEL.LNOMID(3)
C
C     Ne pas recreer le NOMID si deja present
      IF (NOMID.NE.0) THEN
        NBROBL=LESOBL(/2)
        NBRFAC=LESFAC(/2)
        IPNOMC=NOMID
        RETURN
      ENDIF
C
      MELE=IMODEL.NEFMOD
      MFR =NUMMFR(MELE)
      MFR2=NUMFOR(IMODEL)
C     =================================================================
C     FORMULATION THERMOHYDRIQUE
C     =================================================================
      IF (MFR2.EQ.65) THEN
        IF (IFOUR.EQ.2) THEN
          NBROBL=9
          SEGINI,NOMID
          LESOBL(1)='PG,X    '
          LESOBL(2)='PG,Y    '
          LESOBL(3)='PG,Z    '
          LESOBL(4)='PC,X    '
          LESOBL(5)='PC,Y    '
          LESOBL(6)='PC,Z    '
          LESOBL(7)='T,X     '
          LESOBL(8)='T,Y     '
          LESOBL(9)='T,Z     '
        ELSE IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='PG,X    '
          LESOBL(2)='PG,Y    '
          LESOBL(3)='PC,X    '
          LESOBL(4)='PC,Y    '
          LESOBL(5)='T,X     '
          LESOBL(6)='T,Y     '
        ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='PG,R    '
          LESOBL(2)='PG,Z    '
          LESOBL(3)='PC,R    '
          LESOBL(4)='PC,Z    '
          LESOBL(5)='T,R     '
          LESOBL(6)='T,Z     '
        ENDIF
C     =================================================================
C     FORMULATION MECANIQUE/POREUX/MELANGE/NAVIER_STOKES/HHO
C     =================================================================
      ELSE IF (MFR2.EQ.2.OR.MFR2.EQ.33.OR.MFR2.EQ.38.OR.MFR2.EQ.52
     &    .OR. MFR2.EQ.HHO_MFR_ELEMENT) THEN
C
C       Elements MASSIFS (standard,incompressibles,XFEM,Navier_Stokes,HHO)
C       ================
        IF (MFR.EQ.1 .OR. MFR.EQ.31 .OR. MFR.EQ.63 .OR. MFR.EQ.52 .OR.
     &      MFR.EQ.HHO_MFR_ELEMENT) THEN
          IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.
     &        IFOUR.EQ.2) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UX,X    '
            LESOBL(2)='UX,Y    '
            LESOBL(3)='UX,Z    '
            LESOBL(4)='UY,X    '
            LESOBL(5)='UY,Y    '
            LESOBL(6)='UY,Z    '
            LESOBL(7)='UZ,X    '
            LESOBL(8)='UZ,Y    '
            LESOBL(9)='UZ,Z    '
          ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UR,R    '
            LESOBL(2)='UR,Z    '
            LESOBL(3)='UR,T    '
            LESOBL(4)='UZ,R    '
            LESOBL(5)='UZ,Z    '
            LESOBL(6)='UZ,T    '
            LESOBL(7)='UT,R    '
            LESOBL(8)='UT,Z    '
            LESOBL(9)='UT,T    '
          ELSE IF (MFR.EQ.1 .AND. (IFOUR.GE.3.AND.IFOUR.LE.15)) THEN
            NBROBL=9
            SEGINI,NOMID
            IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN
              LESOBL(1)='UX,X    '
              LESOBL(2)='UX,Y    '
              LESOBL(3)='UX,Z    '
              LESOBL(4)='UY,X    '
              LESOBL(5)='UY,Y    '
              LESOBL(6)='UY,Z    '
              LESOBL(7)='UZ,X    '
              LESOBL(8)='UZ,Y    '
              LESOBL(9)='UZ,Z    '
            ELSE
              LESOBL(1)='UR,R    '
              LESOBL(2)='UR,Z    '
              LESOBL(3)='UR,T    '
              LESOBL(4)='UZ,R    '
              LESOBL(5)='UZ,Z    '
              LESOBL(6)='UZ,T    '
              LESOBL(7)='UT,R    '
              LESOBL(8)='UT,Z    '
              LESOBL(9)='UT,T    '
            ENDIF
          ENDIF
C
C       Elements COQUES
C       ===============
        ELSE IF (MFR.EQ.3) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UX,X'
            LESOBL(2)='UX,Y'
            LESOBL(3)='UX,Z'
            LESOBL(4)='UY,X'
            LESOBL(5)='UY,Y'
            LESOBL(6)='UY,Z'
            LESOBL(7)='UZ,X'
            LESOBL(8)='UZ,Y'
            LESOBL(9)='UZ,Z'
          ELSE IF (IFOUR.EQ.1.OR.IFOUR.EQ.0.OR.
     .             IFOUR.EQ.-1.OR.IFOUR.EQ.-2.or.IFOUR.eq.-3) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UU,S'
            LESOBL(2)='UU,T'
            LESOBL(3)='UU,N'
            LESOBL(4)='UV,S'
            LESOBL(5)='UV,T'
            LESOBL(6)='UV,N'
            LESOBL(7)='UW,S'
            LESOBL(8)='UW,T'
            LESOBL(9)='UW,N'
          ENDIF
C
C       Elements COQUES EPAISSES (COQ6,COQ8)
C       ========================
        ELSE IF (MFR.EQ.5) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UX,X'
            LESOBL(2)='UX,Y'
            LESOBL(3)='UX,Z'
            LESOBL(4)='UY,X'
            LESOBL(5)='UY,Y'
            LESOBL(6)='UY,Z'
            LESOBL(7)='UZ,X'
            LESOBL(8)='UZ,Y'
            LESOBL(9)='UZ,Z'
          ELSE IF (IFOUR.EQ.1.OR.IFOUR.EQ.0)  THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UR,R'
            LESOBL(2)='UR,Z'
            LESOBL(3)='UR,T'
            LESOBL(4)='UZ,R'
            LESOBL(5)='UZ,Z'
            LESOBL(6)='UZ,T'
            LESOBL(7)='UT,R'
            LESOBL(8)='UT,Z'
            LESOBL(9)='UT,T'
          ENDIF
C
C       Element de COQUE en CISAILLEMENT TRANSVERSE (tridimensionnel)
C       ===========================================
        ELSE IF (MFR.EQ.9) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UX,X'
            LESOBL(2)='UX,Y'
            LESOBL(3)='UX,Z'
            LESOBL(4)='UY,X'
            LESOBL(5)='UY,Y'
            LESOBL(6)='UY,Z'
            LESOBL(7)='UZ,X'
            LESOBL(8)='UZ,Y'
            LESOBL(9)='UZ,Z'
          ENDIF
C
C       Elements POUTRES ET TUYAUX
C       ==========================
        ELSE IF (MFR.EQ.7 .OR. MFR.EQ.13) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=3
            SEGINI,NOMID
            LESOBL(1)='U,S '
            LESOBL(2)='V,S '
            LESOBL(3)='W,S '
          ELSE
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)='U,S '
            LESOBL(2)='V,S '
          ENDIF
C
C       Element JOINT JOI1
C       =================
        ELSE IF (MFR.EQ.75) THEN
          IF (IFOUR.EQ.2 .OR. IFOUR.EQ.-3) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='UX,X    '
            LESOBL(2)='UX,Y    '
            LESOBL(3)='UX,Z    '
            LESOBL(4)='UY,X    '
            LESOBL(5)='UY,Y    '
            LESOBL(6)='UY,Z    '
            LESOBL(7)='UZ,X    '
            LESOBL(8)='UZ,Y    '
            LESOBL(9)='UZ,Z    '
          ELSE IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='UX,X    '
            LESOBL(2)='UX,Y    '
            LESOBL(3)='UY,X    '
            LESOBL(4)='UY,Y    '
          ELSE IF (IFOUR.EQ.0) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='UR,R    '
            LESOBL(2)='UR,Z    '
            LESOBL(3)='UZ,R    '
            LESOBL(4)='UZ,Z    '
           ENDIF
C
C       Element JOINT JOI2
C       =================
        ELSEIF (MFR.EQ.35) THEN
          IF (IFOUR.EQ.2 .OR. IFOUR.EQ.-3) THEN
            NBROBL=3
            SEGINI,NOMID
            LESOBL(1)='U1,N    '
            LESOBL(2)='U2,N    '
            LESOBL(3)='UN,N    '
          ELSE IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)='US,N    '
            LESOBL(2)='UN,N    '
          ENDIF
C
C       Element XFEM meca rupt
C       ================================
c       ELSE IF (MFR.EQ.63) THEN
CTY       Element XQ4R (xfem meca rupture en 2D)
c         IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) then
c           NBROBL=4
c           NBRFAC=36
c           NBRFAC=20
c           SEGINI,NOMID
c           LESOBL(1)='UX,X'
c           LESOBL(2)='UX,Y'
c           LESOBL(3)='UY,X'
c           LESOBL(4)='UY,Y'
c           LESFAC(1)='AX,X'
c           LESFAC(2)='AX,Y'
c           LESFAC(3)='AY,X'
c           LESFAC(4)='AY,Y'
c           LESFAC(5)='B1XX'
c           LESFAC(6)='B1XY'
c           LESFAC(7)='B1YX'
c           LESFAC(8)='B1YY'
c           LESFAC(9)='C1XX'
c           LESFAC(10)='C1XY'
c           LESFAC(11)='C1YX'
c           LESFAC(12)='C1YY'
c           LESFAC(13)='D1XX'
c           LESFAC(14)='D1XY'
c           LESFAC(15)='D1YX'
c           LESFAC(16)='D1YY'
c           LESFAC(17)='E1XX'
c           LESFAC(18)='E1XY'
c           LESFAC(19)='E1YX'
c           LESFAC(20)='E1YY'
c           LESFAC(21)='B2XX'
c           LESFAC(22)='B2XY'
c           LESFAC(23)='B2YX'
c           LESFAC(24)='B2YY'
c           LESFAC(25)='C2XX'
c           LESFAC(26)='C2XY'
c           LESFAC(27)='C2YX'
c           LESFAC(28)='C2YY'
c           LESFAC(29)='D2XX'
c           LESFAC(30)='D2XY'
c           LESFAC(31)='D2YX'
c           LESFAC(32)='D2YY'
c           LESFAC(33)='E2XX'
c           LESFAC(34)='E2XY'
c           LESFAC(35)='E2YX'
c           LESFAC(36)='E2YY'
CTY       Element XC8R (xfem meca rupture en 3D)
c         ELSE IF (IFOUR.EQ.2) THEN
c           NBROBL=9
c           NBRFAC=45
c           NBRFAC=81
c           SEGINI,NOMID
c           LESOBL(1)='UX,X'
c           LESOBL(2)='UX,Y'
c           LESOBL(3)='UX,Z'
c           LESOBL(4)='UY,X'
c           LESOBL(5)='UY,Y'
c           LESOBL(6)='UY,Z'
c           LESOBL(7)='UZ,X'
c           LESOBL(8)='UZ,Y'
c           LESOBL(9)='UZ,Z'
c           LESFAC(1)='AX,X'
c           LESFAC(2)='AX,Y'
c           LESFAC(3)='AX,Z'
c           LESFAC(4)='AY,X'
c           LESFAC(5)='AY,Y'
c           LESFAC(6)='AY,Z'
c           LESFAC(7)='AZ,X'
c           LESFAC(8)='AZ,Y'
c           LESFAC(9)='AZ,Z'
c           LESFAC(10)='B1XX'
c           LESFAC(11)='B1XY'
c           LESFAC(12)='B1XZ'
c           LESFAC(13)='B1YX'
c           LESFAC(14)='B1YY'
c           LESFAC(15)='B1YZ'
c           LESFAC(16)='B1ZX'
c           LESFAC(17)='B1ZY'
c           LESFAC(18)='B1ZZ'
c           LESFAC(19)='C1XX'
c           LESFAC(20)='C1XY'
c           LESFAC(21)='C1XZ'
c           LESFAC(22)='C1YX'
c           LESFAC(23)='C1YY'
c           LESFAC(24)='C1YZ'
c           LESFAC(25)='C1ZX'
c           LESFAC(26)='C1ZY'
c           LESFAC(27)='C1ZZ'
c           LESFAC(28)='D1XX'
c           LESFAC(29)='D1XY'
c           LESFAC(30)='D1XZ'
c           LESFAC(31)='D1YX'
c           LESFAC(32)='D1YY'
c           LESFAC(33)='D1YZ'
c           LESFAC(34)='D1ZX'
c           LESFAC(35)='D1ZY'
c           LESFAC(36)='D1ZZ'
c           LESFAC(37)='E1XX'
c           LESFAC(38)='E1XY'
c           LESFAC(39)='E1XZ'
c           LESFAC(40)='E1YX'
c           LESFAC(41)='E1YY'
c           LESFAC(42)='E1YZ'
c           LESFAC(43)='E1ZX'
c           LESFAC(44)='E1ZY'
c           LESFAC(45)='E1ZZ'
c           LESFAC(46)='B2XX'
c           LESFAC(47)='B2XY'
c           LESFAC(48)='B2XZ'
c           LESFAC(49)='B2YX'
c           LESFAC(50)='B2YY'
c           LESFAC(51)='B2YZ'
c           LESFAC(52)='B2ZX'
c           LESFAC(53)='B2ZY'
c           LESFAC(54)='B2ZZ'
c           LESFAC(55)='C2XX'
c           LESFAC(56)='C2XY'
c           LESFAC(57)='C2XZ'
c           LESFAC(58)='C2YX'
c           LESFAC(59)='C2YY'
c           LESFAC(60)='C2YZ'
c           LESFAC(61)='C2ZX'
c           LESFAC(62)='C2ZY'
c           LESFAC(63)='C2ZZ'
c           LESFAC(64)='D2XX'
c           LESFAC(65)='D2XY'
c           LESFAC(66)='D2XZ'
c           LESFAC(67)='D2YX'
c           LESFAC(68)='D2YY'
c           LESFAC(69)='D2YZ'
c           LESFAC(70)='D2ZX'
c           LESFAC(71)='D2ZY'
c           LESFAC(72)='D2ZZ'
c           LESFAC(73)='E2XX'
c           LESFAC(74)='E2XY'
c           LESFAC(75)='E2XZ'
c           LESFAC(76)='E2YX'
c           LESFAC(77)='E2YY'
c           LESFAC(78)='E2YZ'
c           LESFAC(79)='E2ZX'
c           LESFAC(80)='E2ZY'
c           LESFAC(81)='E2ZZ'
c         ENDIF
C
C       Elements ZONE COHESIVE
C       ==========================
        ELSE IF (MFR.EQ.77) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='AX,X    '
            LESOBL(2)='AX,Y    '
            LESOBL(3)='AX,Z    '
            LESOBL(4)='AY,X    '
            LESOBL(5)='AY,Y    '
            LESOBL(6)='AY,Z    '
            LESOBL(7)='AZ,X    '
            LESOBL(8)='AZ,Y    '
            LESOBL(9)='AZ,Z    '
          ELSE IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='AX,X    '
            LESOBL(2)='AX,Y    '
            LESOBL(3)='AY,X    '
            LESOBL(4)='AY,Y    '
          ENDIF
        ENDIF
C     =================================================================
C     FORMULATION THERMIQUE
C     =================================================================
      ELSE IF (MFR2.EQ.29) THEN
C
C       Elements MASSIFS
C       ================
        IF (MFR.EQ.1) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=3
            SEGINI,NOMID
            LESOBL(1)='T,X '
            LESOBL(2)='T,Y '
            LESOBL(3)='T,Z '
          ELSE IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)='T,X '
            LESOBL(2)='T,Y '
          ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)='T,R '
            LESOBL(2)='T,Z '
          ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN
            NBROBL=1
            SEGINI,NOMID
            LESOBL(1)='T,X '
          ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN
            NBROBL=1
            SEGINI,NOMID
            LESOBL(1)='T,R '
          ENDIF
C
C       Elements COQUES
C       ===============
        ELSE IF (MFR.EQ.3 .OR. MFR.EQ. 5 .OR. MFR.EQ.9) THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)='T,X '
            LESOBL(2)='T,Y '
          ELSE IF (IFOUR.EQ.0) THEN
            NBROBL=1
            SEGINI,NOMID
            LESOBL(1)='T,Z '
          ENDIF
C
C       Elements TUY2 TUY3 BARR
C       =======================
        ELSE IF ((MFR.EQ.27) .OR. (MFR.EQ.79)) THEN
          NBROBL=1
          SEGINI,NOMID
          LESOBL(1)='T,S     '
        ENDIF
C     =================================================================
C     FORMULATION MAGNETODYNAMIQUE COQUE TRIDIMENSIONNEL
C     =================================================================
      ELSE IF (MFR2.EQ.69) THEN
        NBROBL=3
        SEGINI,NOMID
        LESOBL(1)='FC,X '
        LESOBL(2)='FC,Y '
        LESOBL(3)='FC,Z '
C     =================================================================
C     FORMULATION ELECTROSTATIQUE 
C     =================================================================
      ELSE IF (MFR2.EQ.71) THEN
        IF (IFOUR.EQ.2) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='VE,X    '
          LESOBL(2)='VE,Y    '
          LESOBL(3)='VE,Z    '
        ELSE IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='VE,X    '
          LESOBL(2)='VE,Y    '
        ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='VE,R    '
          LESOBL(2)='VE,Z    '
        ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN
          NBROBL=1
          SEGINI,NOMID
          LESOBL(1)='VE,X    '
        ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN
          NBROBL=1
          SEGINI,NOMID
          LESOBL(1)='VE,R    '
        ENDIF
C     =================================================================
C     FORMULATION DIFFUSION
C     =================================================================
      ELSE IF (MFR2.EQ.73) THEN
        MLMOT1=IMODEL.IVAMOD(1)
        CCOMP =MLMOT1.MOTS(1)
        NBCHAR=LONG(CCOMP)
C
C       Elements MASSIFS
C       ================
        IF (MFR.EQ.1)THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=3
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',X'
            LESOBL(2)=CCOMP(1:NBCHAR)//',Y'
            LESOBL(3)=CCOMP(1:NBCHAR)//',Z'
          ELSE IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',X'
            LESOBL(2)=CCOMP(1:NBCHAR)//',Y'
          ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',R'
            LESOBL(2)=CCOMP(1:NBCHAR)//',Z'
          ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN
            NBROBL=1
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',X'
          ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN
            NBROBL=1
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',R'
          ELSE
            CALL ERREUR(5)
            RETURN
          ENDIF
C
C       Elements TUY2 TUY3 BARR
C       =======================
        ELSE IF (MFR.EQ.27 .OR. MFR.EQ.79)THEN
          NBROBL=1
          SEGINI,NOMID
          LESOBL(1)=CCOMP(1:NBCHAR)//',S'
C
C       Elements COQUES
C       ===============
        ELSE IF (MFR.EQ.3 .OR. MFR.EQ.5 .OR. MFR.EQ.9)THEN
          IF (IFOUR.EQ.2) THEN
            NBROBL=2
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',X '
            LESOBL(2)=CCOMP(1:NBCHAR)//',Y '
          ELSE IF (IFOUR.EQ.0) THEN
            NBROBL=1
            SEGINI,NOMID
            LESOBL(1)=CCOMP(1:NBCHAR)//',Z '
          ELSE
            CALL ERREUR(5)
            RETURN
          ENDIF
 
        ENDIF
C     =================================================================
      ENDIF
 
      IF (NOMID.NE.0) SEGACT,NOMID*NOMOD
      IPNOMC=NOMID
      END
 
 

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