Implement ability to specify structure function calculation by order and type (longitudinal or transverse)

There are a variety of ways this could be implemented, one possibly being passing a method to `calculate_structure_function`