[aspect-devel] The output of nonadiabatic_temperature

[Juliane Dannberg]

On 10/11/2016 02:21 PM, Wolfgang Bangerth wrote:
> On 10/11/2016 11:50 AM, Juliane Dannberg wrote:
>>
>> I haven't found a good way to plot the depth average files in 
>> Paraview, and as
>> far as I know we also don't have a visualization postprocessor that does
>> exactly what you want to do (if anyone has a good idea please speak up!)
>> When I used the "depth average" postprocessor I always used pyplot for
>> visualization, so if you already have all of your output as .csv you 
>> could
>> write a python script to make your plot.
>
> It seems to me that what is needed is a visualization postprocessor 
> that for every node of the mesh simply outputs the "depth average" 
> temperature. This would just produce a temperature field that is 
> constant for every depth. One could then, in paraview or visit, plot 
> the difference between the temperature and the depth average temperature.
>
> This should be fairly trivial to implement as a visualization 
> postprocessor.
Yes, that should certainly work, I was just thinking that you would have 
to be careful to not compute the lateral average more often than you 
actually need (only once per depth slice instead of doing the 
computation once per cell).

Best,
Juliane