面试题：Fortran变量声明与数据类型：基础运用

Question

在Fortran中，声明一个实型数组，数组大小为5，数据类型为双精度，并将其初始化为1.0, 2.0, 3.0, 4.0, 5.0，然后编写一个子例程计算该数组所有元素的和并返回结果。请给出完整代码。

Answer 1

program main
    implicit none
    real(kind=8) :: my_array(5) = [1.0d0, 2.0d0, 3.0d0, 4.0d0, 5.0d0]
    real(kind=8) :: sum_result

    sum_result = calculate_sum(my_array)
    write(*,*) '数组所有元素的和为:', sum_result
contains
    function calculate_sum(arr) result(sum_value)
        real(kind=8), intent(in) :: arr(:)
        real(kind=8) :: sum_value
        integer :: i

        sum_value = 0.0d0
        do i = 1, size(arr)
            sum_value = sum_value + arr(i)
        end do
    end function calculate_sum
end program main