在Fortran中，从数组中随机选择非零元素的最佳方式是什么？

program random_selection
  implicit none
  
  integer, parameter :: n = 10 ! 数组大小
  integer :: array(n) ! 数组
  integer :: non_zero_indices(n) ! 非零元素索引数组
  integer :: non_zero_count ! 非零元素数量
  integer :: random_index ! 随机索引
  integer :: selected_element ! 选择的非零元素
  
  ! 初始化数组
  array = [0, 1, 0, 2, 0, 3, 0, 4, 0, 5]
  
  ! 遍历数组，记录非零元素的索引位置
  non_zero_count = 0
  do i = 1, n
    if (array(i) /= 0) then
      non_zero_count = non_zero_count + 1
      non_zero_indices(non_zero_count) = i
    end if
  end do
  
  ! 生成随机数，范围在1到非零元素的数量之间
  call random_number(random_index)
  random_index = floor(random_index * non_zero_count) + 1
  
  ! 使用随机数作为索引，从非零元素索引位置中选择一个
  selected_element = array(non_zero_indices(random_index))
  
  ! 输出选择的非零元素
  write(*,*) "Selected element:", selected_element
  
end program random_selection