Télécharger simpsn.eso

Retour à la liste

Numérotation des lignes :

  1. C SIMPSN SOURCE CHAT 05/01/13 03:18:58 5004
  2. SUBROUTINE SIMPSN(NNAPPE,NBG,IOP,POIGAU)
  3. *
  4. ************************************************************************
  5. * CALCUL DES POIDS (DES POINTS A TRAVERS L'EPAISSEUR) PAR LA METHODE *
  6. * D'INTEGRATION DE SIMPSON *
  7. * MULTIPLICATION PAR LES POIDS DES POINTS DANS LE PLAN *
  8. ************************************************************************
  9. * *
  10. * ENTREE : *
  11. * *
  12. * NNAPPE = NOMBRE DE NAPPES DANS L'EPAISSEUR *
  13. * NBG = NOMBRE DE POINTS DANS LE PLAN *
  14. * IOP = 1 CALCULE LES POIDS TOTAUX (W(QSI,ETA,ZETA))A PARTIR DES *
  15. * POIDS DE SURFACE(W(QSI,ETA) *
  16. * = 2 CALCULE LES POIDS DE SURFACE (W(QSI,ETA)A PARTIR DES *
  17. * POIDS TOTAUX(W(QSI,ETA,ZETA)) *
  18. * POIGAU = TABLEAU CONTENANT LES POIDS CORRESPONDANT A L'INTEGRATION *
  19. * DANS LE PLAN *
  20. * SORTIE : *
  21. * *
  22. * POIGAU = POIDS TENANT COMPTE DE L'INTEGRATION DANS L'EPAISSEUR *
  23. * *
  24. ************************************************************************
  25. *
  26. IMPLICIT INTEGER(I-N)
  27. IMPLICIT REAL*8 (A-H,O-Z)
  28. *
  29. DIMENSION POIGAU(*),W(100)
  30. *
  31. *
  32. PAS = 2.D0/(3.D0*(DBLE(NNAPPE)-1.D0))
  33. *
  34. DO 1 IZ = 1,NNAPPE
  35. ITS = MOD(IZ,2)
  36. IF (IZ.EQ.1.OR.IZ.EQ.NNAPPE) THEN
  37. COEF = 1.D0
  38. ELSE IF (ITS.NE.0.) THEN
  39. COEF = 2.D0
  40. ELSE
  41. COEF = 4.D0
  42. ENDIF
  43. W(IZ) = COEF * PAS
  44. 1 CONTINUE
  45. *
  46. IF(IOP.EQ.1) THEN
  47. DO 10 IA=1,NNAPPE
  48. NP=(IA-1)*NBG
  49. DO 20 IG=1,NBG
  50. POIGAU(NP+IG)=POIGAU(NP+IG)*W(IA)
  51. 20 CONTINUE
  52. 10 CONTINUE
  53. ELSE
  54. DO 30 IG=1,NBG
  55. POIGAU(IG) = POIGAU(IG)/W(1)
  56. 30 CONTINUE
  57. ENDIF
  58. *
  59. * WRITE(*,*) ' matrice B elem =',IB
  60. * WRITE(*,1191) (POIGAU(IA),IA=1,3)
  61. *1191 FORMAT(12(1X,E8.2))
  62. C
  63. RETURN
  64. END
  65.  
  66.  
  67.  
  68.  

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