erreur lors de la compilation en fortran

**manaiilhem** · 08/03/2013, 10h45

Bonjour j’écris le programme ci-dessous mais au cours de la compilation j'ai eu un retour d'un message d’erreur suivant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
Error: The number of subscripts is incorrect.   [Z]
z(i)=  (dz/(2.0D0))+((i-1)*dz)
 
Error: The number of subscripts is incorrect.   [R]
r(j)= (dr/(2.D0))-(j*dr)+R_e
^

s'il vous plait quelqu'un peut me renseigner l'erreur est du à quoi ,j'ai pas en compris
merci

Module Mod_DONNEES

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
implicit none
 
real(8),parameter :: Pi= 3.1415926535897932D0
real(8),parameter :: q_evap=176E4
real(8),parameter :: T_inf=293.D0
real(8),parameter :: conds=380.0D0
real(8),parameter :: coef=500.0D0
real(8),parameter :: R_e=5E-03
real(8),parameter :: R_i=4E-03
real(8),parameter :: Lt=0.2D0
real(8),parameter :: L_evap=4.E-02
real(8),parameter :: L_adiab=11.8E-02
real(8),parameter :: L_cond=4.2E-02
 
 
integer(8), parameter ::n=3
integer(8), parameter ::m=3
end  Mod_DONNEES
 
PROGRAM MAIN
 
USE Mod_DONNEES
implicit none
 
 
 
REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: r
REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: z
 
REAL*8 :: dr,dz
 
INTEGER::i,j
 
ALLOCATE(r(1:m,1), z(1:n,1) )
 
 
open (10,file='r(j).dat')
open (11,file='z(i).dat')
 
call grille(dr,dz,r,z)
 
!***** sauvegarder r(j)*************
do j=1,m
write(10,*) r(j)
end do
!***********************************
!****** sauvegarder z(i)************
do i=1,n
write(11,*) z(i)
end do
!************************************
CONTAINS
 
!************maillage********************
subroutine grille(dr,dz,r,z)
implicit none
 
 
REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: r
REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: z
integer:: i,j
REAL*8	:: dr,dz
 
ALLOCATE  (r(1:m,1))
ALLOCATE  (z(1:n,1) )
!les pas dans l'espace
dr=(R_e-R_i)/m
dz=Lt/n
 
 
!calcul de z(i)
do i=1,n
z(i)=  (dz/(2.0D0))+((i-1)*dz)
end do
 
!calcul de r(j)
 
do j=1,m
 
r(j)= (dr/(2.D0))-(j*dr)+R_e
end do
return
end subroutine grille
 
close(10)
close(11)
 
end

**__dardanos__** · 08/03/2013, 11h02

Salut.
r et z sont des tableaux de rang 2 :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
ALLOCATE (r(1:m,1))
ALLOCATE (z(1:n,1) )

Tu les utilises comme s'ils étaient de rang 1.

Il faudrait faire par exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

REAL*8, ALLOCATABLE, DIMENSION(:) :: r

puis :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

ALLOCATE (r(n))

Et finir par :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

DEALLOCATE(r)

**manaiilhem** · 08/03/2013, 11h12

Merci pour votre réponse

**manaiilhem** · 08/03/2013, 12h02

s'il vous plait
si je vais l'utiliser dans une autre subroutine.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

call grille(dr,dz,r,z)

Je dois le re-déclarer de la même manière ?

**manaiilhem** · 08/03/2013, 12h05

j'ai écrit une autre subroutine qui nécessite l'appel à cette subroutine .

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
subroutine matrice (r,z,dr,dz,vect,A)
implicit none


REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: A

REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: vect[/COLOR]
REAL*8, ALLOCATABLE, DIMENSION(:)	:: r
REAL*8, ALLOCATABLE, DIMENSION(:)	:: z

[COLOR="Blue"]REAL*8 :: dr,dz,L_test
INTEGER:: l,l1,l2,l3,l4
INTEGER::i,j

ALLOCATE  (A(1:n*m,1:n*m))
ALLOCATE  (vect(1:n*m,1) )
ALLOCATE (r(1:m), z(1:n) )

call grille(dr,dz,r,z)
!do j=1,m
!write(*,*) r(j)
!end do

L_test=L_evap + L_adiab
!******élement1**********************
i=1
j=1
l=(j-1)*n+i
l1=(j-1)*n+i+1
l3=j*n+i

A(l,l)=(r(j)*(dr/dz)) +(r(j)*(dz/dr))-(dz/(2.0D0))
write (*,*) A(l,l)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=q_evap/((2.0D0)*Pi*conds)
!******fin élement1********************
!******élement2************************
j=1


do i=2,n-1
j=1


l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i

A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)

 if (z(i)< L_evap .and. z(i)==L_evap) then


