Cleaned up the thermo routines a bit
Add function to retrieve prognostic variables from thermo
Fix the check for boundary types, all thermo cases can run now with
surface model
Merge branch 'thermo' into anelastic

Conflicts:
	cases/bomex/bomex.ini
	include/thermo.h
	include/thermo_dry.h
	src/thermo_dry.cxx
	src/thermo_moist.cxx
Add the surface pressure to the drycbl les case