Télécharger elpdi1.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDI1 SOURCE CHAT 05/01/12 23:36:35 5004
  2. SUBROUTINE ELPDI1(P0,XN0,XT0
  3. & ,Q,XN,XT,A,B,OAB,TETA,R,XLQ
  4. & ,XD,XNU
  5. & ,NTRAP,CS1,CB,ISTAT)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-B,D-H,O-Z)
  8. IMPLICIT COMPLEX*16(C)
  9. ************************************************************************
  10. *
  11. * PLAQUES PAR EQUATION INTEGRALE:
  12. *
  13. * ON INTEGRE PAR LA METHODE DES TRAPEZES POUR CALCULER
  14. * LES COEFFICIENTS NON DIAGONAUX
  15. *
  16. ************************************************************************
  17. *
  18. *
  19. DIMENSION P0(2)
  20. DIMENSION A (2)
  21. DIMENSION B (2)
  22. DIMENSION OAB(2)
  23. DIMENSION Q (2)
  24. DIMENSION QI(2)
  25. DIMENSION COP(8)
  26. *
  27. DIMENSION XN (2)
  28. DIMENSION XNQI(2)
  29. DIMENSION XN0(2)
  30. *
  31. DIMENSION XT (2)
  32. DIMENSION XTQI(2)
  33. DIMENSION XT0(2)
  34. *
  35. *
  36. DIMENSION CS1 (8)
  37. *
  38. * 1- TAILLE DES TRAPEZES ARC DE CERCLE
  39. *
  40. X0 = XLQ/(2*NTRAP)
  41. X1 = (2*TETA*R) / NTRAP
  42. *
  43. * 2- INTEGRATION
  44. *
  45. DO 10 J=1,8
  46. CS1(J) = 0D0
  47. 10 CONTINUE
  48. DO 100 I=1,NTRAP
  49. TETAI = ( (I-1)*2 + 1) * TETA /NTRAP
  50.  
  51. QI(1) = OAB(1) + COS(TETAI) *( A(1) -OAB(1))
  52. & - SIN(TETAI) *( A(2) -OAB(2))
  53. QI(2) = OAB(2) + SIN(TETAI) *( A(1) -OAB(1))
  54. & + COS(TETAI) *( A(2) -OAB(2))
  55. * -- NORMALE
  56. XNQI(1) = ( QI(1) - OAB(1) ) / R
  57. XNQI(2) = ( QI(2) - OAB(2) ) / R
  58. PS = XNQI(1) * XN(1) + XNQI(2) * XN(2)
  59. XNQI(1) = ( PS / ABS(PS)) * XNQI(1)
  60. XNQI(2) = ( PS / ABS(PS)) * XNQI(2)
  61.  
  62. * -- TANGENTE
  63. XTQI(1) = -1*XNQI(2)
  64. XTQI(2) = XNQI(1)
  65.  
  66. CALL ELPDOP(P0,XN0,QI,XNQI,XTQI,XD,XNU,COP,CB,ISTAT)
  67. DO 110 J=1,8
  68. CS1(J) = CS1(J) + ( X1 * COP(J))
  69. 110 CONTINUE
  70. 100 CONTINUE
  71. *
  72. RETURN
  73. END
  74.  
  75.  

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