- 积分
- 553
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-6-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 mlittlej 于 2017-4-21 09:45 编辑
麻烦帮我看看程序吧。数据都是可以输出到屏幕上面的但是不能写入文件,文件大小一直是0。谢谢谢谢谢谢!!!
program realtime
parameter (k=5)
parameter (MM=6)
parameter (NN=500)
dimension y(NN),x(NN*k),a(MM),b(mm,mm),v(k),result(NN*k),data(NN*k),typhoon(NN*k)
integer ww,n
character cww*3
integer status
open(121,file='/intensity/REG/coefficient01')
write(121,*)'1231111111111111111111111111111111111'
open(131,file='/intensity/REG/reg01')
do ww=111,113,1
write(cww,'(i3.3)')ww
open(ww,file='/intensity/result/intensity01_'//cww)
i=1
do while(.true.)
read(ww,'(i10,1x,i4,15x,f6.1,1x,f6.1,1x,f6.1,1x,f6.1,1x,f6.1,9x,f6.1)',iostat=status)data(i),typhoon(i),x(5*(i-1)+1),x(5*(i-1)+2),x(5*(i-1)+3),x(5*(i-1)+4),x(5*(i-1)+5),y(i)
write(*,'(i10,1x,i4,15x,f6.1,1x,f6.1,1x,f6.1,1x,f6.1,1x,f6.1,1x,f6.1)')data(i),typhoon(i),x(5*(i-1)+1),x(5*(i-1)+2),x(5*(i-1)+3),x(5*(i-1)+4),x(5*(i-1)+5),y(i)
if(status.ne.0)exit
i=i+1
enddo
N=i-1
call isqt2(x,y,k,mm,n,a,q,s,r,v,u,b,dyy)
do j=1,mm
write(121,'(f8.2)') a(j)
write(*,'(i3,1x,f8.2)') n,a(j)
enddo
do i=1,N
result(i)=a(1)+x(5*(i-1)+1)*a(2)+x(5*(i-1)+2)*a(3)+x(5*(i-1)+3)*a(4)+x(5*(i-1)+4)*a(5)+x(5*(i-1)+5)*a(6)
print*,'111111111111',result(i)
write(131,'(i10,1x,i4,1x,f6.1)')date(i),typhoon(i),result(i)
write(*,'(i10,1x,i4,1x,f6.1)')date(i),typhoon(i),result(i)
enddo
enddo
close(121)
close(131)
close(111)
close(112)
close(113)
end |
|