Esope source files

   1 : C UO2TES    SOURCE    CHAT      05/01/13    03:59:30     5004           
   2 :       SUBROUTINE UO2TES(XINVL,W,WMAX,SMAX,SIGMA,XLTR,PRECIZ,
   3 :      &                  BTR,WRUPT,
   4 :      &                  IB,IGAU,ICAS,KERRE)
   5 : *
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)     
   8 : -INC CCOPTIO
   9 : C 
  10 : C
  11 :       DIMENSION SIGMA(6)
  12 : C
  13 :       DIMENSION WMAX(3),SMAX(3),W(3),
  14 :      &    XLTR(3),XINVL(3),WRUPT(3)
  15 : C
  16 :       DIMENSION WREOUV(3)
  17 : 
  18 : C
  19 : C     ICAS = 1   AVANT APPEL A UO2CA
  20 : C     ICAS = 2   APRES APPEL A UO2CA
  21 : C
  22 :       KERRE=0
  23 : C
  24 : C------------------------------------------------------------------------
  25 : C            ON PASSE EN REVUE LES 3 DIRECTIONS POTENTIELLES
  26 : C                          DE FISSURATION
  27 : C------------------------------------------------------------------------
  28 : C
  29 :       DO 88 I=1,3
  30 : 
  31 :          IF(XINVL(I).NE.0.D0) THEN
  32 : 
  33 :            WREOUV(I)=BTR*MIN(WMAX(I),WRUPT(I))
  34 : *
  35 : *          TEST SUR OUVERTURE MINIMALE
  36 : *
  37 :            IF(W(I).LT.WREOUV(I)-1.D-10*WRUPT(I))  THEN
  38 :              KERRE=70
  39 :              WRITE(IOIMP,60000) ICAS
  40 : 60000        FORMAT (/2X, '         PROBLEME   -  ICAS =',I3 //)
  41 :              WRITE(IOIMP,67675) I,W(I), WREOUV(I),WMAX(I)
  42 : 67675        FORMAT(2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$' /
  43 :      &        4X,'ATTENTION   W < WREOUV '/
  44 :      &        2X, ' I=',I4,2X,'W=',1PE12.5,2X,'WREOUV=',1PE12.5,
  45 :      &        2X,'WMAX=',1PE12.5/
  46 :      &        2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'/)
  47 :           ENDIF
  48 : *
  49 : *          TEST SUR DEPASSEMENT DE LA LIMITE EN TRACTION
  50 : *
  51 :            IF(SIGMA(I)-XLTR(I).GT.PRECIZ) THEN
  52 :              KERRE=70
  53 :              WRITE(IOIMP,60000) ICAS
  54 :              WRITE(IOIMP,77058) I,SIGMA(I)
  55 : 77058        FORMAT(2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$' /
  56 :      &        4X,'ATTENTION   DEPASSEMENT DE LTR ' /
  57 :      &        2X, ' I=',I4,2X,'SIGMA=',1PE12.5 /
  58 :      &        2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'/)
  59 :            ENDIF
  60 : 
  61 : *
  62 : *          TEST SUR SECANTE
  63 : *
  64 :           IF(BTR.NE.1.D0) THEN
  65 :             IF(WMAX(I).LT.WRUPT(I).AND.W(I).GT.WREOUV(I)) THEN
  66 :               PPP = SMAX(I)/( WMAX(I)-WREOUV(I))
  67 :               ZZOB = W(I) -WREOUV(I)
  68 :               SSS = ZZOB*PPP
  69 :               IF(SIGMA(I).GT.PRECIZ.AND.
  70 :      &           ABS(SIGMA(I)-SSS).GT.PRECIZ) THEN
  71 : 
  72 :              KERRE=70
  73 :              WRITE(IOIMP,60000) ICAS
  74 :                 WRITE(IOIMP,79912) I,SIGMA(I),SSS,WMAX(I),
  75 :      &                       W(I),WREOUV(I),ZZOB,SMAX(I)
  76 : 79912     FORMAT(2X,' ################################### '/
  77 :      &    2X,'I=',I4,2X,'SIGMA=',1PE12.5,
  78 :      &    2X,'SSS=',1PE12.5/2X,'WMAX=',1PE12.5/
  79 :      &    2X,'W=',1PE12.5,2X,'WREOUV=',1PE12.5,
  80 :      &    2X,'W-WREOUV = ',1PE12.5,2X,'SMAX= ',1PE12.5/
  81 :      &        2X,' ################################### '/)
  82 :            PRINT *,'SIGMA=',SIGMA(I)
  83 :            PRINT *,'SSS  =',SSS
  84 :            PRINT *,'WMAX =',WMAX(I)
  85 :            PRINT *,'W    =',W(I)
  86 :            PRINT *,'WREOUV=',WREOUV(I)
  87 :            PRINT *,'W - WREOUV=',ZZOB
  88 :            PRINT *,'SMAX    =',SMAX(I)
  89 : 
  90 :           ENDIF
  91 :         ENDIF
  92 :         ENDIF
  93 : 
  94 :       ENDIF
  95 : 88    CONTINUE
  96 : 
  97 :       IF(KERRE.NE.0) THEN
  98 :           WRITE(IOIMP,70902) IB,IGAU
  99 : 70902     FORMAT(2X,'SORTIE OTTEST  IB=',I6,2X,'IGAU=',I3////)
 100 :       ENDIF
 101 : C
 102 :       RETURN
 103 :       END     
 104 :  
 105 :  
 106 :  
 107 :  
 108 :  
 109 :  
 110 :  
 111 :  
 112 :  
 113 :  
 114 :  
 115 :  
 116 :  
 117 :