change advectDensity() and advectTemperature()