Fortran 程序：读取风速文件

program read_wind_file
    implicit none
    ! 隐式声明变量
    integer :: i, n
    real :: x, y, z
    real, dimension(:), allocatable :: col1, col2, col3
    ! 声明变量和数组
    open(10, file='wind.txt')
    ! 打开文件
    n = 0   
    ! 初始化计数器
    do while (.true.)
        read(10,*,iostat=i) x, y, z
        ! 读取文件中的数据
        if (i/=0) exit
        ! 如果读取到文件结尾，则退出循环
        n = n + 1
        ! 计数器加1
    end do
    allocate(col1(n), col2(n), col3(n))
    ! 动态分配数组空间
    rewind(10)
    ! 将文件指针重置到文件开头
    do i=1,n
        read(10,*) col1(i), col2(i), col3(i)
        ! 读取文件中的数据并存入数组中
    end do
    close(10)
    ! 关闭文件
    do i=1,n
        write(*,*) col1(i), col2(i), col3(i)
        ! 循环输出数组中的数据
    end do
end program read_wind_file
! 结束程序