Jh output3

This is pull request of changes in output system. All commits in this branch are git cherry-picked from JH-output.
I suggest to make a function 

template <int dim>
Feild<dim, Scalar<int>> make_region_id_field(Mesh *mesh)

similar thing is done in darcy_flow_mh to output partition ids.