- Bugfix:  Throws exception without vertical hyperdiffusion
Cleanup of some commented out code
- BUGFIX: Change reconstruction function order to one higher
- Switched FluxReconstructionType to 2 (mass lumped vertical velocity
matrix)
- Added code for directly computing Jacobian
Merge branch 'master' of https://github.com/paullric/tempestmodel

Conflicts:
	src/atm/VerticalDynamicsFEM.cpp
- Rearranged Jacobian matrix (and F vector) to be a nearly banded
structure
- Added option to call diagonal matrix solver in VerticalDynamicsFEM
- Remedied conflict