Esope source files

   1 : C KAFINT    SOURCE    CHAT      05/01/13    00:50:28     5004           
   2 :       SUBROUTINE KAFINT(R1,Z1,R2,Z2,R,Z,A)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C----------------------------------------------------------------------
   6 : C SP appele par KAINTE
   7 : C     appartenance du point (R,Z) au triangle defini par les
   8 : C     points (R1,Z1) et (R2,Z2)
   9 : C     resultat A : 2 si exterieur
  10 : C                : represente un cosinus sinon
  11 : C----------------------------------------------------------------------
  12 : 
  13 :       DMIN=1.D-3
  14 :       ZMIN=DMIN1(Z1,Z2)
  15 :       ZMAX=DMAX1(Z1,Z2)
  16 :       RMAX=DMAX1(R1,R2)
  17 : 
  18 :       IF(R.GE.RMAX.OR.Z.LE.ZMIN.OR.Z.GE.ZMAX) THEN
  19 :        A=2.
  20 : 
  21 :       ELSEIF(DABS(Z-Z1).LE.DMIN.OR.DABS(Z2-Z).LE.DMIN) THEN
  22 :        A=2.
  23 :       ELSE
  24 :        A=R*R*(Z2-Z1)*(Z2-Z1)-R2*R2*(Z-Z1)*(Z-Z1)-R1*R1*(Z2-Z)*(Z2-Z)
  25 :        A=A/2./R1/R2/(Z-Z1)/(Z2-Z)
  26 :       ENDIF
  27 :       RETURN
  28 :       END
  29 :