| 
 
	积分553贡献 精华在线时间 小时注册时间2012-6-12最后登录1970-1-1 
 | 
 
| 
本帖最后由 mlittlej 于 2017-4-21 09:45 编辑
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  
 麻烦帮我看看程序吧。数据都是可以输出到屏幕上面的但是不能写入文件,文件大小一直是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
 | 
 |