运行计算均值、方差、标准差的Fortran程序无响应该如何解决

• 执行顺序错误：主程序第一行就调用stats子程序，此时数组长度N、变量data都未初始化，全是内存随机值，后续逻辑完全无法正常执行，违背了「先采集数据再处理」的基本逻辑
• 存在多处无限死循环：N_loop、K_loop两个循环都没有设置退出条件，程序进入后会一直卡死运行，无任何输出
• 数组操作完全错误：
  • 变量data被声明为单个整数，而非数组，执行allocate(v(data))时data没有合法值，数组分配长度完全错误
  • 读入数组元素时一直写入v(data)同一个位置，其余数组元素全为内存垃圾值
  • K_loop中循环读取同一个位置的元素计算平方和，计算结果完全错误
• 逻辑结构不符合需求：需求是「除数据采集外所有逻辑放在子程序」，但原代码数据采集逻辑一半在主程序一半在子程序，结构混乱
• 隐藏问题：计算均值时直接用整数除法会导致结果截断，精度损失严重
• 作用域问题：response变量在主程序声明，子程序直接跨作用域访问，虽然contains定义的子程序可以访问主程序变量，但逻辑设计极不规范

调整逻辑：将所有数据采集逻辑保留在主程序，统计计算、输出逻辑完全封装在子程序中，删除无意义的死循环，修正数组操作逻辑，修复后的可运行代码如下：

Program dyn_array
implicit none
integer :: N, i
integer, dimension(:), allocatable :: arr
character(len=1) :: response

do
    ! 仅数据采集逻辑放在主程序
    Print *, "Enter the size of array:"
    read(*,*) N
    allocate(arr(N))
    
    Print *, "Enter the values of the array:"
    do i = 1, N
        read(*,*) arr(i)
    end do
    
    ! 所有统计逻辑调用子程序完成
    call stats(arr)
    
    ! 重试逻辑控制
    Print*, "Do you want to repeat for another array?"
    read*, response
    if (response /= "y" .and. response /= "Y") then
        deallocate(arr)
        exit
    end if
    deallocate(arr)
end do

contains
subroutine stats(v)
    implicit none
    integer, dimension(:), intent(in) :: v
    real :: mean, variance, std, tot, sumsquare, arr_size
    integer :: i
    
    arr_size = real(size(v))
    tot = real(sum(v))
    ! 计算均值
    mean = tot / arr_size
    
    ! 计算平方和
    sumsquare = 0.0
    do i = 1, size(v)
        sumsquare = sumsquare + real(v(i))**2
    end do
    
    ! 计算总体方差、标准差，如果需要样本方差，除数改为arr_size-1即可
    variance = (sumsquare - tot**2 / arr_size) / arr_size
    std = sqrt(variance)
    
    ! 输出结果
    Print*, "Mean                   =", mean
    Print*, "Variance               =", variance
    Print*, "Standard deviation     =", std

End subroutine stats

End Program dyn_array

本次修复核心调整点：

• 调整执行顺序，主程序先完成所有数据采集再调用统计子程序，符合需求要求
• 删除所有无退出条件的死循环，重试逻辑移到主程序统一控制
• 修正数组分配、读写逻辑，所有数组元素可以正常赋值
• 统计逻辑完全封装在子程序内，仅接收输入数组做计算输出，符合需求边界
• 所有整数先转浮点数再做除法运算，避免结果截断
• 变量作用域划分合理，代码结构清晰

内容的提问来源于stack exchange，提问作者ChefP905