站点35年求平均求助

你看那河下的桥

现在有35个站35年逐日的最低温度资料（已经进行过筛选，在3333文档中存放的格式是非低温日tmin(i,j,k)=0，低温日有数据）现在想对35个站分别求35年平均，但是不求最后的值，只是对低温日进行统计再除以35年，例如第一个站对35年进行逐日统计，统计出来的低温日的天数是300，那么第一个站35年平均为：300/35，这就是我最后要输出的值
但是怎么都算不对，感觉是统计那里出了问题，求大神帮忙解决


program ex02_02
   implicit none
   !定义变量名
   integer,parameter::sta_num=35,ns=35,day=365    ！站数  年份  天数
   integer   yr(sta_num,ns),mon(sta_num,ns),da(sta_num,ns),sta(sta_num,ns)
  integer i,j,k,a
real tmin(sta_num,ns,day),ave2(sta_num),tmin1(sta_num)
open(1111,file='d:\lunwen\zhandian\diwenri\35yearpingjun.txt')

open(3333,file='d:\lunwen\zhandian\diwenri\huadongshaixuanwendu.txt')
do i=1,35
do j=1,ns
do k=1,365
read(3333,*) sta(i,j),yr(i,j),mon(i,j),da(i,j),tmin(i,j,k)
enddo
enddo
enddo
close(3333)
a=0
do i=1,sta_num
     do j=1,ns
         do k=1,365
if(tmin(i,j,k)/=0)then
a=a+1
tmin1(i)=a
endif
  enddo
enddo
enddo
do i=1,sta_num
    do j=1,ns

ave2(i)=tmin1(i)/ns
            write(1111,*) sta(i,j),ave2(i)
         enddo

         enddo
end

你看那河下的桥

求大神帮助

胡图图

想要问你是怎么筛选的。我现在处理数据，不懂怎么筛选，求赐教