Several things:

Now the convection coefficient is defined as an operator variable,
so there should be no need to update it manually.

The equation is been rewritten to its dimensionless form.

The unused code has been deleted.
Merge branch 'master' of https://github.com/garmilla/circumbinary