$$$$ PAS_EPTH
* PAS_EPTH  PROCEDUR  MB234859  25/03/13    21:15:03     12175          
'DEBP' PAS_EPTH   PRECED*TABLE MODEVAL*MMODEL MATEVAL*MCHAML TEVAL*MCHAML ;
************************************************************************
* Cette procedure calcule les contraintes thermiques pour les modeles
* poreux
*
* PRECED : La TABLE donnee en entree de PASAPAS. On a besoin des indices
*           - WTAB.'TEMPERATURE_REFERENCE'
* MODEVAL: Modele    d'entree
* MATEVAL: Materiaux d'entree
*
* TEVAL : Temperature (MCHAML) a laquelle est evalué EPSILON_THERMIQUE
*         Type MCHAML Obligatoire pour les JOINTS
*         Necessaire  pour certains éléments linéiques
*
* Le resultat est :
* MSRTHER2 : contraintes thermiques de dilatation
************************************************************************

   WTAB = PRECED. 'WTABLE';
  
   MOPOR  = WTAB.'MOD_POR' ;
   MAIPOR = 'EXTR' MOPOR  'MAIL' ;
   MAPOR  = WTAB.'MAT_POR' ;
   ZALM   = ('EXCO' 'ALPM' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
   ZBI    = ('EXCO' 'COB'  MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
   ZALPH  = ('EXCO' 'ALPH' MAPOR 'MSR0') 'CHAN' 'TYPE' 'SCALAIRE' ;
  
   ZFAC = 3.D0 '*' ( ZALM '-' ( ZBI '*' ZALPH ) );
   ZFAC = ZFAC  'CHAN' 'TYPE' 'SCALAIRE';
   ZFAC = 'CHAN' 'STRESSES' MOPOR ZFAC  ;
  
   ZDT  = 'REDU' (TEVAL '-' WTAB.'TEMPERATURE_REFERENCE') MAIPOR;
   ZDDT = 'CHAN' 'CHAM' ZDT MOPOR 'STRESSES' 'SCALAIRE';
  
   MSRTHER2 = ZFAC  '*' ZDDT ;

'FINP' MSRTHER2;