Extracting elements for given layer from 3D model

lawrence-bird

Hi folks:

I'd like to extract data that are defined on elements for a given layer in a 3D model. This seems to be easy enough to do for data on vertices because md.mesh.numberofvertices is equal to md.mesh.numberofvertices2d * md.mesh.numberoflayers. However, md.mesh.numberofelements is not equal to mdd.mesh.numberofelements2d * md.mesh.numberoflayers.

I see there are md.mesh.upperelements and md.mesh.lowerelements fields, but I am not sure how to extract, for example, layer #3.

Thanks!

Lawrence

heleneseroussi

Hi Lawrence,

Thanks for asking. In the 3d models, there are always one more layer of vertices than elements, so md.mesh.numberofelements should actually be equal to md.mesh.numberofelements2d * (md.mesh.numberoflayers - 1). And to plot the values, you can easily use the 'layer' option of plotmodel and specify which layer you are interested in, so for example:
plotmodel(md,'data',md.materials.rheology_n,'layer',2)

Best
Helene

lawrence-bird

heleneseroussi

Wow, it's easy when you know how, hey?! 🙂

Thanks Helene!