代码好长，我有阅读Fortran恐惧症，建议是你首先在逻辑上研究一下有没有并行写，就是最外层循环的两个下标有对数组同一个位置的更新。
另外你刚开始那个代码我觉得珠海要是循环顺序的问题。考虑试一下：
DO L=1,4
c$acc parallel loop
DO K=ifirst1,ilast1
c$acc data copyout(dux(:, k, l))
DO J=ifirst0,ilast0
DUX(J,K,L)= 0.0
ENDDO
c$acc end data
ENDDO
c$acc end parallel loop
ENDDO
