读取1024 ^ 3个数据点的文件

时间:2013-07-17 12:28:56

标签: fortran77

我正在尝试读取具有1024 ^ 3个单精度数据点的文件。代码看起来像这样

  real, allocatable :: stored_data(:)

  character(255):: in_file
  integer n
  integer*8 length, length4

  length = n*n*n            ! number of data points
  length4 = 4*length
  in_file = "filename.dat"

  allocate(stored_data(length))

  open(10,file=trim(in_file),access="direct",
 &     form="unformatted",recl=length4)
  read(10,rec=1) stored_data
  close(10)

  deallocate(stored_data)

但是当我运行代码时,我收到以下错误:

Fortran runtime error: RECL parameter is non-positive in OPEN statement

我猜它是因为变量“length4”的值很大。 任何解决方法?

0 个答案:

没有答案