Panel a: integrated intensity image (NOEMA+30 m) of the ¹³CO 1 → 0 line. Contour levels are 15 (≈150 ×σ), 20, and 25 K km s^-1. Crosses indicate the positions of B1b-N and S. The blue and red contours correspond to the blue and red lobes of the B1b-N and S bipolar outflows as traced by the H₂CO emission published by Gerin et al. (2015). The beam is plotted in the left bottom corner. Panel b: the same as a) for the C¹⁸O 1 → 0 line. Contour levels are 4 (≈40 ×σ), 5, 6, 7, 7.5, 8, and 8.5 km s^-1. Panel c: the same as a) for the SO 3₂→ 2₁ line. Contour levels are 1 (≈10 ×σ) to 8 in steps of 0.5 K km s^-1. Panel d: the same as a) for the integrated intensity of all the hyperfine components of the NH₂D 1_1,1a→ 1_0,1s line. Contour levels are 1 (≈10 ×σ) to 8 in steps of 0.5 K km s^-1. Panel e: I(NH₂D)/I(SO) integrated intensity ratio (ratio between the images in panels d and c). The spectra toward the positions labeled in this panel are shown in Fig. A.1.