- 积分
- 1971
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-4-7
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 咩咩的羊 于 2018-9-27 20:18 编辑
fortran程序运行不出错,调试出错。主要程序段是读取三个grd文件处理之后形成新的grd文件。求助!
program main
integer nx,ny,i,j,t
parameter(nx=480,ny=78,nt=31,nt2=28)
real t2m1(nx,ny,nt),t2m2(nx,ny,nt),t2m3(nx,ny,nt),num(nx,ny),num1(nx,ny),num2(nx,ny),num3(nx,ny)
num=0
open(2,file='F:\academy\team\200701.grd',form='binary')
read(2)(((t2m1(i,j,t),i=1,nx),j=1,ny),t=1,nt)
do j=1,ny
num1=0
do i=1,nx
do t=1,nt
if(t2m1(i,j,t)>=36)then
num1(i,j)=num1(i,j)+1
endif
enddo
end do
enddo
open(4,file='F:\academy\team\200702.grd',form='binary')
read(4)(((t2m2(i,j,t),i=1,nx),j=1,ny),t=1,nt2)
do j=1,ny
num2=0
do i=1,nx
do t=1,nt2
if(t2m2(i,j,t)>=36)then
num2(i,j)=num2(i,j)+1
endif
enddo
end do
enddo
open(6,file='F:\academy\team\200712.grd',form='binary')
read(6)(((t2m3(i,j,t),i=1,nx),j=1,ny),t=1,nt)
do j=1,ny
num3=0
do i=1,nx
do t=1,nt
if(t2m3(i,j,t)>=36)then
num3(i,j)=num3(i,j)+1
endif
enddo
end do
enddo
open(3,file='F:\academy\team\2007d.grd',form='binary')
num(i,j)=num1(i,j)+num2(i,j)+num3(i,j)
write(3)((num(i,j),i=1,nx),j=1,ny)
close(2)
close(3)
close(4)
close(6)
end
|
-
|