How to calculate the thickness of a specific compositional field?

[jbnaliboff]

Hi Xie,

See the following post for importing ASPECT VTU files into numpy arrays:

Loading VTK data into Python

Hi all, This topic has come up a number of times in recent conversations, so I thought I would share the general workflow for importing VTK data directly into python (numpy arrays). Below is a snippet of python code that should largely work out of the box for 3D cartesian models, where all one has to do is provide the full path and name of the .pvtu in a string. Please note that lines in bold are comments, and you would need to add a # symbol before them in python. I excluded them as they prod…

Cheers,
John