lecture de plusieurs paramètres, fortran 77

Version imprimable

05/07/2012, 11h47
Boris09

lecture de plusieurs paramètres, fortran 77

Bonjour à tous,

je suis débutant en programmation et je me retrouve avec un programme en fortran 77, qui m'a été fourni. En fait j'ai 2 fichiers: un premier contenant mon programme et un second contenant mes données.
Pour décrire un peu, mon fichier données contient une première ligne de descriptif, une seconde avec différents paramètres et plus d'une centaine contenant les données.
Mon problème se situe au niveau de la lecture de la seconde ligne, à travers mon programme. Le programme ne reconnait pas les différents paramètres et n'arrivent pas à les lire.

Je ne sais pas si j'explique bien et je vais prendre un exemple:

Voici mon fichier de données:

Pau 31 Mai 1986
-81.16 52.30 160 5
-90.60 53.20 6
-80.20 45.95 6
-79.62 65.50 6
...etc..

Je veux lire la seconde ligne contenant les paramètres, donc après description des variables, input, output, constantes ... je lis la première ligne:

read(iunit,'(a)') eqid

jusque là ça marche, mais ensuite le problème arrive. Je veux lire les 4 paramètres de la seconde ligne et que je veux appeler xlon, xlat, spread et del:

read(iunit,*) xlon, xlat, spread, del

et là ça plante.

Je n'arrive pas à trouver d'informations concernant la prise en compte de plusieurs paramètres avec la fonction read et je ne comprends pourquoi ça ne fonctionne pas. Le programme sait-il que je fais référence à cette ligne là?

Ce programme m'a été fourni par quelqu'un de très sérieux, m'ayant certifié qu'il fonctionnait parfaitement.

Merci pour toutes informations ou conseils que vous pouvez me donner
05/07/2012, 13h54
FR119492

Salut!
Essaie:

Code:

read (iunit,'(2F6.2,I4,I2)') xlon,xlat,spread,del

Jean-Marc Blanc
05/07/2012, 15h34
Grame

Bonjour

A ce qu'a écrit Jean-Marc, j'ajouterais :

As-tu bien déclaré tes variables xlon et xlat réelles et spread et del entières ?
05/07/2012, 16h24
Boris09

j'ai bien déclaré les variables comme entiers et réels. J'ai changé le format comme proposé par Jean-Marc, mais ça ne marche toujours pas. J'ai aussi changé le nom de la variable spread car c'est le nom d'une fonction fortran.

Petite question supplémentaire: j'utilise la ligne print '(2a)', 'marche' pour savoir à quelle ligne mon programme s'arrête. Est ce une bonne solution? Ça a fonctionné jusqu'à cette ligne

merci pour votre aide
05/07/2012, 16h33
Grame

Peux-tu nous poster ton programme (entre balises pour faciliter la lecture) ?

S'il est trop long (> 50-60 lignes), poste juste le début.

Voici le début du programme

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
61
62
63
64
      program pro
      real x,y
      parameter (ndatdim=3000)
      parameter (ndim=161)
      character*28 eqid,infile,outfile
      real zlat(ndatdim), zlon(ndatdim),  zmmi(ndatdim)
      real xlat(ndatdim), xlon(ndatdim),  xmmi(ndatdim)
      real wts(ndatdim)
      real xloneq, xlateq
      integer sprd, dellat
c  xlat(i), xlon(i) are the coords for intensity xmmi(i) 
      real  xmii3(ndatdim),r3(26500), s3(26500)
      real  glat(26500),glon(26500)
c  xmii1, xmii3, xmii4 are scratch arrays
      integer iunit,ounit1,ounit2
      parameter (iunit=21, ounit1=22, ounit2=23)
      pio180=4.0*atan(1.0)/180.
      Distmax=150.
c Open infile...
 1    print '(a,$)', ' Infile:  '
      read '(a)', infile
      open (iunit, file=infile, form='formatted', status='old',
     &iostat=ios)
      if (ios.ne.0) then
         print *, ' Something wrong with infile = ', infile
         print *, ' ** Try again**'
        go to 1
      endif
 
c     Open output files...
    2 print '(a,$)', ' Outfile for MI3:  '
      read '(a)', outfile
      open (ounit1, file=outfile, form='formatted', status='new', 
     &  iostat=ios)
      if (ios.ne.0) then
         print *, ' Something wrong with outfile = ', outfile
         print *, ' ** Try again **'
        goto 2
      endif
    3 print '(a,$)', ' Outfile for RI3:  '
      read '(a)', outfile
      open (ounit2, file=outfile, form='formatted', status='new', 
     &  iostat=ios)
      if (ios.ne.0) then
         print *, ' Something wrong with outfile = ', outfile
         print *, ' ** Try again **'
        goto 3
      endif
 
