Télécharger kafint.eso

Retour à la liste

Numérotation des lignes :

  1. C KAFINT SOURCE CB215821 17/11/30 21:16:32 9639
  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=MIN(Z1,Z2)
  15. ZMAX= MAX(Z1,Z2)
  16. RMAX= MAX(R1,R2)
  17.  
  18. IF(R.GE.RMAX.OR.Z.LE.ZMIN.OR.Z.GE.ZMAX) THEN
  19. A=2.
  20.  
  21. ELSEIF(ABS(Z-Z1).LE.DMIN.OR.ABS(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.  
  30.  
  31.  
  32.  
  33.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales