next up previous contents
Next: 29.2.4 Island equations Up: 29.2 stream_function Previous: 29.2.2 The coefficient matrices

   
29.2.3 Solving the equation

After choosing whether five point numerics enabled by option sf_5_point or nine point numerics enabled by option sf_9_point is to be used with Equation (29.14), the elliptic equation is inverted by the method of conjugate gradients. Note that the conjugate gradient solver may fail to converge for large values of the implicit Coriolis parameter acor, which de-symmetrize the equations.



RC Pacanowski and SM Griffies, GFDL, Jan 2000