Numérotation des lignes :

eurodi
1. C EURODI SOURCE CB215821 16/04/21 21:16:42 8920
2. SUBROUTINE EURODI(HO,ROH,FCM,S,TMAX,MAILLE,TEMPS
3. 1 ,TREL,EVO,YOUNG,JDIM)
4.
5. *
6. IMPLICIT INTEGER(I-N)
7. IMPLICIT REAL*8(A-H,O-Z)
8.
9. -INC PPARAM
10. -INC CCOPTIO
11. *
12. REAL*8 TEMPS(0:1000),REL(0:1000),YOUNG(0:1000)
13. REAL*8 TEMP0(0:100),TR1(0:4),GAMMA(0:4),TR(0:8)
14. REAL*8 MODULE(0:10),TABMOD(0:100,0:10)
15. REAL*8 TABGAM(0:10,0:4),TREL(0:10)
16. REAL*8 EVO(0:1000,0:10)
17. *
18. TINI=1.
19. TFIN=TMAX
20. DO 10 I=0,100
21. DO 11 J=0,10
22. TABMOD(I,J)=0.
23. 11 CONTINUE
24. 10 CONTINUE
25. DO 12 I=0,10
26. DO 13 J=0,4
27. TABGAM(I,J)=0.
28. 13 CONTINUE
29. 12 CONTINUE
30. *
31. * liste des t0 a calculer
32. *
33. TEMP0(0)=TINI
34. TMP=LOG((TFIN-TINI)/.25)/.5
35. NMAX=INT(TMP)
36. IF (NMAX.GT.100) THEN
37. NMAX=100
38. ENDIF
39. DO 20 I=1,NMAX
40. TEMP0(I)=.25*EXP(.5*FLOAT(I))+TINI
41. 20 CONTINUE
42. *
43. * boucle sur les t0
44. *
45. DO 30 L=0,NMAX
46. T0=TEMP0(L)
47. TINI=TEMP0(L)
48. CALL TABTMP(HO,ROH,FCM,S,T0,TINI,TFIN,TEMPS,JDIM)
49. CALL FLUREL(HO,ROH,FCM,S,T0,TEMPS,JDIM,REL)
50. *
51. CALL CALMOD(TEMPS,JDIM,REL,MAILLE,MODULE,TR)
52. IF(IERR.NE.0) RETURN
53. *
54. DO 31 I=0,MAILLE
55. TABMOD(L,I)=MODULE(I)
56. TREL(I)=TR(I)
57. 31 CONTINUE
58. 30 CONTINUE
59. *
60. DO 40 INDICE=0,MAILLE
61. CALL CALEVO(TEMP0,NMAX,TABMOD,INDICE,GAMMA,TR1)
62. IF(IERR.NE.0) RETURN
63. DO 41 J=0,4
64. TABGAM(INDICE,J)=GAMMA(J)
65. 41 CONTINUE
66. 40 CONTINUE
67. *
68. * on recupere la liste des temps dans TEMPS
69. *
70. T0=TEMP0(0)
71. TINI=TEMP0(0)
72. CALL TABTMP(HO,ROH,FCM,S,T0,TINI,TFIN,TEMPS,JDIM)
73. *
74. * on recupere la liste des gamma(t)
75. *
76. DO 50 I=0,MAILLE
77. EVO(0,I)=0.
78. DO 51 J=0,JDIM
79. EVO(J,I)=TABGAM(I,0)
80. 1 +TABGAM(I,1)*EXP(-TR1(1)*TEMPS(J))
81. 2 +TABGAM(I,2)*EXP(-TR1(2)*TEMPS(J))
82. 3 +TABGAM(I,3)*EXP(-TR1(3)*TEMPS(J))
83. 4 +TABGAM(I,4)*EXP(-TR1(4)*TEMPS(J))
84.
85. 51 CONTINUE
86. 50 CONTINUE
87. *
88. * on recupere le module d'young
89. *
90. YOUNG(0)=0.
91. DO 60 J=0,JDIM
92. YOUNG(J)=0.
93. DO 61 I=0,MAILLE
94. YOUNG(J)=YOUNG(J)+EVO(J,I)
95. 61 CONTINUE
96. 60 CONTINUE
97.
98.
99. END
100.
101.
102.
103.
104.
105.
106.

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