读数出错

胖胖

下载了欧洲中心的四维格点相对湿度数据，x,y,z,t分别为95,55,37,365。将nc转为GrADS能读的dat，用fortran将dat写为txt，结果在检验数据的时候（输出该txt各个点的值），fortran提示forrtl:severe(24):end-of-file during read,在GrADS中打开该dat时，分明有值。实在找不到哪里出错，各位帮忙看看。下面是dat转txt的程序，运行没出错。为什么读txt会出错？
program main
implicit none
!---------------------------------------------------
integer,parameter ::mx=95,my=55,mz=37,mt=365
integer                ::i,j,k,l,nt,m
real                     ::r(mx,my,mz,mt)
!--------------------------------------------------------------------
open(11,file='j:\shuju\r\2001.dat',access='direct',form='binary',recl=mx*my*mz*mt*4)
read(11,rec=1) ((((r(m,j,k,l),m=1,mx),j=1,my),k=1,mz),l=1,mt)
close(11)
  
!--------------------------------------------------------------------
open(12,file='j:\shuju\r\2001.txt')
do l=1,mt
    do k=1,mz
          do j=1,my
          write(12,'(95f15.7)')(r(m,j,k,l),m=1,mx)
          enddo
    enddo
enddo
close(12)
end program

river

首先说nc是grads可以直接读取画图的，不需要其他的处理。第二，那就不存在楼主后面说的问题了

胖胖

river 发表于 2013-10-16 20:23

                               
登录/注册后可看大图

首先说nc是grads可以直接读取画图的，不需要其他的处理。第二，那就不存在楼主后面说的问题了

我不要图，想要里面的数据拟合廓线。所以只能转数据

river

胖胖 发表于 2013-10-16 20:26

                               
登录/注册后可看大图

我不要图，想要里面的数据拟合廓线。所以只能转数据

没明白你的意思，这个程序里没有读txt文件的语句

胖胖

river 发表于 2013-10-16 20:35

                               
登录/注册后可看大图

没明白你的意思，这个程序里没有读txt文件的语句

这是dat转txt 的程序，运行没有出错，生成的txt很大，1G多，打不开。想看看里面的数据，就写了另一个读这个txt的程序，结果显示上面的错误。

mofangbao

二进制就是为了压缩数据，加速存取的，做计算干嘛要用txt？
你自己算算就知道为什么输出来会这么大了，你的输出过程没什么问题，但是不一定是你期望的，因为你可能对输入输出还不太理解，按照你的输入输出方法计算：
总共至少有：55*37*365=742775行，
每行95个小块，每个小块是15个字符，也就是95*15=1425列，加上最后的换行（假设为2个字节），那就是1427个字节的大小，再乘上行数，也就是1059939925字节，转换一下，也就是1010.8mb,也就是你的txt至少有1个多G，何必呢，二进制存储的话，只需要270mb左右。
所以，数据不一定要肉眼能看到是多少，你知道里面的结构能用来计算就行了，就算你看到数据是多少，你也不能凭肉眼算出结果来~

言深深

错误的意思是没有那么长的数据，你检查一下看看吧，应该是数据量不足所致

胖胖

mofangbao 发表于 2013-10-16 23:05

                               
登录/注册后可看大图

二进制就是为了压缩数据，加速存取的，做计算干嘛要用txt？
你自己算算就知道为什么输出来会这么大了，你的 ...

我需要各点每天的37层数据去拟合廓线（三次方函数），得到它的各项系数，再接着做下面的工作。现在只知道使用fortran拟合，自己也不知道其他方法。确实很苦恼！

胖胖

言深深 发表于 2013-10-17 08:57

                               
登录/注册后可看大图

错误的意思是没有那么长的数据，你检查一下看看吧，应该是数据量不足所致

就觉得理论上，不该出错的。我再看看，哪里不对。

胖胖

mofangbao 发表于 2013-10-16 23:05

                               
登录/注册后可看大图

二进制就是为了压缩数据，加速存取的，做计算干嘛要用txt？
你自己算算就知道为什么输出来会这么大了，你的 ...

谢谢清风，换了个方法做。原来老想着拟合的时候看不到数据做不了，是自己想错咯。
PS:灰色的界面，心情很沉痛。