Bonjour à tous;
j'ai essayer d'executer un programme en fortran qui me fait des icrementations en 3 directions càd un maillage dans un cylindre creux,les résultats que j'obtient sont toujours des 0 je n'arrive pas à trouver où est ce l'erreur pouvez vous m'aider?

voici mon prog:
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
 
	DIMENSION  X(100,100,100),Y(100,100,100),Z(100,100,100)
 
 
		pi=4.0*atan(1.0)
		 ni=80;       nj=80;      nk=80 
	m=ni+1
	n=nj+1	   
	zp=1.5	            
 
	 zl=5.
	r1=8.
	r0=4.
	dth=2*pi/float(nj)
 
	do 100 k=1,nk
 
	   do 200 j=1,n
 
	   x(1,j,k)=(r0+(r1-r0)*(j-1)/float(nj))
	y(1,j,k)=0.0
	x(m,j,k)=x(1,j,k)
	y(m,j,k)=y(1,j,k)
  200 continue
      do 300 i=1,m
	theta=(1-i)*dth+2.0*pi
	x(i,1,k)=r0*cos(theta)
	y(i,1,k)=r0*sin(theta)
	x(i,n,k)=r1*cos(theta)
	y(i,n,k)=r1*sin(theta)
  300 
 
	continue
	z(1,1,1)=0.0
	z(i,j,k)=zl*((k-1)/nk)**zp
 
  100	continue
 
 
	    do 110 i=1,m
  110   write(*,*)x(i,j,k)
       do 120 j=1,n
  120  write(*,*) y(i,j,k)
        do 130 k=1,nk
  130  write(*,*)z(i,j,k) 
    	WRITE(*,*)x(i,j,k),y(i,j,k),z(i,j,k)