Discussion :

[kevinl]

Bonjour,
J'ai un problème lors de l'ouverture d'un fichier de données binaire.

Mon fichier a les caractéristiques suivante :
Integers, signed
Little-endian ordering
taille du fichier : 10800,4800
10800*4800=51840000

J'ai essayer de l'ouvrier en faisant :

Code :

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

1
2
open (10,file='/home/kevin/data.dir/topography.dir/DEM1kGlobalInPart.dir/a10g',&
form='unformatted',access='direct',recl=51840000)

l'erreur que j'obtiens est:
Fortran runtime error: Format present for UNFORMATTED data transfer

Merci

[Marlan]

Salut,

Je me demande juste au niveau de la lecture de quelle façon la fais tu?

Le fichier binaire a t'il été produit par fortran?

[Grame]

Bonjour

Autre piste : la machine sur laquelle tu essaies de lire ton fichier non-formatté est-elle du même type (i.e même architecture) que celle sur laquelle le fichier a été créé. Les normes Fortran 90/95 précisent bien dans le paragraphe 9.1.2 "Unformatted record"

An unformatted record consists of a sequence of values in a processor-dependent form.

[kevinl]

Merci pour vos réponse.

Il y a bien un problème lors de la lecture du fichier (si je commente la ligne ou je lis le fichier je n'ai plus d'erreur a la compilation). Je pensais lire le fichier en faisant :

Code :

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

read (10,*) hb

Il semblerai que se ne soit pas aussi simple.

Je ne sais pas comment les données ont été produite (je l'ai ai récupérées sur le site de la NOAA, http://www.ngdc.noaa.gov/mgg/topo/gltiles.html).
Cependant je sais qu'elles ont été produites pour : IBM PC-compatible byte order, http://www.ngdc.noaa.gov/mgg/topo/re...1/s11Giii.html

Dans un premier temps, j'avais du utiliser matlab et avais pu ouvrir ce fichier. Cependant cela est un peu plus difficile pour moi a faire en fortran étant donné que je suis débutant.

[Marlan]

Salut,

Le problème avec la commande que tu utilises ici (read(10,*)) est que tu demandes à ton programme de lire un fichier formatté (symbolisé par l'étoile).

Essayes plutôt avec read(18) hb, là tu lui diras de lire un fichier correspondant à l'ouverture demandée (FORM='unformatted').

Bonne journée,

Marlan

[kevinl]

Merci,
quand j'essaie read(10) hb, je suppose que le 18 était une faute de frappe, j'obtiens:
Direct access data transfer requires record number

Si j'essaie : read (10,rec=51840000) hb ça me donne :
Non-existing record number

Une suggestion?

[Marlan]

Salut,

Effectivement, le 18 était une faute de frappe.

Alors, concernant le 'rec', il ne correspond pas à la taille de ton fichier, mais au nombre de nombres enregistrés dedans (ie : le nombre de ligne).

Admettons donc que ton fichier fasse 51840000 octets, sauvegardé en réels 8 (ou double précision), alors rec devrait avoir pour valeur 810000 (51840000/64), le 64 étant 8 octets codés sur 8 bits chacun (si je ne me trompes pas).

A mon avis, l'erreur que tu fais là est d'indiquer la taille du fichier au lieu du nombre de nombres à prendre en considération.

En espérant avoir pu t'aider,

Marlan

[kevinl]

Salut et merci de ton aide.
RECL dans open correspond t'il a la meme chose que REC dans read?

Merci

[Marlan]

Alors, sauf erreur de ma pars, le 'recl' dans l'ouverture indique le nombre de nombres, celui dans le 'read' indique directement la ligne à laquelle aller lire le nombre.

ie : read(10,rec=2) VAL te dirait que VAL aura pour valeur la seconde donnée du fichier.

A plus,

Nota : regardes ce post pour plus de précision :

http://www.developpez.net/forums/d21...-acc-s-direct/