31 |
#include "SIZE.h" |
#include "SIZE.h" |
32 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
33 |
#include "PARAMS.h" |
#include "PARAMS.h" |
34 |
|
#include "GRID.h" |
35 |
#include "CG2D.h" |
#include "CG2D.h" |
36 |
|
|
37 |
C === Routine arguments === |
C === Routine arguments === |
123 |
& -aW2d(I ,J ,bi,bj)*cg2d_x(I ,J ,bi,bj) |
& -aW2d(I ,J ,bi,bj)*cg2d_x(I ,J ,bi,bj) |
124 |
& -aW2d(I+1,J ,bi,bj)*cg2d_x(I ,J ,bi,bj) |
& -aW2d(I+1,J ,bi,bj)*cg2d_x(I ,J ,bi,bj) |
125 |
& -aS2d(I ,J ,bi,bj)*cg2d_x(I ,J ,bi,bj) |
& -aS2d(I ,J ,bi,bj)*cg2d_x(I ,J ,bi,bj) |
126 |
& -aS2d(I ,J+1,bi,bj)*cg2d_x(I ,J ,bi,bj)) |
& -aS2d(I ,J+1,bi,bj)*cg2d_x(I ,J ,bi,bj) |
127 |
|
& -freeSurfFac*zA(i,j,bi,bj)* |
128 |
|
& cg2d_x(I ,J ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm |
129 |
|
& ) |
130 |
err = err + |
err = err + |
131 |
& cg2d_r(I,J,bi,bj)*cg2d_r(I,J,bi,bj) |
& cg2d_r(I,J,bi,bj)*cg2d_r(I,J,bi,bj) |
132 |
sumRHS = sumRHS + |
sumRHS = sumRHS + |
173 |
& +pW(I+1,J ,bi,bj)*cg2d_r(I+1,J ,bi,bj) |
& +pW(I+1,J ,bi,bj)*cg2d_r(I+1,J ,bi,bj) |
174 |
& +pS(I ,J ,bi,bj)*cg2d_r(I ,J-1,bi,bj) |
& +pS(I ,J ,bi,bj)*cg2d_r(I ,J-1,bi,bj) |
175 |
& +pS(I ,J+1,bi,bj)*cg2d_r(I ,J+1,bi,bj) |
& +pS(I ,J+1,bi,bj)*cg2d_r(I ,J+1,bi,bj) |
176 |
C & cg2d_r(I ,J ,bi,bj) |
CcnhDebugStarts |
177 |
|
C cg2d_q(I,J,bi,bj) = cg2d_r(I ,J ,bi,bj) |
178 |
|
CcnhDebugEnds |
179 |
etaN = etaN |
etaN = etaN |
180 |
& +cg2d_q(I,J,bi,bj)*cg2d_r(I,J,bi,bj) |
& +cg2d_q(I,J,bi,bj)*cg2d_r(I,J,bi,bj) |
181 |
ENDDO |
ENDDO |
225 |
& -aW2d(I+1,J ,bi,bj)*cg2d_s(I ,J ,bi,bj) |
& -aW2d(I+1,J ,bi,bj)*cg2d_s(I ,J ,bi,bj) |
226 |
& -aS2d(I ,J ,bi,bj)*cg2d_s(I ,J ,bi,bj) |
& -aS2d(I ,J ,bi,bj)*cg2d_s(I ,J ,bi,bj) |
227 |
& -aS2d(I ,J+1,bi,bj)*cg2d_s(I ,J ,bi,bj) |
& -aS2d(I ,J+1,bi,bj)*cg2d_s(I ,J ,bi,bj) |
228 |
|
& -freeSurfFac*zA(i,j,bi,bj)* |
229 |
|
& cg2d_s(I ,J ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm |
230 |
alpha = alpha+cg2d_s(I,J,bi,bj)*cg2d_q(I,J,bi,bj) |
alpha = alpha+cg2d_s(I,J,bi,bj)*cg2d_q(I,J,bi,bj) |
231 |
ENDDO |
ENDDO |
232 |
ENDDO |
ENDDO |