四维场读取和写入问题

nickbsb

program main
implicit none
integer i,j,k,t,it,nt
real uwnd(144,29,17,121)
real x(121),y(144,29,17,121)
open(10,file='E:\uwndfield\uwnd.grd',form='binary')
do k=1,17
  do j=1,29
   do i=1,144
   do t=1,121
read(10) uwnd(i,j,k,t)
x(t)= uwnd(i,j,k,t)
   
end do
end do
end do
end do
do k=1,17
  do j=1,29
   do i=1,144
   do it=1,121
y(i,j,k,it)= x(it)
write(*,*) y(i,j,k,it)
   
end do
end do
end do
end do
end


我先读取uwnd这个四维场写入一个一维时间序列里， 再将这个时间序列转回到一个四维场，进行读取发现fortran运行的黑框出现的都是0.0000000E+00这个程序是不是格式有什么问题，望指教
如果将121改成更大的数3630，fortran运行的黑框干脆就是瞬间闪现一下就消失了，这是肿么回事啊

做个霸气的木头

grd数据的顺序不是你这样的！用循环的话，从外到内依次是：时间、层次、y、x：
do t=1,it
   do z=1,iz
　　do y=1,iy
              read(100)（ss(x,y,z,t),x=1,ix）
enddo;enddo;enddo;
手机打字不容易，你看着对应改吧！

言深深

记得前几次貌似和楼主说过一次，另外一个你发的帖子里面。http://bbs.06climate.com/forum.php?mod=viewthread&tid=8114请阅读第一点，关于二进制数据的读取问题，第二点是纠正你程序错误。第三点是你的思路可能并正确，对滤波的了解可能不对。

再给你贴到这边来
*******************************************

你就是对每个格点逐一滤波吧？这样是可以的，但是程序上面有一些漏洞！
1、二进制的数据在读取方式上面是有要求的，类似read(*,*)((((x(i,j,k),i=1,m),j=1,n),k=1,o)
2、在call调用子过程得到结果的时候你的输出是不妥的，你这样做的输出会导致文件中只有一个相同的结果，就是最后一个个点的滤波结果，正确的做法是在执行循环的外端打开文件，而不是每执行一次循环打开一次
3、输出的时候是不妥的，至少保证输出的时候是一条序列或者多条序列，应该在输出语句里面执行一次或多次循环。

【建议】楼主参考一下论坛里发过的不少二进制文件的输入输出，并弄清楚滤波的输入输出文件以及设计一下输出方案。

*******************************************

topmad

言深深 发表于 2012-6-23 12:16

                               
登录/注册后可看大图

记得前几次貌似和楼主说过一次，另外一个你发的帖子里面。http://bbs.06climate.com/forum.php?mod=viewthr ...

哇........
深深的回答好专业 看不懂了

言深深

topmad 发表于 2012-6-23 20:51

                               
登录/注册后可看大图

哇........
深深的回答好专业 看不懂了

么有表述好。；楼主前几次发过一个帖子啦，之前我回复了一下，这次又贴出来了，应该楼主能看明白的吧。就是fortran语句的一些应用嘛