[firedrake] building a function space for slice vector fields

[David Ham]

Have you tried it? I think it might Just Work. The first thing which might
break is the HDiv operator, but you never know....
On Mon, 15 Jun 2015 at 11:58 Cotter, Colin J <colin.cotter at imperial.ac.uk>
wrote:

>    Dear all,
>    A lot of what we do in the vertical slice model would be a lot easier
> if we could combine together our u-w components of our velocity field
> (stored as an RT field in the slice) with our v component (stored as a DG
> field in the slice) so that the basis functions are 3D vectors on a 2D x-z
> plane. In particular, then I could assemble the inverse of the combined
> u-w-v velocity mass matrix combined with Coriolis, and make things a lot
> more implicit in the Coriolis term, leading to bigger timesteps when we
> sent the Coriolis parameter to infinity (something we have to do in our
> test case unfortunately).
>
>  Is this supported by current code, if not, how much work would it take
> for us to have this?
>
>  all the best
>  --cjc
>
-------------- next part --------------
HTML attachment scrubbed and removed