- 积分
- 318
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-3-21
- 最后登录
- 1970-1-1
![[KYV_] 粉丝数:198 微博数:926 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我想对pc时间序做每年的DJF合成,我的pc文件是8个模态,时间是40年*12月=480个,做合成的时候提醒我数组越界了,有没有大佬指导我修改一下?
以下是我的fortran程序:
program djf_pc_compst
implicit none
integer, parameter :: nt = 480, nm = 40, kvt=8 !
integer :: i, k !k是时间,i是模态
real, allocatable,dimension(:,:) :: F, pc_djf
open(11,file='D:\DATA\teleconnection\data\anomaly(tropical)\mveof\pc.dat',access='direct',recl=kvt)
do k=1, nt
read(11,rec=k) (F(i,k),i=1,kvt)
enddo
close(11)
do i=1, kvt
do k=1, nm-1
pc_djf(i,k)=(F(i,12*(k-1)+12)+F(i,12*(k-1)+13)+F(i,12*(k-1)+14))/3
enddo
enddo
open(13,file='D:\DATA\teleconnection\data\anomaly(tropical)\djf_pc_compst.dat',form='unformatted',access='direct',recl=kvt)
do i=1, kvt
write(11,rec=i) (pc_djf(i,k),k=1,nm-1)
enddo
close(13)
end
报的错是这样的:
|
|