A(l,l)=((2.0D0)*r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/2.0D0)
!A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=q_evap/((2.0D0)*Pi*conds)

else if ( z(i) > L_evap .and. z(i) < L_test .and. z(i)==L_test)  then
A(l,l)=((2.0D0)*r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/2.0D0)
!A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)= 0.0D0

else

A(l,l)=((2.0D0)*r(j)*(dr/dz))+(coef/(2.0D0*Pi*conds))+(r(j)*(dz/dr))-(dz/2.0D0)
!A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=((-1.0D0)*(coef*T_inf))/(2.0D0*Pi*conds)
end if
!end if

end do


!******fin élement2**********************
!******élement3**************************
i=n
j=1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i

A(l,l)=(r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/2.0D0)-(coef/(2.0D0*Pi*conds))
A(l,l2)=(-1.0D0)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=((-1.0D0)*(coef*T_inf))/(2.0D0*Pi*conds)
!******fin élement3***********************
!******élement4***************************
i=1

do j=2,m-1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(2.0D0*r(j)*(dz/dr))
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
!******fin élement4***********************
!******élement5***************************
do i=2,n-1
 do j=2,m-1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i
l4=(j-2)*n+i

A(l,l)=(2.0D0*r(j)*(dr/dz))+(2.0D0*r(j)*(dz/dr))
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
 end do
!******fin élément5*************************
!******élement6*****************************
i=n

do j=2,m-1

l=(j-1)*n+i
l2=(j-1)*n+i-1
l3=j*n+i
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(2.0D0*r(j)*(dz/dr))
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
!******fin élement6*************************
!******élement7*****************************
i=1
j=m

l=(j-1)*n+i
l1=(j-1)*n+i+1
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(r(j)*(dz/dr))+(dz/2.0D0)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0
!******fin élement7*************************
!******élement8*****************************
j=m
do i=2,n-1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l4=(j-2)*n+i

A(l,l)=(2.0D0*r(j)*(dr/dz))+(r(j)*(dz/dr))+(dz/2.0D0)
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
!******fin élement8**************************
!******élement9******************************
i=n
j=m

l=(j-1)*n+i
l2=(j-1)*n+i-1
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(r(j)*(dz/dr))+(dz/2.0D0)
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0
!******fin élement9************************

return
end subroutine matrice

ce qui est marqué en rouge il faut re-déclarer le r et z de cette manière ?

**__dardanos__** · 08/03/2013, 13h22

Par exemple, tu peux mettre la subroutine dans un autre module (module_grille).
Les procédures l'utilisant l'appelle en ajoutant l'instruction use module_grille en début de programme (comme pour ton module paramètres).
Il faut impérativement compiler le module avant de compiler la procédure l'utilisant.

L'allocation que tu as faite est correcte.
On peut regrouper les déclarations :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

REAL*8, ALLOCATABLE, DIMENSION(:) :: r, z

Par défaut, l'indice minimal vaut 1, on peut dans ce cas ne pas l'écrire :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

ALLOCATE (r(m), z(n) )

Et ne pas oublier :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

DEALLOCATE(r,z)

**manaiilhem** · 08/03/2013, 13h52

dans ci-dessous comment je fait mon programme et j'ai eu un retour de l'erreur suivante lors de l'exécution

mon programme est le suivant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
PROGRAM MAIN

USE Mod_DONNEES
implicit none



REAL*8, ALLOCATABLE, DIMENSION(:)	:: r
REAL*8, ALLOCATABLE, DIMENSION(:)	:: z

REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: A
REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: vect




REAL*8 :: dr,dz

INTEGER::i,j

ALLOCATE(r(m))
ALLOCATE ( z(n) )
ALLOCATE  (A(1:n*m,1:n*m))
ALLOCATE  (vect(1:n*m,1))
open (10,file='r(j).dat')
open (11,file='z(i).dat')
open (12,file='A(i,j).dat')
open (13,file='vect(j).dat')


call grille(dr,dz,r,z)


call matrice (r,z,dr,dz,vect,A)
!***** sauvegarder r(j)*************
do j=1,m
write(10,*) r(j)
end do

!***********************************
!****** sauvegarder z(i)************
do i=1,n
write(11,*) z(i)
end do
!************************************

do i=1,n*m
do j=1,n*m


write(12,*) A(i,j)
 
end do
write(12,*) '    ' 
end do
!************************************
do j=1,n*m
write(13,*) vect(j,1)

	END DO


!************************************
close(10)
close(11)
close(12)
close(13)


END PROGRAM
Module Mod_DONNEES
implicit none

real(8),parameter :: Pi= 3.1415926535897932D0
real(8),parameter :: q_evap=176E4
real(8),parameter :: T_inf=293.D0
real(8),parameter :: conds=380.0D0
real(8),parameter :: coef=500.0D0
real(8),parameter :: r_max=5E-03
real(8),parameter :: r_min=4E-03
real(8),parameter :: Lt=0.2D0
real(8),parameter :: L_evap=4.E-02
real(8),parameter :: L_adiab=11.8E-02
real(8),parameter :: L_cond=4.2E-02


