Esope source files

   1 : C KAFTAN    SOURCE    CHAT      05/01/13    00:50:31     5004           
   2 :       SUBROUTINE KAFTAN(R1,Z1,RA,ZA,RB,ZB,A)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C----------------------------------------------------------------------
   6 : C SP appele par KAINTE
   7 : C     existence d un point tangent associe au segment obstructeur AB
   8 : C     points (RA,ZA) et (RB,ZB)
   9 : C     resultat A : 2 si exterieur
  10 : C                : represente un cosinus sinon
  11 : C----------------------------------------------------------------------
  12 : 
  13 :       DMIN=1.D-3
  14 :       IF(DABS(ZB-ZA).LE.DMIN) THEN
  15 :         A=2.
  16 :       ELSE
  17 :         A=Z1*(RB-RA)+(RA*ZB-RB*ZA)
  18 :         A=A/R1/(ZB-ZA)
  19 :       ENDIF
  20 :       RETURN
  21 :       END
  22 :