Hướng dẫn python read fortran formatted file - python đọc tệp được định dạng fortran

Một đối tượng tệp cho các tệp tuần tự không được định dạng từ mã Fortran.

Mở đối tượng tệp hoặc tên tệp.

Chế độ đọc-viết, mặc định là ‘r.

Kiểu dữ liệu của tiêu đề. Kích thước và độ cuối phải phù hợp với tệp đầu vào/đầu ra.

Ghi chú

Các tệp này được chia thành các bản ghi của các loại không xác định. Kích thước của mỗi bản ghi được đưa ra khi bắt đầu (mặc dù kích thước của tiêu đề này không phải là tiêu chuẩn) và dữ liệu được ghi vào đĩa mà không có bất kỳ định dạng nào. Các trình biên dịch Fortran hỗ trợ câu lệnh Backspace sẽ viết một bản sao thứ hai có kích thước để tạo điều kiện tìm kiếm ngược.

Lớp này chỉ hỗ trợ các tệp được viết với cả hai kích thước cho bản ghi. Nó cũng không hỗ trợ các subrecord được sử dụng trong các trình biên dịch Intel và GFORTRAN cho các bản ghi lớn hơn 2GB với tiêu đề 4 byte.

Một ví dụ về một tệp tuần tự không được định dạng trong Fortran sẽ được viết là:

OPEN(1, FILE=myfilename, FORM='unformatted')

WRITE(1) myvariable

Vì đây là định dạng tệp không chuẩn, có nội dung phụ thuộc vào trình biên dịch và tính nội dung của máy, nên thận trọng. Các tệp từ GFORTRAN 4.8.0 và GFORTRAN 4.1.2 trên x86_64 được biết là hoạt động.

Xem xét sử dụng các tệp hoặc tệp truy cập trực tiếp Fortran từ luồng I/O mới hơn, có thể dễ dàng đọc bởi

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

Ví dụ

Để tạo một tệp Fortran tuần tự không được định dạng:

>>> from scipy.io import FortranFile
>>> f = FortranFile('test.unf', 'w')
>>> f.write_record(np.array([1,2,3,4,5], dtype=np.int32))
>>> f.write_record(np.linspace(0,1,20).reshape((5,4)).T)
>>> f.close()

Để đọc tệp này:

>>> f = FortranFile('test.unf', 'r')
>>> print(f.read_ints(np.int32))
[1 2 3 4 5]
>>> print(f.read_reals(float).reshape((5,4), order="F"))
[[0.         0.05263158 0.10526316 0.15789474]
 [0.21052632 0.26315789 0.31578947 0.36842105]
 [0.42105263 0.47368421 0.52631579 0.57894737]
 [0.63157895 0.68421053 0.73684211 0.78947368]
 [0.84210526 0.89473684 0.94736842 1.        ]]
>>> f.close()

Hoặc, trong Fortran:

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

Phương pháp

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

>>> float('4.50000+1')
Traceback (most recent call last):
  File "", line 1, in 
ValueError: invalid literal for float(): 4.50000+1

>>> import fortranformat as ff
>>> reader = ff.FortranRecordReader('(6F13.7)')
>>> reader.read(' 4.500000+1 1.894719-3 4.600000+1 8.196721-3 4.700000+1 2.869539-3')
[45.0, 0.001894719, 46.0, 0.008196721, 47.0, 0.002869539]

>>> import fortranformat as ff
>>> reader = ff.FortranRecordReader('(6F13.7)')
>>> reader.read(' 4.500000+1 1.894719-3 4.600000+1 8.196721-3 4.700000+1 2.869539-3')
[45.0, 0.001894719, 46.0, 0.008196721, 47.0, 0.002869539]