| 103 |
recip_deltaT = 1. _d 0 / SEAICE_deltaTdyn |
recip_deltaT = 1. _d 0 / SEAICE_deltaTdyn |
| 104 |
|
|
| 105 |
iOutFGMRES=0 |
iOutFGMRES=0 |
| 106 |
C iOutFgmres=1 gives a little bit of output |
C with iOutFgmres=1, seaice_fgmres prints the residual at each iteration |
| 107 |
IF ( debugLevel.GE.debLevA .AND. |
IF ( debugLevel.GE.debLevC .AND. |
| 108 |
& DIFFERENT_MULTIPLE( SEAICE_monFreq, myTime, deltaTClock ) ) |
& DIFFERENT_MULTIPLE( SEAICE_monFreq, myTime, deltaTClock ) ) |
| 109 |
& iOutFGMRES=1 |
& iOutFGMRES=1 |
| 110 |
|
|
| 144 |
I uIce, vIce, |
I uIce, vIce, |
| 145 |
O uIceRes, vIceRes, |
O uIceRes, vIceRes, |
| 146 |
I newtonIter, 0, myTime, myIter, myThid ) |
I newtonIter, 0, myTime, myIter, myThid ) |
| 147 |
|
C probably not necessary, will be removed later: |
| 148 |
CALL EXCH_UV_XY_RL( uIceRes, vIceRes,.TRUE.,myThid) |
CALL EXCH_UV_XY_RL( uIceRes, vIceRes,.TRUE.,myThid) |
| 149 |
C local copies of precomputed coefficients that are to stay |
C local copies of precomputed coefficients that are to stay |
| 150 |
C constant for the preconditioner |
C constant for the preconditioner |
| 358 |
&' // =======================================================' |
&' // =======================================================' |
| 359 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
| 360 |
& SQUEEZE_RIGHT, myThid ) |
& SQUEEZE_RIGHT, myThid ) |
| 361 |
WRITE(msgBuf,'(A)') ' // Begin JFNK statistics' |
WRITE(msgBuf,'(A)') ' // End JFNK statistics' |
| 362 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
| 363 |
& SQUEEZE_RIGHT, myThid ) |
& SQUEEZE_RIGHT, myThid ) |
| 364 |
WRITE(msgBuf,'(A)') |
WRITE(msgBuf,'(A)') |