Hi Jemma, I think I figured it out now. Since the function I want to use is more complicated, I actually wrote a direct par_loop which modifies mesh.coordinates in place, I guess that’s also ok. When I plot the mesh it looks like what I want. Thanks, Eike -- Dr Eike Hermann Mueller Research Associate (PostDoc) Department of Mathematical Sciences University of Bath Bath BA2 7AY, United Kingdom +44 1225 38 5803 e.mueller@bath.ac.uk http://people.bath.ac.uk/em459/
On 27 Feb 2015, at 11:51, Shipton, Jemma <j.shipton@imperial.ac.uk> wrote:
Hi Eike,
Take a look at line 304 of this:
https://bitbucket.org/colinjcotter/slicemodels/src/b6d5eac9dbf28cbcf24aef5be...
An example of how to use this is in examples/nonlinear_mountain.py
I should be around in the chat for a couple of hours if you want to catch me…
Jemma ________________________________________ From: firedrake-bounces@imperial.ac.uk [firedrake-bounces@imperial.ac.uk] on behalf of Eike Mueller [E.Mueller@bath.ac.uk] Sent: 27 February 2015 09:01 To: firedrake Subject: [firedrake] Perturbed extruded meshes
Dear firdrakers,
how can I perturb an extruded mesh? I want to add 'mountains' (i.e. bumps) to my spherical icosahedral mesh, so I want to change the radius according to something like
r = r(theta) = R_earth+H*exp(-0.5*theta^2/D^2)
Can I do this by moving around the underlying 2d host grid and then extruding that? How can I perturb the 2d grid?
Thanks a lot,
Eike
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake <winmail.dat>_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake