Fortran Coder

查看: 139|回复: 1

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

[复制链接]

3

帖子

1

主题

0

精华

入门

早日渡过科研难关

F 币
129 元
贡献
5 点

规矩勋章

QQ
发表于 2024-1-31 17:25:13 | 显示全部楼层 |阅读模式
1F 币

    !---------------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)


3

帖子

1

主题

0

精华

入门

早日渡过科研难关

F 币
129 元
贡献
5 点

规矩勋章

QQ
 楼主| 发表于 2024-1-31 17:53:48 | 显示全部楼层
解决了,谢谢大家帮忙,感谢感谢
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-2-25 12:33

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表