RT with T0 = T1 top, free bottom, balance
density
temperature