search

firedrakeproject / gusto

Three dimensional atmospheric dynamical core using the Gung Ho numerics.
http://firedrakeproject.org/gusto/
MIT License
10 stars 10 forks source link

Expand the used to MixedFSOptions to one time discretisation used on the whole equation set #549

Open atb1995 opened 6 days ago

atb1995 commented 6 days ago

Currently, if you try to used MixedFSOptions on the whole equations set it fails. Example code is:

suboptions = {'b': EmbeddedDGOptions()}
opts = MixedFSOptions(suboptions=suboptions)
transport_methods = [DGUpwind(eqns, "u"),
                     DGUpwind(eqns, "p"),
                     DGUpwind(eqns, "b")]
scheme = SSPRK3(domain, options=opts)
! Time stepper
stepper = Timestepper(eqns, scheme, io, transport_methods)

Error is: raise RuntimeError("%r defined on incompatible FunctionSpace!" % r)