Télécharger calg3.eso

Retour à la liste

Numérotation des lignes :

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

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