Structure of the numerical grid used in the model. Top panel: slice of one quadrant of the xz plane. To enhance clarity and due to symmetry, only half of the domain is shown. Bottom panel: slice in the xy (equatorial) plane. The dark red dot represents the planet. To enhance clarity, the black grid is four times coarser than the adopted resolution (i.e., each box contains 4 × 4 grid points).