integer(8), parameter ::n=3
integer(8), parameter ::m=3

CONTAINS

!************maillage********************
subroutine grille(dr,dz,r,z)
implicit none


REAL*8, ALLOCATABLE, DIMENSION(:) :: r
REAL*8, ALLOCATABLE, DIMENSION(:)	:: z

integer:: i,j
REAL*8	:: dr,dz

ALLOCATE  (r(m))
ALLOCATE  (z(n) )
!les pas dans l'espace
dr=(r_max-r_min)/m
dz=Lt/n


!calcul de z(i)
do i=1,n
z(i)=  (dz/(2.0D0))+((i-1)*dz)
end do

DEALLOCATE(z)

!calcul de r(j)

do j=1,m
r(j)= (dr*0.5D0)-(j*dr)+r_max

end do

DEALLOCATE(r)

return
end subroutine grille


!************remplissage de la matrice*************

subroutine matrice (r,z,dr,dz,vect,A)
implicit none



REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: A

REAL*8, ALLOCATABLE, DIMENSION(:,:)	:: vect
REAL*8, ALLOCATABLE, DIMENSION(:)	:: r
REAL*8, ALLOCATABLE, DIMENSION(:)	:: z

REAL*8 :: dr,dz,L_test
INTEGER:: l,l1,l2,l3,l4
INTEGER::i,j

ALLOCATE  (A(1:n*m,1:n*m))
ALLOCATE  (vect(1:n*m,1) )
ALLOCATE (r(1:m), z(1:n) )

call grille(dr,dz,r,z)
do j=1,m
write(*,*) r(j)
end do

L_test=L_evap + L_adiab
!******élement1**********************
i=1
j=1
l=(j-1)*n+i
l1=(j-1)*n+i+1
l3=j*n+i

A(l,l)=(r(j)*(dr/dz)) +(r(j)*(dz/dr))-(dz/(2.0D0))
write (*,*) A(l,l)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=q_evap/((2.0D0)*Pi*conds)
!******fin élement1********************
!******élement2************************
j=1


do i=2,n-1
j=1


l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i

A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)

 if (z(i)< L_evap .and. z(i)==L_evap) then


A(l,l)=((2.0D0)*r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/2.0D0)
!A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=q_evap/((2.0D0)*Pi*conds)

else if ( z(i) > L_evap .and. z(i) < L_test .and. z(i)==L_test)  then
A(l,l)=((2.0D0)*r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/2.0D0)
!A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)= 0.0D0

else

A(l,l)=((2.0D0)*r(j)*(dr/dz))+(coef/(2.0D0*Pi*conds))+(r(j)*(dz/dr))-(dz/2.0D0)
!A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
!A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=((-1.0D0)*(coef*T_inf))/(2.0D0*Pi*conds)
end if
!end if

end do


!******fin élement2**********************
!******élement3**************************
i=n
j=1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i

A(l,l)=(r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/2.0D0)-(coef/(2.0D0*Pi*conds))
A(l,l2)=(-1.0D0)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=((-1.0D0)*(coef*T_inf))/(2.0D0*Pi*conds)
write (*,*) vect(l,1)
!******fin élement3***********************
!******élement4***************************
i=1

do j=2,m-1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(2.0D0*r(j)*(dz/dr))
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
!******fin élement4***********************
!******élement5***************************
do i=2,n-1
 do j=2,m-1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i
l4=(j-2)*n+i

A(l,l)=(2.0D0*r(j)*(dr/dz))+(2.0D0*r(j)*(dz/dr))
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
 end do
!******fin élément5*************************
!******élement6*****************************
i=n

do j=2,m-1

l=(j-1)*n+i
l2=(j-1)*n+i-1
l3=j*n+i
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(2.0D0*r(j)*(dz/dr))
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
!******fin élement6*************************
!******élement7*****************************
i=1
j=m

l=(j-1)*n+i
l1=(j-1)*n+i+1
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(r(j)*(dz/dr))+(dz/2.0D0)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0
!******fin élement7*************************
!******élement8*****************************
j=m
do i=2,n-1

l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l4=(j-2)*n+i

A(l,l)=(2.0D0*r(j)*(dr/dz))+(r(j)*(dz/dr))+(dz/2.0D0)
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0

end do
!******fin élement8**************************
!******élement9******************************
i=n
j=m

l=(j-1)*n+i
l2=(j-1)*n+i-1
l4=(j-2)*n+i

A(l,l)=(r(j)*(dr/dz))+(r(j)*(dz/dr))+(dz/2.0D0)
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l4)=(-1.0D0)*(r(j)+(dr/2.0D0))*(dz/dr)
vect(l,1)=0.0D0
!******fin élement9************************
DEALLOCATE(r)
DEALLOCATE(z)
return
end subroutine matrice

