Dear Will, Make a Function for outputting; e.g. energy_fn = Function(R) [if you want to output in the function space R] Then do energy_fn.interpolate(dot(u, u) + rho*rho) [for interpolation] or energy_fn.project(dot(u, u) + rho*rho) [for L^2 projection] followed by localenergyfile.write(energy_fn, time=t) Apologies for any typos in the above. Thanks, Andrew On 3 May 2016 at 14:00, William Booker <scwb@leeds.ac.uk> wrote:

Dear Firedrakers,

For the model,

https://github.com/booker213/compressible_acoustic/blob/master/Firedrake/1D/...

I'm interested in finding out how the energy in my system is distributed at the end of my model.

My naive approach, below, was to try and create a function, so I could visualise it.

localenergyfile = File('./Results/local_energy.pvd') E_local = Function(R, name="Local Energy") E_local = 0.5*(inner(u,u) + rho**2 localenergyfile.write(E_local , time =t)

But this does not work. How would I go about implementing this?

Regards,

Will