逐分钟资料用fortran求日平均

q863510557

lqouc 发表于 2013-11-2 15:24

                               
登录/注册后可看大图

write（2，*）
有格式的要声明输出格式，即使是默认格式。

在前辈的指导下已经搞定了  谢谢前辈！！！

q863510557

andrewsoong 发表于 2013-11-2 14:39

                               
登录/注册后可看大图

版主好人做到底~~~给改改程序~~~

嘿嘿  已经搞定了  thks

andrewsoong

q863510557 发表于 2013-11-2 19:48

                               
登录/注册后可看大图

嘿嘿  已经搞定了  thks

搞定了，那就把程序贴一下，让后辈们也学习一下~~~

q863510557

andrewsoong 发表于 2013-11-2 19:52

                               
登录/注册后可看大图

搞定了，那就把程序贴一下，让后辈们也学习一下~~~

ok
程序如下：
program main
implicit none
integer n,nn
integer j,jj,i,d
integer k,sum
parameter(n=1578240,nn=1096)
real  cs1(n),cs2(n),cs3(n),cs4(n)
real  dat(nn),ave(nn)
character cha(n),hh
open(1,file='d:\gbno2.txt',status='old',form='formatted')
do j=1,n
read(1,"(I7,A10,f8.5,A8,f8.5,A8,f8.5)") cs1(j),cha(j),cs2(j),hh,cs3(j),hh,cs4(j)
end do
close(1)
i=1
d=0
sum=1440
k=0
do j=1,n
d=d+cs2(j)
k=k+1
if(cs2(j)==0.0) then
sum=sum-1
end if
if(k==1440.and.sum==0) then
ave(i)=0
dat(i)=cs1(j)
i=i+1

d=0
k=0
sum=1440
else if(k==1440) then
ave(i)=d/sum
dat(i)=cs1(j)

i=i+1

d=0
k=0
sum=1440
endif
end do
open(2,file='d:\ave.txt',status='new')
do jj=1,nn
write(2,'(I10,A5,f9.5)') dat(jj),hh,ave(jj)

end do
close(2)
end
求批评指正！

flowers940815

谢楼主分享，共同学习~~