[Débutant] Les tableaux en Fortran

Bonjour,
Je veux faire entrer un tableau (d'une dimension) par le clavier, j'ai écrit le programme mais ça n'a pas marché; pouvez-vous me dire où est l'erreur ou où faut il modifier ? Merci
Voilà le programme :

Code:

---

1 2 3 4 5 6 7 8 9

PROGRAM mat IMPLICIT NONE integer i INTEGER, PARAMETER :: n=10 real, DIMENSION(n) :: A   do i = 1, n write(*,*) A(i) enddo

---

end

Bonjour,

write(*,*) pour écrire et read(*,*) pour lire...

Bonjour,
J ai une autre question sur les tableaux,
J arrive pas a voire la signification des indices negatifs dans les tableaux par exemple : real, dimension(-4:0,0:3) ::A
merci

c'est juste une manière d'avoir des indice négatif...

c'est utilie pour de transformation de fourier par exemple...
mais lleur utilité principale est de simplifier les formule d'indice, plutôt que de faire:

Code:

---

1 2 3

do i=1,n a(i)=i-8 enddo

---

tu peut le faire comme

Code:

---

1 2 3

do i=-8,n-8 a(i)=i enddo

---

certes sur cet exemple c'est pas flagrant mais quand tu a des recherche dans un tableau tu peux fixé le point milieu au 0 et ainsi simplifier les recherches... (par exemple)

Bonjour,

Je travaille en Fortran 90 (Visual Studio), mon problème est simpliste mais déconcertant:

Voici le contexte: Je lis un fichier ZoneK.dat contenant 7920 lignes et 2 colonnes de réels.

Mon objectif : Créer un tableau de même dimension pouvant contenir les valeurs du fichier.dat

Mon second objectif : travailler sur les valeurs au sein de ce tableau pour plus de facilité. Est-ce justifié? Je le suppose...le fichier ZoneK.dat n'est pas considéré comme un tableau une fois ouvert.

Voici un extrait de mon fichier

Code:

---

1 2 3 4 5 6 7

7.7450000e+003 2.0000000e+000 7.7460000e+003 2.0000000e+000 7.7470000e+003 2.0000000e+000 . . . . . .

---

Voici ce que je fais, vous passant les étapes inintéressantes ici :
Nb : je n'ouvre pas encore ici le fichier ZoneK.dat. Je ne fais qu'initialiser mon tableau.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

! déclaration de mon TABLEAU real, dimension (7920,2):: zone   ! Ouverture du fichier qui va contenir les valeurs de ZoneK.dat open(unit=15,file='test.dat', access='sequential', form='formatted', action = 'write', status='new', iostat=err)   ! Initialisation du TABLEAU zone(:,:) = 0   ! Et hop, on affiche le résultat do i=1,7920 write(15,'(f16.7))')(zone(i,j),j=1,2) end do

---

Mon résultat (à l'écran ou dans un fichier) : un colonne de 15820 éléments 0, soit mes deux colonnes l'une à la suite de l'autre.

J'imagine que mon erreur vous sauteras aux yeux, mais à le différence de MATLAB, je trouve Fortran moins intuitif et je ne parviens pas à voir où se situe l'erreur. Aussi, le fichier qui est crée l'est souvent en complément à un fichier créé par défaut 'fort' avec l'index du unit correspondant. Comment se limiter au fichier nommé dans le Open?

J'ose espérer que vous pourrez m'aider.

Merci d'avance.

S.

Problème #1 - Pourquoi une colonne de 15xxx éléments et non 2 colonnes de 7 9xx éléments : Parce que bien que le write produise 2 nombres, le format n'en prévoit qu'un. Le programme doit donc utiliser 2 fois le format. Et à chaque fois, il y a saut de ligne (c'est une règle de Fortran). Tu dois donc prévoir 2 entrées dans le format 2fxxx ou fxxx,fxxx.

Problème #2 - Création de fichier fort15.dat ou similaire. Dans le open, tu as mis iostat=err et status ='new'. La clause iotstat fait que s'il y a erreur, la librairie sera silencieuse et il n'y aura pas d'arrêt automatique. Le contrôle est retourné au programme. Mais dans le programme, tu ne testes pas la variable err. S'il y a erreur, le fichier n'est pas ouvert, mais tout est silencieux. Vient alors le write(15. Mais il n'y a pas de fichier associé au unit 15. Que fait la librairie ? Elle est tellement gentille, elle t'ouvre un beau fichier fort15.... La source de l'erreur, c'est que tu exige un nouveau fichier (clause status new), mais le fichier existe probablement déjà, ce qui provoque l'erreur. Solution ? Enlève la clause status si ça ne te dérange pas d'écraser le fichier à chaque fois (mais si ça te dérange, tu dois préciser ce que tu veux faire, et ça risque d'être plus compliqué...)

Citation:

---

Envoyé par Sylvain Bergeron

Problème #1 - Pourquoi une colonne de 15xxx éléments et non 2 colonnes de 7 9xx éléments : Parce que bien que le write produise 2 nombres, le format n'en prévoit qu'un. Le programme doit donc utiliser 2 fois le format. Et à chaque fois, il y a saut de ligne (c'est une règle de Fortran). Tu dois donc prévoir 2 entrées dans le format 2fxxx ou fxxx,fxxx.

Problème #2 - Création de fichier fort15.dat ou similaire. Dans le open, tu as mis iostat=err et status ='new'. La clause iotstat fait que s'il y a erreur, la librairie sera silencieuse et il n'y aura pas d'arrêt automatique. Le contrôle est retourné au programme. Mais dans le programme, tu ne testes pas la variable err. S'il y a erreur, le fichier n'est pas ouvert, mais tout est silencieux. Vient alors le write(15. Mais il n'y a pas de fichier associé au unit 15. Que fait la librairie ? Elle est tellement gentille, elle t'ouvre un beau fichier fort15.... La source de l'erreur, c'est que tu exige un nouveau fichier (clause status new), mais le fichier existe probablement déjà, ce qui provoque l'erreur. Solution ? Enlève la clause status si ça ne te dérange pas d'écraser le fichier à chaque fois (mais si ça te dérange, tu dois préciser ce que tu veux faire, et ça risque d'être plus compliqué...)

---

Je réitère mon Merci.

Problèmes #1 et #2 résolus.

Merci pour votre aide et votre assistance.

S.