Hi once more, one remaining question is about the kinetic energy and the equation of state. The specific internal energy should be e = e_int + 1/2 (vr^2 + vz^2 + vphi^2), So I seems neccessary to include the vphi^2-term "everywhere". But where is "everywhere" in the FLASH-code? Thank you, Peter Woitke