Télécharger pouro2.eso

Retour à la liste

Numérotation des lignes :

  1. C POURO2 SOURCE CHAT 05/01/13 02:19:46 5004
  2. C POUROT SOURCE CASTEM 85/06/13 12:34:25
  3. SUBROUTINE POURO2(REL,LRE,P,RP)
  4. C-----------------------------------------------------------------------
  5. C ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES
  6. C GLOBAUX EN 2D
  7. C D APRES TEDEL
  8. C
  9. C ENTREE
  10. C LRE NOMBRE DE LIGNES DU TABLEAU REL
  11. C P MATRICE DE PASSAGE
  12. C RP MATRICE ELEMENTAIRE 12X12 EN AXES LOCAUX
  13. C
  14. C SORTIE
  15. C REL MATRICE ELEMENTAIRE EN AXES GLOBAUX
  16. C
  17. C VIENT DE TEDEL AOUT 84
  18. C
  19. C D'APRES POUROT DC 98
  20. C-----------------------------------------------------------------------
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. DIMENSION REL(LRE,1),P(2,2),RP(6,6)
  24. DIMENSION NI(6),NJ(6)
  25. DATA ZER,UN/0.D0,1.0D0/
  26. C
  27. DATA NI/1,1,1,4,4,4/
  28. DATA NJ/1,2,3,1,2,3/
  29. DO 10 K=1,6
  30. I3=NI(K)
  31. I4=I3+2
  32. LK=NJ(K)
  33. DO 11 L=K,6
  34. REL(K,L)= 0.D00
  35. I1=NI(L)
  36. I2=I1+2
  37. LL=NJ(L)
  38. LM=0
  39. DO 12 M=I3,I4
  40. LM=LM+1
  41. S=0.
  42. II=0
  43. DO 13 N=I1,I2
  44. II=II+1
  45. C
  46. IF (II.EQ.3) THEN
  47. IF (LL.EQ.3) THEN
  48. XFACT1 =UN
  49. ELSE
  50. XFACT1 = ZER
  51. ENDIF
  52. ELSEIF (LL.EQ.3) THEN
  53. IF (II.EQ.3) THEN
  54. XFACT1 =UN
  55. ELSE
  56. XFACT1 = ZER
  57. ENDIF
  58. ELSE
  59. XFACT1 = P(II,LL)
  60. ENDIF
  61. C
  62. IF (LM.EQ.3) THEN
  63. IF (LK.EQ.3) THEN
  64. XFACT2 =UN
  65. ELSE
  66. XFACT2 = ZER
  67. ENDIF
  68. ELSEIF (LK.EQ.3) THEN
  69. IF (LM.EQ.3) THEN
  70. XFACT2 =UN
  71. ELSE
  72. XFACT2 = ZER
  73. ENDIF
  74. ELSE
  75. XFACT2 = P(LM,LK)
  76. ENDIF
  77. C
  78. 13 S=S+ RP(M,N)*XFACT1
  79. 12 REL(K,L)= REL(K,L)+ S*XFACT2
  80. 11 REL(L,K)= REL(K,L)
  81. 10 CONTINUE
  82. RETURN
  83. END
  84.  
  85.  
  86.  
  87.  

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