Fluxes are lost

SNES Object: 3 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=27
  total number of function evaluations=2
  SNESLineSearch Object:   3 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:   3 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:   3 MPI processes
    type: asm
      Additive Schwarz: total subdomain blocks = 3, 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 7.55587
          Factored matrix follows:
            Matrix Object:             1 MPI processes
              type: seqaij
              rows=16064, cols=16064
              package used to perform factorization: petsc
              total: nonzeros=3347552, allocated nonzeros=3347552
              total number of mallocs used during MatSetValues calls =0
                using I-node routines: found 4016 nodes, limit used is 5
      linear system matrix = precond matrix:
      Matrix Object:       1 MPI processes
        type: seqaij
        rows=16064, cols=16064
        total: nonzeros=443040, allocated nonzeros=443040
        total number of mallocs used during MatSetValues calls =0
          using I-node routines: found 4016 nodes, limit used is 5
    linear system matrix = precond matrix:
    Matrix Object:     3 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 3495 nodes, limit used is 5
  0 SNES Function norm 3.156023513435e+04 
    0 KSP Residual norm 2.706691690707e-03 
    1 KSP Residual norm 1.646464508792e-03 
    2 KSP Residual norm 9.170781030862e-04 
    3 KSP Residual norm 6.009481698118e-04 
    4 KSP Residual norm 3.417126712781e-04 
    5 KSP Residual norm 4.787681320432e-05 
    6 KSP Residual norm 4.560454472526e-05 
    7 KSP Residual norm 3.859597601641e-05 
    8 KSP Residual norm 3.673079306728e-05 
    9 KSP Residual norm 3.445974271820e-05 
   10 KSP Residual norm 2.159415537911e-05 
   11 KSP Residual norm 7.543016336165e-06 
   12 KSP Residual norm 7.416944825904e-06 
   13 KSP Residual norm 6.403254025759e-06 
   14 KSP Residual norm 2.644176505099e-06 
   15 KSP Residual norm 2.496415668030e-06 
   16 KSP Residual norm 2.045313649056e-06 
   17 KSP Residual norm 1.722288627483e-06 
   18 KSP Residual norm 1.461651056621e-06 
   19 KSP Residual norm 7.618542869548e-07 
   20 KSP Residual norm 6.221008888244e-07 
   21 KSP Residual norm 3.435033013676e-07 
   22 KSP Residual norm 2.878937008503e-07 
   23 KSP Residual norm 2.595440945017e-07 
   24 KSP Residual norm 1.259062360140e-07 
   25 KSP Residual norm 7.623152596428e-08 
   26 KSP Residual norm 6.607573140865e-08 
   27 KSP Residual norm 4.066679914334e-08 
   28 KSP Residual norm 2.876876789151e-08 
   29 KSP Residual norm 2.034792834785e-08 
  Linear solve converged due to CONVERGED_RTOL iterations 29
  1 SNES Function norm 2.027633283204e+01