C TUFIC1    SOURCE    CHAT      05/01/13    03:55:18     5004           
      SUBROUTINE TUFIC1(XM,XP,THETA,IRETOU,PRECIS,XM0,XP0)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C    Include contenant quelques constantes dont XPI :
-INC CCREEL
C=======================================================================
C      POUR SAVOIR SI ON EST ELASTIQUE OU PAS
C     
C         ENTREES:XM,XP      ETAT DE CONTRAINTES
C                 THETA    ANGLE DE LA FISSURE
C     
C         SORTIE:IRETOU    0 SI ELASTIQUE,1 SINON
C=======================================================================
      IRETOU=0
C     
C      CALCUL DE LA VALEUR DE LA FONCTION DE CHARGE EN OUVERTURE
C     
      CRIT=SIN(THETA)/2.D0
      CRIT=CRIT-COS((THETA+XPI*XP/XP0)/2.D0)
      CRIT=CRIT+XM/XM0
      IF(CRIT.GE.PRECIS) THEN
         IRETOU=1
         RETURN
      ENDIF
C     
C      CALCUL DE LA VALEUR DE LA FONCTION DE CHARGE EN FERMETURE
C     
      CRIT=SIN(THETA)/2.D0
      CRIT=CRIT-COS((THETA-XPI*XP/XP0)/2.D0)
      CRIT=CRIT-XM/XM0
      IF(CRIT.GE.PRECIS) THEN
         IRETOU=1
         RETURN
      ENDIF
      RETURN
      END