Esope source files

   1 : C TUFINO    SOURCE    CHAT      05/01/13    03:55:59     5004           
   2 :       SUBROUTINE TUFINO(XM,XP,THETA,IRET,XM0,XP0)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : C    Include contenant quelques constantes dont XPI :
   6 : -INC CCREEL
   7 : C=======================================================================
   8 : C     SOUS-PROGRAMME PERMETTANT DE SAVOIR DANS QUELLE REGION
   9 : C     DU PLAN M-P ON SE TROUVE
  10 : C     
  11 : C     ENTREES : XM     :VALEUR DU MOMENT
  12 : C               XP     :VALEUR DE L EFFORT NORMAL
  13 : C               THETA  :ANGLE DE LA FISSURE
  14 : C     
  15 : C     SORTIE : IRET    :ENTIER DEFINISSANT LA REGION DU PLAN
  16 : C=======================================================================
  17 :       CRI=2.D0*XP0*(1.D0-THETA/XPI)*XM
  18 :       CRI=CRI+XM0*SIN(THETA)*XP
  19 :       IF(CRI.GE.XZERO) THEN
  20 :          GO TO 1
  21 :       ELSE
  22 :          GO TO 2
  23 :       ENDIF
  24 :  1    CONTINUE
  25 :       CRIT1=SIN(THETA)/2.D0
  26 :       CRIT1=CRIT1-COS((THETA-XPI*XP/XP0)/2.D0)
  27 :       CRIT1=CRIT1-XM/XM0
  28 :       IF(CRIT1.LT.XZERO) IRET=1
  29 :       IF(CRIT1.GE.XZERO.AND.XP.GT.0.D0) IRET=2
  30 :       IF(CRIT1.GE.XZERO.AND.XP.LE.0.D0) IRET=4
  31 :       RETURN
  32 :  2    CONTINUE
  33 :       CRIT1=SIN(THETA)/2.D0
  34 :       CRIT1=CRIT1-COS((THETA+XPI*XP/XP0)/2.D0)
  35 :       CRIT1=CRIT1+XM/XM0
  36 :       IF(CRIT1.LT.XZERO) IRET=3
  37 :       IF(CRIT1.GE.XZERO.AND.XP.GT.0.D0) IRET=2
  38 :       IF(CRIT1.GE.XZERO.AND.XP.LE.0.D0) IRET=4
  39 :       RETURN
  40 :       END
  41 : 
  42 : 
  43 :  
  44 :  
  45 :