Discussion :

[feynman]

Bonjour,
si j'ai un fichier data qui contient 3 colonnes : la première, c'est (i*h) qui représente la discrétisation de l'axe des x, la deuxieme c'est (j*h) qui représente la discrétisation de l'axe des y, et la troisième c'est A(i,j) qui représente la matrice d'une courbe donnée.

Ma question est : si je veux réutiliser la matrice A(i,j) dans un programme f77 donné, comment je la lis à partir de ce fichier .dat ?
Merci

[genteur slayer]

est_ce que tu as autant de lignes sur chaque colonne, parce que si c ça c'est trrop facil...

cela dit si c'est des nombre différent...

bon je suppose que ton fichier c'est:

Code :

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

1
2
3
4
5
6
7
 
...
i j-1   A(i,j-1)
i j      A(i,j)
i j+1   A(i,j+1)
...
i+1 j  A(i+1,j)

et bien tu fait un read tout bete:

Code :

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

1
2
3
4
5
6
 
do i=1,Nx
  do j=1,NY
    read(11,*)x(i),y(j),A(i,j)
  end do
end do

les x et y seront lus plusieur fois mais en principe ces valeurs seront les même, sinon il faut les indicer en i,j

[feynman]

pour ta premiere question oui,

[feynman]

c est a dire je fait

Code :

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

1
2
3
4
5
6
7
8
9
10
 
program erreur
integer i,j,N
real A(0:N,0:N),x(0:N),y(0:N)
open(1,file='s.dat',status='unknown')
do i=0,N
      do j=0,N
           read(1,*) x(i),y(j),A(i,j)
       enddo
enddo

?

[genteur slayer]

en plus tu as autan de ligne que de colonne??? lol

sinon, bin cela dépend fortement de comment est écrit ton fichier dat...

mais c'est ça l'idée, oui!

[feynman]

et ensuite je met close(1)?
sinon, comment je quites du fichier 1?

[genteur slayer]

heu oui, évidement...
je te conseil aussi de mettre 'old' au lieu de 'unknown' dans le open: comme on veuit le lire, c'est qu'il existe déjà...