C ENERGI    SOURCE    FANDEUR   10/12/14    21:16:04     6812           
      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 CCOPTIO                                                                    
-INC SMCOORD                                                                    
*                                                                               
*   LECTURE DU MODELE                                                           
*                                                                               
      CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)                                     
      IF (IERR.NE.0) RETURN                                                     
*                                                                               
*   LECTURE DU PREMIER MCHAML                                                   
*                                                                               
      CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU)                                     
      IF (IERR.NE.0) RETURN                                                     
*                                                                               
*   LECTURE DU DEUXIEME MCHAML                                                  
*                                                                               
      CALL LIROBJ('MCHAML',IPCHE2,1,IRETOU)                                     
      IF (IERR.NE.0) RETURN                                                     
*                                                                               
      CALL ENERCA(IPMODL,IPCHE1,IPCHE2,IPCHR)                                   
*                                                                               
      IF (IPCHR.NE.0) CALL ECROBJ('MCHAML',IPCHR)                               
*                                                                               
      RETURN                                                                    
      END