Fortran模块文件读取优化：实现仅首次调用时读取数据

先直接解答你的第一个疑问：默认情况下，如果你把文件读取逻辑放在模块的顶层代码段（不在任何子程序/函数内部），模块只会在首次被主程序引用时执行这段读取逻辑——而且只会执行一次。因为Fortran模块的初始化代码是在模块被加载时运行的，不管主程序多少次调用该模块，模块都只会被初始化一次。但如果你的读取逻辑是放在某个子程序里，每次调用子程序都触发读取，那肯定会重复读写文件。

下面给你两种实用的实现方案，帮你实现「仅首次读取、全局存储数据」的需求：

方法一：利用模块自动初始化特性（最简单） #

Fortran 2003及以后的标准支持在模块顶层直接执行代码，这部分代码会在模块首次加载时自动运行，之后不会再执行。你可以把文件读取逻辑放在这里，把数据存在模块的全局变量中：

module data_storage
    implicit none
    ! 定义全局存储变量（设为private避免主程序意外修改）
    real, allocatable, private :: my_data(:,:)
    integer, private :: rows, cols

    ! 模块初始化代码：仅执行一次
    call load_data()
    ! 对外暴露需要访问的变量/接口
    public :: get_data_shape, get_data_element

contains

    subroutine load_data()
        implicit none
        integer :: iunit, iostat

        ! 打开数据文件
        open(newunit=iunit, file='data.txt', status='old', action='read', iostat=iostat)
        if (iostat /= 0) error stop "Error: Failed to open data.txt"

        ! 假设文件第一行是数据的行列数，后续是数据内容
        read(iunit, *) rows, cols
        allocate(my_data(rows, cols))
        read(iunit, *) my_data

        close(iunit)
        print *, "Data loaded successfully (only once!)"
    end subroutine load_data

    ! 对外提供获取数据的接口（避免直接暴露全局变量）
    subroutine get_data_shape(r, c)
        integer, intent(out) :: r, c
        r = rows
        c = cols
    end subroutine get_data_shape

    function get_data_element(i, j) result(val)
        integer, intent(in) :: i, j
        real :: val
        val = my_data(i, j)
    end function get_data_element

end module data_storage

主程序使用示例：

program main
    use data_storage
    implicit none
    integer :: r, c

    ! 直接调用接口获取数据，无需手动触发读取
    call get_data_shape(r, c)
    print *, "Data shape: ", r, "x", c
    print *, "First element: ", get_data_element(1,1)

    ! 就算多次引用模块，也不会重复读取文件
    print *, "Another element: ", get_data_element(2,2)
end program main

这种方式完全自动，主程序不需要额外操作，模块加载时就完成了一次性读取，非常省心。

方法二：用标志位手动控制（更灵活） #

如果你的编译器版本较旧（不支持模块顶层call语句），或者你希望在主程序的特定时机触发读取，可以加一个initialized标志位，控制读取逻辑仅执行一次：

module data_storage
    implicit none
    real, allocatable, private :: my_data(:,:)
    integer, private :: rows, cols
    logical, private :: initialized = .false.  ! 初始化状态标志

    public :: load_data_once, get_data_shape, get_data_element

contains

    subroutine load_data_once()
        implicit none
        integer :: iunit, iostat

        ! 已经初始化过，直接返回
        if (initialized) return

        ! 读取文件逻辑和方法一一致
        open(newunit=iunit, file='data.txt', status='old', action='read', iostat=iostat)
        if (iostat /= 0) error stop "Error: Failed to open data.txt"

        read(iunit, *) rows, cols
        allocate(my_data(rows, cols))
        read(iunit, *) my_data

        close(iunit)
        initialized = .true.
        print *, "Data loaded successfully (only once!)"
    end subroutine load_data_once

    ! 以下接口和方法一相同
    subroutine get_data_shape(r, c)
        integer, intent(out) :: r, c
        call load_data_once()  ! 确保数据已加载
        r = rows
        c = cols
    end subroutine get_data_shape

    function get_data_element(i, j) result(val)
        integer, intent(in) :: i, j
        real :: val
        call load_data_once()  ! 确保数据已加载
        val = my_data(i, j)
    end function get_data_element

end module data_storage

主程序使用示例：

program main
    use data_storage
    implicit none
    integer :: r, c

    ! 第一次调用会触发读取，后续调用无操作
    call load_data_once()
    call get_data_shape(r, c)
    print *, "Data shape: ", r, "x", c

    ! 再次调用也不会重读文件
    call load_data_once()
end program main

额外注意事项 #

• 尽量将模块的全局变量设为private，通过公共接口提供访问，避免主程序意外修改数据导致错误。
• 测试时可以在读取子程序里加打印语句，验证是否只执行一次。
• 如果数据文件路径需要灵活配置，可以把路径设为模块的可配置参数，或者通过子程序传入。

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