fortran做eof分析

1、!此程序?yàn)?EOF 程序 ! 運(yùn)行時(shí)要改動(dòng)前面的空間、時(shí)間格點(diǎn)以及文件路徑，ks 和 kvt 根據(jù)自己的需要進(jìn)行改動(dòng)!程序中自動(dòng)去除缺省值并將其寫回生成數(shù)據(jù)(生成數(shù)據(jù)中缺省值為-9999.0)!對(duì)程序中 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 為時(shí)間長(zhǎng)度INTEGER,PARAMETER ::M=nt,KS=0,KVT=8 !kvt 為輸出的模態(tài)數(shù)! KS 的設(shè)置： ks>0 計(jì)算前先將數(shù)據(jù)標(biāo)準(zhǔn)化 ，! ks=0 時(shí)取距平，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 !時(shí)間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