Profiles from Run 7 at times t_O (i.e. the moment at which the primary protostar forms), t_O+ 5.9 kyr, t_O+ 12.2 kyr, t_O+ 16.1 kyr, t_O+ 20.1 kyr, t_O+ 24.0 kyr, t_O+ 28.0 kyr and t_O+ 31.9 kyr (the end of the simulation). (a) The radial-density profile on the equatorial plane, ρ(r, z= 0). The axes are logarithmic, and lines with slopes −1.5 and −2.0 are shown to mimic the run of density in the disc from 10 to 100 au, and from 100 to 1000 au, respectively. The key gives the line-styles used to represent the different times. (b) The radial-temperature profile on the equatorial plane, T(r, z= 0). The axes are logarithmic, and a line with slope − 0.5 is shown to mimic the run of temperature in the disc between 10 and 1000 au. (c) The radial-velocity profile on the equatorial plane, v_r(r, z= 0). The velocity axis is linear, and velocity is given in units of the local sound speed. The radius axis is logarithmic. (d) The azimuthal velocity profile on the equatorial plane, v_φ(r, z= 0). The axes are logarithmic, and lines with slopes of −1/2 and −1/4 are included to show that in the disc v_φ∝r^−1/4, corresponding to non-Keplerian differential rotation. (e) The radial profile of the disc aspect ratio, ⁠. The axes are linear.