Fortran Coder

标题: 数组传递 [打印本页]

作者: caimaxwell 时间: 2023-12-13 20:33
标题: 数组传递
本帖最后由 caimaxwell 于 2023-12-13 20:35 编辑

大佬们问个问题，TmpiSubblocks数组有值，为什么call之后sbtmp没有被赋值呢?这个数组的值为什么没有被传递过去呢？

3.png (5.95 KB, 下载次数: 64)

2.png (13.32 KB, 下载次数: 70)

1.png (4.83 KB, 下载次数: 61)

作者: fcode 时间: 2023-12-14 08:40
这个不好说，需要更多的代码才好判断。

你可以检查一下 TmpiSubblocks 数组的大小，是否和子程序中 sbtmp 的一致？
尝试用假定形状传递方式？

作者: caimaxwell 时间: 2023-12-14 20:26
十分感谢您的回答！对于接口处的数组，应该设置为传入数组的最大维度这样应该没错吧？问题在于，我在call subroutine前，输出

do j = 1, NLayer
do i = 1,sbnum(j)
write(*,*) j,i,TmpiSubblocks(1,i,j),TmpiSubblocks(2,i,j)
end do
end do

是有值的，然后

subroutine FiLLMatrix(sbtmp, untmp, atmp, amax, unmax, sbmax, cts, sbnum)
  !为了将远近区块填充至所构建的矩阵中
  implicit none
  integer*4  :: unmax, amax, sbmax, cts
  integer*4  :: untmp(2, unmax), atmp(2, amax)
  integer*4  :: sbtmp(2, sbmax, NLayer)
  integer*4  :: i, j, k, j1, j2, k1, k2, t
  integer*4  :: BParam1(sbmax),sbnum(NLayer)
  integer*4  :: row, coL

  allocate(iSubblocks(sbmax,NLayer))

  write(*,*) "FiLLMatrix",sbmax,NLayer

  do i = 1,NLayer
do j = 1,sbnum(i)
      write(*,*) i,j,sbtmp(1,j,i),sbtmp(2,j,i)
end do
  end do
end subroutine !省略版

到这里sbtmp只能输出sbtmp(1,1,1)和sbtmp(2,1,1)的值，其余值为0，我实在想不明白这样为什么会出错。之前的代码是sbtmp(i,j)，有值，为什么改成sbtmp(i,j,k)就出错呢？
我正尝试用假定形状传递。

作者: fcode 时间: 2023-12-14 21:40
试试传入真实的数组大小。
认真检查实参数组的大小，虚参数组的大小，是否一致。比如 write(*,*) shape(tmpiSubblocks)

作者: caimaxwell 时间: 2023-12-26 20:29
确实是这样的问题！是我没有认真细致的检查导致的，谢谢大佬们的解答！十分感谢！

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)