Correct fluxes

SNES Object: 2 MPI processes
  type: newtonls
  maximum iterations=30, maximum function evaluations=2000
  tolerances: relative=0.99, absolute=1e-05, solution=1e-10
  total number of linear solver iterations=9
  total number of function evaluations=2
  SNESLineSearch Object:   2 MPI processes
    type: basic
    maxstep=1.000000e+08, minlambda=1.000000e-12
    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
    maximum iterations=1
  KSP Object:   2 MPI processes
    type: gmres
      GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      GMRES: happy breakdown tolerance 1e-30
    maximum iterations=10000, initial guess is zero
    tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
    left preconditioning
    using PRECONDITIONED norm type for convergence test
  PC Object:   2 MPI processes
    type: asm
      Additive Schwarz: total subdomain blocks = 2, amount of overlap = 5
      Additive Schwarz: restriction/interpolation type - NONE
      Local solve is same for all blocks, in the following KSP and PC objects:
    KSP Object:    (sub_)     1 MPI processes
      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:    (sub_)     1 MPI processes
      type: lu
        LU: out-of-place factorization
        tolerance for zero pivot 2.22045e-14
        matrix ordering: nd
        factor fill ratio given 5, needed 8.07489
          Factored matrix follows:
            Matrix Object:             1 MPI processes
              type: seqaij
              rows=22212, cols=22212
              package used to perform factorization: petsc
              total: nonzeros=4959792, allocated nonzeros=4959792
              total number of mallocs used during MatSetValues calls =0
                using I-node routines: found 5553 nodes, limit used is 5
      linear system matrix = precond matrix:
      Matrix Object:       1 MPI processes
        type: seqaij
        rows=22212, cols=22212
        total: nonzeros=614224, allocated nonzeros=614224
        total number of mallocs used during MatSetValues calls =0
          using I-node routines: found 5553 nodes, limit used is 5
    linear system matrix = precond matrix:
    Matrix Object:     2 MPI processes
      type: mpiaij
      rows=41820, cols=41820, bs=4
      total: nonzeros=1161136, allocated nonzeros=1161136
      total number of mallocs used during MatSetValues calls =0
        using I-node (on process 0) routines: found 5009 nodes, limit used is 5
  0 SNES Function norm 2.418625311837e+01 
    0 KSP Residual norm 5.492008726991e-03 
    1 KSP Residual norm 4.662294753714e-06 
    2 KSP Residual norm 1.641881871011e-06 
    3 KSP Residual norm 1.325240892268e-07 
    4 KSP Residual norm 1.021559349793e-07 
    5 KSP Residual norm 9.411877567164e-08 
    6 KSP Residual norm 9.410247634626e-08 
    7 KSP Residual norm 8.716883682711e-08 
    8 KSP Residual norm 7.167112060299e-08 
    9 KSP Residual norm 3.466561784610e-08 
  Linear solve converged due to CONVERGED_RTOL iterations 9
  1 SNES Function norm 6.413093110621e-03 
Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1