换种方法，对于[3, 0, 0]，先找到其中非0数组，然后对非0数组的每一个找到其在gammaA中的位置，并定义stoichh=selio(index)

如果你想要通过这种方式处理数据，首先需要创建一个辅助函数来找到非零元素的位置，然后再遍历这些位置去取selio的相应元素。以下是步骤的示例代码：

pure function find_nonzero_positions(arr) result(positions)
implicit none
integer(kind=int32), intent(in) :: arr(:) ! 输入的整数向量
integer(kind=int32), allocatable :: positions(:) ! 存储非零元素位置的数组

positions = integer(kind=int32) :: [] ! 初始化空数组
integer :: i

do i = 1, size(arr)
   if (arr(i) /= 0) then
      call append(positions, i) ! 将非零元素的索引添加到positions
   end if
end do

if (size(positions) > 0) then
   allocate(positions(size(positions))) ! 确保数组已分配空间
else
   write(*,*) 'All elements are zero.'
   return
end if

find_nonzero_positions = positions
end function find_nonzero_positions

subroutine append(array, value)
implicit none
integer(kind=int32), allocatable, target :: array(:)
integer(kind=int32), intent(in) :: value

if (.not. allocated(array)) then
   allocate(array(1))
   array(1) = value
else
   allocate(new_array(size(array) + 1))
   new_array(1:size(array)) = array
   new_array(size(array) + 1) = value
   deallocate(array)
   array = new_array
end if
end subroutine append

... (在主程序里)

integer(kind=int32), dimension(3) :: c = [3, 0, 0]
integer(kind=int32), dimension(5) :: gammaA = ... ! 定义gammaA
real(kind=real64), dimension(size(c)) :: selio = ... ! 定义selio
integer(kind=int32), allocatable :: stoichh(:)

stoichh = selio(find_nonzero_positions(c))

! 输出 stoichh 的值
if (allocated(stoichh)) then
   print *, "Stoichiometry values:", stoichh
else
   print *, "No non-zero elements in c to get stoichiometry."
end if

这段代码首先定义了一个`find_nonzero_positions`函数来找出`c`中非零元素的索引，然后将这些索引传递给`selio`，生成`stoichh`数组。