Fixed index bug in the cross sections lngrad function
Make the thermal case running with the new buoyancy model
Merge branch 'master' into pebudget