fortran算月平均值

q863510557

用fortran算月平均值，结果出现有一行没有值（就是说等于0），其他行都有值，不知道那一步错了，看了半天也没看出来，望哪位大神指教一下下，如下图为输出文件的截图（第一列为月份，080100代表08年一月；第二列为所求的平均值，最后一行的值为什么是0）
程序如下（有点长，不太会编）：
program main
implicit none
integer n,nn
parameter(n=1578240,nn=36)
integer j,jj,f
integer cs1(n),su(nn),sum,dat(nn),s
real ave(nn)
real cs2(n)
real cs3(n)
real cs4(n)
real d
character hh,cha(n)
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)
j=1
do jj=1,nn
if(cs1(j)>80200.and.cs1(j)<80300) then
su(jj)=29
dat(jj)=cs1(j)-1
j=j+29*1440
else  if(cs1(j)>90200.and.cs1(j)<90300) then
su(jj)=28
dat(jj)=cs1(j)-1
j=j+28*1440
else if(cs1(j)>100200.and.cs1(j)<100300) then
su(jj)=28
dat(jj)=cs1(j)-1
j=j+28*1440
else if(cs1(j)>80400.and.cs1(j)<80500) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>80600.and.cs1(j)<80700) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>80900.and.cs1(j)<81000) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>81100.and.cs1(j)<81200) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>90400.and.cs1(j)<90500) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>90600.and.cs1(j)<90700) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>90900.and.cs1(j)<91000) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>91100.and.cs1(j)<91200) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>100400.and.cs1(j)<100500) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>100600.and.cs1(j)<100700) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>100900.and.cs1(j)<101000) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)>91100.and.cs1(j)<91200) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else
su(jj)=31
dat(jj)=cs1(j)-1
j=j+31*1440
end if
end do
d=0
s=su(1)*1440
jj=1
f=0
do j=1,n
d=d+cs2(j)
if(cs2(j)==0) then
f=f-1
end if
if(j==s) then
sum=su(jj)*1440
ave(jj)=d/(sum+f)
jj=jj+1
s=s+su(jj)*1440
d=0
f=0
end if
end do
open(2,file='d:\avemonth.txt',form='formatted',status='new')
do jj=1,36
write(2,'(I10,A5,f20.15)') dat(jj),hh,ave(jj)
end do
close(2)
end

腐朽的我Zero

andrewsoong

这个我好像遇到过~~~需要在你读取的文件最后一行加上一行，加上这一行任意~~~只要是在最后一行就行~~~你试试~~~

mofangbao

请自己先逐步调试  一点点的输出各个变量

oneworker

是否读取数据有问题，建议把读取的数据存储在文本文件中，对照原始数据，检验一下。

q863510557

腐朽的我Zero 发表于 2013-11-4 08:36

                               
登录/注册后可看大图

                                

q863510557

andrewsoong 发表于 2013-11-4 08:47

                               
登录/注册后可看大图

这个我好像遇到过~~~需要在你读取的文件最后一行加上一行，加上这一行任意~~~只要是在最后一行就行~~~你试试 ...

我发现我有一行写错了，就是中间求dat（jj）和su（jj）那一段，我复制粘贴之前的内容，有个东西搞忘改了，就导致最后j多了一天的量，所以最后一个出来是0，谢谢了

q863510557

mofangbao 发表于 2013-11-4 08:49

                               
登录/注册后可看大图

请自己先逐步调试  一点点的输出各个变量

问题解决了  谢谢！！

q863510557

oneworker 发表于 2013-11-4 09:37

                               
登录/注册后可看大图

是否读取数据有问题，建议把读取的数据存储在文本文件中，对照原始数据，检验一下。

问题解决了，谢谢了
改过以后的程序：
program main
implicit none
integer n,nn
parameter(n=1578240,nn=36)
integer j,jj,f
integer cs1(n),su(nn),sum,dat(nn),s
real ave(nn)
real cs2(n)
real cs3(n)
real cs4(n)
real d
character hh,cha(n)
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)
j=1
do jj=1,nn
if(cs1(j)==80201) then
su(jj)=29
dat(jj)=cs1(j)-1
j=j+29*1440
else  if(cs1(j)==90201) then
su(jj)=28
dat(jj)=cs1(j)-1
j=j+28*1440
else if(cs1(j)==100201) then
su(jj)=28
dat(jj)=cs1(j)-1
j=j+28*1440
else if(cs1(j)==80401) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==80601) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==80901) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==81101) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==90401) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==90601) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==90901) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==91101) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==100401) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==100601) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==100901) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else if(cs1(j)==101101) then
su(jj)=30
dat(jj)=cs1(j)-1
j=j+30*1440
else
su(jj)=31
dat(jj)=cs1(j)-1
j=j+31*1440
end if
end do
d=0
s=su(1)*1440
jj=1
f=0
do j=1,n
d=d+cs2(j)
if(cs2(j)==0) then
f=f-1
end if
if(j==s) then
sum=su(jj)*1440
ave(jj)=d/(sum+f)
jj=jj+1
s=s+su(jj)*1440
d=0
f=0
end if
end do
open(2,file='d:\avemonth.txt',form='formatted',status='new')
do jj=1,36
write(2,'(I10,A5,f20.15)') dat(jj),hh,ave(jj)
end do
close(2)
end

andrewsoong

q863510557 发表于 2013-11-4 10:28

                               
登录/注册后可看大图

问题解决了，谢谢了
改过以后的程序：
program main

恭喜~~~