Divergence of B. In order to preserve to machine precision we need to use H as the fundamental quantity. @tobson has previously suggested that we start out by seeing how bad it is too use B as the fundamental quantity in order to always have incremental translations.
Order of interpolation necessary for the translation operations. We already have fft in place but how well does that work with high frequency noise?
How to handle different timestaggerings without a ton of if-loops.
We could start to think about how to implement the G operators for use in the shearing sheet.