Fixed index bug in the cross sections lngrad function
Make the thermal case running with the new buoyancy model
Make sure the corners of the grid are written in case of a 2D. Otherwise
some interpolations might pull in bad values