- 积分
- 275
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-10-31
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
运行fortran程序出现这个结果,是怎么回事?
cccccccccc prepare data for svd analysis ccccccccc
cccccccccc the program is to normalize sea surface temperature(SST) cccccccccc
cccccccccc mt: the length of time series;
cccccccccc mo: the month numbers; my: the year numbers;
cccccccccc sst: sea surface temperature data;
cccccccccc sst3: the work array; avf: the average of SST;
cccccccccc df: the variance of SST;
program main
parameter(mo=12,my=43,nx=18,ny=12,mt=516)
dimension sst2(nx,ny,mo,my),avf(mo,nx,ny),df(mo,nx,ny)
dimension sst(nx,ny,mt),sst3(nx,ny,mt),II(NX,NY),
& jp(nx,ny,mo,my)
open(1,file='g:sstpx.grd',form='unformatted',
& access='direct',recl=nx*ny)
irec=1
do 5 it=1,mt
read(1,rec=irec)((sst(i,j,it),i=1,nx),j=1,ny)
irec=irec+1
5 continue
do i=1,nx
it1=i
do j=1,ny
do it=1,mt
do m=1,mo
do myear=1,my
sst2(i,j,m,myear)=sst(i,j,it1)
it1=it1+12
end do
end do
end do
end do
end do
do i=1,nx
do j=1,ny
do m=1,12
avf(i,j,m)=0.0
do myear=1,my
avf(i,j,m)=avf(i,j,m)+sst2(i,j,m,myear)/my !求平均值
end do
end do
end do
end do
do i=1,nx
do j=1,ny
do m=1,mo
do myear=1,my
jp(i,j,m,myear)=sst2(i,j,m,myear)-avf(i,j,m) !求距平
end do
end do
end do
end do
open(2,file='g:\jp.grd',form='binary')
do i=1,nx
do j=1,ny
do myear=1,my
do m=1,mo
kk=(myear-1)*12+m
sst3(i,j,kk)=jp(i,j,m,myear)
end do
end do
end do
end do
do it=1,mt
write(2)((sst3(i,j,it),i=1,nx),j=1,ny)
end do
ccccccccccccccccccccccccccccccccccccccccccc
close (1)
close (2)
stop
end
这是我编的一段程序,需要将这段程序和另一段程序结合,构成eof分析的程序,没错,但是在电脑运行不出来,想请大神帮忙看看
|
|