Sources Esope | Cast3M

eximod
C EXIMOD    SOURCE    OF166741  24/12/18    21:15:10     12092          
 
C----------------------------------------------------------------------
C  ENTREE :
C         ICOLAC  POINTEUR SUR LE CHAPEAU DES PILES
C         IMODEL  OBJET IMODEL A EXAMINER
C         IIICHA  =1 POUR CHANGER LES POINTEURS
C-----------------------------------------------------------------------
C  REMARQUE : ICOLAC EST UN SEGMENT ACTIF(*MOD)  EN ENTREE ET SORTIE
C             PAS DE CHANGEMENT DE STATUT AU COURS DU SP 
C=======================================================================
C  TABLEAU KCOLA : VOIR LE SOUS-PROGRAMME TYPFIL
C=======================================================================
 
      SUBROUTINE EXIMOD(ICOLAC,IMODEL,IIICHA,NIVEAU)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
C==DEB= FORMULATION HHO == INCLUDE =====================================
-INC CCHHOPA
C==FIN= FORMULATION HHO ================================================
 
-INC SMMODEL
 
-INC TMCOLAC
 
      CHARACTER*(8) ITYPE
 
      i_1 = 1
 
      ILISSE = icolac.ILISSF
c*      SEGACT,ILISSE*MOD
      ILISSE = icolac.ILISSG
c*      SEGACT,ILISSE*MOD
 
      ICO1  = icolac.KCOLA( 1)
      ICO10 = icolac.KCOLA(10)
      ICO29 = icolac.KCOLA(29)
      ICO40 = icolac.KCOLA(40)
 
      SEGACT,IMODEL*MOD
 
      IVA = imodel.IMAMOD
      IF (IVA.GT.0) THEN
        CALL AJOUN(ICO1,IVA,ILISSE,i_1)
        IF (IIICHA.EQ.1) imodel.IMAMOD = -IVA
      ENDIF
      IVA = imodel.IPDPGE
      IF (IVA.GT.0) THEN
        CALL AJOUN(ICO1,IVA,ILISSE,i_1)
        IF (IIICHA.EQ.1) imodel.IPDPGE = -IVA
      ENDIF
      NFOR = imodel.FORMOD(/2)
      NM3  = imodel.INFMOD(/1)
C Cas 'NAVIER_STOKES' : INFMOD(2) contient une table
      IF (NFOR.GT.0 .AND. NM3.GT.1) THEN
        IF ((imodel.FORMOD(1).EQ.'NAVIER_STOKES   ') .OR.
     &      (imodel.FORMOD(1).EQ.'DARCY           ') .OR.
     &      (imodel.FORMOD(1).EQ.'EULER           ')) THEN
          IVA = imodel.INFMOD(2)
          IF (IVA.GT.0) THEN
            CALL AJOUN(ICO10,IVA,ILISSE,i_1)
            IF (IIICHA.EQ.1) imodel.INFMOD(2) = -IVA
          ENDIF
        ENDIF
      ENDIF
      DO io = 3, NM3
        IVA = imodel.INFMOD(io)
        IF (IVA.GT.0) THEN
          IF (io.EQ.14) THEN
            CALL AJOUN(ICO29,IVA,ilisse,i_1)
          ELSE
            CALL AJOUN(ICO40,IVA,ilisse,i_1)
          ENDIF
          IF (IIICHA.EQ.1) imodel.INFMOD(io) = -IVA
        ENDIF
      ENDDO
      ntym = imodel.TYMODE(/2)
      DO io = 1, ntym
        IVA = imodel.IVAMOD(io)
        IF (IVA.LT.0) GOTO 17
        ITYPE = imodel.TYMODE(io)
        j = 0
        CALL TYPFIL(ITYPE,j)
        IF (j.LE.0) THEN 
          MOTERR(1:8) = ITYPE
          CALL ERREUR(336)                   
          GOTO 17
        ELSE IF (j.EQ.32) THEN
          GOTO 17
        ENDIF
        icoj = icolac.KCOLA(j)
        NUMLIS = 1
        ilissd = icolac.ILISSG
        IF (j.EQ.24) NUMLIS = 6
        IF (j.EQ.25) THEN
          NUMLIS = 3
          ilissd = icolac.ILISSF
        ENDIF
        IF (j.EQ.26) NUMLIS = 2
        IF (j.EQ.27) NUMLIS = 5
c-        IF (j.EQ.32) THEN
c-          NUMLIS = 3
c-          ilissd = icolac.ILISSP
c-        ENDIF
        IF (j.EQ.36) NUMLIS = 7
        IF (j.EQ.45) NUMLIS = 5
c*dbg    if (j.eq.51 .or. j.eq.38)
c*dbg   &  write(ioimp,*) 'eximod (51)',iel,imodel,'tymode',io,itype,j
        CALL AJOUN(icoj,IVA,ilissd,NUMLIS)
        IF (IIICHA.EQ.1) imodel.IVAMOD(io) = -IVA
 
 17     CONTINUE
      ENDDO
 
C==DEB= FORMULATION HHO == Sauvegarde des MAILLAGES HHO ================
      IF (imodel.NEFMOD .EQ. HHO_NUM_ELEMENT) THEN
c-dbg        CALL HHOPIL(2,i_1,i_1)
        CALL HHOPIL(3,ICOLAC,IIICHA)
      ENDIF
C==FIN= FORMULATION HHO ================================================
 
      SEGDES,IMODEL
 
      RETURN
      END