generic fields - subdomain, region_id

Static functions to create fields returning subdomain (partition) and region_id
of elements for given mesh.
Unit test, fully test only region_id.

Use them in darcy_flow_mh_output to output subdomain.

Add FieldFlag::equation_external_output.
Fix some documentation.