Esope source files

   1 : C ENERGI    SOURCE    FANDEUR   10/12/14    21:16:04     6812           
   2 :       SUBROUTINE ENERGI                                                         
   3 : *_______________________________________________________________________        
   4 : *                                                                               
   5 : *               OPERATEUR DENSITE D'ENERGIE                                     
   6 : *                                                                               
   7 : *     CALCULE LE CHAMELEM PRODUIT   SIGMA * EPS                                 
   8 : *                                                                               
   9 : *       CHEL3 = ENER MODL CHEL1 CHEL2 ;                                         
  10 : *                                                                               
  11 : *          MODL   MODELE DE CALCUL ,TYPE MMODEL                                 
  12 : *          CHEL1  CHAMP PAR ELEMENT DE SOUS TYPE CONTRAINTE OU                  
  13 : *                 DEFORMATION ( TYPE MCHAML )                                   
  14 : *          CHEL2  CHAMP PAR ELEMENT DE SOUS TYPE CONTRAINTE OU                  
  15 : *                 DEFORMATION ( TYPE MCHAML )                                   
  16 : *          CHEL3  CHAMP PAR ELEMENT DE SOUS TYPE SCALAIRE ( TYPE                
  17 : *                 MCHAML )                                                      
  18 : *                                                                               
  19 : *  PASSAGE AUX NOUVEAUX MCHAML PAR JM CAMPENON LE 01/91                         
  20 : *_______________________________________________________________________        
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8 (A-H,O-Z)
  23 : 
  24 : -INC CCOPTIO                                                                    
  25 : -INC SMCOORD                                                                    
  26 : *                                                                               
  27 : *   LECTURE DU MODELE                                                           
  28 : *                                                                               
  28 :       CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)                                     
  29 :       IF (IERR.NE.0) RETURN                                                     
  30 : *                                                                               
  31 : *   LECTURE DU PREMIER MCHAML                                                   
  32 : *                                                                               
  32 :       CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU)                                     
  33 :       IF (IERR.NE.0) RETURN                                                     
  34 : *                                                                               
  35 : *   LECTURE DU DEUXIEME MCHAML                                                  
  36 : *                                                                               
  36 :       CALL LIROBJ('MCHAML',IPCHE2,1,IRETOU)                                     
  37 :       IF (IERR.NE.0) RETURN                                                     
  38 : *                                                                               
  38 :       CALL ENERCA(IPMODL,IPCHE1,IPCHE2,IPCHR)                                   
  39 : *                                                                               
  39 :       IF (IPCHR.NE.0) CALL ECROBJ('MCHAML',IPCHR)                               
  40 : *                                                                               
  41 :       RETURN                                                                    
  42 :       END                                                                       
  43 : 
  44 :