[aspect-devel] Postprocessing scripts

[Marine Lasbleis]

Yes, sure Wolgang :-)

So far, I have done test runs, so basics info like average velocities 
and direct visualisation with paraview were enough.

But to be able to do a good comparisons with existing results, I'll need 
things like horizontal averages (of composition, temperature, 
velocities) and difference between the field and the horizontal average 
(of composition, temperature). But the horizontal here is average over a 
sphere, so I am guessing the post-processors are not so easy. (to get 
average radial profiles)

I'm definitely not used to paraview, and I'm more used to do all my data 
processing by hand... So a good tutorial on how to implement that with 
paraview or equivalent, with a scripting (I usually have several dozens 
of runs to process at a time) would be nice :-) But so far, I think that 
a way to obtain a list of the points and values in the volume (without 
duplicates), and maybe something to do the interpolation at required 
points would be nice.


Best,

Marine


On 2018年08月06日 03:58, Wolfgang Bangerth wrote:
> On 08/04/2018 10:54 PM, Marine Lasbleis wrote:
>> Which, so far, I have not been able to do with ASPECT. I'm kind of 
>> guessing I'm the only one using the full sphere geometry here, and I 
>> don't understand exactly how the post-processors in ASPECT can deal 
>> with what I need.
>
> Marine -- there's nothing inherently wrong to doing these things in 
> Python, but out of curiosity: can you describe what you want to do and 
> what it is you couldn't figure out so far? There are people here on 
> this mailing list who could help! :-)
>
> Best
>  Wolfgang
>