报错为 error #5082: Syntax error, found END-OF-STATEMENT when expecting o...
!---------------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)
解决了,谢谢大家帮忙,感谢感谢
页:
[1]