Esope source files

   1 : C FUCHXR    SOURCE    CHAT      11/03/16    21:24:00     6902           
   2 : 
   3 :       SUBROUTINE FUCHXR (IMODLA,IMODLB, IMODLS)
   4 : 
   5 : C-----------------------------------------------------------------------
   6 : C     SOUS-PROGRAMME REALISANT LA FUSION DES MCHAMLS D'ENRICHISSEMENT 
   7 : C     ISSUS DE DEUX MODELES CONTENANT DES ELEMENTS DE TYPE XFEM/RUPTURE
   8 : C-----------------------------------------------------------------------
   9 : C     IMODLA   1er modele elementaire contenant des elements XFEM
  10 : C              Segment IMODEL : ACTIF en Entree/Sortie
  11 : C     IMODLB   2e modele elementaire contenant des elements XFEM
  12 : C              Segment IMODEL : ACTIF en Entree/Sortie
  13 : C     IMODLS   Modele elementaire contenant le MCHAML d'ENRICHISSMENT
  14 : C              resultant de la fusion des MCHAMLs de IMODLA et IMODLB
  15 : C              Segment IMODEL : ACTIF en Sortie
  16 : C-----------------------------------------------------------------------
  17 : C     ON SUPPOSE QUE LES MAILLAGES ASSOCIES AUX SOUS-MODELES ONT ETE 
  18 : C     FUSIONNES ET QUE LE RESULTAT EST STOCKE DANS IMODLA.IMAMOD.
  19 : C-----------------------------------------------------------------------
  20 : 
  21 :       IMPLICIT INTEGER (I-N)
  22 :       IMPLICIT REAL*8 (A-H,O-Z)
  23 : 
  24 : -INC CCOPTIO
  25 : 
  26 : C*-INC SMELEME
  27 : -INC SMMODEL
  28 : 
  29 :       IMODE1 = IMODLA
  30 :       IMODE2 = IMODLB
  31 : C*    SEGACT,IMODE1,IMODE2
  32 : C- Recuperation du MCHAML d enrichissement dans IMODE1
  32 :       CALL ICHAXR(IMODE1, ICHAX1,MCHAM1)
  33 :       IF (ICHAX1 .EQ. 0) GOTO 900
  34 : C- Recuperation du MCHAML d enrichissement dans IMODE2
  34 :       CALL ICHAXR(IMODE2, ICHAX2,MCHAM2)
  35 :       IF (ICHAX2 .EQ. 0) GOTO 900
  36 : 
  37 :       IF (IMODLS.GT.0) THEN
  38 :         IMODEL = IMODLS
  39 :         SEGACT,IMODEL*MOD
  40 : C- Recuperation du MCHAML d enrichissement dans IMODEL
  40 :         CALL ICHAXR(IMODEL, ICHAXS,MCHAMS)
  41 :         IF (ICHAXS .EQ. 0) GOTO 900
  42 :       ELSE
  43 :         IMODEL = IMODLA
  44 :         ICHAXS = ICHAX1
  45 :       ENDIF
  46 : 
  47 : C- Fusion des MCHAMLs d'enrichissement
  47 :       CALL FUSCHL(MCHAM1,MCHAM2, MCHAMS)
  48 :       IF (IERR.NE.0) GOTO 900
  49 :       IMODEL.IVAMOD(ICHAXS) = MCHAMS
  50 : C*C- Reactivation du maillage support du modele suite a fusion
  51 : C*      MELEME = IMODEL.IMAMOD
  52 : C*      SEGACT,MELEME
  53 : C* Pour l'instant inutile de reactiver le maillage support car il 
  54 : C* correspond deja au maillage fusionne (equivalent a celui du
  55 : C* MCHAML fusionne ci-dessus mais de pointeur different).
  56 : C* Le MCHAML d'enrichissement fusionne sera mis en correspondance 
  57 : C* avec le maillage fusionne a la fin (appel a PARTXR).
  58 : 
  59 :  900  CONTINUE
  60 : C*    SEGDES,IMODE1,IMODE2
  61 : C*    IF (IMODLS.GT.0) SEGDES,IMODEL
  62 : 
  63 :       RETURN
  64 :       END
  65 : 
  66 :  
  67 :