Télécharger elpdk2.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDK2 SOURCE CHAT 05/01/12 23:36:54 5004
  2. SUBROUTINE ELPDK2(XLP0
  3. & ,XD,XNU,CB,ISTAT,NTRAP2
  4. & ,CS1)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-B,D-H,O-Z)
  7. IMPLICIT COMPLEX*16(C)
  8. ************************************************************************
  9. *
  10. * PLAQUES EN DYN PAR EQUATION INTEGRALE:
  11. *
  12. * CALCUL DES TERMES "DIAGONAUX"
  13. * V .. V,N0 0
  14. * VN 0 VN,N0 ..
  15. * MN .. MN,N0 0
  16. * KN 0 KN,N0 ****** voir elpdm1
  17. *
  18. ************************************************************************
  19. DIMENSION CS0 (8)
  20. DIMENSION CS1 (8)
  21. DIMENSION S1 (8)
  22. *
  23. DO 10 I=1,8
  24. S1(I) = 0D0
  25. CS0(I) = 0D0
  26. 10 CONTINUE
  27. *
  28. * 1-- CALCUL ANALYTIQUE DES TERMES STATIQUES sauf knn0
  29. *
  30. CALL ELPI02(XLP0,XD,XNU,S1)
  31. *
  32. * 2-- CALCUL PAR TRAPEZE SUR op(W - WS)
  33. *
  34. IF ( ISTAT.NE.1) THEN
  35. CALL ELPDK1(XLP0,XD,XNU,NTRAP2,CS0,CB)
  36. ENDIF
  37. *
  38. * 3-- ON AJOUTE
  39. *
  40. CS1(1) = CS0(1) + S1(1)
  41. CS1(2) = 0D0
  42. CS1(3) = CS0(3) + S1(3)
  43. CS1(4) = 0D0
  44. CS1(5) = 0D0
  45. CS1(6) = CS0(6) + S1(6)
  46. CS1(7) = 0D0
  47. CS1(8) = CS0(8) + S1(8)
  48. *
  49. RETURN
  50. END
  51. *
  52.  
  53.  

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