Hello, I am trying to fix up some plots based on Colin's Camassa-Holm demo. http://firedrakeproject.org/demos/camassaholm.py.html I have two questions on plotting. 1. How do we set the y-limits in the context of the plot function? In matplotlib I would do ylim([0,1]) for example but that doesn't seem to work in this context. 2. If I wanted to make a Hovmoeller plot instead of a move I would like to be able to store the snapshots into an array and then plot those. Does someone have a suggestion on how one can do this easily? Cheers, Francis ------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo email: fpoulin@uwaterloo.ca Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637