Add a thermo module with buoyancy as the prognostic variable
Fixed the dns statistics such that they take the proper buoyancy switch
into account