请教一下诸位 关于周期边界

hy304hy315

下面是计算泊肖叶流动 和  coutte流动  处理左右边界的代码  请教一下诸位 下面代码有问题吗？ 先谢谢了
j是y方向循环变量 Q是离散的速度个数 f数组是分布函数 F数组是临时存储分布函数的
for(j=1; j<NY; j++)     //左右边界
      for(k=0; k<Q; k++)
      {
         
                  
        f[0][j][1]=F[NX][j][1];
                f[0][j][5]=F[NX][j-1][5];
                f[0][j][8]=F[NX][j+1][8];
                f[0][j][2]=F[0][j-1][2];
                f[0][j][4]=F[0][j+1][4];
                f[0][j][3]=F[1][j][3];
                f[0][j][6]=F[1][j-1][6];
                f[0][j][7]=F[1][j+1][7];
                f[NX][j][3]=F[0][j][3];
                f[NX][j][6]=F[0][j-1][6];
                f[NX][j][7]=F[0][j+1][7];
                f[NX][j][1]=F[NX-1][j][1];
                f[NX][j][2]=F[NX][j-1][2];
                f[NX][j][4]=F[NX][j+1][4];
                f[NX][j][5]=F[NX-1][j-1][5];
                f[NX][j][8]=F[NX-1][j+1][8];
}