Fortran已知数组索引但访问at%mag(1,i)报段错误的排查求助

• 派生类型引用逻辑错误：如果at本身是派生类型数组，而非包含二维mag数组的单个派生类型实例，就会触发该问题。Fortran 2003及以后版本支持对派生类型数组的批量成员引用，你写at%mag(:,i)时，实际是批量取所有at元素下mag数组的第i个值，编译器会自动生成临时数组存储结果，所以能正常打印；但你写at%mag(1,i)时，编译器会默认at是单个派生类型实例、mag是二维数组，实际寻址时会访问非法内存，触发段错误。
• 子程序隐式接口导致的参数传递异常：如果你的print子程序没有显式接口（既没有放在module中定义，也没有写interface块声明），传递数组切片at%mag(:,i)时，编译器会自动生成连续的临时数组传入，哪怕原数组地址非法，临时数组刚好拷贝到了栈上的残留有效数据，所以打印正常；但传递单个元素at%mag(1,i)时，编译器直接传入原数组的元素地址，原地址本身非法就会触发段错误。
• 数组维度校验不完整：你只查询了切片at%mag(:,i)的第一维上下界，没有校验原数组at%mag第二维的上下界是否包含i，部分编译器对切片越界的检测逻辑宽松，可能返回残留数据不报错，但单元素访问时会严格触发内存错误。

• 首先校验at的类型属性：执行print *, lbound(at), ubound(at)，确认at是单个派生类型实例还是派生类型数组。如果是派生类型数组，正确的单元素访问写法应为at(k)%mag(1)而非at%mag(1,k)。
• 校验原数组完整上下界：执行print *, lbound(at%mag, dim=1), ubound(at%mag, dim=1), lbound(at%mag, dim=2), ubound(at%mag, dim=2)，确认二维数组两个维度的实际分配范围是否包含你用到的索引1和i。
• 开启编译器调试选项重编运行：使用GCC Fortran编译时添加-g -fcheck=all -fbacktrace参数，使用Intel Fortran编译时添加-g -check all -traceback参数，运行后会直接输出越界访问的具体位置和变量信息，快速定位问题。
• 补全子程序显式接口：将print子程序放到module中定义，或者在调用位置添加interface块声明接口，避免隐式接口导致的参数传递异常掩盖实际问题。

内容的提问来源于stack exchange，提问作者高玉强