35 |
#endif |
#endif |
36 |
#include "SOLVE_FOR_PRESSURE.h" |
#include "SOLVE_FOR_PRESSURE.h" |
37 |
|
|
38 |
|
C === Functions ==== |
39 |
|
LOGICAL DIFFERENT_MULTIPLE |
40 |
|
EXTERNAL DIFFERENT_MULTIPLE |
41 |
|
|
42 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
43 |
C == Routine arguments == |
C == Routine arguments == |
44 |
C myTime - Current time in simulation |
C myTime - Current time in simulation |
202 |
ENDIF |
ENDIF |
203 |
#endif |
#endif |
204 |
|
|
205 |
_BEGIN_MASTER( myThid ) |
C- dump CG2D output at monitorFreq (to reduce size of STD-OUTPUT files) : |
206 |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_init_res =',firstResidual |
IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime, |
207 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
& myTime-deltaTClock) ) THEN |
208 |
WRITE(msgBuf,'(A34,I6)') 'cg2d_iters =',numIters |
_BEGIN_MASTER( myThid ) |
209 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_init_res =',firstResidual |
210 |
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_res =',lastResidual |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
211 |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
WRITE(msgBuf,'(A34,I6)') 'cg2d_iters =',numIters |
212 |
_END_MASTER( ) |
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
213 |
|
WRITE(msgBuf,'(A34,1PE24.14)') 'cg2d_res =',lastResidual |
214 |
|
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
215 |
|
_END_MASTER( ) |
216 |
|
ENDIF |
217 |
|
|
218 |
C-- Transfert the 2D-solution to "etaN" : |
C-- Transfert the 2D-solution to "etaN" : |
219 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |