请问fortran 怎么写入,读取大数组(可以二进制吗?)数据过大,发现读取不了。谢谢

program array
 implicit none

 integer(8)::i,j,m,n,mm,nn
 integer(8),parameter::zz=777
 real(8):: dt
 real(8):: RR(800,76981),IR(800,76981),RdR(800,76981),IdR(800,76981),f(12,76981),ff(2,800)
 real(8)::t(76981),a(76981),H(76981),k(800),pk(800),c1(76981),c2(76981),c3(76981),c4(76981),c5(76981),&
 c6(76981),c7(76981),c8(76981),aa1(76981),qs(76981),AA2(76981),AA3(76981),AA4(76981),AA5(76981),AA6(76981)&
,AA7(76981),AA8(76981)
 integer(8)::tc(800),te(800),ts(800)
 complex(16)::sum1(800),sum2(800),sum3(800),sum4(800),sum5(800),sum6(800),sum71(800),sum72(800),&
sum81(800),sum82(800)
 complex(16)::fac1(800),fac2(800),fac3(800),fac4(800),fac5(800),fac6(800),fac71(800),fac81(800),&
fac72(800),fac82(800),sum00(800),sum(800)
 
 complex(16),dimension(:,:),allocatable::R,DR,RRR
 allocate( R(800,76981),DR(800,76981),RRR(800,76981))
 
  Open(0 , File ='/home/cly/桌面/look(复件)/result1-4.log')
  Open(11 , File ='/home/cly/桌面/look(复件)/Re_Rll.log' )
  Open(22 , File ='/home/cly/桌面/look(复件)/Im_Rll.log' )
  Open(33 , File ='/home/cly/桌面/look(复件)/Re_dRll.log' )
  Open(44 , File ='/home/cly/桌面/look(复件)/Im_dRll.log' )
  Open(55 , File ='/home/cly/桌面/look(复件)/ts.log' )
  Open(66 , File ='/home/cly/桌面/look(复件)/tc.log' )
  Open(77 , File ='/home/cly/桌面/look(复件)/te.log' )
  Open(99, File ='/home/cly/桌面/look(复件)/spectra1-4.log' )
! Open(100, File ='/home/jm/桌面/test/equi/equi_1.log' )

do i=1,76981
  read(0,*)  f(:,i)
  read(11,*) RR(:,i)
  read(22,*) IR(:,i)
  read(33,*) Rdr(:,i)
  read(44,*) Idr(:,i)
 end do
do I=1,800
  read(99,*) (ff(j,i),j=1,2)
end do
do i=1,800
 read(55,*) ts(i)
end do
do i=1,800
 read(66,*) tc(i)
end do
do i=1,777
 read (77,*) te(i)
end do

write(*,*)  'reading is ok'
 

end program

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值