-
Notifications
You must be signed in to change notification settings - Fork 19
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Forcing for ISDAC #3025
Forcing for ISDAC #3025
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'll take a closer look at the AtmosphericProfiles shortly and comment here once that is done, but this LGTM overall.
acbc977
to
08dc509
Compare
@akshaysridhar could you have a look at the config file and see if there's something I should change to make the test pass? I'm probably just misconfiguring resolution / time step / other knobs relative to a robust LES setup. |
Where is |
When I ran this, the simulation crashed in the first few iterations until I reduced the timestep to satisfy the cfl condition (thanks @akshaysridhar). |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
device: CUDA
is invalid. Use device: CUDADevice
@@ -189,3 +189,58 @@ function external_forcing_tendency!(Yₜ, Y, p, t, ::GCMForcing) | |||
|
|||
return nothing | |||
end | |||
|
|||
# ISDAC external forcing (i.e. nudging) | |||
external_forcing_cache(Y, external_forcing::ISDACForcing) = (;) # Don't need to cache anything |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@Sbozzolo
Line is missing parameters argument.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Resolved in latest push.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
GCMForcing{DType}([...])
is in the wrong elseif block. Moved up into appropriate section for GCMForcing.
Purpose
See CliMA/AtmosphericProfilesLibrary.jl#41 for forcing
To-do
Content