'''General steps required to define materials with the Lagrangian cloud'''

''Pre-modification thoughts''
 1. Add flag to select between defining materials with material volume from level set functions or from the Lagrangian cloud.
 2. Add material number to cloud particles defined by the initial surface geometries. Provide option for removing selected surface. Surface 1 must never be removed or divFEM will break. Should this even be an option? Should this occur in {{{nliter}}} 1 or using a separate subroutine following the surface geometry definition???
 3. At the time of elemental material definition, count number of each material from cloud in a given element. Provide option for majority rule or others (can we do divFEM with Lag particles??)??? How do we resolve issues with equal numbers of materials? Use old material number or something else??? Note that the general criterion for using the Lag cloud for material definition should check to be sure to use only completely uncut elements, not uncut subelements.
 4. Take decision from count/rule above and give element new material number. For certain properties, be aware we may want to override the default from the material number using something stored on the cloud particles. The material number stored on the cloud should not be changed based on the elemental material number. However, if we want material changes available, we should be able to update those numbers.

''Record of code changes by routine''

'''Dynamic refinement of the surfaces'''