报错为 error #5082: Syntax error, found END-OF-STATEMENT when expecting o...

bingo · 发表于 2024-1-31 17:25:13

!---------------initialize--------
subroutine initialize
implicit none
include "head.inc"

! 流体变量声明
integer i, j, k, q, x, y, z, obst(lx, ly, lz)
real*8 rho(lx, ly, lz), u_x(lx, ly, lz), u_y(lx, ly, lz), u_z(lx, ly, lz), ff(0:18, lx, ly, lz),ff_old(0:18, lx, ly, lz), ff_eq(0:18, lx, ly, lz),t_s(0:18)

! 初始化流体变量
do X = 1, lx
      do Y = 1, ly
         do Z = 1, lz
            rho(X, Y, Z) = 1.0d0
            u_x(X, Y, Z) = 0.0d0
            u_y(X, Y, Z) = 0.0d0
            u_z(X, Y, Z) = 0.0d0
            obst(X, Y, Z) = 2 !
         end do
      end do
end do
!分布函数
do X = 1, lx
      do Y = 1, ly
         do Z = 1, lz
            do q = 0, 18
                  ff(q, lx, ly, lz) = 0
                  ff_old(q, lx, ly, lz) = 0
            end do
         end do
      end do
end do

! 初始化分布函数
!这里加并行 omp parallel do private schedule(daynamic)
do X = 1, lx
      do Y = 1, ly
         do Z = 1, lz
            do k = 0, 18
                  ff_eq(k, x, y, z) = t_s(k) * rho(X, Y, Z) * &
                     (1.0d0 + 3.0d0 * (u_x(X, Y, Z) * ex(k) + u_y(X, Y, Z) * ey(k) + u_z(X, Y, Z) * ez(k)) + &
                     4.5d0 * (u_x(X, Y, Z) * ex(k) + u_y(X, Y, Z) * ey(k) + u_z(X, Y, Z) * ez(k))**2 - &
                     1.5d0 * (u_x(X, Y, Z)**2 + u_y(X, Y, Z)**2 + u_z(X, Y, Z)**2)
            end do

            do q = 0, 18
                  ff_old(q, X, Y, Z) = ff_eq(q, X, Y, Z)
                  ff(q, X, Y, Z) = ff_eq(q, X, Y, Z)
            end do
         end do
      end do
end do
报错为 error #5082: Syntax error, found END-OF-STATEMENT when expecting one of: ) ,
具体位置在原文这一行1.5d0 * (u_x(X, Y, Z)**2 + u_y(X, Y, Z)**2 + u_z(X, Y, Z)**2)

bingo · 发表于 2024-1-31 17:53:48

解决了，谢谢大家帮忙，感谢感谢

[求助] 报错为 error #5082: Syntax error, found END-OF-STATEMENT when expecting o...

规矩勋章