怎样对一个有四五万条记录的txt文件进行统计运算？？？困扰很久了。求助呀%>_<%

锦城の烟雨

《四川地面气候日值数据集》1951-2011年四川省各测站各气象要素的逐日观测数据，现在想求56492（宜宾）站1951-2011年逐年11月上旬日平均气温的旬平均值。
每个txt文件里有四五万行，不知道怎样求旬平均值。
还是将每个txt文件里的数据存到一个数组里再对数组进行运算吗？那这个数组会不会太大了？？？
这个问题困扰我好久了。。。不解决论文就没法动笔。。下周就要中期检查了。。。。求助各位大神。。。。

lqouc

同志，才四五万行不要说得那么悲壮吧，循环读取加两层判断就好了，数组怎么会太大？我算全球的六小时间隔站点也没溢出啊，楼主先自己动手吧。不会有人直接编好给你的。

冰点

我刚好做了类似的工作，第一步我计算了每个站每个月的月平均温度
我的数据是60年的全国756站的降水，最大、最小及平均日温度数据，程序中我只用到了温度，程序运行后可以得到60年的每个站的月平均温度数据。 years.txt中存放的是60年的文件名，作批处理用，temp文件夹中放有60个文件，最后一张图是我随意打开的一年的数据，数据格式如图，也有缺测值，平均下来每年的数据也是几万行，这儿没用数组，希望你能看明白。现把程序贴上来： monthly.f90 (1.09 KB, 下载次数: 30)

2013-4-20 12:05 上传

点击文件名下载附件

冰点

另外，我看你的数据文件名也都是连续的，作批处理的时候也很方便……

topmad

4w+不算大，对fortran来说 小case
再说你都有时间了，判断一下 循环一个  都出来了

做个霸气的木头

你要相信fortran的数组能力，楼上说的都对~

river

补充一点，楼上的说的都对
从楼主说的那个数组会不会太大来看，楼主自己还没有具体思路，而且也不理解fortran的一些东西，所以还是自己先努力一下吧

锦城の烟雨

river 发表于 2013-4-20 15:38

                               
登录/注册后可看大图

补充一点，楼上的说的都对
从楼主说的那个数组会不会太大来看，楼主自己还没有具体思路，而且也不 ...

我的想法就是把文件里的数据读取到数组中，再对数组进行运算。下面是我的程序。编译没有错误。但是运行不出来，能不能帮我看看。。

program main

character filei,fileii
integer stn(100000),year(100000),mon(100000),day(100000),temp(100000),n
real sum

filei='SURF_CLI_CHN_MUL_DAY_1348578049292_01.txt'
fileii='xpj'

sum=0

do i=1,100000
stn(k)=0
year(k)=0
mon(k)=0
day(k)=0
temp(k)=0
enddo

print*,stn

do j=1951,2011
!write(fileii(1:4),'(i4)') j
open(2,file=fileii)

do i=1,10
write(filei(36:37),'(i2)') i
open(1,file=filei)

do k=1,100000
read(1,100) stn(k),year(k),mon(k),day(k),temp(k)
100 format(I5,I7,I5,I7,I12)
close(1)
enddo

do k=1,100000
if(stn(k).eq.56492.and.mon(k).eq.11.and.day(k)<=10.and.temp(k).ne.32766.and.year(k).ne.j)then
sum=sum+temp(k)
n=n+1
endif
sum=sum/(n*10)
write(2,*)         j,sum
enddo
enddo

close(2)
enddo

end
      

锦城の烟雨

topmad 发表于 2013-4-20 12:53

                               
登录/注册后可看大图

4w+不算大，对fortran来说 小case
再说你都有时间了，判断一下 循环一个  都出来了

我的想法就是把文件里的数据读取到数组中，再对数组进行运算。下面是我的程序。编译没有错误。但是运行不出来，能不能帮我看看。。

program main

character filei,fileii
integer stn(100000),year(100000),mon(100000),day(100000),temp(100000),n
real sum

filei='SURF_CLI_CHN_MUL_DAY_1348578049292_01.txt'
fileii='xpj'

sum=0

do i=1,100000
stn(k)=0
year(k)=0
mon(k)=0
day(k)=0
temp(k)=0
enddo

print*,stn

do j=1951,2011
!write(fileii(1:4),'(i4)') j
open(2,file=fileii)

do i=1,10
write(filei(36:37),'(i2)') i
open(1,file=filei)

do k=1,100000
read(1,100) stn(k),year(k),mon(k),day(k),temp(k)
100 format(I5,I7,I5,I7,I12)
close(1)
enddo

do k=1,100000
if(stn(k).eq.56492.and.mon(k).eq.11.and.day(k)<=10.and.temp(k).ne.32766.and.year(k).ne.j)then
sum=sum+temp(k)
n=n+1
endif
sum=sum/(n*10)
write(2,*)         j,sum
enddo
enddo

close(2)
enddo

end

锦城の烟雨

lqouc 发表于 2013-4-20 11:20

                               
登录/注册后可看大图

同志，才四五万行不要说得那么悲壮吧，循环读取加两层判断就好了，数组怎么会太大？我算全球的六小时间隔站 ...

我的想法就是把文件里的数据读取到数组中，再对数组进行运算。下面是我的程序。编译没有错误。但是运行不出来，能不能帮我看看。。

program main

character filei,fileii
integer stn(100000),year(100000),mon(100000),day(100000),temp(100000),n
real sum

filei='SURF_CLI_CHN_MUL_DAY_1348578049292_01.txt'
fileii='xpj'

sum=0

do i=1,100000
stn(k)=0
year(k)=0
mon(k)=0
day(k)=0
temp(k)=0
enddo

print*,stn

do j=1951,2011
!write(fileii(1:4),'(i4)') j
open(2,file=fileii)

do i=1,10
write(filei(36:37),'(i2)') i
open(1,file=filei)

do k=1,100000
read(1,100) stn(k),year(k),mon(k),day(k),temp(k)
100 format(I5,I7,I5,I7,I12)
close(1)
enddo

do k=1,100000
if(stn(k).eq.56492.and.mon(k).eq.11.and.day(k)<=10.and.temp(k).ne.32766.and.year(k).ne.j)then
sum=sum+temp(k)
n=n+1
endif
sum=sum/(n*10)
write(2,*)         j,sum
enddo
enddo

close(2)
enddo

end