Consider factoring out sigma coordinate code into a standalone module

Currently, it's associated with the class Model. However, it makes little or no use of class attributes and could be implemented as standalone code. This would make testing it easier. The Model class could maintain wrappers for functions relating to sigma coordinate transformations meaning client code would not need to change.