Télécharger elpdgf.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDGF SOURCE CHAT 05/01/12 23:36:29 5004
  2. SUBROUTINE ELPDGF(R,CB,CG)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-B,D-H,O-Z)
  5. IMPLICIT COMPLEX*16(C)
  6. C
  7. C =====================================================================
  8. C ELFE
  9. C CALCUL DE LA FONCTION DE GREEN ET DE DES DERIVEES
  10. C POUR UNE PLAQUE EN DYNAMIQUE
  11. C
  12. C =====================================================================
  13. C
  14. C
  15. C
  16. DIMENSION CG(5)
  17. DIMENSION CZ(1)
  18. DIMENSION CGF0(1)
  19. DIMENSION CGF1(1)
  20. DIMENSION CGF2(1)
  21. DIMENSION CGF3(1)
  22. DIMENSION CGF4(1)
  23. C
  24. C ON UTILISE LES CGFI COMME TABLEAUX DE TRAVAIL
  25. C POUR CALCULER LES H10,H11,H10(IBR),H11(IBR)
  26. C
  27. np = 10
  28. xm = 5d0
  29. CI = (0d0,1d0)
  30. CZ(1) = CB* R
  31. N=1
  32. CALL ELPDBE(CZ ,N,np,CGF0,CGF1,xm)
  33. DO 100 I = 1,N
  34. CZ(I) = CI*CZ(I)
  35. 100 CONTINUE
  36. CALL ELPDBE(CZ ,N,np,CGF2,CGF3,xm)
  37. DO 110 I = 1,N
  38. CZ(I) = CZ(I)/CI
  39. 110 CONTINUE
  40.  
  41. DO 150 I = 1,N
  42. CH0M= CGF2(I) - CGF0(I)
  43. CH0P= CGF2(I) + CGF0(I)
  44. CH1M= CI * CGF3(I) - CGF1(I)
  45. CH1P= CI * CGF3(I) + CGF1(I)
  46. CGF0(I) = CH0M
  47. CGF1(I) =-1.* CH1M
  48. CGF2(I) = CH0P + ( CH1M / CZ(I) )
  49. CGF3(I) = -1/CZ(I) *CH0P
  50. & -1 *CH1P
  51. & -2 /(CZ(I)** 2) *CH1M
  52. CGF4(I) = CH0M
  53. & +3 /(CZ(I)** 2) *CH0P
  54. & +2 / CZ(I) *CH1P
  55. & +6 /(CZ(I)** 3) *CH1M
  56. 150 CONTINUE
  57. C
  58. CG(1)= -1*( CGF0(1))*CI/(8*CB*CB)
  59. CG(2)= -1*CGF1(1)*CI/(8*CB)
  60. CG(3)= -1*CGF2(1)*CI/8
  61. CG(4)= -1*CGF3(1)*CI*(CB)/8
  62. CG(5)= -1*CGF4(1)*CI*(CB*CB)/8
  63. C
  64. RETURN
  65. END
  66.  
  67.  

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