* Fix problem where h wasn't initialised properly in GBH
* Add output fn to attach a label to models
* Use H0 instead of fH0*h where possible
* Add some temp. code for testing the offcentre geodesics