Output of p-T-t paths


#1

Hi all,

I would like to create p-T-t paths (pressure-temperature-time paths) for certain points in my models, and I was wondering if there is an easy way to produce those.

For example, is there an easy way to output values of certain particles over time? I am using particles in my models, so an easy way could be to simply follow 1 or more particles as they move through the computational domain through time. Temperature then needs to be interpolated on the particles, and the depth could then be taken as a proxy for (lithostatic) pressure.

If something isn’t readily available, what would be the easiest or most logical way to implement this?

Many thanks for any suggestions.

Jeroen


#2

Hi Jereon,

There is in fact a particle property to track P-T paths via particles. This particle property is utilized in the cookbook composition-passive-particles-properties.prm.

Please let us know if you have any questions about the specific functionality!

Cheers,
John


#3

Hi Jeroen,

there is a particle property called “pT path” you may be able to use for this. It is explained in the cookbook “5.2.4 Using passive and active compositional fields” in the section on Using particle properties. In essence, you just have to add that property to your “List of particle properties” in the input file.

Every time you write particle output, it will now include the pressure and temperature of that particle, so if you want to track a particular particle, you can identify it by its particle id, which is also in the output. I don’t know what would be the best way to plot the pT path though; you can use the ascii data format for the particle output to make it human readable, but each time step would still be a separate file. So you could write a python script (or something like this) to go though each file and extract the values for the particle with the id you want; but alternatively it might be easier to just write a new postprocessor for ASPECT that writes out these values for a particle with a given id into a separate file (similar to the point values postprocessor, which just writes the solution at a given point into a separate file).

Best,
Juliane