Zéro résultat à l'éxecution

Version imprimable

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?
:bug::bug::bug:
voici mon prog:

Code:

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)

29/11/2007, 08h25
yogitetradim
bonjour,

une piste (pas vérifiée...désolé) : ta déclaration de tableau semble ne pas indiquer le type de X,Y et Z. Ainsi, si le tableau est enregistré comme un tableau d'entiers, il se pourrait que tes 0 correspondent à la partie entière de ton résultat.

Quelquechose comme celà devrait être plus précis :
Code:

1 2 real , dimension(1:100,1:100,1:100) :: x,y,z
De plus, un conseil utilise le" Implicit none" , ça permet de bien vérifier que toutes tes variables ont bien été déclarées, ça pourrait permettre d'identifier des causes d'erreurs dans ton programme.
En espérant que ceci ait pu t'aider ;-)
29/11/2007, 09h06
gregory2526

il me semble que le type implite de ces tableaux est REAL
(si je ne me trompe pas, c'est les variables dont la première lettre est entre i et n qui sont des INTEGER)

mais effectivement, l'utilisation de "IMPLICIT NONE" permet d'une part de bien contrôler le type de toutes les variables, mais aussi d'éviter les erreurs de frappe dans le nom d'une variable.

Pour ton problème, tu devrai afficher les résultats intermédiaires (en réduisant les ni, nj et nk pour que ce soit plus simple à débugger). Tu verrai ainsi a quel moment le calcul devient faux.
09/12/2007, 17h09
solo12

Bonjour à tous;
merci yogitetradim et gregory2526 de votre aide oui j'ai déclaré mes variables comme real ca a marché pour les deux premières composantes mais la troisième n'est pas calculé
09/12/2007, 17h11
solo12

Citation:

Envoyé par solo12

Bonjour à tous;
merci yogitetradim et gregory2526 de votre aide oui j'ai déclaré mes variables comme real ca a marché pour les deux premières composantes mais la troisième n'est pas calculé (le prog est toujours executable)
09/12/2007, 18h00
FR119492

Salut !

Lorsque tu testes un programme, affiche les résultats intermédiaires en format E pour distinguer les valeurs très petites de zéro.

D'autre part, on obtient souvent de bons résultats avec la méthode de la bissection: Tu affiches des valeurs intermédiaires à peu près au milieu de ton programme. Si ces résultats sont faux, l'erreur se trouve dans la première moitié; s'ils sont justes,c'est dans la seconde. Puis tu recommences avec la moitié qui contient l'erreur, et ainsi de suite; ça ne marche pas toujours, mais assez souvent et c'est assez rapide.

Jean-Marc Blanc
10/12/2007, 18h06
solo12

Bonjour
merci FR119492 de votre aide mais je ne vois pas comment utiliser le format E ainsi que la methode de bissection ???

Salut !

Connais-tu au moins la notion de format de lecture/écriture?

Exemple:

Code:

1
2
3
4
5
6
7
8
      X = ...
      Write (*,'(A4,E12.3)') 'X : ',X
      Y = ...
      Write (*,'(A4,E12.3)') 'Y : ',Y
      Do k=1,4
        A(k) = ...
      End Do
      Write (*,'(A4,4E12.3)') 'A : ',(A(k),k=1,4)

Jean-Marc Blanc