We have developed an self-consistent PDR model to synthesize warm CO lines of Herschel/PACS observations more accurately. The PDR model solves the FUV continuum radiative transfer, gas energetics, and chemistry simultaneously. A local FUV radiation flux is calculated by using a Monte Carlo method taking anisotropic scattering into account. A new (r, ${delta}$) coordinate system was used, where the r is the distance from the origin and the ${delta}$ is z/$R^2$ in the cylindrical coordinate of (R,z). This is an adequate coordinate system to represent a power-law density of an envelope and a high spatial resolution near the outflow wall. The gas enegetics and chemistry are solved locally and considered $10^4K$ blackbody radiation field instead of the interstellar radiation filed. This newly developed model can be used to analyze quantitatively the effect of UV-heated outflow walls on the warm molecular lines in the embedded proto-stellar objects.