面试题：Fortran在石油勘探数据处理中的数组操作

Question

在石油勘探数据处理中，经常需要处理大量的数组数据。假设你有一个二维数组，代表地震波数据，行表示不同的采样点，列表示不同的时间点。请用Fortran代码实现对该数组每一列的均值计算，并将结果存储在一个一维数组中。

Answer 1

program calculate_column_mean
    implicit none
    integer, parameter :: num_samples = 100  ! 假设采样点数
    integer, parameter :: num_time_points = 50  ! 假设时间点数
    real :: seismic_data(num_samples, num_time_points)
    real :: column_mean(num_time_points)
    integer :: i, j

   ! 初始化地震波数据，这里假设数据已有值，实际应用中应从文件读取等方式获取
    do i = 1, num_samples
        do j = 1, num_time_points
            seismic_data(i, j) = real(i * j)
        end do
    end do

    do j = 1, num_time_points
        column_mean(j) = sum(seismic_data(:, j)) / real(num_samples)
    end do

    ! 输出结果
    do j = 1, num_time_points
        write(*,*) 'Column ', j,'mean: ', column_mean(j)
    end do
end program calculate_column_mean