利用Fortran做单一站点多时次txt转换grd时候出错！

二子搞研究

之前参考的清风大神的贴子，但是因为我是单一站点多时次的（见下图），stid部分是我自己补上的，因为没有这个好像grads就不承认是站点数据（个人理解）。数据总共是338行，即总时次为338，之前数据的时间序列是不连续的，-32766部分是我自己补上的缺测值方便数据时间上连续，好处理一些

Fortran部分：
  Program sta2grd
    parameter (itt=338)
    character*8 stid
    character lat,lon,value(itt),tim   #这个地方我写real老出错，所以改成了character
    integer nlev,flag,j
    open(1,file='F:\cygwindata\lagcor\fortran\value.txt',status='old')
      do j=1,itt
       read(1,*)   
#这个地方我不知道到底要不要加上，因为第一行确实没数据啊，从第二行开始的，但看论坛里其他人的有些是没有这些的，但加上这个就会出错！提示"end of file encountered"
       read(1,*)stid,lat,lon,value(j)
       end do
       close(1)
    open(2,file='F:\cygwindata\lagcor\fortran\value.grd',status='replace',form='binary')
    do j=1,itt
    tim=0.0
    nlev=1
    flag=1
    write(2)stid,lat,lon,tim,nlev,flag,value(j)
    enddo
    nlev=0
    write(2)stid,lat,lon,tim,nlev,flag
    close(2)
      end
ctl部分

最后提示错误：

求大神指点！

nunu18

你这个还设置要贡献，谁会帮你解决问题啊

二子搞研究

nunu18 发表于 2016-3-1 11:22
你这个还设置要贡献，谁会帮你解决问题啊

不好意思啊，我的本意是想给解决问题的人2个贡献呢，设置错了！现在改过来了

Lancelot

还没尝试运行，也很久没写站点资料了，但是很明显，character的数据属性肯定是错的，而且很明显fortran读入的时候lon和lat是 不对应的，所以你先改掉第一处错误试试看，还有，检查一下你的数据文件是不是只用【一个】tab或者空格分割数据，以上

wongsubao

do j=1,itt
       read(1,*)   
#这个地方我不知道到底要不要加上，因为第一行确实没数据啊，从第二行开始的，但看论坛里其他人的有些是没有这些的，但加上这个就会出错！提示"end of file encountered"

原因，read(1,*)  应该在循环外，否则读数据会跳行读。

nunu18

  character lat,lon,value(itt),tim ,是REAL，如果出错你把VALUE(ITT)单独拿出来放一行。
提示"end of file encountered'是因为你的数据溢出了，所以你循环的时候(itt=338+1),就行了。加上implicte none 避免程序出错。这样试试

二子搞研究

Lancelot 发表于 2016-3-1 11:43
还没尝试运行，也很久没写站点资料了，但是很明显，character的数据属性肯定是错的，而且很明显fortran读入 ...

亲亲，多谢，确实txt文件里有两列数据是两个空格的，后来删掉了，再结合其他人的意见，改成如下的版本
                  Program sta2grd
                 parameter (itt=338)
                 character*8 stid
                real lat,lon,tim
               real iso(itt)  #把它单独拿出来定义了
                integer nlev,flag,j
     open(1,file='F:\cygwindata\lagcor\fortran\iso2.txt',status='old')
                 read(1,*) #放在了循环的外面
               do j=1,itt
                 read(1,*)stid,lat,lon,iso(j)
               end do
               close(1)
      open(2,file='F:\cygwindata\lagcor\fortran\value.grd',status='replace',form='binary')
                 tim=0.0
                nlev=1
                flag=1
              do j=1,itt
                  write(2)stid,lat,lon,tim,nlev,flag,iso(j)
              enddo
                nlev=0
                write(2)stid,lat,lon,tim,nlev,flag
                close(2)
               end
还是用的之前的ctl文件，出来的错误，看起来比之前好了一些，但还是有错误，说possible causes：premature EOF

二子搞研究

nunu18 发表于 2016-3-1 12:04
character lat,lon,value(itt),tim ,是REAL，如果出错你把VALUE(ITT)单独拿出来放一行。
提示"end of fi ...

已经按照亲你说的改了，回复在楼下，还是出错了

二子搞研究

@river 呼叫大神！

lqouc

二子搞研究 发表于 2016-3-1 14:05
亲亲，多谢，确实txt文件里有两列数据是两个空格的，后来删掉了，再结合其他人的意见，改成如下的版本
  ...

再认真看下清风的帖子哈。
你想要的是1个站点338个时刻，现在你写的方法表示的是338个站点1个时刻。改改循环的位置就行了。