Esope source files

   1 : C CLCISA    SOURCE    CHAT      05/01/12    22:04:30     5004           
   2 :       SUBROUTINE CLCISA (ET,YUNG,XNU,BETR,XLTR,CISA)                            
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C                                                                               
   6 :       DIMENSION ET(*),CISA(*),XLTR(*),ES(3)
   7 : C                                                                               
   8 : C***********************************************************************        
   9 : C**************************** INITIALISATION ***************************        
  10 : C***********************************************************************        
  11 : C                                                                               
  12 :       RFE=1.D-8
  13 :       RFS=YUNG*RFE
  14 : C
  15 :       ES(1)=ET(1)
  16 :       ES(2)=ET(2)
  17 :       ES(3)=ET(3)
  18 :       DUNPNU=2.D0*(1.D0+XNU)                                                    
  19 :       CISARS=YUNG*BETR/DUNPNU                                                   
  20 : C                                                                               
  21 :       XLTRI=XLTR(1)+XLTR(2)+XLTR(3)                                             
  22 : C                                                                               
  23 : C***********************************************************************        
  24 : C*********************** DIFFERENTS CAS PREVUS *************************        
  25 : C***********************************************************************        
  26 : C                                                                               
  27 :       IF(XLTRI.GT.0.D0) THEN
  28 : C                                                                               
  29 :         IF(ES(1).LT.0.D0.AND.ABS(ES(1)).GT.RFS) ES(1)=YUNG
  30 :         IF(ES(2).LT.0.D0.AND.ABS(ES(2)).GT.RFS) ES(2)=YUNG
  31 :         IF(ES(3).LT.0.D0.AND.ABS(ES(3)).GT.RFS) ES(3)=YUNG
  32 : C
  33 :         CISA(1)=MIN(ES(1),ES(2))/DUNPNU
  34 :         CISA(2)=MIN(ES(2),ES(3))/DUNPNU
  35 :         CISA(3)=MIN(ES(3),ES(1))/DUNPNU
  36 : C
  37 :         IF(CISA(1).LT.CISARS) CISA(1)=CISARS
  38 :         IF(CISA(2).LT.CISARS) CISA(2)=CISARS
  39 :         IF(CISA(3).LT.CISARS) CISA(3)=CISARS
  40 : C                                                                               
  41 :       ELSE                                                                      
  42 : C                                                                               
  43 :         CISA(1)=CISARS                                                          
  44 :         CISA(2)=CISARS
  45 :         CISA(3)=CISARS
  46 : C                                                                               
  47 :       ENDIF                                                                     
  48 : C
  49 :       RETURN                                                                    
  50 :       END                                                                       
  51 :