Einlesen einer FORTRAN Binärdatei mit Delphi

[Uwe Raabe]

Zitat von Tenobaal:

Problem gelöst. Ich habe "Stream.Position:=1;" vor den ersten Read-Befehl meiner Delphi Prozedur gehangen, damit erst ab dem 2 Byte gelesen wird

Das sieht mir aber nicht nach einer generellen Lösung aus. Zeig doch mal den Hexdump der gesamten Fortran-Datei oder hänge sie hier an.

[Tenobaal]

War auch zu früh gefreut. Bei mehreren Datentypen in der Binärdatei funktioniert die Lösung nicht mehr. Scheinbar werden weitere EOT-Zeichen wahllos eingestreut

[Klaus01]

.. kann es sein, dass jedes Integer*4 Datenfeld mit 04 beginnt oder ist es ein Feldtrenner?

Code:

04 DF 02 00 00 04

Anhand der wenigen Daten ist das schwer zu erkennen.

Grüße
Klaus

[Uwe Raabe]

Zitat von Tenobaal:

Scheinbar werden weitere EOT-Zeichen wahllos eingestreut

Na eben nicht wahllos! Aber wir können anhand der gezeigten sechs Bytes natürlich auch nur was ahnen.

[p80286]

Zitat von Uwe Raabe:

Das sieht mir aber nicht nach einer generellen Lösung aus. Zeig doch mal den Hexdump der gesamten Fortran-Datei oder hänge sie hier an.

Gruß
K-H

[Tenobaal]

Hier eine Funktion die eine Double, zwei Integer Zahlen und ein Array der Länge 3 vom Datentyp Char in eine Binärdatei schreibt.
Inhalt der Binärdatei aus dem Hexeditor:

Code:

13000000000000F03F020000000200000041424313

Programmcode:

Code:

      SUBROUTINE OUTNET
C**********************************************************************
C
C
C
C IMPLICIT
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
      COMMON /STEUER/ IERR,LIN,KPER,NPRINT,MPRINT,MAXP,MAXT,IRECH,IB2,
     ,                NSPG,NROT1,NROT2,NSTROM
      COMMON /TEXTE/  TITEL
      COMMON /NETZP/  NPMAX,NPUNKT,NPVAR,NZMAX,NZELLE,NSMAX,NSTR,MMAX
      COMMON /NETZD/  X(10000),Y(10000),A(10000),NUE(20000),SD(20000),
     ,                GAM(20000),
     ,                D(20000),GK(3,20000),NACHBZ(60000),ENDZ(10000),
     ,                NPER(10000),MGEB(20000),
     ,                ECKP(3,20000),RANDZ(3,20000),RANDS(3,20000),
     ,                AP(30000),EP(30000),GL(30000),GR(30000),
     ,                RAND(30000),ZL(30000),ZR(30000)
      COMMON /KONTUR/ GECKP(3,4000),IGECKP
      COMMON /IODEF / VERSIO,NTEST,IKART,NLIST,ND1,INET,NNET,NPLOT,
     ,                NUEPL

      DOUBLE PRECISION    NUE,GAM
      COMPLEX*16 A,SD,AComplex
      CHARACTER*3 TITEL
      REAL*8 VERSIO
      INTEGER*4 NACHBZ,ENDZ,NPER,MGEB,ECKP,RANDZ,RANDS
      INTEGER*4 AP,EP,GL,GR,RAND,ZL,ZR

      ADOUBLE=1.0 !Double 64-Bit
      NINTEGER=2  !Integer 32-Bit
      TITEL='ABC' ! Array of Char
      
      WRITE (NNET)  ADOUBLE,NINTEGER,NINTEGER,TITEL !Schreibe in Binärdatei
      RETURN
      END

[Uwe Raabe]

Die $13 sind offensichtlich die Länge des Records (19 Zeichen). Warum die am Ende nochmal kommt, weiß ich auch nicht. Mach doch mal eine Datei mit mehreren Datensätzen. (Mit Fortran hatte ich das letzte Mal vor 40 Jahren zu tun...)

[BUG]

Zitat von Uwe Raabe:

Die $13 sind offensichtlich die Länge des Records (19 Zeichen). Warum die am Ende nochmal kommt, weiß ich auch nicht.

Ein Suchergebnis hat Rückwärtsspringen über Records vermutet ... würde Sinn machen.

@tenobaal: Du scheinst den implementierungsabhängigen Unformatted-Mode zu benutzen, insofern wirst du keine Lösung finden, die mit allen Fortran-Compilern kompatibel ist.

[Uwe Raabe]

Zitat von BUG:

Zitat von Uwe Raabe:

Die $13 sind offensichtlich die Länge des Records (19 Zeichen). Warum die am Ende nochmal kommt, weiß ich auch nicht.

Ein Suchergebnis hat Rückwärtsspringen über Records vermutet ... würde Sinn machen.

Das wäre denkbar. Ein ähnliches Format hatte ich schon mal bei Daten für Magnetbänder.

[Mavarik]

Also für mich sieht das ganz und gar nicht nach eine Binärdatei aus...

Schau Dir mal das hier an...

Mavarik