[firedrake] Perturbed extruded meshes

[Eike Mueller]

Dear all,

is there a 2d square mesh with periodic boundary conditions or is this only supported in 1d? I realised that to get steep orography in my setup I would either have to work at very high resolution or shrink the size of the earth.

So I thought instead it would be easier if I work in a 2+1 dimensional box which represents part of the atmosphere, but I do not want to have to worry about boundary conditions in the horizontal condition.

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

> On 2 Mar 2015, at 10:00, Eike Mueller <E.Mueller at bath.ac.uk> wrote:
> 
> 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/ <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 <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 <https://mailman.ic.ac.uk/mailman/listinfo/firedrake>
>> <winmail.dat>_______________________________________________
>> firedrake mailing list
>> firedrake at imperial.ac.uk <mailto:firedrake at imperial.ac.uk>
>> https://mailman.ic.ac.uk/mailman/listinfo/firedrake <https://mailman.ic.ac.uk/mailman/listinfo/firedrake>
> 
> _______________________________________________
> firedrake mailing list
> firedrake at imperial.ac.uk <mailto:firedrake at imperial.ac.uk>
> https://mailman.ic.ac.uk/mailman/listinfo/firedrake <https://mailman.ic.ac.uk/mailman/listinfo/firedrake>
-------------- next part --------------
HTML attachment scrubbed and removed