Maillage à 3D

Version imprimable

Maillage à 3D

Bonjour,
j'ai à écrire un programme qui fait le maillage dans une geométrie cylindrique à 3dimensions, mais lors de l'exécution j'obtiens les messages suivants :

Citation:

error LNK2001: unresolved external symbol _ZK@12
fatal error LNK1120: 1 unresolved externals

Je suis débutant en Fortran et mon programme est le suivant :

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
DIMENSION  XU(20,20,20),YV(20,20,20),ZW(20,20,20)
 
		pi=3.1416
		 ni=10;       nj=5;      nk=5             
       OPEN( 6, FILE='aa.DAT')
	   xl=2.*pi; yl=0.5;  zl=2.;    zpwr=1.
	   do 100 k=1,nk+1
	   do 200 j=1,nj+1
	   do 300 i=1,ni+1
	   yr=yl*(j-1)/nj
	   xx=xl*(i-1)/ni
	   xu(i,j,k)=yr*cos(xx)
	   yv(i,j,k)=xx*sin(xx)
	   zw(i,j,k)=zl*((k-1)/nk)**zpwr
  300 	continue
  200	   continue
  100  continue
	   do 110 i=1,ni+1
  110  write(*,*)xu(i,j,k)
       do 120 j=1,nj+1
   120 write(*,*) yv(i,j,k)
       do 130 k=1,nk+1
   130  write(*,*)zw(i,j,k)
 WRITE(40,405) xu(i,j,k),yv(i,j,k),zk(i,j,k)
 405      FORMAT(3(3X,1E15.7))   
          stop   
 
        end

Merci d'avance

03/08/2007, 10h21
DaZumba
Ben tu as un zk qui tombe du ciel
Code:

1 2 WRITE(40,405) xu(i,j,k),yv(i,j,k),zk(i,j,k)
J'imagine que tu voulais ecrire zw.

Comme la variable zk n'est pas declaree, le compilateur suppose que c'est un appel de fonction, que l'editeur de lien ne trouve nul part. D'ou le unresolved external. C'est une maladresse de la syntaxe Fortran, malheureusement : l'acces a un element de tableau a la meme forme qu'un appel de fonction.
03/08/2007, 13h20
genteur slayer

et je suppose que tu as oublié le 'implicit none' au départ, du coup le compilateur ne revoie par d'erreur car il interprète le zk comme un réel, cependant comme il y a une parenthèse devant et que tu n'a pas déclarré cela comme un tableau il interprète cela comme une fonction à valeur réelle. à l'edition de lien (LNK2001) il ne trouve nul part une fonction zk donc ERREUR...

je pense aussi que tu as simplement fait une faute de frappe
04/08/2007, 07h42
solo12

oui je crois que j'ai fait une faute de frappe
merci