Esope source files

   1 : C GAMCAF    SOURCE    CHAT      05/01/13    00:15:23     5004           
   2 :       SUBROUTINE GAMCAF(SI,DSI,GAM)                                             
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C                                                                               
   6 : C--------------------------------------------                                   
   7 : C     CETTE ROUTINE  PERMET DE CALCULER     |                                   
   8 : C     GAMMA POUR ARRIVER SUR LE CRITERE     |                                   
   9 : C     DE COMPRESSION                        |                                   
  10 : C--------------------------------------------                                   
  11 : C                                                                               
  12 :       GAM=0.D0                                                                  
  13 :       IF(DSI.EQ.0.D0) GOTO 10                                                   
  14 :       GAM=-SI/DSI                                                               
  15 :    10 CONTINUE                                                                  
  16 :       IF(GAM.LT.0.D0) GAM=0.D0                                                  
  17 :       IF(GAM.GT.1.D0) GAM=1.D0                                                  
  18 :       RETURN                                                                    
  19 :       END                                                                       
  20 :