[firedrake] Perturbed extruded meshes

Shipton, Jemma j.shipton at imperial.ac.uk
Mon Mar 2 10:03:59 GMT 2015


Hi Eike,

Yes, that should work too - we do something like that for periodic coordinates.

Glad it's working for you!

Jemma
________________________________
From: firedrake-bounces at imperial.ac.uk [firedrake-bounces at imperial.ac.uk] on behalf of Eike Mueller [e.mueller at bath.ac.uk]
Sent: 02 March 2015 10:00
To: firedrake
Subject: Re: [firedrake] Perturbed extruded meshes

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 at bath.ac.uk<mailto:e.mueller at bath.ac.uk>
http://people.bath.ac.uk/em459/

On 27 Feb 2015, at 11:51, Shipton, Jemma <j.shipton at imperial.ac.uk<mailto:j.shipton at imperial.ac.uk>> wrote:

Hi Eike,

Take a look at line 304 of this:

https://bitbucket.org/colinjcotter/slicemodels/src/b6d5eac9dbf28cbcf24aef5be60e2c5756d0b677/slicemodels.py?at=master

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 at imperial.ac.uk [firedrake-bounces at imperial.ac.uk] on behalf of Eike Mueller [E.Mueller at 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 at imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake
<winmail.dat>_______________________________________________
firedrake mailing list
firedrake at imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake

-------------- next part --------------
A non-text attachment was scrubbed...
Name: winmail.dat
Type: application/ms-tnef
Size: 7886 bytes
Desc: not available
URL: <http://mailman.ic.ac.uk/pipermail/firedrake/attachments/20150302/8b774538/attachment.bin>


More information about the firedrake mailing list