Boussinesq Term in computeMomentumEquation eingefuegt und main um
Temperaturgrid T (initialisiert mit TI) und Waermequellengrid
(initialisiert mit 0) ergaenzt