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