[WIP] Reference density state

This is a first idea of how a reference density implementation for different approximations could work.
I think I don't understand this part, Rene. Why do you only set the temperature if the model is compressible? I would understand that you don't set the *pressure* if the model is not compressible, but why the temperature?