Question

如何用fortran读出csv格式的一列数据，算出平均值和标准差？

一直从1955年2月到2005年5月的数据 一共604个月 求具体

Answer 1

楼上的只算了平均值。

如果是这样，那就不必那么麻烦了。

Integer , parameter :: N = 604
Real :: t = 0.0 , d , r
Open( 11 , File = 'xxx.csv' )
Do i = 1 , 4
Read( 11 , * ) line
End Do
Do i = 1 , N
Read( 11 , * ) d , d , r
t = t + r / N
End Do
!// t 为平均值

如果你要算标准差和平均值，那就用一个数组就可以了。
Program BaiduThlws
Implicit None
Integer , parameter :: N = 604
Real :: t = 0.0 , sd = 0.0 , d , wind( N )
Integer :: i
Character( Len = 1 ) :: line
Open( 11 , File = 'xxx.csv' )
Do i = 1 , 4
Read( 11 , * ) line
End Do
Do i = 1 , N
Read( 11 , * ) d , d , wind( i )
t = t + wind(i) / N
End Do
Do i = 1 , N
sd = sd + ( wind(i) - t )**2
End Do
sd = sqrt( sd / N )
!// t 为平均值，sd 为标准差
End Program BaiduThlws

Answer 2

open(10,file='xxx.csv',action='read')
read(10,*)charline1 ! 读取表头
read(10,*)charline2 ! 读取表头
read(10,*)charline3 ! 读取表头
read(10,*)charyear,charmonth, charwind ! 读取表头
do i =1,604
read(10,*) year(i),month(i),wind(i)
end do
totalwind = 0
do i = 1,604
totalwind = totalwind + wind(i)
end do
avewind = totalwind/604
write(*,*)avewind
end