-
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
Implicit rain fall always #3515
Conversation
I'll start the |
47fdbfb
to
75b0e5a
Compare
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.
Looks correct to me, thanks! Let's see if the test works.
3133708
to
67da7d6
Compare
MatrixFields.unrolled_map( | ||
name -> (name, @name(c.ρ)) => similar(Y.c, TridiagonalRow), | ||
available_tracer_names_precip, | ||
)..., |
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.
Instead of initializing these blocks and setting them to 0, let's just remove them altogether.
dtγ * -(ᶜprecipdivᵥ_matrix()) ⋅ | ||
DiagonalMatrixRow(Geometry.WVector(1) * ᶠwinterp(ᶜJ, ᶜρ)) ⋅ | ||
ᶠright_bias_matrix() ⋅ DiagonalMatrixRow(-(ᶜwₚ) / ᶜρ) - (I,) | ||
∂ᶜρqₚ_err_∂ᶜρ .= (zero(eltype(∂ᶜρqₚ_err_∂ᶜρ)),) |
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.
Since we're not planning to add these blocks later, we can just get rid of this. Hopefully it doesn't change anything.
67da7d6
to
cf03bf9
Compare
Fixes issue #3511