Provides the base case for running a simulation. New cases are
implemented by extending this to specify the initial and boundary
conditions, forcing terms and case-specific postprocessing and analysis.
"Application storage" stores spatial data with a directionality for better cache locality
This set of functions converts indices from this application storage (_dir) to cartesian indices (_ijk)