253 |
firstResidual=0. |
firstResidual=0. |
254 |
lastResidual=0. |
lastResidual=0. |
255 |
numIters=cg2dMaxIters |
numIters=cg2dMaxIters |
256 |
|
c CALL TIMER_START('CG2D [SOLVE_FOR_PRESSURE]',myThid) |
257 |
CALL CG2D( |
CALL CG2D( |
258 |
U cg2d_b, |
U cg2d_b, |
259 |
U cg2d_x, |
U cg2d_x, |
262 |
U numIters, |
U numIters, |
263 |
I myThid ) |
I myThid ) |
264 |
_EXCH_XY_R8(cg2d_x, myThid ) |
_EXCH_XY_R8(cg2d_x, myThid ) |
265 |
|
c CALL TIMER_STOP ('CG2D [SOLVE_FOR_PRESSURE]',myThid) |
266 |
|
|
267 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
268 |
IF ( debugLevel .GE. debLevB ) THEN |
IF ( debugLevel .GE. debLevB ) THEN |
413 |
firstResidual=0. |
firstResidual=0. |
414 |
lastResidual=0. |
lastResidual=0. |
415 |
numIters=cg3dMaxIters |
numIters=cg3dMaxIters |
416 |
|
CALL TIMER_START('CG3D [SOLVE_FOR_PRESSURE]',myThid) |
417 |
CALL CG3D( |
CALL CG3D( |
418 |
U cg3d_b, |
U cg3d_b, |
419 |
U phi_nh, |
U phi_nh, |
422 |
U numIters, |
U numIters, |
423 |
I myThid ) |
I myThid ) |
424 |
_EXCH_XYZ_R8(phi_nh, myThid ) |
_EXCH_XYZ_R8(phi_nh, myThid ) |
425 |
|
CALL TIMER_STOP ('CG3D [SOLVE_FOR_PRESSURE]',myThid) |
426 |
|
|
427 |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime,deltaTClock) |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime,deltaTClock) |
428 |
& ) THEN |
& ) THEN |