2024年4月3日发(作者：)

!此程序为EOF程序

! 运行时要改动前面的空间、时间格点以及文件路径，ks和kvt根据自己的需要进行

改动

!程序中自动去除缺省值并将其写回生成数据(生成数据中缺省值为-9999.0)

!对程序中data_in到F的传递 进行调整后 此程序也可用于s-eof和mv-eof

!

PROGRAM EOF

IMPLICIT NONE

INTEGER,PARAMETER :: nt=12,nx=23,ny=34 ! you need change,NT

为时间长度

INTEGER,PARAMETER ::M=nt,KS=0,KVT=8 !kvt为输出的模态数

! KS的设置： ks>0 计算前

先将数据标准化 ，

! ks=0时取距平，ks<0时不进行这一步处理

INTEGER :: i,j ,MNH,N ,K,IM , m1

REAL, allocatable,dimension(:,:,:)::DATA_IN

REAL, allocatable,dimension(:,:)::F,S,ER,A,S1,F1

CHARACTER(LEN=20) :: NOW , TRACK

REAL :: land(nx,ny), D,AVE,PT(NX,NY,kvt) ,ran1

TRACK='E:aatEA' !输出的目标文

件夹，默认为程序所在文件夹

call time(now)

print*, now

!!1111111111读入数据并去掉缺省值111

ALLOCATE(DATA_IN(NX,NY,NT))

OPEN(1,file='E:',access='direct',recl=nx*ny*nt) !****修改

路径

READ(1,rec=1) (((data_in(I,J,K),I=1,nx),J=1,ny),K=1,nt)