Commit 45221cef authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Fix the calculation of the equilibrium pressure

parent a6760b65
......@@ -108,7 +108,7 @@ double get_thermal_equilibrium_pressure(
/* still cooling, need to decrease u */
if (Lambda_current < 0) {
u_current = u_prev / du;
u_current = u_current / du;
/* heating: interpolate to find eq */
} else {
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment