Applying a coordinate system retroactively changes all points, but when calculating new points they are not in the newly applied coordinate system. C
Currently the user must do the following to get all their calculated points in the same coordinate system:
h = shoulder.Humerus(stl_bone)
h.canal.axis([0.5, 0.8]) #csys applied
h.trans_epiconylar.axis() #csys applied
h.bicipital_groove.axis() #csys applied
h.apply_csys_canal_transepiconylar()
h = shoulder.Humerus(stl_bone)
h.canal.axis([0.5, 0.8]) #csys applied
h.trans_epiconylar.axis() #csys applied
h.apply_csys_canal_transepiconylar()
h.bicipital_groove.axis() #csys NOT applied