Discussion :

[yoyo24]

Bonjour

Je code en fortran 77.
Je possède un fichier de données de plus de 700000 lignes, composés d’une vingtaine de colonne (x,y,z,P,V........)
Les lignes sont regroupées par x constant.pour chaque x il y a entre 400 et 5000 lignes.

L’utilisateur doit choisir un ou plusisieurs x, et le programme doit permettre de retrouver ces valeurs, et de recopier toutes les lignes correspondantes pour ce x dans un autre fichier.

Mon problème est que les différentes valeurs de x sont séparées par quelques lignes de textes.

Je fais une boucle do itérant du début du programme jusqu'à la fin pour rechercher la valeur de x correspondante. La valeur que je recherche est donc un REAL, mais des que mon programme rencontre les lignes de textes entre les différentes sections, il plante.

J’aimerai donc savoir si il est possible dans cette boucle de créer une condition afin de savoir si ma valeur est bien un réel et non une chaîne de caractère, et si c’est le cas passer à l’itération suivante.

Voici la partie de code auquel j’avais pensé, mais je ne sais pas si il existe une instruction ou une méthode permettant de d’éviter cette chaîne de caractère.

Code :

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

1
2
3
4
5
6
7
8
9
do 777 j=1,700000
         if (xc ????????)  then
            goto 777
         endif
         read(101,*) xc  
            if ((x-xplan(i)).lt.0.01)then
             numligne=j
            endif
777       continue

xplan(i) sont les valeurs choisies par l'utilisateur que je souhaite réécrire dans un nouveau fichier.

J’espère avoir été assez clair.

[Ehouarn]

Bonjour,

Pour ce que tu souhaites faire, le plus simple est de controler que la lecture de 'xc' s'est bien déroulée où non. Il existe un argument (optionnel) de 'read' fait pour ça, 'iostat', qui renvoie la valeur 0 (zéro) si la lecture s'est correctement déroulée.

Comme un exemple vaut généralement une longue dissertation, le petit code suivant illustre comment utiliser 'iostat' pour tester si la lecture d'un nombre réel s'est correctement déroulée; la transposition à ton cas est immédiate:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
      program lecture
      implicit none
      integer ierr
      real val
 
      write(*,*) "Entrez une valeur"
      read(5,*,iostat=ierr) val
      if (ierr.eq.0) then
        write(*,*) "OK, c'est bien un reel"
      else
        write(*,*) "Non, ce n'est pas un reel"
      endif
 
      end

Bonne continuation.

[yoyo24]

merci de ton aide