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