Esope source files

   1 : C GAMPO     SOURCE    CHAT      05/01/13    00:15:35     5004           
   2 :       FUNCTION GAMPO(S,DS,R)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 : C
   7 :       DIMENSION S(*),DS(*)
   8 : C
   9 :       TRDSIG=DS(1)+DS(2)+DS(3)
  10 :       IF(TRDSIG.NE.0.D0) GO TO 1
  11 :       GAMPO=100.D0
  12 :       RETURN
  13 :     1 A=-0.333333333333333333333333D0
  14 :       TRSIG0=S(1)+S(2)+S(3)
  15 :       GAMPO=(R-A*TRSIG0)/(A*TRDSIG)
  16 : C
  17 :       RETURN
  18 :       END
  19 :  
  20 :