Esope source files

   1 : C SEUIL     SOURCE    CHAT      06/03/16    21:26:20     5333           
   2 :         SUBROUTINE SEUIL(Q1,Q2,Q3,XI1,XJ2,SIGM,ft,VALc,Fi)
   3 : C
   4 : C CALCUL DE LA SURFACE SEUIL DE GURSON MODIFIE
   5 : C
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C
   8 :       PARAMETER (XUnDemi=0.5D0,XZERO=0.D0,UN=1.D0,DEUX=2.D0)
   9 :       PARAMETER (XPETIT=1.D-12)
  10 : C
  11 :         VAL=Q2*XI1/(DEUX*SIGM)
  12 : C
  13 :         IF (abs(VAL).LT.valc) THEN
  14 :           Fi=XJ2/(SIGM**2)
  15 :      1     +DEUX*Q1*fT*COSH(VAL)-(UN+(Q3*fT)**2)
  16 :         ELSE
  17 :           IF (val.GT.XZERO) THEN
  18 :             Fi=XJ2/(SIGM**2)
  19 :      1        +DEUX*Q1*fT*(cosh(valc)+sinh(valc)*(val-valc)
  20 :      2        +XUnDemi*cosh(valc)*(val-valc)**2)/(DEUX*SIGM)
  21 :      3        -(UN+(Q3*fT)**2)
  22 :           ELSE
  23 :             Fi=XJ2/(SIGM**2)
  24 :      1        +DEUX*Q1*fT*(cosh(valc)+sinh(valc)*(val+valc)
  25 :      2        +XUnDemi*cosh(valc)*(val+valc)**2)/(DEUX*SIGM)
  26 :      3        -(UN+(Q3*fT)**2)
  27 :           ENDIF
  28 :         ENDIF
  29 :         IF(ABS(Fi).LT.XPETIT) THEN 
  30 :          Fi=XZERO
  31 :         ENDIF
  32 : C
  33 :         END
  34 : 
  35 :