- 积分
- 648
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-9-18
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
parameter(ii=37,jj=17,zz=12,nt=48)
real var(ii,jj,nt)
real:: ave(ii,jj,12)=0
real ::have(ii,jj,nt)=0
real ::sav(ii,jj,nt) =0
open(5,file='d:\h500.dat')
DO 200 it=1,48
read(5,1000)
read(5,3000) ((var(i,j,it),i=1,ii),j=1,jj)
200 continue
open(6,file='d:\trans.grd',form='unformatted',access='direct',recl=ii*jj*4)
irec=0
do it=1,nt
irec=irec+1
write(6,rec=irec) ((var(i,j,it),i=1,ii),j=1,jj)
enddo
1000 format(2i7)
2000 FORMAT(37F6.2)
3000 format(37f8.1)
4000 format(37f7.2)
!求平均
do im=1,12
do it=im,48,12
do j=1,jj
do i=1,ii
ave(i,j,it)=ave(i,j,it)+var(i,j,it)/4
enddo
enddo
enddo
enddo
open(6,file='d:\ave.grd',form='unformatted',access='direct',recl=ii*jj*4)
irec=0
do im=1,12
irec=irec+1
write(6,rec=irec) ((pj(i,j,i),i=1,ii),j=1,jj)
enddo
!求距平,方差
do it=1,12
do j=1,jj
do i=1,ii
have(i,j,it)= have(i,j,it)+var(i,j,it)-ave(i,j,it)
sav(i,j,it)=sav(i,j,it)+(have(i,j,it)**2)/4
enddo
enddo
enddo
open(7,file='d:\have.grd',form='unformatted',access='direct',recl=ii*jj*4)
open(8,file='d:\sav.grd',form='unformatted',access='direct',recl=ii*jj*4)
irec=0
do it=1,nt
irec=irec+1
write(7,rec=irec) ((ave(i,j,i),i=1,ii),j=1,jj)
enddo
irec=0
do it=1,nt
irec=irec+1
write(8,rec=irec) ((sav(i,j,i),i=1,ii),j=1,jj)
enddo
end
下面是错误提示,调试无错,执行时有错Text1.obj : error LNK2005: _main already defined in Text3.obj
Text1.obj : error LNK2001: unresolved external symbol _PJ@12
Text1.exe : fatal error LNK1120: 1 unresolved externals
Error executing link.exe.
Text1.exe - 3 error(s), 0 warning(s)
|
|