-
Notifications
You must be signed in to change notification settings - Fork 13
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
use correct beta_u in shallow water linear solvers #544
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.
This puts things in line with the rest of the code ... but the awkward thought I've had is that maybe we don't want the beta
values to be constant at all?
In theory we could update dt
, but if we did then we would want the beta
values to be updated too. At the moment they won't be, so do we actually want:
beta = self.dt*Constant(tau)
instead of
beta_ = self.dt*tau
beta = Constant(beta_)
gusto/solvers/linear_solvers.py
Outdated
beta_u = dt*self.tau_values.get("u", self.alpha) | ||
beta_t = dt*self.tau_values.get("theta", self.alpha) | ||
beta_r = dt*self.tau_values.get("rho", self.alpha) | ||
beta_u_cp = Constant(beta_u * cp) |
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.
Can we remove this variable and just have beta_u*cp in the code below?
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.
Should've known I couldn't get away with being lazy!
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.
Thanks!
Sorry @Witt-D I didn't notice this in your PR.