Compute a boundary flux after solving the PDE

See Read before posting: How do I get my question answered?