135 |
ENDDO |
ENDDO |
136 |
ENDDO |
ENDDO |
137 |
|
|
138 |
|
#ifndef EXCLUDE_DEBUGMODE |
139 |
|
CALL DEBUG_STATS_RL(1,cg2d_b,'cg2d_b (SOLVE_FOR_PRESSURE)', |
140 |
|
& myThid) |
141 |
|
#endif |
142 |
|
|
143 |
C-- Find the surface pressure using a two-dimensional conjugate |
C-- Find the surface pressure using a two-dimensional conjugate |
144 |
C-- gradient solver. |
C-- gradient solver. |
155 |
I myThid ) |
I myThid ) |
156 |
_EXCH_XY_R8(cg2d_x, myThid ) |
_EXCH_XY_R8(cg2d_x, myThid ) |
157 |
|
|
158 |
|
#ifndef EXCLUDE_DEBUGMODE |
159 |
|
CALL DEBUG_STATS_RL(1,cg2d_x,'cg2d_x (SOLVE_FOR_PRESSURE)', |
160 |
|
& myThid) |
161 |
|
#endif |
162 |
|
|
163 |
_BEGIN_MASTER( myThid ) |
_BEGIN_MASTER( myThid ) |
164 |
WRITE(*,'(A,I6,1PE30.14)') ' CG2D iters, err = ', |
WRITE(*,'(A,I6,1PE30.14)') ' CG2D iters, err = ', |
165 |
& 0, firstResidual |
& 0, firstResidual |