RT with T0 = T1 top, free bottom
density
temperature