fortran子程序怎样输出数组

Fortran子程序可以使用DO循环语句遍历数组中的元素，并使用WRITE语句将每个元素输出到控制台。下面是一个简单的示例：

SUBROUTINE PRINT_ARRAY(A)
INTEGER, DIMENSION(:) :: A
INTEGER :: I

DO I = 1, SIZE(A)
WRITE(*,*) A(I)
END DO
END SUBROUTINE

在这个子程序中，我们首先定义了一个名为“A”的整数数组，然后使用DO循环遍历数组中的元素，并使用WRITE语句将每个元素输出到控制台。在调用这个子程序时，我们只需要将要输出的数组作为实参传递给该子程序即可：CALL PRINT_ARRAY(MY_ARRAY)。

相关问题

fortran子程序传入可变数组要在module里实现

在 Fortran 中，如果要在一个子程序中传入可变数组，你需要使用 Fortran 中的 `POINTER` 数据类型。同时，为了在模块中实现这个功能，你需要将这个数组声明为 `POINTER` 类型的变量，并将其在模块中声明。在子程序中，你需要使用 `ALLOCATABLE` 属性来为这个数组分配内存，并将其绑定到 `POINTER` 变量上。这样就可以在子程序中对这个可变数组进行操作了。下面是一个示例代码：

module example_module
  implicit none
  
  interface
    subroutine example_subroutine(arr, n)
      integer, dimension(:), pointer :: arr
      integer :: n
    end subroutine example_subroutine
  end interface
  
contains
  
  subroutine example_subroutine(arr, n)
    integer, dimension(:), pointer :: arr
    integer :: n
    allocate(arr(n))
    
    ! 对数组进行操作
    
    deallocate(arr)
  end subroutine example_subroutine
  
end module example_module

在这个示例代码中，我们在模块中声明了一个名为 `example_subroutine` 的子程序，并将其声明为 `POINTER` 类型的 `arr` 数组和一个整型变量 `n` 的输入参数。在子程序中，我们使用 `ALLOCATABLE` 属性为 `arr` 数组分配了 `n` 个整型空间，并在操作完成后使用 `DEALLOCATE` 释放了这些空间。

fortran 子程序矩阵求逆

下面是一个简单的 Fortran 子程序，用于计算矩阵的逆：

subroutine invert(n, a, ainv)
! 计算矩阵 a 的逆矩阵 ainv
! n 是矩阵维数
! a 是原始矩阵
! ainv 是逆矩阵
implicit none
integer, intent(in) :: n
real*8, intent(in) :: a(n,n)
real*8, intent(out) :: ainv(n,n)
integer :: i, j, k
real*8 :: factor, temp(n,n)
! 初始化逆矩阵为单位矩阵
ainv = 0.0
do i = 1, n
  ainv(i,i) = 1.0
end do
! 高斯-约旦消元法计算逆矩阵
do k = 1, n
  factor = a(k,k)
  do j = 1, n
    a(k,j) = a(k,j) / factor
    ainv(k,j) = ainv(k,j) / factor
  end do
  do i = 1, n
    if (i /= k) then
      factor = a(i,k)
      do j = 1, n
        a(i,j) = a(i,j) - factor * a(k,j)
        ainv(i,j) = ainv(i,j) - factor * ainv(k,j)
      end do
    end if
  end do
end do
end subroutine

该程序使用高斯-约旦消元法计算逆矩阵。在调用该子程序之前，您需要将矩阵 `a` 和逆矩阵 `ainv` 声明为 `real*8` 类型的二维数组，并将矩阵大小 `n` 传递给子程序。