Initial T_d = T_g for our 2D self-shadowing test Sect. (3.3), based on dg3t100 from Melon Fuksman & Klahr (2022), shown above, and the same for our 3D planet-driven spiral test Sect. (3.4), based on dg3t0.1, shown below. We overplot the contours for the radial optical depth, τ_r, in red (computed using opacities Planck-averaged over the stellar spectrum) and cooling parameter of β = t_cΩ. Our calculations of these quantities assume an MRN distribution of small silicon-graphite grains and a 1M_⊕ classical T Tauri star as the source of light and gravity. In the lower plot, a white dot indicates the position of the embedded 1 M_J planet, while grey shading covers the region outside the r = [8,100] au boundaries of our 3D simulation.