/[MITgcm]/MITgcm/pkg/seaice/seaice_jfnk.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_jfnk.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by mlosch, Wed Oct 17 14:53:51 2012 UTC revision 1.3 by mlosch, Tue Nov 6 12:53:14 2012 UTC
# Line 132  C     variables DWATN, zeta, and eta) Line 132  C     variables DWATN, zeta, and eta)
132       I      uIce, vIce,       I      uIce, vIce,
133       O      uIceRes, vIceRes,       O      uIceRes, vIceRes,
134       I      newtonIter, 0, myTime, myIter, myThid )       I      newtonIter, 0, myTime, myIter, myThid )
135           CALL EXCH_UV_XY_RL( uIceRes, vIceRes,.TRUE.,myThid)
136  C     local copies of precomputed coefficients that are to stay  C     local copies of precomputed coefficients that are to stay
137  C     constant for the preconditioner  C     constant for the preconditioner
138         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
# Line 196  C     in that routine Line 197  C     in that routine
197          WRITE(msgBuf,'(2A,2(1XI6),2E12.5)')          WRITE(msgBuf,'(2A,2(1XI6),2E12.5)')
198       &       ' S/R SEAICE_JFNK: newtonIter,',       &       ' S/R SEAICE_JFNK: newtonIter,',
199       &       ' total newtonIter, JFNKgamma_lin, initial norm = ',       &       ' total newtonIter, JFNKgamma_lin, initial norm = ',
200       &       newtonIter, SEAICEnewtonIterMax*myIter+newtonIter,       &       newtonIter,SEAICEnewtonIterMax*(myIter-nIter0)+newtonIter,
201       &       JFNKgamma_lin, JFNKresidual       &       JFNKgamma_lin, JFNKresidual
202          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
203       &       SQUEEZE_RIGHT, myThid )       &       SQUEEZE_RIGHT, myThid )
# Line 217  C     solution vector sol = du/vIce Line 218  C     solution vector sol = du/vIce
218  C     residual vector (rhs) Fu = u/vIceRes  C     residual vector (rhs) Fu = u/vIceRes
219  C     output work vectors wk1, -> input work vector wk2  C     output work vectors wk1, -> input work vector wk2
220  C      C    
          CALL EXCH_UV_XY_RL( uIceRes, vIceRes,.TRUE.,myThid)  
          CALL EXCH_UV_XY_RL( duIce, dvIce,.TRUE.,myThid)  
221           CALL SEAICE_FGMRES_DRIVER(           CALL SEAICE_FGMRES_DRIVER(
222       I        uIceRes, vIceRes,       I        uIceRes, vIceRes,
223       U        duIce, dvIce, iCode,       U        duIce, dvIce, iCode,
# Line 248  C     some output diagnostics Line 247  C     some output diagnostics
247          IF ( debugLevel.GE.debLevA ) THEN          IF ( debugLevel.GE.debLevA ) THEN
248           WRITE(msgBuf,'(3(A,I6))')           WRITE(msgBuf,'(3(A,I6))')
249       &        ' S/R SEAICE_JFNK: Newton iterate / total = ', newtonIter,       &        ' S/R SEAICE_JFNK: Newton iterate / total = ', newtonIter,
250       &        ' / ', SEAICEnewtonIterMax*myIter+newtonIter,       &        ' / ', SEAICEnewtonIterMax*(myIter-nIter0)+newtonIter,
251       &        ', Nb. of FGMRES iterations = ', krylovIter       &        ', Nb. of FGMRES iterations = ', krylovIter
252           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
253       &        SQUEEZE_RIGHT, myThid )       &        SQUEEZE_RIGHT, myThid )
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22