如何在Fortran中为嵌套do循环使用多处理？

Question

我有一个包含3个计数器（i，j和k）的嵌套循环。因为我在具有8核多处理器CPU的计算机上运行此代码。我打算删除内部的do循环（do k =），并使代码并行运行，因此，如果k = N（N = 1,2，...，8），则每个内核将分别计算f（i，j，k）。有人可以帮忙吗？

do i=1,nx
    do j=1,ny
        do k=1,8
            f(i,j,k)= omegaP*f(i,j,k)+omega*feq(i,j,k)+fi(i,j,k)*dt
        end do
    end do
end do