EOF分解程序,附测试数据,图

言深深

不相信 发表于 2013-7-4 21:21

                               
登录/注册后可看大图

谢谢分享，正好我需要的程序，我想问一下，如果我要用站点降水量数据的话，在下面的surfer画图准备数据的程 ...

烦请认真些哈，多看一会应该能看懂的。你解释的貌似都有问题···汗···

koperz

很厉害啊，学习

又是那隻貓

谢谢分享！！

xiaohundanpipi

深深，我用这个程序怎么一直也运行不了，说是array bounds exceeded我是二进制的数据，用fortran转成txt格式，程序如下
program iiimplicit none
integer::i,j,t,irec,k,q,p
real,dimension(10512,65)::aa
real,dimension(144,73,65)::bb
irec=0
open(1,file='e:/eof/summer/ave.dat',status='old',access='direct',form='unformatted',recl=73*144)
do t=1,65
irec=irec+1
read(1,rec=irec)((bb(i,j,t),i=1,144),j=1,73)
end do
close(1)
print*,bb(1,:,1)
print*, 'a'

j=0
k=0
do t=1,65
do i=1,10512,144
  do p=0,143
   q=i+p
   j=int((i-1)/144)+1
   k=q-(j-1)*144
   aa(q,t)=bb(k,j,t)
  end do
end do
end do
print*,aa(10507:10512,65)

open(1,file='ave.txt',status='replace')
do i=1,10512
write(1,'(65f30.8)')(aa(i,t),t=1,65)
end do
close(1)

生成的txt怎么也运行不了，好奇怪，你能帮我看看么？
我还想问你你的测试数据，是71*144*62的，那每一列是144个71存的，还是71个144存的？

xiaohundanpipi

xiaohundanpipi 发表于 2013-7-15 21:35

                               
登录/注册后可看大图

深深，我用这个程序怎么一直也运行不了，说是array bounds exceeded我是二进制的数据，用fortran转成txt格式 ...

找到问题了，这个程序里 mnl代表m,n中小的那个数，求出的矩阵交叉积的维数是(mnl,mnl)。如果给mnl赋了别的值，就会报错。

王跃峰

楼主，帮忙看看这是咋回事。

宋斌

大赞一个

WEN

厉害！

WEN

。。。。我还在想怎么去挣两个贡献下文件呢。。。。。

言深深

WEN 发表于 2013-7-30 22:01

                               
登录/注册后可看大图

。。。。我还在想怎么去挣两个贡献下文件呢。。。。。

额，这悲剧的···写点儿帖子给你加分去啊，嘻嘻···