Télécharger pb302.eso

Retour à la liste

Numérotation des lignes :

pb302
  1. C PB302 SOURCE MAGN 10/05/18 21:16:33 6675
  2. SUBROUTINE PB302(XREF,X,Y,PG,FN,GR,FM,GM,ND,NP,MP,NPG,NOM2,ITYPI)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C************************************************************************
  6. C
  7. C CALCULE LES FONCTIONS DE FORME D'UN : TRI3
  8. C
  9. C ^ eta
  10. C |
  11. C a |n3
  12. C |\
  13. C | \ a=sqrt(2)
  14. C | \
  15. C |___\ _____>ksi
  16. C 0 a
  17. C n1 n2
  18. C
  19. C
  20. C************************************************************************
  21.  
  22. REAL*8 XREF(ND,NP),X(NPG),Y(NPG)
  23. DIMENSION FN(NP,NPG),GR(ND,NP,NPG),PG(NPG)
  24. DIMENSION FM(MP,NPG),GM(ND,MP,NPG)
  25. REAL*8 A,B,C,D,U(5),H(5)
  26. CHARACTER*4 NOM2
  27. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  28.  
  29. NP2=NP*ND
  30.  
  31. R2=SQRT(2.D0)
  32.  
  33. XREF(1,1)=0.D0
  34. XREF(2,1)=0.D0
  35. XREF(1,2)=R2
  36. XREF(2,2)=0.D0
  37. XREF(1,3)=0.D0
  38. XREF(2,3)=R2
  39.  
  40. IF(ITYPI.EQ.2)THEN
  41. X(1)=0.D0
  42. Y(1)=0.D0
  43. X(2)=R2
  44. Y(2)=0.D0
  45. X(3)=0.D0
  46. Y(3)=R2
  47. DO 2 L=1,3
  48. PG(L)=1.D0/3.D0
  49. 2 CONTINUE
  50. ELSE
  51. CALL CALUHH(X,Y,PG,NPG)
  52. ENDIF
  53. DO 1 L=1,NPG
  54. C
  55. FN(1,L)=-R2/2.D0*(X(L)+Y(L)-R2)
  56. FN(2,L)=R2/2.D0*X(L)
  57. FN(3,L)=R2/2.D0*Y(L)
  58. C
  59. GR(1,1,L)=-R2/2.D0
  60. GR(2,1,L)=-R2/2.D0
  61. GR(1,2,L)=R2/2.D0
  62. GR(2,2,L)=0.D0
  63. GR(1,3,L)=0.D0
  64. GR(2,3,L)=R2/2.D0
  65. C
  66.  
  67. IF(NOM2.EQ.'P1P1')THEN
  68. FM(1,L)=FN(1,L)
  69. FM(2,L)=FN(2,L)
  70. FM(3,L)=FN(3,L)
  71. GM(1,1,L)=-R2/2.D0
  72. GM(2,1,L)=-R2/2.D0
  73. GM(1,2,L)=R2/2.D0
  74. GM(2,2,L)=0.D0
  75. GM(1,3,L)=0.D0
  76. GM(2,3,L)=R2/2.D0
  77. ELSE
  78. FM(1,L)=1.D0
  79. GM(1,1,L)=0.D0
  80. GM(2,1,L)=0.D0
  81. ENDIF
  82.  
  83. 1 CONTINUE
  84.  
  85. C WRITE(6,101)
  86. C WRITE(6,1002)FM
  87. C WRITE(6,1002)GM
  88. C WRITE(6,1002)FN
  89. C WRITE(6,1002)GR
  90. C WRITE(6,101)
  91. RETURN
  92. 1002 FORMAT(10(1X,1PD11.4))
  93. 1001 FORMAT(20(1X,I5))
  94. 101 FORMAT(1X,'... SUBPB302 ... FM,GM,FN,GR ',9(10H..........)/)
  95. C
  96. END
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  

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