END Module Mod_DONNEES

la ligne 94 correspond à celle en bleu
la ligne 34 correspond à celle en rose
s'il vous plait et prière me renseigner de cet erreur dont je le comprend pas
aussi un autre erreur pour celle coloré en violet c'est comme il ne reconnait pas les valeur de r(j) pour calculer les éléments de la matrice

**__dardanos__** · 08/03/2013, 14h21

Il faut que tu prennes l'habitude de mettre le code de tes posts entre balises (le boutton #).

r et z sont déclarés deux fois : dans le module et dans le programme principal. C'est une fois de trop.

Pourquoi ne pas déclarer les variables dans un module ? (comme pour Mod_DONNEES)

**manaiilhem** · 08/03/2013, 14h39

citation:
"Il faut que tu prennes l'habitude de mettre le code de tes posts entre balises (le boutton #).
merci pour cette information
la déclaration de r et z dans le programme principale n'est pas une faute car si je l’enleve il me retourne un message d'erreur pour les déclarer

**__dardanos__** · 08/03/2013, 19h34

Bonsoir Ilhem.

Ce qui pose problème, c'est que les tableaux dynamiques ont plusieurs allocations.
Les tableaux dynamiques ne doivent être alloués qu'une seule fois sur l'ensemble du programme, et pas dans chaque procédure. De la même façon, une seule déallocation suffit.

J'ai dû réorganiser le programme en modules pour y voir plus clair. Mes modules contiennent tes procédures; ils peuvent être mieux conçus. Le sous-programme matrice m'est pas propre, je te conseille de le reprendre.

**manaiilhem** · 11/03/2013, 10h10

Bonjour,
Je vous remercie pour votre organisation de mon programme.
s'il te plait j'ai pas compris la déclaration que tu l'avais fait pour r,z,dr,dz dans la subroutine maillage

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
REAL*8, intent(out) :: dr, dz
        REAL*8, DIMENSION(:), intent(out) :: r
        REAL*8, DIMENSION(:), intent(out) :: z

pourquoi t'as ajouté intent(out)aussi dans la subroutine construction de la matrice t'as mis

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 REAL*8, DIMENSION(:), intent(in) :: r
        REAL*8, DIMENSION(:), intent(in) :: z

t'as écrit (in)
je te remercie encore

**manaiilhem** · 11/03/2013, 10h12

une autre question:
est ce qu'il ya quelqu'un qui a une idée comment on pourrait tracer des résultats par matlab à travers un program écrit en fortran?
merci

**__dardanos__** · 11/03/2013, 13h21

Salut Ilhem.

J'ai ajouté des vocations (intent) aux arguments des sous-programmes. Ils servent à l'interfaçage; des erreurs sont levées si les règle suivantes ne sont pas respectées :
Avec la vocation in, les arguments doivent avoir une valeur à l'entrée de la procédure, et la valeurs ne doit pas être modifiée dans le programme.
Avec la vocation out, c'est l'inverse.

Je pense qu'il y a une erreur dans la matrice du système.

Je ne connais pas ta mise en équations, mais pour avoir déjà fait ce type de calcul, je me souviens que les éléments de matrice étaient constants : ils ne dépendaient ni de r ni de z.
Or ta matrice dépends de z (et bizarrement pas de r). C'est suspect.

Je te conseille de commencer par exprimer les conditions aux limites en fonction de r et de z dans les 9 volumes.
Puis d'écrire et de simplifier les 9 équations bilans, et ensuite de construire la matrice.

Bon courage.

**manaiilhem** · 11/03/2013, 14h23

les conditions aux limites ont été prises en considération dans les équations.
je te communiquerai un fichier pour expliquer pour un volume comment on tient compte des conditions aux limites lors de la discrétisation et je le fais pas séparément.
merci

**manaiilhem** · 11/03/2013, 16h50

s'il vous plait dans votre réorganisation j'ai trouvé les instructions suivantes pour stoker les élement de la matrice mais j'arrive pas à en comprendre.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
print '(a)', "A ="
    do i = lbound(A,1), ubound(A,1)
        print '(500e14.6)', A(i,:)
    enddo

tu veux dire quoi par do i = lbound(A,1), ubound(A,1)

merci

**__dardanos__** · 11/03/2013, 18h02

lbound (lower boundary) et ubound (upper boundary) donnent les indices mini et maxi d'un tableau. Le second argument sert pour les tableaux à plusieurs indices.

Dans notre cas, lbound(A,1) vaut 1, et ubound est équivalent à size(A,1), c'est-à-dire la taille de A suivant son premier indice.

erreur lors de la compilation en fortran

Fortran

Discussions similaires

Partager

Partager