Variation d'une condition dans un IF

Version imprimable

Bonjour
j'ai commencé à utiliser Fortran y a quelques semaines seulement, et étant donné qu'on l'utilise que faire du calcul de physique stats, on a vu que les choses de bases.

Voici mon programme :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 program analyse implicit none integer i double precision somme somme=0 open(11,file='valeur_vitesse') do i=1,188 read(11,*)vitesse(i) if(vitesse(i).lt.-1000)then somme=somme+1 endif enddo print *,'La cumulee pour Vx=-1000 est',somme/188 end
Mon fichier valeur_vitesse est composé d'un unique colonne avec des valeurs de vitesses.

Ce que je voudrais, c'est faire varier la valeur -1000 dans if(vitesse(i).lt.-1000) entre -1000 et 1000.

J'avais essayé de rajouter une boucle pour le faire varier,comme ci-dessous, mais ça a jamais fonctionné.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 program analyse implicit none integer i integer vx double precision somme somme=0 open(11,file='valeur_vitesse') do vx=-1000,1000 do i=1,188 read(11,*)vitesse(i) if(vitesse(i).lt.vx)then somme=somme+1 endif enddo print *,'La cumulee pour Vx=',vx,somme/188 enddo end
J'imagine que ça doit pas être compliqué à faire mais j'ai pas assez de vécu en Fortran pour voir où est le problème :lol:

J'espère que vous pourrez me trouver cette solution.
Merci d'avance, Frédérick.

19/11/2011, 13h05
rogue-spectre

Salut,

Le problème principal, que tu rencontre est que tu fais k boucles pour lire un fichier de k lignes (jusque là tout va bien), mais le truc c'est que ces k boucles sont dans n autres boucles, donc à la fin de la première boucle sur les n, tu arrives à la fin de ton fichier, au début de la seconde boucle sur n tu essayes donc de lire la k+1 ème ligne de ton fichier qui n'en a que k.

+ tu as oublié de déclarer le tableau vitesse (integer, dimension(ta_dimension) :: vitesse)
+ tu n'as pas fermé ton fichier (close(unit_de_ton_fichier))

Le mieux je pense est que tu fasses une première boucle pour lire ton fichier et sauvegarder les données dans un tableau puis dans une seconde boucle tu feras le traitement des données.

En espérant avoir éclairé ta lanterne, bon week end et bon travail ;)

J'ai tenté une modif avec ce que tu m'as dit mais j'ai pas réussi :mouarf:
Voilà ce que j'ai essayé :

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
 
program analyse
implicit none
 
integer i
double precision somme
integer vx
double precision vitesse(1000)
 
open(11,file='valeur_vitesse')
do i=1,188
read(11,*)vitesse(i)
enddo
close(11)
 
somme=0
open(12,file='sortie')
do vx=-1000,1000
if(vitesse(i).lt.vx)then
somme=somme+1
endif
enddo
write(12,*)vx,somme
enddo
close(12)
end

Pour avoir une bonne idée de ce que je veux faire, j'ai un programme qui fonctionne bien mais est trop long :

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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
 
program analyse
implicit none
 
integer i
integer x
double precision somme
double precision somme1
double precision somme2
double precision somme3
double precision somme4
double precision somme5
double precision somme6
double precision somme7
double precision somme8
double precision vitesse(1000)
 
somme=0
open(11,file='valeur_vitesse')
do i=1,188
read(11,*)vitesse(i)
if(vitesse(i).lt.-1000)then
somme8=somme8+1
endif
if(vitesse(i).lt.-800)then
somme7=somme7+1
endif
if(vitesse(i).lt.-500)then
somme6=somme6+1
endif
if(vitesse(i).lt.-300)then
somme=somme+1
endif
if(vitesse(i).lt.0)then
somme1=somme1+1
endif
if(vitesse(i).lt.233)then
somme2=somme2+1
endif
if(vitesse(i).lt.541)then
somme3=somme3+1
endif
if(vitesse(i).lt.700)then
somme4=somme4+1
endif
if(vitesse(i).lt.911)then
somme5=somme5+1
endif
enddo
 
print *,'La cumulee pour Vx=-1000 est',somme8/188
print *,'La cumulee pour Vx=-800 est',somme7/188
print *,'La cumulee pour Vx=-500 est',somme6/188
print *,'La cumulee pour Vx=-300 est',somme/188
print *,'La cumulee pour Vx=0 est',somme1/188
print *,'La cumulee pour Vx=233 est',somme2/188
print *,'La cumulee pour Vx=541 est',somme3/188
print *,'La cumulee pour Vx=700 est',somme4/188
print *,'La cumulee pour Vx=911 est',somme5/188
end

Si quelqu'un peut m'aider ! :cry: :lol:
Merci d'avance.

19/11/2011, 21h39
rogue-spectre

… tu as tous les ingrédients, il te manque juste… un peu de sommeil ;) somme1, somme2, somme3 et si on mettait tout ça dans un tableau on pourrait ensuite faire une boucle et réduire le code ! + un autre truc qui peut être sympathique est de faire un tri de ta liste de vitesse pour aller encore plus vite puisque tu n'aurais plus qu'un seul balayage de ton tableau à faire pour savoir combien de vitesses de ta liste sont inférieures à une vitesse donnée.

Bonne nuit !

Problème résolu ! :lol:

Je poste mon programme, si ça peut aider quelqu'un plus tard, qui sait.
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 program analyse implicit none integer i integer vx double precision somme double precision vitesse(1000) do vx=-2500,1500 somme=0 open(11,file='valeur_vitesse') do i=1,188 read(11,*)vitesse(i) if(vitesse(i).lt.vx)then somme=somme+1 endif enddo open(12,file='resultat') write (12,*),vx,somme/188 rewind(11) enddo close(11) close(12) end
Je ne sais pas si c'est bien écrit mais dans tous les cas, ça fonctionne comme je le souhaite :D
Merci à toi de m'avoir éclairer ;)

Bon dimanche, Frédérick.