Updated this routine for using element families, leading to reduced indexing.  At present, everything is set up so that an element family simply consists of a single material type, and it is no longer possible to mix elements.  This means that there are no longer any arrays with final dimension netypes, and strides through all elements in a family are constant.  It will probably not be too hard to add back in the option of mixing element types, but it is not a priority at present.
Updated this routine for using element families, leading to
reduced indexing.  At present, everything is set up so that an element
family simply consists of a single material type, and it is no longer
possible to mix elements.  This means that there are no longer any
arrays with final dimension netypes, and strides through all elements
in a family are constant.  It will probably not be too hard to add back
in the option of mixing element types, but it is not a priority at
present.


geodynamics/lithomop,geodynamics/lithomop,geodynamics/lithomop,geodynamics/lithomop