c     Read header information.
      read(iunit,'(a)') eqid
 
c     Print header info at terminal:
      print '(2a)', ' id = ', eqid
 
c     Write header to outfile.
      write(ounit1,'(a)') eqid
      eqid(21:28)='RMS MI-3'
      write(ounit2,'(a)') eqid
      read (iunit,'(2F6.2)',I4,I2)  xloneq, xlateq,sprd,dellat
      print '(2a)', 'marche'
 
      write(ounit1,'(1x,4f10.2)') xloneq,xlateq,sprd,dellat
      write(ounit2,'(1x,4f10.2)') xloneq,xlateq,sprd,dellat

et le début du fichier de données

Pau 31 Jan 1986
-81.16 41.65 160 5
-81.35 41.4 6
-83.65 41.37 6
-81.14 41.61 6
-80.94 41.8 6
-81.28 41.73 6
-81.05 41.53 6
-81.36 41.62 6
-81.14 41.68 6
-81.06 41.78 6
-81.35 41.66 6
-81.28 41.43 6
-81.07 41.46 6
-81.26 41.73 6

05/07/2012, 16h59
Boris09

pardon pour les balises.

J'ai oublié de vous dire le message d'erreur:

read start: end of file
apparent state: unit 21 named OH.1
last format: (2F6.2,I4,I2)
lately reading sequential formatted external IO
Aborted (core dumped)
05/07/2012, 23h43
FR119492

Salut!
Je vois en tout cas une erreur: tu écris les variables sprd et dellat en format f alors qu'elles sont déclarées comme entières. Il est possible qu'il y en ait une autre.
Jean-Marc Blanc
06/07/2012, 11h18
Grame
Bonjour

Quelques erreurs à corriger.

- ligne 60 : tu as mal écrit le format
Remplace
Code:

1 2 read (iunit,'(2F6.2)',I4,I2) xloneq, xlateq,sprd,dellat
par
Code:

1 2 read (iunit,'(2F6.2,I4,I2)') xloneq, xlateq,sprd,dellat
- lignes 63 et 64
Remplace
Code:

1 2 3 write(ounit1,'(1x,4f10.2)') xloneq,xlateq,sprd,dellat write(ounit2,'(1x,4f10.2)') xloneq,xlateq,sprd,dellat
par
Code:

1 2 3 write(ounit1,'(1x,2f10.2,i4,i2)') xloneq,xlateq,sprd,dellat write(ounit2,'(1x,2f10.2,i4,i2)') xloneq,xlateq,sprd,dellat
sprd et dellat sont des entiers que tu essaies d'écrire en format réel.

Par ailleurs, je ne sais pas quel OS et compilateur tu utilises mais les constructions du type
Code:

1 2 print '(a,$)'
ne sont pas standards : le $ qui évite le saut de ligne est accepté par un certain nombre de compilateurs comme une extension mais ne fait pas partie de la norme FORTRAN77. Il est donc plus correct et plus sûr d'écrire :
Code:

1 2 print '(a)'
09/07/2012, 11h38
Boris09

merci pour vos réponses et votre aide. J'ai changé mon programme avec vos suggestions. Le programme ne tourne toujours pas et je ne vois vraiment pas d'où vient l'erreur. J'ai passé le week-end avec un livre de fortran 77 mais je n'ai pas trouvé d'où ça pourrait venir...

Je ne sais pas ce que tu entends par "ne tourne pas" mais en simplifiant et surtout en regardant que le fichier de données fournit 2 réels et 2 entiers dans sa première ligne puis 1 seul entier ensuite :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 59 read (iunit,*) xloneq, xlateq,sprd,dellat 60 write(ounit1,*) xloneq,xlateq,sprd,dellat 61 write(ounit2,*) xloneq,xlateq,sprd,dellat 62 10 read (iunit,*,end=20) xloneq, xlateq,sprd 63 write(ounit1,*) xloneq,xlateq,sprd 64 write(ounit2,*) xloneq,xlateq,sprd 65 goto 10 66 20 continue 67 close (iunit) 68 close (ounit1) 69 close (ounit2) 70 end
j'obtiens les sorties suivantes

Citation:

bash-4.2$ ./a.out
Infile:
data
Outfile for MI3:
o1
Outfile for RI3:
o2
id = Pau 31 Jan 1986
bash-4.2$ cat o1
Pau 31 Jan 1986
-81.160004 41.650002 160 5
-81.349998 41.400002 6
-83.650002 41.369999 6
-81.139999 41.610001 6
-80.940002 41.799999 6
-81.279999 41.730000 6
-81.050003 41.529999 6
-81.360001 41.619999 6
-81.139999 41.680000 6
-81.059998 41.779999 6
-81.349998 41.660000 6
-81.279999 41.430000 6
-81.070000 41.459999 6
-81.260002 41.730000 6
bash-4.2$ cat o2
Pau 31 Jan 1986 RMS MI-3
-81.160004 41.650002 160 5
-81.349998 41.400002 6
-83.650002 41.369999 6
-81.139999 41.610001 6
-80.940002 41.799999 6
-81.279999 41.730000 6
-81.050003 41.529999 6
-81.360001 41.619999 6
-81.139999 41.680000 6
-81.059998 41.779999 6
-81.349998 41.660000 6
-81.279999 41.430000 6
-81.070000 41.459999 6
-81.260002 41.730000 6
bash-4.2$

