The formation mechanism of CO clouds observed with the NANTEN2 and Mopra telescopes toward the stellar cluster Westerlund 2 is studied by 3D magnetohydrodynamic simulations, taking into account the interstellar cooling. These molecular clouds show a peculiar shape composed of an arc-shaped cloud on one side of the TeV gamma-ray source HESS J1023-575 and a linear distribution of clouds ( jet clouds) on the other side. We propose that these clouds are formed by the interaction of a jet with clumps of interstellar neutral hydrogen (H I). By studying the dependence of the shape of dense cold clouds formed by shock compression and cooling on the filling factor of H I clumps, we found that the density distribution of H I clumps determines the shape of molecular clouds formed by the jet-cloud interaction: arc clouds are formed when the filling factor is large. On the other hand, when the filling factor is small, molecular clouds align with the jet. The jet propagates faster in models with small filling factors.