- 积分
- 71
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-23
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2013-7-27 19:44:06
|
显示全部楼层
生成数据的程序
program main
implicit none
real,parameter::nn=8,nx=144,ny=45
real::a(nx,ny),b(nx,ny),c(nx,ny),sum1(nx,ny),ave(nx,ny)
integer::n,i,j
character*2::nam1(nn)=(/'1','2','3','4','5','6','7','8'/)
do n=1,nn
open(1,file='d:\MJO\v200'//trim(nam1(n))//'-3.grd',form='binary')
open(2,file='d:\MJO\v200'//trim(nam1(n))//'-4.grd',form='binary')
open(3,file='d:\MJO\v200'//trim(nam1(n))//'-5.grd',form='binary')
read(1) ((a(i,j),i=1,nx),j=1,ny)
close(1)
read(2) ((b(i,j),i=1,nx),j=1,ny)
close(2)
read(3) ((c(i,j),i=1,nx),j=1,ny)
close(3)
open(10,file='d:\MJO\v200'//trim(nam1(n))//'-345.grd',form='binary')
open(11,file='d:\MJO\v200'//trim(nam1(n))//'-345.txt')
do j=1,ny
do i=1,nx
sum1(i,j)=a(i,j)+b(i,j)+c(i,j)
ave(i,j)=sum1(i,j)/3
write(10) ave(i,j)
write(11,*) ave(i,j)
enddo
enddo
close(10)
enddo
end
|
|