在GDB中调试Fortran代码时，如何将大数组输出到TXT文件？

嘿，我完全懂你这种调试Fortran时的头疼事——要把几百上千个元素的数组按指定格式导出到文本文件，又因为没C基础搞不定printf，自带的dump命令又是二进制格式完全不符合需求。别慌，这里有几个纯GDB+Fortran原生的方法，不用碰C代码就能搞定！

相信我，这个方法绝对是Fortran用户的福音——直接在GDB里调用你本来就熟悉的Fortran文件读写语句，完全不用切换到C的思维模式。

步骤很简单：

1. 当程序在断点处暂停后，先调用open打开一个输出文件：

   call open(10, file='array_output.txt', status='replace', action='write')
   

   这里的10是文件单元号，选一个你的程序没在使用的就行（比如10、20这类常用备用号）；status='replace'表示覆盖已有文件，action='write'指定打开用于写入。

2. 用write按你需要的格式输出数组。比如你有一维数组my_arr(1:2000)，想每行输出10个浮点数、每个占8位宽度保留2位小数，就用：

   call write(10, '(10F8.2)') my_arr
   

   这里的(10F8.2)就是标准Fortran格式符，和你平时写程序用的完全一致，上手毫无难度。

3. 输出完成后别忘了关闭文件：

   call close(10)
   

如果是多维数组（比如二维数组my_mat(150, 300)），想按行输出避免Fortran默认的列优先顺序，可以用GDB循环实现：

call open(10, file='matrix_output.txt', status='replace', action='write')
set $row = 1
while $row <= 150
  call write(10, '(300F6.1)') my_mat($row,:)
  set $row = $row + 1
end
call close(10)

如果你只是想快速输出数组、格式要求不复杂，可以用GDB的日志功能把printf输出定向到文件：

1. 先设置日志文件：

   set logging file quick_output.txt
   set logging on
   

2. 用GDB的printf命令输出数组。注意这里要按GDB的内存索引访问（从0开始，和Fortran的1起始不同），arr[0]@N表示从第一个元素开始的N个元素：

   printf "%8.2f ", my_arr[0]@2000
   

3. 输出完后关闭日志：

   set logging off
   

   这个方法胜在快捷，但格式控制不如Fortran的WRITE灵活，比如很难实现自动换行。

如果经常需要导出数组，可以把常用命令写成GDB脚本，比如新建save_fortran_arr.gdb文件：

define save_f_arr
  # 参数：文件名、格式符、数组名
  call open(10, file=$arg0, status='replace', action='write')
  call write(10, $arg1) $arg2
  call close(10)
  printf "数组已保存到 %s\n", $arg0
end

在GDB里加载脚本：source save_fortran_arr.gdb，之后就能一键保存数组：

save_f_arr "my_arr.txt" "(10F8.2)" my_arr

注意事项 #

• 编译Fortran程序时一定要加调试选项，比如gfortran -g -O0，不然GDB没法正确识别数组结构和内容。
• 如果数组在模块里定义，调用时要加上模块名，比如call write(10, '(10F8.2)') my_module::my_arr。
• 选文件单元号时，避开程序正在使用的单元，避免文件操作冲突。

内容的提问来源于stack exchange，提问作者A. Selimov