Télécharger calg2.eso

Retour à la liste

Numérotation des lignes :

  1. C CALG2 SOURCE CHAT 05/01/12 21:46:02 5004
  2. SUBROUTINE CALG2(A,B,C,D,NG,H,U,X,Y,PG)
  3. C***********************************************************************
  4. C CALCUL DES COORDONNEES ET DES POIDS DES POINTS DE GAUSS SUR LE PAVE
  5. C (A,B) X (C,D) LES CONNAISSANT SUR LE SEGMENT (-1,1)
  6. C
  7. C NG NOMBRE DE POINTS DE GAUSS SUR LE SEGMENT (-1,1)
  8. C U(NG) ABCISSES DES POINTS DE GAUSS
  9. C H(NG) POIDS DES POINTS DE GAUSS
  10. C
  11. C X,Y COORDONNEES DES POINTS DE GAUSS SUR LE PAVE (A,B) X (C,D)
  12. C PG POIDS DES POINTS DE GAUSS SUR LE PAVE (A,B) X (C,D)
  13. C
  14. C***********************************************************************
  15. IMPLICIT INTEGER(I-N)
  16. IMPLICIT REAL*8 (A-H,O-Z)
  17. REAL*8 A,B,C,D,H(NG),U(NG),X(NG,NG),Y(NG,NG)
  18. REAL*8 PG(NG,NG)
  19.  
  20. C***
  21.  
  22. BMA=(B-A)/2.D0
  23. BPA=(B+A)/2.D0
  24. DMC=(D-C)/2.D0
  25. DPC=(D+C)/2.D0
  26. C
  27. DO 1 LY=1,NG
  28. DO 1 LX=1,NG
  29. X(LX,LY)=BMA*U(LX)+BPA
  30. Y(LX,LY)=DMC*U(LY)+DPC
  31. PG(LX,LY)=BMA*DMC*H(LX)*H(LY)
  32. 1 CONTINUE
  33.  
  34. C WRITE(6,*)' SUB CALG2 : POIDS ET POSITIONS DES PTS DE GAUSS '
  35. C WRITE(6,*)' SUB CALG2 : H,U,X,Y ET PG '
  36. C WRITE(6,1002)H
  37. C WRITE(6,1002)U
  38. C WRITE(6,1002)X
  39. C WRITE(6,1002)Y
  40. C WRITE(6,1002)PG
  41. 1002 FORMAT(10(1X,1PD11.4))
  42. RETURN
  43. END
  44.  
  45.  
  46.  

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