Télécharger caluhh.eso

Retour à la liste

Numérotation des lignes :

  1. C CALUHH SOURCE CHAT 05/01/12 21:48:09 5004
  2. SUBROUTINE CALUHH(X,Y,PG,NG)
  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 HAMMER SUR LE
  7. C TRIANGLE DE REFERENCE
  8. C
  9. C NG NOMBRE DE POINTS D'INTEGRATION
  10. C
  11. C X,Y COORDONNEES DES POINTS DE HAMMER SUR LE TRIANGLE DE REFERENCE
  12. C DE SOMMETS (0.,0.) , (SQRT(2),0) , (0,SQRT(2))
  13. C PG POIDS DES POINTS DE HAMMER
  14. C
  15. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  16. REAL*8 X(NG),Y(NG)
  17. DIMENSION PG(NG)
  18. REAL*8 A1,A2,B1,B2,C1,P1,P2,P3,R2
  19. C***
  20. R2=SQRT(2.D0)
  21.  
  22. IF(NG.NE.7) GO TO 10
  23. A1=(9.D0-2.D0*SQRT(15.D0))/21.D0
  24. A2=(9.D0+2.D0*SQRT(15.D0))/21.D0
  25. B1=(6.D0+SQRT(15.D0))/21.D0
  26. B2=(6.D0-SQRT(15.D0))/21.D0
  27. C1=1.D0/3.D0
  28. C
  29. P1=9.D0/40.D0
  30. P2=(155.D0+SQRT(15.D0))/1200.D0
  31. P3=(155.D0-SQRT(15.D0))/1200.D0
  32. C
  33. X(1)=C1
  34. Y(1)=C1
  35. X(2)=B1
  36. Y(2)=B1
  37. X(3)=A1
  38. Y(3)=B1
  39. X(4)=B1
  40. Y(4)=A1
  41. X(5)=B2
  42. Y(5)=B2
  43. X(6)=A2
  44. Y(6)=B2
  45. X(7)=B2
  46. Y(7)=A2
  47. DO 71 L=1,NG
  48. X(L)=X(L)*R2
  49. Y(L)=Y(L)*R2
  50. 71 CONTINUE
  51. C
  52. PG(1)=P1
  53. PG(2)=P2
  54. PG(3)=P2
  55. PG(4)=P2
  56. PG(5)=P3
  57. PG(6)=P3
  58. PG(7)=P3
  59. C
  60. RETURN
  61. 10 CONTINUE
  62. IF(NG.NE.1) GO TO 30
  63. X(1)=1.D0/3.D0*R2
  64. Y(1)=1.D0/3.D0*R2
  65. PG(1)=1.D0
  66. RETURN
  67. 30 CONTINUE
  68. IF(NG.NE.3) CALL ARRET(0)
  69. X(1)=0.D0
  70. Y(1)=0.5D0
  71. X(2)=0.5D0
  72. Y(2)=0.5D0
  73. X(3)=0.5D0
  74. Y(3)=0.D0
  75. DO 31 L=1,NG
  76. X(L)=X(L)*R2
  77. Y(L)=Y(L)*R2
  78. 31 CONTINUE
  79. PG(1)=1.D0/3.D0
  80. PG(2)=1.D0/3.D0
  81. PG(3)=1.D0/3.D0
  82. RETURN
  83. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  84. END
  85.  
  86.  
  87.  
  88.  

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