C MUTU      SOURCE    OF166741  24/10/23    21:15:04     12046          

*_______________________________________________________________________
*
*            OPERATEUR MUTU
*
*  ELEMENTS FINIS :
*  ==============
*
*  SYNTAXE   :
*  -----------
*
*           RIG1=MUTU  MOD1  CAR1 GEO2
*
*                  MOD1 MODELE DE CALCUL, TYPE MMODEL
*                  CAR1 CHAMP PAR ELEMENT DE CARACTERISTIQUES
*                       GEOMETRIQUES ET MATERIELLES, TYPE MCHAML
*                  GEO2 MAILLAGE INDUCTEUR SUPPORT DU COURANT, TYPE MELEME
*
*-----------------------------------------------------------------------
*    Y. STEPHAN   (FEVRIER 97)
*_______________________________________________________________________
      SUBROUTINE MUTU

      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)

-INC PPARAM
-INC CCOPTIO
-INC SMRIGID

      IPMODL=0
      IPCHEL=0
      IPGEOM=0
      IPMUTU=0
C_______________________________________________________________________
C     LECTURE DU MODELE
C
      MOTERR(1:8)=' MODELE '
      CALL MESLIR(-137)
      CALL LIROBJ('MMODEL',IPMODL,1,iret)
      IF (IERR.NE.0) RETURN
      CALL ACTOBJ('MMODEL  ',IPMODL,1)
      IF (IERR.NE.0) RETURN
C_______________________________________________________________________
C     LECTURE DU 1 IER MCHAML
C
      MOTERR(1:8)=' MCHAML '
      CALL MESLIR(-137)
      CALL LIROBJ('MCHAML',IPIN,1,iret)
      IF (IERR.NE.0) RETURN
      CALL ACTOBJ('MCHAML  ',IPIN,1)
      IF (IERR.NE.0) RETURN

      CALL REDUAF(IPIN,IPMODL,IPCHEL,0,IR,KER)
      IF (IR  .NE.1) CALL ERREUR(KER)
      IF (IERR.NE.0) RETURN
C_______________________________________________________________________
C     LECTURE DU MAILLAGE
C
      MOTERR(1:8)=' MAILLAGE '
      CALL MESLIR(-137)
      CALL LIROBJ('MAILLAGE',IPGEOM,0,iret)
      IF (IERR.NE.0) RETURN
      IF (IPGEOM.NE.0) THEN
        CALL ACTOBJ('MAILLAGE',IPGEOM,1)
        IF (IERR.NE.0) RETURN
      ENDIF
C_______________________________________________________________________
C     CALCUL DE LA MATRICE DE MUTUELLE
C
      CALL MUTU1(IPMODL,IPCHEL,IPGEOM,IPMUTU)
      IF (IERR.NE.0) RETURN

C_______________________________________________________________________
C     ECRITURE DU RESULTAT
      CALL ECROBJ('RIGIDITE',IPMUTU)

C      RETURN
      END