Density xz-midplane stratifications superimposed on velocity field vectors (left panels) and an xz-map of the Mach number (right panels), M = v / a (a is the local sound speed) for models M2a (top) and M2b (bottom), at t = 200 yr. The snapshots are a closer view of the star-disk system being photoevaporated. The distance (in AU) to the origin of the coordinate system is depicted for the x and z axes in the left panels. The right-hand side bars on each panel indicate the density scale (in g cm^-3; left panels) and the Mach number (right panels). The labels in the top-rigth panel indicates the region where we have a supersonic outflow (1 and 3), a subsonic outflow (2), and a subsonic inflow (4). We also note that region (1) is neutral, at a temperature T₂ = 10³ K, whereas (3) is ionized, at a temperature T₁ = 10⁴ K.