lid velocity = 0.0204082, prandtl # = 1., grashof # = 1.
  0 SNES Function norm 0.146194
    0 KSP Residual norm 0.146194
    1 KSP Residual norm 0.0554997
    2 KSP Residual norm 0.0103331
    3 KSP Residual norm 0.00349351
    4 KSP Residual norm 0.00162771
    5 KSP Residual norm 0.000420559
    6 KSP Residual norm 0.000131423
    7 KSP Residual norm 4.9738e-05
    8 KSP Residual norm 1.91653e-05
    9 KSP Residual norm 5.05572e-06
   10 KSP Residual norm 1.51319e-06
   11 KSP Residual norm 5.43836e-07
  1 SNES Function norm 2.50119e-05
    0 KSP Residual norm 2.50119e-05
    1 KSP Residual norm 5.42803e-06
    2 KSP Residual norm 1.19801e-06
    3 KSP Residual norm 7.19919e-07
    4 KSP Residual norm 3.23248e-07
    5 KSP Residual norm 8.65301e-08
    6 KSP Residual norm 3.40266e-08
    7 KSP Residual norm 1.01576e-08
    8 KSP Residual norm 1.77734e-09
    9 KSP Residual norm 3.403e-10
   10 KSP Residual norm 1.094e-10
  2 SNES Function norm 1.093e-10
SNES Object: 1 MPI process
  type: newtonls
  maximum iterations=50, maximum function evaluations=10000
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=21
  total number of function evaluations=3
  norm schedule ALWAYS
  Jacobian is built using colored finite differences on a DMDA
  SNESLineSearch Object: 1 MPI process
    type: bt
      interpolation: cubic
      alpha=1.000000e-04
    maxstep=1.000000e+08, minlambda=1.000000e-12
    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
    maximum iterations=40
  KSP Object: 1 MPI process
    type: fgmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    maximum iterations=10000, initial guess is zero
    tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 1 MPI process
    type: fieldsplit
      FieldSplit with MULTIPLICATIVE composition: total splits = 4
      Solver info for each split is in the following KSP objects:
    Split number 0 Defined by IS
    KSP Object: (fieldsplit_x_velocity_) 1 MPI process
      type: preonly
      maximum iterations=10000, initial guess is zero
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using NONE norm type for convergence test
    PC Object: (fieldsplit_x_velocity_) 1 MPI process
      type: ilu
        out-of-place factorization
        0 levels of fill
        tolerance for zero pivot 2.22045e-14
        matrix ordering: natural
        factor fill ratio given 1., needed 1.
          Factored matrix follows:
            Mat Object: (fieldsplit_x_velocity_) 1 MPI process
              type: seqaij
              rows=49, cols=49
              package used to perform factorization: petsc
              total: nonzeros=217, allocated nonzeros=217
                not using I-node routines
      linear system matrix = precond matrix:
      Mat Object: (fieldsplit_x_velocity_) 1 MPI process
        type: seqaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node routines
    Split number 1 Defined by IS
    KSP Object: (fieldsplit_y_velocity_) 1 MPI process
      type: preonly
      maximum iterations=10000, initial guess is zero
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using NONE norm type for convergence test
    PC Object: (fieldsplit_y_velocity_) 1 MPI process
      type: ilu
        out-of-place factorization
        0 levels of fill
        tolerance for zero pivot 2.22045e-14
        matrix ordering: natural
        factor fill ratio given 1., needed 1.
          Factored matrix follows:
            Mat Object: (fieldsplit_y_velocity_) 1 MPI process
              type: seqaij
              rows=49, cols=49
              package used to perform factorization: petsc
              total: nonzeros=217, allocated nonzeros=217
                not using I-node routines
      linear system matrix = precond matrix:
      Mat Object: (fieldsplit_y_velocity_) 1 MPI process
        type: seqaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node routines
    Split number 2 Defined by IS
    KSP Object: (fieldsplit_Omega_) 1 MPI process
      type: preonly
      maximum iterations=10000, initial guess is zero
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using NONE norm type for convergence test
    PC Object: (fieldsplit_Omega_) 1 MPI process
      type: ilu
        out-of-place factorization
        0 levels of fill
        tolerance for zero pivot 2.22045e-14
        matrix ordering: natural
        factor fill ratio given 1., needed 1.
          Factored matrix follows:
            Mat Object: (fieldsplit_Omega_) 1 MPI process
              type: seqaij
              rows=49, cols=49
              package used to perform factorization: petsc
              total: nonzeros=217, allocated nonzeros=217
                not using I-node routines
      linear system matrix = precond matrix:
      Mat Object: (fieldsplit_Omega_) 1 MPI process
        type: seqaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node routines
    Split number 3 Defined by IS
    KSP Object: (fieldsplit_temperature_) 1 MPI process
      type: preonly
      maximum iterations=10000, initial guess is zero
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using NONE norm type for convergence test
    PC Object: (fieldsplit_temperature_) 1 MPI process
      type: ilu
        out-of-place factorization
        0 levels of fill
        tolerance for zero pivot 2.22045e-14
        matrix ordering: natural
        factor fill ratio given 1., needed 1.
          Factored matrix follows:
            Mat Object: (fieldsplit_temperature_) 1 MPI process
              type: seqaij
              rows=49, cols=49
              package used to perform factorization: petsc
              total: nonzeros=217, allocated nonzeros=217
                not using I-node routines
      linear system matrix = precond matrix:
      Mat Object: (fieldsplit_temperature_) 1 MPI process
        type: seqaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node routines
    linear system matrix = precond matrix:
    Mat Object: 1 MPI process
      type: seqaij
      rows=196, cols=196, bs=4
      total: nonzeros=3472, allocated nonzeros=3472
      total number of mallocs used during MatSetValues calls=0
        using I-node routines: found 49 nodes, limit used is 5
Number of SNES iterations = 2
