fortran做eof分析

1、!此程序為 EOF 程序 ! 運行時要改動前面的空間、時間格點以及文件路徑，ks 和 kvt 根據(jù)自己的需要進(jìn)行改動!程序中自動去除缺省值并將其寫回生成數(shù)據(jù)(生成數(shù)據(jù)中缺省值為-9999.0)!對程序中 data_in 到 F 的傳遞 進(jìn)行調(diào)整后 此程序也可用于 s-eof 和 mv-eof!PROGRAM EOFIMPLICIT NONEINTEGER,PARAMETER :: nt=12,nx=23,ny=34

2、! you need change,NT 為時間長度INTEGER,PARAMETER ::M=nt,KS=0,KVT=8 !kvt 為輸出的模態(tài)數(shù)! KS 的設(shè)置： ks>0 計算前先將數(shù)據(jù)標(biāo)準(zhǔn)化 ，! ks=0 時取距平，ks=N) THEN DO K=1,KVTAVE=SUM(F(K,1:M))/REAL(M)F1(K,1:M)=F(K,1:M)-AVED= SQRT(SUM(F1(

3、K,1:M)*F1(K,1:M))/REAL(M))F(K,1:M)= F(K,1:M)/D !時間S(K,1:N)= S(K,1:N)*D !空間ENDDOm1=0DO i=1,nxDO j=1,nyIF(land(i,j).eq.0.0)thenm1=m1+1PT(i,j,1:kvt)=S(1:kvt,m1)ELSEPT(i,j,1:kvt)=-9999.00ENDIFENDDOENDDOOPEN (1,

4、file=TRIM(track)//'pt.dat',access='direct',recl=NX*NY)DO k=1,KVTWRITE(1,rec=k)((PT(i,j,k),i=1,nx),j=1,ny)ENDDOCLOSE(1)OPEN (2,FILE=TRIM(track)//'PC.DAT',ACCESS='DIRECT',RECL=M)DO K=1,KVTWR