Sources Esope | Cast3M

kaftan
C KAFTAN    SOURCE    CB215821  17/11/30    21:16:33     9639           
      SUBROUTINE KAFTAN(R1,Z1,RA,ZA,RB,ZB,A)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C----------------------------------------------------------------------
C SP appele par KAINTE
C     existence d un point tangent associe au segment obstructeur AB
C     points (RA,ZA) et (RB,ZB)
C     resultat A : 2 si exterieur
C                : represente un cosinus sinon
C----------------------------------------------------------------------
 
      DMIN=1.D-3
      IF(ABS(ZB-ZA).LE.DMIN) THEN
        A=2.
      ELSE
        A=Z1*(RB-RA)+(RA*ZB-RB*ZA)
        A=A/R1/(ZB-ZA)
      ENDIF
      RETURN
      END