代码错误Incompatible ranks 0 and 2 in assignment

fcode

matmul 的结果是二维数组，而 ta4 是单变量。无法存储matmul的结果

fcode

是 1*1 的二维数组。就好像 real :: a(1,1) 一样
你可以用
ta4 = sum(REAL(MATMUL(MATMUL(ta2,vel(1+(i-1)*num_wan:i*num_wan,:)),ta1)))
得到单变量。
这个 sum 是假装的。对一个1*1的二维数组求和，就能得到它里面唯一的那个元素。