Question

求教，fortran程序求月平均值

我有58年的每日温度数据 求每月平均气温的fortran改如何编写
如
年 月 日 T
1951 1 1 10
1951 1 1 9
.
.
2008 12 31 8

Answer 1

这个问题可以很简单，也可以很复杂。这取决于：

1. 年月日是否按顺序来的？即，同一个月的数据是否连续在一起？

2. 每天的数据都有么？也就是，是否有缺失？

3. 每天的数据是否只有一行？也就是，是否有重复？
   

追问

是按照顺序来的 缺失的数据有特定的数字来表达 每天的数据只有一行 我主要是不知道计算的循环怎样弄 就是怎样判断是否为同一个月的数据

是按照顺序来的 缺失的数据有特定的数字来表达 每天的数据只有一行 我主要是不知道计算的循环怎样弄 就是怎样判断是否为同一个月的数据

追答

那就简单了。

Program Baidu_Thlws
Implicit None
Integer :: year , mon , day , curryear , currmon , iCount , iErr
Real :: total , t
Logical :: bOver
Open( 13 , File = '输出文件.txt' )
Write( 13 , * ) '年 月 平均值'
Open( 12 , File = '文件名.txt' )
Read( 12 , * ) !// 这句是为了跳过表头
Read( 12 , * ) year , mon , day , t
curryear = year
currmon = mon
iCount = 1
total = t
bOver = .false.
Do
Read( 12 , * , ioStat = iErr ) year , mon , day , t
if ( iErr /= 0 ) bOver = .true. !// 如果读错了，就表示文件结束了
if ( ( year /= curryear ) .or. ( mon /= currmon ) .or. ( bOver ) ) then ! // 如果年或月改变了，或者文件结束了
write( 13 , * ) curryear , currmon , total/iCount
total = t
iCount = 1
curryear = year
currmon = mon
else !// 否则是同一个月
if ( abs( t - 特定值 ) > 0.0001 ) then !// 如果 t 不等于特定值（即缺失代表的值）
total = total + t !// 算上总和
iCount = iCount + 1 !// 计数 + 1
end if
end if
if ( bOver ) Exit
End Do
Close( 12 )
End Program Baidu_Thlws