Extract boundary region from a DG field
The current code would fail to extract the boundary field using "FieldConvert -m extract ...", if the field is compressible 3D and there is boundary layer mesh.
The error could vary, as it showed "Segmentation fault" in a real case while it just freezed in the attached test case.
According to Spencer, it seems that since by default the FieldConvert assumes it is a CG field it also defaults to a DirectMultiLevelStaticCond for the GlobalSysSoln which then tries to set up a fancy numbering scheme on this mesh which is failing. test_case.zip