C CRIT1     SOURCE    CHAT      05/01/12    22:30:38     5004           
      FUNCTION CRIT1(XM,XP,THETA,XM0,XP0)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C    Include contenant quelques constantes dont XPI :
-INC CCREEL
C=======================================================================
C     CALCULE LA VALEUR DE LA FONCTION DE CHARGE
C     
C     ENTREES:XM      VALEUR DU MOMENT
C             XP      VALEUR DE L EFFORT NORMAL
C             THETA   ANGLE DE LA FISSURE
C     
C     SORTIE:CRIT1    VALEUR DE LA FONCTION
C=======================================================================
      CRI=2.D0*XP0*(1.D0-THETA/XPI)*XM
      CRI=CRI+XM0*SIN(THETA)*XP
      IF(CRI.GE.XZERO) THEN
         GO TO 1
      ELSE
         GO TO 2
      ENDIF
 1    CONTINUE
      CRIT1=SIN(THETA)/2.D0
      CRIT1=CRIT1-COS((THETA+XPI*XP/XP0)/2.D0)
      CRIT1=CRIT1+XM/XM0
      RETURN
 2    CONTINUE
      CRIT1=SIN(THETA)/2.D0
      CRIT1=CRIT1-COS((THETA-XPI*XP/XP0)/2.D0)
      CRIT1=CRIT1-XM/XM0
      RETURN
      END