175 |
ENDDO |
ENDDO |
176 |
|
|
177 |
#ifndef DISABLE_DEBUGMODE |
#ifndef DISABLE_DEBUGMODE |
178 |
IF (debugMode) THEN |
IF ( debugLevel .GE. debLevB ) THEN |
179 |
CALL DEBUG_STATS_RL(1,cg2d_b,'cg2d_b (SOLVE_FOR_PRESSURE)', |
CALL DEBUG_STATS_RL(1,cg2d_b,'cg2d_b (SOLVE_FOR_PRESSURE)', |
180 |
& myThid) |
& myThid) |
181 |
ENDIF |
ENDIF |
197 |
_EXCH_XY_R8(cg2d_x, myThid ) |
_EXCH_XY_R8(cg2d_x, myThid ) |
198 |
|
|
199 |
#ifndef DISABLE_DEBUGMODE |
#ifndef DISABLE_DEBUGMODE |
200 |
IF (debugMode) THEN |
IF ( debugLevel .GE. debLevB ) THEN |
201 |
CALL DEBUG_STATS_RL(1,cg2d_x,'cg2d_x (SOLVE_FOR_PRESSURE)', |
CALL DEBUG_STATS_RL(1,cg2d_x,'cg2d_x (SOLVE_FOR_PRESSURE)', |
202 |
& myThid) |
& myThid) |
203 |
ENDIF |
ENDIF |
206 |
C- dump CG2D output at monitorFreq (to reduce size of STD-OUTPUT files) : |
C- dump CG2D output at monitorFreq (to reduce size of STD-OUTPUT files) : |
207 |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime, |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime, |
208 |
& myTime-deltaTClock) ) THEN |
& myTime-deltaTClock) ) THEN |
209 |
_BEGIN_MASTER( myThid ) |
IF ( debugLevel .GE. debLevA ) THEN |
210 |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_init_res =',firstResidual |
_BEGIN_MASTER( myThid ) |
211 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_init_res =',firstResidual |
212 |
WRITE(msgBuf,'(A34,I6)') 'cg2d_iters =',numIters |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
213 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,I6)') 'cg2d_iters =',numIters |
214 |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_res =',lastResidual |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
215 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_res =',lastResidual |
216 |
_END_MASTER( ) |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
217 |
|
_END_MASTER( ) |
218 |
|
ENDIF |
219 |
ENDIF |
ENDIF |
220 |
|
|
221 |
C-- Transfert the 2D-solution to "etaN" : |
C-- Transfert the 2D-solution to "etaN" : |
356 |
|
|
357 |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime, |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime, |
358 |
& myTime-deltaTClock) ) THEN |
& myTime-deltaTClock) ) THEN |
359 |
_BEGIN_MASTER( myThid ) |
IF ( debugLevel .GE. debLevA ) THEN |
360 |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg3d_init_res =',firstResidual |
_BEGIN_MASTER( myThid ) |
361 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg3d_init_res =',firstResidual |
362 |
WRITE(msgBuf,'(A34,I6)') 'cg3d_iters =',numIters |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
363 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,I6)') 'cg3d_iters =',numIters |
364 |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg3d_res =',lastResidual |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
365 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg3d_res =',lastResidual |
366 |
_END_MASTER( ) |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
367 |
|
_END_MASTER( ) |
368 |
|
ENDIF |
369 |
ENDIF |
ENDIF |
370 |
|
|
371 |
ENDIF |
ENDIF |