C ENERGI    SOURCE    CB215821  19/08/01    21:15:45     10279          
      SUBROUTINE ENERGI
*_______________________________________________________________________
*
*               OPERATEUR DENSITE D'ENERGIE
*
*     CALCULE LE CHAMELEM PRODUIT   SIGMA * EPS
*
*       CHEL3 = ENER MODL CHEL1 CHEL2 ;
*
*          MODL   MODELE DE CALCUL ,TYPE MMODEL
*          CHEL1  CHAMP PAR ELEMENT DE SOUS TYPE CONTRAINTE OU
*                 DEFORMATION ( TYPE MCHAML )
*          CHEL2  CHAMP PAR ELEMENT DE SOUS TYPE CONTRAINTE OU
*                 DEFORMATION ( TYPE MCHAML )
*          CHEL3  CHAMP PAR ELEMENT DE SOUS TYPE SCALAIRE ( TYPE
*                 MCHAML )
*
*  PASSAGE AUX NOUVEAUX MCHAML PAR JM CAMPENON LE 01/91
*_______________________________________________________________________
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)


-INC PPARAM
-INC CCOPTIO
      IPMODL = 0
      IPCHE1 = 0
      IPCHE2 = 0
      IPCHR  = 0

*
*   LECTURE DU MODELE
*
      CALL LIROBJ('MMODEL  ',IPMODL,1,IRETOU)
      CALL ACTOBJ('MMODEL  ',IPMODL,1)
      IF (IERR.NE.0) RETURN
*
*   LECTURE DU PREMIER MCHAML
*
      CALL LIROBJ('MCHAML  ',IPIN,1,IRETOU)
      CALL ACTOBJ('MCHAML  ',IPIN,1)
      IF (IERR.NE.0) RETURN
      CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
      IF(IR   .NE. 1) CALL ERREUR(KER)
      IF(IERR .NE. 0) RETURN
*
*   LECTURE DU DEUXIEME MCHAML
*
      CALL LIROBJ('MCHAML  ',IPIN,1,IRETOU)
      CALL ACTOBJ('MCHAML  ',IPIN,1)
      IF (IERR.NE.0) RETURN
      CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
      IF(IR   .NE. 1) CALL ERREUR(KER)
      IF(IERR .NE. 0) RETURN
*
      CALL ENERCA(IPMODL,IPCHE1,IPCHE2,IPCHR)
      IF(IERR .NE. 0) RETURN
*
      CALL ACTOBJ('MCHAML',IPCHR,1)
      CALL ECROBJ('MCHAML',IPCHR)

      END