"文件结束"将ASCII转换为GRD时出错

时间:2017-12-12 12:14:46

标签: unix fortran ascii gfortran

我正在尝试运行FORTRAN程序(使用gfortran),我将ASCII文件转换为GRD文件。我希望程序读取文本文件,然后输出grd文件。

运行程序后,我收到以下错误消息:

Fortran运行时错误:文件结束

以下是我正在使用的文件:

终端命令:

gfortran ascii2grd.f90
./a.out

FORTRAN计划:

!PROGRAM TO CONVERT ASCII TO GRD 
program ascii2grd
parameter(nx=25,ny=23,np=1)
real u(nx,ny,np),temp1(nx,ny)

open(2,file='/home/weather/test.grd',
form='unformatted',access='direct',recl=nx*ny*4)  
open(1,file='/home/weather/out.dat')

 do k=1,np
 read(1,*)((u(j,i,k),j=1,nx),i=1,ny)
10  continue    
 enddo

 rec=1
 do kk=1,np
 write(2,rec=irec)((u(j,i,kk),j=1,nx),i=1,ny)
 write(*,*)'Processing...'
irec=irec+1
 enddo
 write(*,*)'Done!'

stop
end

来自out.dat的示例

   999.000000    
   11.1000004    
   999.000000    
   999.000000    
   16.7000008    
   999.000000    
   14.6000004    
   999.000000    
   999.000000    
   999.000000      

0 个答案:

没有答案