0
点赞
收藏
分享

微信扫一扫

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

两岁时就很帅 2022-01-10 阅读 59

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 条评论