C PERMEA    SOURCE    CHAT      06/08/24    21:57:30     5529           
      SUBROUTINE PERMEA                                                         
*______________________________________________________________________         
*                                                                               
*      OPERATEUR PERMEABILITE   ( MILIEUX POREUX )                              
*                                                                               
*  syntaxe :
* ----------
*                                                                               
*           CHPER =  PERM MODL1 CHMAT ;
*                                                                               
*       MODL1  Modele de calcul (type MMODEL)                                   
*       CHMAT  champ par element de materiau                                    
*                                                                               
*_______________________________________________________________________        
*                                                                               
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
*                                                                               
-INC CCOPTIO                                                                    
*                                                                               
*     Lecture du modele                                                         
*                                                                               
      CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)
      IF (IERR.NE.0) RETURN                                                     
*                                                                               
      CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU)
      IF(IERR.NE.0) RETURN
*                                                                               
      CALL PERMAB(IPMODL,IPCHE1,IPRIGI,IRET)                                    
*                                                                               
      IF (IRET.EQ.1) CALL ECROBJ('RIGIDITE',IPRIGI)                             
*                                                                               
      RETURN                                                                    
      END