Hi Colin, oh, that's even better! The other test I want to do is do the outer iteration over the velocity, pressure, buoyancy system, i.e. not eliminating the buoyancy pointwise in the linear solve, and switch on orography. I know that it probably does not really make sense to do this in the linear solve (which has not to be done to high accuracy either), and you really want to use the non-linear Newton solve with pointwise buoyancy elimination in the linear solve. But looking at the linear solve could at least give us some ideas, in particular since John seems to be saying that in their paper the model does not work if the orography terms are not treated correctly in the linear solve. Or have you already played around with this and the PETSc solver? Eike On 27/02/15 09:05, Cotter, Colin J wrote:
Hi Eike, It's best to edit the coordinate field of the extruded mesh, particularly as you probably want to use a linear transform in r so that you don't have mountains at the top of the atmosphere too.
Jemma has done this, she can fill you in with the details.
all the best --Colin ________________________________________ 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
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake