[Numpy-discussion] example reading binary Fortran file

[David Froger]

Sorry, I still don't understand how to use FortranFile ...

============
The fortran code
============

program writeArray

implicit none
integer,parameter:: nx=2,ny=5
real(4),dimension(nx,ny):: ux,uy,p
integer :: i,j

do i = 1,nx
do j = 1,ny
   ux(i,j) = 100.  + j+(i-1.)*10.
   uy(i,j) = 200. + j+(i-1.)*10.
   p(i,j) = 300. + j+(i-1.)*10.
enddo
enddo

open(11,file='uxuyp.bin',form='unformatted')
write(11) ux,uy
write(11) p
close(11)

end program writeArray

=============
The Python script
=============

from fortranfile import FortranFile

f = FortranFile('uxuyp.bin')

x = f.readReals()

=============
The output
=============

Traceback (most recent call last):
  File "readArray.py", line 5, in <module>
    x = f.readReals()
  File
"/home/users/redone/file2/froger/travail/codes/lib/Tests/fortranread/fortranfile.py",
line 181, in readReals
    data_str = self.readRecord()
  File
"/home/users/redone/file2/froger/travail/codes/lib/Tests/fortranread/fortranfile.py",
line 128, in readRecord
    raise IOError('Could not read enough data')
IOError: Could not read enough dat


=> How to read the file 'uxuyp.bin' ?


2009/5/28 David Froger <david.froger.info at gmail.com>

> Thank you very much :-)
>
> 2009/5/28 Neil Martinsen-Burrell <nmb at wartburg.edu>
>
>> On 2009-05-28 09:32 , David Froger wrote:
>>
>>> Hy Neil Martinsen-Burrell,
>>>
>>> I'm trying the FortranFile class,
>>> http://www.scipy.org/Cookbook/FortranIO/FortranFile
>>>
>>> It looks like there are some bug in the last revision (7):
>>>
>>>    * There are errors cause by lines 60,61,63 in
>>>
>>>    * There are indentation errors on lines 97 and 113.
>>>
>>
>> There seem to have been some problems in putting the file on the wiki
>> ("Proxy-Connection: keep-alive\nCache-Control: max-age=0" seems to come from
>> an HTML communication).  I've attached my current version of the file to
>> this email.  Let me know if you have problems with this.  I will try to get
>> the working version up on the wiki.  Peace,
>>
>> -Neil
>>
>> _______________________________________________
>> Numpy-discussion mailing list
>> Numpy-discussion at scipy.org
>> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20090528/b58db76e/attachment.html>