- Added theta diffusivity to direct diagonal solver
- Changed default vertical solver to direct diagonal solver