[数值问题] 数组越界

159 帖子	2 主题	1 精华

大师

Vim

F 币: 961 元
贡献: 469 点

发消息

楼主

发表于 2023-9-8 18:23:32 | 显示全部楼层

[Fortran] 纯文本查看 复制代码

  function flux(rho, u, p, i) result(f)
    real, intent(in) :: rho(0:), u(0:), p(0:)
    integer, intent(in) :: i
    real :: f
   
    f = rho(i) * u(i)
   
  end function flux

这种传递数组的方式如果没有写下界，默认是1，你这里需要改为默认下界是0

这个最好写成integer

[Fortran] 纯文本查看 复制代码

integer, parameter :: nx = 20         ! 空间网格数

回复支持反对

使用道具举报

返回列表发新帖