Fortran计算任意两个阶段的平均值之差

今天帮朋友写个程序,给出的是一个44年,每年300个站点的月平均数据,排放顺序是44行,每行300列,要求出任意两个阶段平均值的差,并且输出为GrADS的站点文件。数据列数太多,就不贴了。
由于要计算好几组这样的数据,比如1970-1980的均值减去1980-1990的均值这种,前后两个阶段会不断的变化,这也是这个程序里面比较有意思的一个地方,其实整个计算是很简单的,关键是怎样让别人用起来方便,只修改阶段的开始和结束值就能重新计算,而不需要修改程序的主体。
所以这里就要设计一个专门的子程序来计算均值,传入开始和结束的标记,然后调用两次这个子程序,就能求出两个阶段的平均值,再相减即可。
还有一个有趣的就是根据年份来设计输出的文件名,这里用到了内部文件。
程序如下:

Fortran计算任意两个阶段的平均值之差:等您坐沙发呢!

发表评论

点击正确的图片来提交评论:
快捷键:Ctrl+Enter