Sorry, I forgot mentioning that my test run shell_simple_3d.prm with the output option nonadiabatic temperature. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161028/283271e6/attachment.html>