164 |
DO J=1,sNy |
DO J=1,sNy |
165 |
DO I=1,sNx |
DO I=1,sNx |
166 |
cg2d_q(I,J,bi,bj) = |
cg2d_q(I,J,bi,bj) = |
167 |
C & pW(I ,J ,bi,bj)*cg2d_r(I-1,J ,bi,bj) |
& pC(I ,J ,bi,bj)*cg2d_r(I ,J ,bi,bj) |
168 |
C & +pW(I+1,J ,bi,bj)*cg2d_r(I+1,J ,bi,bj) |
& +pW(I ,J ,bi,bj)*cg2d_r(I-1,J ,bi,bj) |
169 |
C & +pS(I ,J ,bi,bj)*cg2d_r(I ,J-1,bi,bj) |
& +pW(I+1,J ,bi,bj)*cg2d_r(I+1,J ,bi,bj) |
170 |
C & +pS(I ,J+1,bi,bj)*cg2d_r(I ,J+1,bi,bj) |
& +pS(I ,J ,bi,bj)*cg2d_r(I ,J-1,bi,bj) |
171 |
C & +pC(I ,J ,bi,bj)*cg2d_r(I ,J ,bi,bj) |
& +pS(I ,J+1,bi,bj)*cg2d_r(I ,J+1,bi,bj) |
172 |
& cg2d_r(I ,J ,bi,bj) |
C & cg2d_r(I ,J ,bi,bj) |
173 |
etaN = etaN |
etaN = etaN |
174 |
& +cg2d_q(I,J,bi,bj)*cg2d_r(I,J,bi,bj) |
& +cg2d_q(I,J,bi,bj)*cg2d_r(I,J,bi,bj) |
175 |
ENDDO |
ENDDO |