用fortran计算站点距平

好的，下面是一个简单的Fortran程序，可以计算站点距平。

假设你已经获取了站点的气温数据，数据按照时间顺序存储在一个文本文件中，每行包含站点编号、日期和气温值，用空格或逗号分隔。例如：

1, 1951-01-01, 10.1
1, 1951-01-02, 9.8
...

接下来是Fortran程序的代码：

program calc_anomaly
  implicit none
  
  integer, parameter :: nyears = 61 ! 总共61年
  integer, parameter :: nstations = 160 ! 站点数
  integer :: i, j, k, m, n
  
  real :: temp(nyears, nstations) ! 存储气温数据
  real :: avg_temp(nstations) ! 存储每个站点的平均气温
  real :: anomaly(nyears, nstations) ! 存储每个站点每年的距平
  
  character(10) :: datestr ! 存储日期字符串
  integer :: year, month, day ! 存储年月日
  
  real :: sum_temp, sum_anomaly ! 用于计算平均气温和距平
  
  ! 读取气温数据
  open(10, file='temperature.txt')
  do i = 1, nyears
    do j = 1, nstations
      read(10, *) datestr, temp(i,j)
      ! 解析日期字符串
      read(datestr, '(I4,1X,I2,1X,I2)') year, month, day
      ! 计算每个站点的平均气温
      avg_temp(j) = avg_temp(j) + temp(i,j) / nyears
    end do
  end do
  close(10)
  
  ! 计算每个站点每年的距平
  do i = 1, nyears
    do j = 1, nstations
      anomaly(i,j) = temp(i,j) - avg_temp(j)
    end do
  end do
  
  ! 输出结果
  open(20, file='anomaly.txt')
  do i = 1, nyears
    do j = 1, nstations
      write(datestr, '(I4,1X,I2,1X,I2)') 1951+i-1, 1, 1 ! 构造日期字符串
      write(20, '(I3,1X,A10,1X,F8.2)') j, trim(datestr), anomaly(i,j)
    end do
  end do
  close(20)
  
end program calc_anomaly

这个程序的大致思路是，首先读取气温数据，计算每个站点的平均气温，然后计算每个站点每年的距平，并将结果输出到文件中。

程序中使用了三个数组：`temp` 存储气温数据，`avg_temp` 存储每个站点的平均气温，`anomaly` 存储每个站点每年的距平。其中，`temp` 和 `anomaly` 都是 61 行 160 列的二维数组，分别表示每年和每个站点的数据。`avg_temp` 是 160 个元素的一维数组，表示每个站点的平均气温。

程序中使用了两个文件：`temperature.txt` 存储气温数据，`anomaly.txt` 存储距平数据。在读取气温数据时，程序会依次读取每行数据，并解析出日期和气温值。在计算距平时，程序会依次遍历每个站点和每年的数据，计算出每个站点的平均气温，并将距平存储到 `anomaly` 数组中。最后，程序将距平数据输出到文件中。

请注意，这只是一个简单的示例程序，实际应用中需要根据具体情况进行修改和完善。同时，程序中没有包含任何输入数据的检查和错误处理，需要自行添加。