χ² converted probability distribution of the 10 000 parameter set of RADMC-3D models for all clumps. The orange point indicates the best-fit model. Blue crosses give the positions of the 30 best-fit models, so there could be overlaps between the different parameter sets due to the binning. For sources G31, G13, G08a, and G19 the results before re-adjusting T(r) based on SED are shown as light blue contours.