- 积分
- 38420
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-20
- 最后登录
- 1970-1-1
![[言之深深] 粉丝数:128 微博数:2073 新浪微博达人勋](source/plugin/sina_login/img/light.png)
成长值: 0
|
发表于 2012-5-23 15:58:45
|
显示全部楼层
本帖最后由 言深深 于 2012-5-23 16:00 编辑
http://bbs.06climate.com/forum.php?mod=viewthread&tid=7133
主要是第二部分,我算的是年,去掉头上两个经纬度,从3--14表示算了12个月的,你要是单算那一个月的也是可以类似的方法。不过你的数据和程序匹配的话也可以不这么做。
parameter(m=?,n=6) !m是你数据的总行数
real x(m,n)
open(1,file="")
do i=1,m
read (1,*) (x(i,j),j=1,n)
enddo
close(1)
do i=1,m
s4=0
s5=0
s6=0
k=0
if(x(i,1).eq.53691) then
if(x(i,2).eq.2005.or.x(i,2).eq.2006)then
if(x(i,3).eq.7)then
s4=s4+x(i,4)
s5=s5+x(i,5)
s6=s6+x(i,6)
k=k+1
endif
endif
endif
enddo
s4=s4/k
s5=s5/k
s6=s6/k
print*,s4,s5,s6end
这是在回复里面临时敲的(只求了一个站点7月的),没有编译一下,你自己试试看吧,能看懂多少就看你自己的了。
|
|