- 积分
- 9004
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-3-17
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-10-30 14:15:47
|
显示全部楼层
本帖最后由 绯宜物语 于 2015-10-30 15:30 编辑
program cmip5
integer*4,parameter:: x=360,y=181,t=105,m=17,n=10
integer*4 :: ix,iy,im,in,it,irec
real*4,parameter :: undef = -9.99E+8
real*4 pra(x,y,t),pr(m,n,t)
open(1,file='f:\pr\h\cru\prcsu.grd',form='binary')
do it=1,t
do iy=1,y
do ix=1,x
read(1) pra(ix,iy,it)
enddo
enddo
enddo
close(1)
do it=1,t
do im=1,m
do in=1,n
ix=4*im+69
iy=4*in+15
pr(im,in,it)=pra(ix,iy,it)+pra(ix,iy+1,it)+pra(ix,iy+2,it)+pra(ix,iy+3,it)+pra(ix+1,iy,it)+pra(ix+1,iy+1,it)+pra(ix+1,iy+2,it)+pra(ix+1,iy+3,it)+pra(ix+2,iy,it)+pra(ix+2,iy+1,it)+pra(ix+2,iy+2,it)+pra(ix+2,iy+3,it)+pra(ix+3,iy,it)+pra(ix+3,iy+1,it)+pra(ix+3,iy+2,it)+pra(ix+3,iy+3,it)
enddo
enddo
enddo
open(3,file='f:\prmm.grd',form='unformatted',access='direct',recl=m*n*4)
irec = 0
do it = 1,t
irec = irec + 1
write(3,rec=irec) ((pr(im,in,it),im=1,m),in=1,n)
enddo
close(3)
end
写了Fortran程序做,但是输出的值基本都是0.
有人可以帮我看看哪有问题吗?谢谢。 |
|