/[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.16 by mlosch, Thu Jan 17 08:51:15 2013 UTC revision 1.18 by mlosch, Fri Feb 15 15:19:17 2013 UTC
# Line 170  C     constant for the preconditioner Line 170  C     constant for the preconditioner
170         ENDDO         ENDDO
171  C     compute convergence criterion for linear preconditioned FGMRES  C     compute convergence criterion for linear preconditioned FGMRES
172         JFNKgamma_lin = JFNKgamma_lin_max         JFNKgamma_lin = JFNKgamma_lin_max
173         IF ( newtonIter.GT.1.AND.newtonIter.LE.100         IF ( newtonIter.GT.1.AND.newtonIter.LE.SEAICE_JFNK_tolIter
174       &      .AND.JFNKresidual.LT.JFNKres_t ) THEN       &      .AND.JFNKresidual.LT.JFNKres_t ) THEN
175  C     Eisenstat, 1996, equ.(2.6)        C     Eisenstat, 1996, equ.(2.6)      
176          phi_e = 1. _d 0          phi_e = 1. _d 0
# Line 254  C     some output diagnostics Line 254  C     some output diagnostics
254          IF ( krylovIter.EQ.SEAICEkrylovIterMax ) THEN          IF ( krylovIter.EQ.SEAICEkrylovIterMax ) THEN
255           krylovFails = krylovFails + 1           krylovFails = krylovFails + 1
256          ENDIF          ENDIF
257  C     Set the stopping criterion for the Newton iteration  C     Set the stopping criterion for the Newton iteration and the
258          IF ( newtonIter .EQ. 1 ) JFNKtol=JFNKgamma_nonlin*JFNKresidual  C     criterion for the transition from accurate to approximate FGMRES
259            IF ( newtonIter .EQ. 1 ) THEN
260             JFNKtol=JFNKgamma_nonlin*JFNKresidual
261             IF ( JFNKres_tFac .NE. UNSET_RL )
262         &        JFNKres_t = JFNKresidual * JFNKres_tFac
263            ENDIF
264  C     Update linear solution vector and return to Newton iteration  C     Update linear solution vector and return to Newton iteration
265  C     Do a linesearch if necessary, and compute a new residual.  C     Do a linesearch if necessary, and compute a new residual.
266  C     Note that it should be possible to do the following operations  C     Note that it should be possible to do the following operations
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22