C SEUIL     SOURCE    CHAT      06/03/16    21:26:20     5333           
        SUBROUTINE SEUIL(Q1,Q2,Q3,XI1,XJ2,SIGM,ft,VALc,Fi)
C
C CALCUL DE LA SURFACE SEUIL DE GURSON MODIFIE
C
      IMPLICIT REAL*8(A-H,O-Z)
C
      PARAMETER (XUnDemi=0.5D0,XZERO=0.D0,UN=1.D0,DEUX=2.D0)
      PARAMETER (XPETIT=1.D-12)
C
        VAL=Q2*XI1/(DEUX*SIGM)
C
        IF (abs(VAL).LT.valc) THEN
          Fi=XJ2/(SIGM**2)
     1     +DEUX*Q1*fT*COSH(VAL)-(UN+(Q3*fT)**2)
        ELSE
          IF (val.GT.XZERO) THEN
            Fi=XJ2/(SIGM**2)
     1        +DEUX*Q1*fT*(cosh(valc)+sinh(valc)*(val-valc)
     2        +XUnDemi*cosh(valc)*(val-valc)**2)/(DEUX*SIGM)
     3        -(UN+(Q3*fT)**2)
          ELSE
            Fi=XJ2/(SIGM**2)
     1        +DEUX*Q1*fT*(cosh(valc)+sinh(valc)*(val+valc)
     2        +XUnDemi*cosh(valc)*(val+valc)**2)/(DEUX*SIGM)
     3        -(UN+(Q3*fT)**2)
          ENDIF
        ENDIF
        IF(ABS(Fi).LT.XPETIT) THEN 
         Fi=XZERO
        ENDIF
C
        END