gfortran 9.3.0运行CFD代码出现不定位置Segmentation Fault的原因是什么？

核心成因 #

• 数组越界未被即时捕获：当前编译选项未开启边界检查，数组越界后只会破坏后续使用的内存空间，不会立刻触发报错，因此表现为不同时间步、不同代码位置随机触发段错误。你当前代码中数组维度ID=66,JD=134,KD=66为硬编码，若运行时传入的NI/NJ/NK大于对应硬编码值，或循环边界参数K11/KNK为非法值，都会触发越界。
• Common块不匹配：你使用了-fno-align-commons关闭了common块内存对齐，若不同子程序中声明的同命名common块的变量顺序、类型、数组长度存在差异，会直接导致变量内存错位，读取到的网格参数为垃圾值，引发随机非法内存访问。
• 未初始化变量：网格参数NI/NJ/NK/K11/KNK/I23等如果未在首次使用前被正确赋值，会取内存随机值，部分情况下会触发非法数组下标访问。
• 栈溢出：若代码其他位置存在大尺寸自动数组，系统默认栈大小限制会引发随机段错误。

解决方案 #

1. 新增编译检查选项重新编译：

gfortran -ffree-line-length-none -fno-align-commons -fdefault-real-8 -mcmodel=medium -g -fcheck=all -Wall -Wextra znb3.f90 -o nb3

运行后如果直接抛出数组越界报错，即可定位到非法下标来源，确认是硬编码的ID/JD/KD不足，还是循环边界参数赋值错误。
2. 核验所有common块一致性：逐一核对所有用到/IGRID/、/GRIDB/、/LEVFN/三个common块的子程序的声明内容，确保变量顺序、类型、数组维度完全一致，无缺漏、错写的情况。如果确认common块完全匹配，可以尝试去掉-fno-align-commons参数，避免内存对齐相关的潜在问题。
3. 解除栈大小限制：运行程序前执行命令：

ulimit -s unlimited

排除栈空间不足导致的段错误。
4. 检查网格参数初始化逻辑：确认所有网格相关参数在首次被调用前已经被正确赋值，无未初始化直接使用的情况。

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