Velocity-integrated intensity maps (color bars in units of Jy beam⁻¹ km s⁻¹) of para-H₂CO 3₀₃–2₀₂ (left panels),3₂₂–2₂₁ (central panels), and 3₂₁–2₂₀ (right panels) of N113 (integrated local standard of rest velocity range: 229–240 km s⁻¹; top panels) and N159W (integrated velocity range: 232–244 km s⁻¹; bottom panels). The centers of the fields are those also given in Sect. 2.2, i.e., α₂₀₀₀ = 05:13:18.2, δ₂₀₀₀ = –69:22:25.0 and α₂₀₀₀ = 05:39:36.0, δ₂₀₀₀ =– 69:45:25.0. The contours provide, with respect to the peak intensity, levels from 10 to 100% with steps of 10% (~10σ and ~7σ) for para-H₂CO 3₀₃–2₀₂ in N113 and N159W, and from 10% to 100% with steps of 10% (~5σ) for para-H₂CO 3₂₂-2₂₁ and 3₂₁ –2₂₀. The para-H₂CO 3₀₃–2₀₂ peak intensities are 1.2 and 0.5 Jy beam⁻¹ km s⁻¹, the para-H₂CO 3₂₂–2₂₁ peak intensities are 0.5 and 0.2 Jy beam⁻¹ km s⁻¹, and the para-H₂CO 3₂₁–2₂₀ peak intensities are 0.5 and 0.2 Jy beam⁻¹ km s⁻¹ in N113 and N159W, respectively. The pixel size of each image is 0.′′2×0.′′2. The beam size of each image is shown in the lower left corner. Gray circles show the half power beam sizes of the N113 and N159W APEX-12m data, reported by Tang et al. (2017b). The shown fields of view encompass slightly more than our ALMA primary half power beam width (Sect. 2.2). No primary beam correction has been applied. Angular scales of 4″ correspond to a linear scale of ~1 pc.