Pour moi, "ca marche"

09/07/2012, 15h28
Boris09
En exécutant le même programme, j'obtiens les fichiers de sortie suivants:
Code:

1 2 3 4 5 bash-4.2$ cat o1 -81.16 4Jan 1986 bash-4.2$ cat o2 RMS MI-3Jan 1986
Là je ne comprends plus du tout.

Je t'envoie mon code complet, il t'en manque peut être un bout puisque tu n'as pas l'air d'avoir de boucle de lecture.

Quel compilateur utilises-tu et (éventuellement) avec quelles options ?

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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
 
     program pro
 
      real x,y
      parameter (ndatdim=3000)
      parameter (ndim=161)
      character*28 eqid,infile,outfile
      real zlat(ndatdim), zlon(ndatdim),  zmmi(ndatdim)
      real xlat(ndatdim), xlon(ndatdim),  xmmi(ndatdim)
      real wts(ndatdim)
      real xloneq, xlateq
      integer sprd, dellat
c  xlat(i), xlon(i) are the coords for intensity xmmi(i) 
      real  xmii3(ndatdim),r3(26500), s3(26500)
      real  glat(26500),glon(26500)
c  xmii1, xmii3, xmii4 are scratch arrays
      integer iunit,ounit1,ounit2
      parameter (iunit=21, ounit1=22, ounit2=23)
      pio180=4.0*atan(1.0)/180.
      Distmax=150.
c Open infile...
 1    print '(a)', ' Infile:  '
      read '(a)', infile
      open (iunit, file=infile, form='formatted', status='old',
     &iostat=ios)
      if (ios.ne.0) then
         print *, ' Something wrong with infile = ', infile
         print *, ' ** Try again**'
        go to 1
      endif
 
c     Open output files...
    2 print '(a)', ' Outfile for MI3:  '
      read '(a)', outfile
      open (ounit1, file=outfile, form='formatted', status='new', 
     &  iostat=ios)
      if (ios.ne.0) then
         print *, ' Something wrong with outfile = ', outfile
         print *, ' ** Try again **'
        goto 2
      endif
    3 print '(a)', ' Outfile for RI3:  '
      read '(a)', outfile
      open (ounit2, file=outfile, form='formatted', status='new', 
     &  iostat=ios)
      if (ios.ne.0) then
         print *, ' Something wrong with outfile = ', outfile
         print *, ' ** Try again **'
        goto 3
      endif
 
c     Read header information.
      read(iunit,'(a)') eqid
 
c     Print header info at terminal:
      print '(2a)', ' id = ', eqid
 
c     Write header to outfile.
      write(ounit1,'(a)') eqid
      eqid(21:28)='RMS MI-3'
      write(ounit2,'(a)') eqid
 
      read (iunit,*)  xloneq, xlateq,sprd,dellat
      write(ounit1,*) xloneq,xlateq,sprd,dellat
      write(ounit2,*) xloneq,xlateq,sprd,dellat
10    read (iunit,*,end=20)  xloneq, xlateq,sprd
      write(ounit1,*) xloneq,xlateq,sprd
      write(ounit2,*) xloneq,xlateq,sprd
      goto 10
20    continue
      close (iunit)
      close (ounit1)
      close (ounit2)
      end

09/07/2012, 16h04
Boris09

j'ai fait un copié collé du code et j'obtiens la même chose que précédemment... le problème peut être au niveau du compilateur? j'utilise FORT77 sous Linux avec les options - g -v et -cpp.
09/07/2012, 16h37
Grame

Citation:

j'ai fait un copié collé du code et j'obtiens la même chose que précédemment... le problème peut être au niveau du compilateur?

Du compilateur, non, mais peut être des options.

Citation:

j'utilise FORT77 sous Linux avec les options - g -v et -cpp.

Je ne connais pas FORT77. J'ai utilisé simplement gfortran sans options.
09/07/2012, 17h12
Boris09

ça marche! Bon, il y a encore des erreurs dans la suite du programme mais j'ai franchi une étape!
J'ai juste changé de compilateur, sans options, et ça marche.

Un très grand merci pour vos disponibilités et patience. Un remerciement spécial pour Grame qui m'a sorti une grosse épine du pied :D