[aspect-devel] Evenly spaced vectors, dynamic pressure

[Adam Holt]

Hi all!

I am a relatively new user of ASPECT, and have ran into some tricky things related to visualizing/post-processing ASPECT output using Paraview. I wondered if anyone had experience with the following:

First, plotting evenly spaced velocity vectors for models with mesh refinement. As the element sizes vary dramatically, plotting vectors using "Every Nth Point" (Paraview Glyph option) produces very uneven vector coverage. I thought I had solved this by interpolating the data onto a plane (using Paraview function "Resample with dataset"), but the result of this function (an evenly spaced velocity field) does not time evolve with the simulation (at least for Paraview 5.2.0). Has anybody ran into a similar issue?

Second, I am interested in the dynamic pressure field and wondered how best to retrieve it from my (incompressible) models. For such models, I assume it can be computed by subtracting the horizontally-constant static pressure from the pressure outputted (the 'nonadiabatic pressure' output variable). Is this something that should be done in Paraview, or by writing a new post-processor plugin?

Thanks in advance for any input!
Adam Holt