如何使用Fortran命令移除矩阵指定列？附代码场景咨询

Fortran标准库里并没有专门的"移除列"单一命令，但有几种简洁高效的方式能实现你想要的操作，完全适配你现有代码的逻辑。

方法1：直接用数组切片拼接 #

既然你要移除的列位置是固定的（第97和113列），可以直接通过数组切片把需要保留的列分段赋值到AA2中。结合你的代码来看，AA1是120×120的矩阵，移除两列后AA2是120×118的矩阵，刚好匹配你定义的RN=118：

! 赋值前96列
AA2(:, 1:96) = AA1(:, 1:96)
! 赋值第98到112列（跳过第97列），对应AA2的第97到111列
AA2(:, 97:111) = AA1(:, 98:112)
! 赋值第114到120列（跳过第113列），对应AA2的第112到118列
AA2(:, 112:118) = AA1(:, 114:120)

这种方法的优势是性能最优，Fortran编译器会对数组切片操作做高度优化，比手动写嵌套循环效率高得多。

方法2：使用逻辑掩码（更灵活） #

如果后续需要移除的列数不确定或者位置动态变化，用逻辑掩码的方式会更简洁易维护：

Logical, allocatable :: keep_cols(:)
! 创建与AA1列数相同的逻辑数组
Allocate(keep_cols(size(AA1, 2)), stat=ierr)
If (ierr .ne. 0) Then
  print *, 'allocate mask failed'
  pause
End If

! 默认保留所有列，标记要移除的列
keep_cols = .True.
keep_cols([97, 113]) = .False.

! 直接提取保留的列到AA2
AA2 = AA1(:, keep_cols)

这种方法不需要计算切片的位置偏移，代码可读性更强，适合复杂的列筛选场景。

注意事项 #

• 你已经定义RN=3*NBE-2=118，刚好是120-2=118，所以AA2的维度分配是完全正确的；
• Fortran是列优先存储的，数组切片操作会直接利用连续内存块，效率比行优先语言更高。

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