Télécharger ricoq.eso

Retour à la liste

Numérotation des lignes :

  1. C RICOQ SOURCE CHAT 05/01/13 02:56:17 5004
  2. SUBROUTINE RICOQ(RE,XA,HM,E,RIGL,REP,REM,P,XE,YE,ZE,XX,YY,RAI)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. DIMENSION E(3,3),RIGL(18,18),REP(9,9),REM(6,6),
  6. 1P(3,3),XE(3),YE(3),ZE(3),XX(3,3),YY(3,3)
  7. DIMENSION RE(18,18),XA(3,3)
  8. DATA EPSI/1.E-9/
  9. *
  10. DO9I=1,18
  11. DO9J=1,18
  12. RE(I,J)=0.D0
  13. 9 RIGL(I,J)=0.D0
  14. DO 1 I=1,3
  15. XE(I)=XA(1,I)
  16. YE(I)=XA(2,I)
  17. 1 ZE(I)=XA(3,I)
  18. CALL PASSA(XE,YE,ZE, P,X13,X23,Y13,Y23)
  19. S=0.5*Y13*(X13-X23)
  20. * SPECIAL RIGI PSI Z
  21. RKZ= RAI*HM*S*EPSI
  22. CALL RICOQ0(RIGL,X13,X23,Y13,Y23,RKZ)
  23. *
  24. CALL RIMT6(REM,E,RAI,X13,X23,Y13,HM,S)
  25. DO100I=1,6
  26. DO100J=I,6
  27. 100 REM(J,I)=REM(I,J)
  28. CALL RIPP9(REP,E,RAI,X13,X23,Y13,HM,S)
  29. DO101I=1,9
  30. DO101J=I,9
  31. 101 REP(J,I)=REP(I,J)
  32. DO2K=1,3
  33. K1=K-1
  34. DO2L=K,3
  35. L1=L-1
  36. DO 3 I=1,2
  37. M=6*K1+I
  38. ML=2*K1+I
  39. DO 3 J=1,2
  40. N=6*L1+J
  41. NL=2*L1+J
  42. 3 RIGL(M,N)=RIGL(M,N)+REM(ML,NL)
  43. DO 4 I=1,3
  44. M=6*K1+2+I
  45. ML=3*K1+I
  46. DO 4 J=1,3
  47. N=6*L1+J+2
  48. NL=3*L1+J
  49. 4 RIGL(M,N)=RIGL(M,N)+REP(ML,NL)
  50. 2 CONTINUE
  51. DO 6 K=1,6
  52. K1=3*(K-1)
  53. DO6L=K,6
  54. L1=3*(L-1)
  55. DO 7 I=1,3
  56. IM=I+K1
  57. DO 7 J=1,3
  58. JM=J+L1
  59. 7 XX(I,J)=RIGL(IM,JM)
  60. CALL PRODT(YY,XX,P,3,3)
  61. DO 8 I=1,3
  62. IM=I+K1
  63. DO 8 J=1,3
  64. JM=J+L1
  65. 8 RE(IM,JM)=YY(I,J)
  66. 6 CONTINUE
  67. DO 10 I=1,18
  68. DO 10 J=1,I
  69. 10 RE(I,J)=RE(J,I)
  70. * DO 555 I=1,18
  71. * WRITE(6,77882) I
  72. *77882 FORMAT('0 LIGNE ',I4//)
  73. * WRITE(6,77881) ( RE(I,J),J=1,18)
  74. *77881 FORMAT(1X,6(1X,1PE12.5)/)
  75. *555 CONTINUE
  76. RETURN
  77. END
  78.  
  79.  

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