为什么读出的数据都是常数呢?

wygohappy

nc格式的资料，通过gs文件转换为grd，如下：
'reinit'
'sdfopen d:\q\vwnd.mon.mean.nc'
'set gxout fwrite'
'set fwrite d:\q\vwnd850.grd'
tt=421
while(tt<=780)
'set t 'tt
'set z 3'
'set x 25 73'
'set y 33 57'
'd vwnd'
tt=tt+1
endwhile
'disable fwrite'
;
要计算编程，因此用该grd文件，因画不出图，考虑到可能是数据的问题，我把该grd转换为dat，发现均为常数：
open(11,file='d:\q\vwnd850.grd',form='binary')
do t=1,nt
   do j=1,ny
   read(11)  (v(nx,ny,nt),i=1,nx)
   enddo
enddo
open(15,file='d:\q\vwnd850.dat',form='formatted')
do t=1,nt
   do j=1,ny
   write(15,*) (v(nx,ny,nt),i=1,nx)
   enddo
enddo
dat文件：
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998        2.649998        2.649998        2.649998
        2.649998
为什么数据始终为常数呢，特别困惑，求助各位了！！！谢谢各位的解答了！！！

lqouc

gs看上去没问题
fortran部分尽量不要用binary
open(11,file='d:\q\vwnd850.grd',form='unformatted',recl=nx*ny)
do t=1,nt
    read(11,rec=t)  ((v(nx,ny,nt),i=1,nx),j=1,ny)
enddo
但是估计问题不在这里，你看看你的数组声明的时候类型和原来数据一不一致？

wygohappy

lqouc 发表于 2013-4-27 11:27

                               
登录/注册后可看大图

gs看上去没问题
fortran部分尽量不要用binary
open(11,file='d:\q\vwnd850.grd',form='unformatted',recl ...

非常非常感谢你热心的回答！
出现了新的问题

程序里定义的数组类型均是real型，我不太知道如何查看原来数据的类型~~

之前做u风场的时候，采用的方法跟这个一样，数据没有问题
可是这次v风场不知为何就是常数了？

lqouc

real的话应该没问题的啊，话说open里面可能要加上access='direct'
嗯.....你说你的导出的grd画不出图？用grads画不出来么？
那样的话可能是导出来的时候就有问题了。

wygohappy

lqouc 发表于 2013-4-27 13:26

                               
登录/注册后可看大图

real的话应该没问题的啊，话说open里面可能要加上access='direct'
嗯.....你说你的导出的grd画不出图？用g ...

加上了以后运行的错误变为：too many bytes read from unformatted record
画出的图就是一个常数场
首先读出的数据就都是常数了，后面的算法应该可以暂时不考虑吧。

我觉得很神奇，真是不知道哪里出错了。。。。

lqouc

wygohappy 发表于 2013-4-27 13:39

                               
登录/注册后可看大图

加上了以后运行的错误变为：too many bytes read from unformatted record
画出的图就是一个常数场
首先 ...

确实很神奇，那直接用nc数据化一下看看呢？应该不是原始数据的错误吧？

wygohappy

lqouc 发表于 2013-4-27 14:05

                               
登录/注册后可看大图

确实很神奇，那直接用nc数据化一下看看呢？应该不是原始数据的错误吧？

确定不是原始数据的错误，直接绘了几张图，还是有数值的~~~
虽然问题还是没有解决，但是非常感谢你的答复~~

lqouc

wygohappy 发表于 2013-4-27 14:11

                               
登录/注册后可看大图

确定不是原始数据的错误，直接绘了几张图，还是有数值的~~~
虽然问题还是没有解决，但是非常感谢你的答复 ...

唉，惭愧啊，搞了这么久也没弄出来，话说今天斑竹们走不在啊，其实可以召唤一下管理员求助。

wygohappy

lqouc 发表于 2013-4-27 14:05

                               
登录/注册后可看大图

确实很神奇，那直接用nc数据化一下看看呢？应该不是原始数据的错误吧？

啊啊啊，经过询问我班的大神，他给我找到了错误
是程序中的数组，应该是v(i,j,t)，都没变化，怪不得是常数，我太弱智了。。。

lqouc

我去.......