针对你提出的问题——想要以可移植方式输出双精度实数到ASCII文件，确保其他程序能读回完全一致的值，这里有几个可靠的方案：

核心思路：输出足够的有效数字 #

双精度（REAL64）浮点数的尾数是53位二进制，对应大约15-17位十进制有效数字。只要输出至少17位有效数字，就能完整保留所有二进制信息，保证读取时可以精确还原原始值。

1. 直接使用G格式指定17位有效数字（双精度专用） #

G格式会自动适配数字的大小选择科学计数法或普通小数格式，G0则让编译器自动调整输出宽度，避免手动指定的麻烦。这是最简洁的可移植写法：

program main
    use, intrinsic :: ISO_FORTRAN_ENV, only : dp=>REAL64, stdout=>OUTPUT_UNIT
    implicit none
    real(kind=dp), parameter :: pi = 3.141592653589793238462643383279502884197_dp
    write(stdout, '(G0.17)') pi  ! 可移植的全精度输出
end program main

在ifort和gfortran中，这个语句都会输出17位有效数字，足够让MATLAB、Python等程序精确读回原始的双精度值。

2. 实现kind无关的通用写法 #

如果你想让代码适配任意实数kind（单精度、双精度甚至更高精度），可以利用ISO_FORTRAN_ENV提供的REAL_DIGITS常量，它会返回对应kind类型的十进制有效数字位数上限：

program main
    use, intrinsic :: ISO_FORTRAN_ENV, only : dp=>REAL64, stdout=>OUTPUT_UNIT, REAL_DIGITS
    implicit none
    real(kind=dp), parameter :: pi = 3.141592653589793238462643383279502884197_dp
    ! 动态生成格式字符串，适配任意real kind
    write(stdout, '(G0.'//trim(int(REAL_DIGITS(dp)))//')') pi
end program main

这种写法完全不依赖特定编译器的默认行为，跨编译器的可移植性拉满。

3. 避开不可移植的格式写法 #

你代码中使用的'(f)'无参数格式是典型的不可移植写法——ifort能接受，但gfortran要求必须指定宽度和小数位数。这类依赖编译器扩展的写法一定要避免，改用前面提到的明确格式说明。

补充说明 #

你测试中用到的(f30.20)属于过度精度，虽然不会出错，但输出的多余数字并没有实际意义（因为双精度本身无法提供那么多有效信息）。用17位有效数字的G格式既简洁又能保证精度。

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