Esope source files

   1 : C PERMEA    SOURCE    CHAT      06/08/24    21:57:30     5529           
   2 :       SUBROUTINE PERMEA                                                         
   3 : *______________________________________________________________________         
   4 : *                                                                               
   5 : *      OPERATEUR PERMEABILITE   ( MILIEUX POREUX )                              
   6 : *                                                                               
   7 : *  syntaxe :
   8 : * ----------
   9 : *                                                                               
  10 : *           CHPER =  PERM MODL1 CHMAT ;
  11 : *                                                                               
  12 : *       MODL1  Modele de calcul (type MMODEL)                                   
  13 : *       CHMAT  champ par element de materiau                                    
  14 : *                                                                               
  15 : *_______________________________________________________________________        
  16 : *                                                                               
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 : *                                                                               
  20 : -INC CCOPTIO                                                                    
  21 : *                                                                               
  22 : *     Lecture du modele                                                         
  23 : *                                                                               
  23 :       CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)
  24 :       IF (IERR.NE.0) RETURN                                                     
  25 : *                                                                               
  25 :       CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU)
  26 :       IF(IERR.NE.0) RETURN
  27 : *                                                                               
  27 :       CALL PERMAB(IPMODL,IPCHE1,IPRIGI,IRET)                                    
  28 : *                                                                               
  28 :       IF (IRET.EQ.1) CALL ECROBJ('RIGIDITE',IPRIGI)                             
  29 : *                                                                               
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :  
  33 :