Discussion :

[manaiilhem]

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__]

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]

Merci pour votre réponse

[manaiilhem]

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]

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__]

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)

[__dardanos__]

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]

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]

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__]

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]

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