Comparison between the analytical models and the simulation after eight outer orbits and at z = 20 au. Left panel: density field. Right panel: Pressure gradient term of the rotation curve. The red dots represent the simulation data, while the blue and orange line show the thermally stratified model and the isothermal one, respectively. The model with thermal stratification matches very well the simulation.