- Cleanup of VerticalDynamicsFEM
- Added 2nd order vertical hyperdiffusion to theta
- Added access to vertical hyperdiffusion to ThermalBubbleCartesianTest
- Added fourth order vertical hyperviscosity on theta
- Bugfix:  Throws exception without vertical hyperdiffusion
Merge remote-tracking branch 'origin/WorkingJEG'
Schar mountain case runs but shows instability. May involve the
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
Did some work testing the memory issue. Hyperdiffusion seems to be an
issue among the 3 existing cartesian grid test cases. Working on
implementing more test cases now.
- Rearranged Jacobian matrix (and F vector) to be a nearly banded
structure
- Added option to call diagonal matrix solver in VerticalDynamicsFEM
- Remedied conflict
- Removed erroneously retained version marker
- Added theta diffusivity to direct diagonal solver
- Changed default vertical solver to direct diagonal solver
- A little bit of code cleanup / maintenance in
InertialGravityCartesianXZTest
- Added model top Rayleigh damping to HorizontalDynamicaFEM
- Moved differentiation functionality from VerticalDynamicsFEM to
GridGLL
- GridCartesianGLL::Initialize() and GridCSGLL::Initialize() now call
GridGLL::Initialize()
- Moved call to Grid::InitializeVerticalCoordinate() from
GridCartesianGLL::Initialize() and GridCSGLL::Initialize() to
GridGLL::Initialize()
- Added new pressure variable in Grid to store pointwise pressure data;
pressure is now pre-computed in HorizontalDynamicsFEM
- Removed Orthonormalization objects from GridPatch since they’re not
being used
- Learned how to use the “mutable” keyword
- BUGFIX: Incorrect set of quadrature points in GridGLL
- Random cleanup