Pass thv as argument to bu function
Subtract horizontal average pressure field
Some optimization of moist thermo; row running at 125% of dry thermo case
Fixed bug in calcbuoyancy(), replaced 2D arrays with tmp array
Merge branch 'bart' into masterupstream