Pablo Marchant

Not a proper fix, but here is a workaround I found. I just create a twin axis in each direction, match the x and y boundaries and remove labels using...

This ok for me as well, though it would be ideal to do a run with fpe checks. In the way rotation is split up, with various things computed at...

This is also a good opportunity to make it much clearer what is and what is not kept constant through a timestep. This way we can have `s% start% D_mix`...

> Remember that we can pass these sub-structs directly to subroutines, so we could directly pass the iter struct and then call iter%lnR. I know that, was bringing it up...

Think we're converging too. Don't mind keeping `iter` as the name even if it is updated at other places. A call to update the values on iter would happen after...

Can provide a few examples of cases that produce supersonic convection with Cox next week. It's not just crazy exploding stuff, happens as well in much more run-of-the-mill stars.

The calls to cash_karp in hydro_rotation come from the old way to compute the rotation corrections that is no longer the default. Kept those as an alternative when we introduced...

A bit of a tangential thing, but how would people feel about changing `secyer` to `secyear`? Guess it would break a bit of external code, but in a pretty trivial...

Thanks for pointing this out. This is an issue that's been hanging around since a couple of MESA versions where the indexing of the mixing regions changed. As a workaround...

I can reproduce locally, checking now