fortarn学渣求教读取EC文件，保存在数组中时会出现forrtl:serve(24):end-of-file d...

须佐能乎

fortarn学渣求教读取EC文件，保存在数组中时会出现forrtl:severe(24):end-of-file during read

就是一个简单的程序 ，但是当数组第二层循环j=1，1改为其他数时就会出现end-of-file during read的错误，我该怎么办啊
读取的是ec数据 数据是从第四行开始的，前三行是描述，主要问题出在程序中感叹号标注的区域
求助论坛的大大们帮助

    program text

    implicit none


    character(len = 79) :: filename
    character(len = 79) :: temp
    character(len = 79) :: dy,dx,sy,ey,sx,ex
    integer(4) :: xnum,ynum
    integer :: i,j,ii,jj
    logical :: alive
    real,allocatable :: EC_data(:,:)

    inquire(file="16061108.012",exist=alive)
    if(alive)then
        open(10,file="16061108.012")
        print*, 'file exist'
    else
        print*, 'file dosen`t exist'
        stop
    endif

    allocate(EC_data(xnum,ynum))
   
    EC_data=0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    do i=1,3
        read(10,*)temp
        print*, temp
    end do

    do i=1,xnum
        do j=1,ynum
            read(10,*) EC_data(i,j)
            write(*,"(1X,F5.2)") EC_data(i,j)
        end do
    end do
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

    close(10)
    open(20,file="text1.txt")
    do
    write(20,"(1X,F5.2)") EC_data
    stop
    end do

    end program text

HSY

数据还没读完就结束了，应该是数组长度长了，比如你循环到i=100，但文件里面只到99

HSY

额，说错了，提前读完了

hengkai1993

想请教楼主EC文件是啥？

zhouling

fortarn？是什么鬼玩意

zhimin0122

厉害啊。。。。。。。。。。。

问题少女

我不知道你EC数据的排列格式是怎样的，如果是xnum*ynum格式的话，我觉得应该是这样读取的：
do i=1,xnum
    read(10,*) (EC_data(i,j),j=1,ynum)
end do
这种情况就是要这样读取的，你可以试一下。

aas314

学习学习  赞

须佐能乎

hengkai1993 发表于 2016-6-19 08:45
想请教楼主EC文件是啥？

是欧洲中心的预报文件啊

须佐能乎

hengkai1993 发表于 2016-6-19 08:45
想请教楼主EC文件是啥？

是欧洲中心的预报文件啊