Hello Jemma, Thanks for your reply! First of all, what does moist physics imply? do you mean a convection scheme, or that the equations will contain an additional equation for an condensible tracer? Back to what I was talking about, to start with, if I want to run a Held-Suarez benchmark, how would I do it? Even if I begin without having to use climt, how would I go about this? Going through the examples, I see the following code at the end: stepper = Timestepper(state, advection_dict, linear_solver, forcing) ######################################################################## ###### # Run! ############################################################ ################## stepper.run(t=0, tmax=tmax) and the forcing seems to be a class of some sort: forcing = CompressibleForcing(state) for example. Now, I have a couple of questions: 1) what is involved in implementing a Forcing? 2) Suppose I want to implement a GCM with a nominal time step model_time_step = t so I do stepper.run(t=T, tmax=T+model_time_step) Then, I would like to access U,V,T (required for H-S forcing), create a forcing for the next t seconds, and excecute stepper again. Will I be able to do something like this? Pardon me if these issues seem self-evident, but I have no experience working with finite-element frameworks! Thanks, Joy On Sun, Apr 2, 2017 at 5:35 PM, Shipton, Jemma <j.shipton@imperial.ac.uk> wrote:
OK David, thanks - I was just trying not to clutter the list.
So Joy, what kind of problems would you like to run with Gusto?
Jemma ------------------------------ *From:* firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> on behalf of David Ham <David.Ham@imperial.ac.uk> *Sent:* 02 April 2017 16:18:08 *To:* firedrake; joy.merwin@gmail.com *Subject:* Re: [firedrake] Interacting with Gusto
Hi all,
please keep discussions on list to spread the support load and ensure everyone has the chance to learn from the answers.
Regards,
David On Sun, 2 Apr 2017 at 15:58, Shipton, Jemma <j.shipton@imperial.ac.uk> wrote:
Hi Joy,
Yes, sorry, I'm afraid that there isn't much documentation on Gusto at the moment - I hope to write some soon.
We don't have a straightforward way to interact with "physics" just yet, although we are working on implementing moist physics.
I'll email you off-list to chat more about what you need...
Thanks!
Jemma ------------------------------ *From:* firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial. ac.uk> on behalf of Joy merwin monteiro <joy.merwin@gmail.com> *Sent:* 02 April 2017 15:44:21 *To:* firedrake *Subject:* [firedrake] Interacting with Gusto
Hello,
I have recently heard about Firedrake and Gusto, and am looking forward to using it!
I am working on a project named climt (http://climt.readthedocs.io/ en/latest/readme.html), and I was wondering whether I could get Gusto/Firedrake to work with the components already available there.
my question was what interface does Gusto provide for interaction with "physics", like say, a radiative transfer code?
In other words, is there a straightforward way to access model arrays, calculate physics tendencies, and provide them back to the dynamical core?
I tried to look for more documentation about Gusto, but could not find any. I would be grateful if you could point me in the right direction.
TIA, Joy
-- The best ruler, when he finishes his tasks and completes his affairs, the people say “It all happened naturally”
- Te Tao Ch'ing
-- Dr David Ham Department of Mathematics Imperial College London
-- The best ruler, when he finishes his tasks and completes his affairs, the people say “It all happened naturally” - Te Tao Ch'ing