Skip to content

Conversation

@trontrytel
Copy link
Member

This PR fixes the latent heat bug discovered by @jbphyswx in our buoyancy gradients.

It would be great to rewrite parts of that file to get rid of the getter functions at the top and the bg_model, and to reduce the number of functions we define. But maybe we can do it in the next PRs.

get_qi_sat(thermo_params, bg_model) *
TD.latent_heat_sublim(thermo_params, t_sat)
) /
(get_ql_sat(thermo_params, bg_model) + get_qi_sat(thermo_params, bg_model))
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

add an eps to make sure this is not divided by zero?

Copy link
Member

@szy21 szy21 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good, just one comment.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants