We have developed a PDR code to reproduce the high rotational transitions of CO observed with Herschel-PACS. Part of these high-J CO line emission is produced by UV heated outflow walls around protostars. The local FUV radiation flux is calculated by using Monte Carlo method in (${gamma}$, ${alpha}$) grid taking anisotropic scattering into account. Kinetic temperature and Abundance of molecules were computed self-consistently. CO Line fluxes are calculated using RIG. We compare our PDR model with the results by Visser et al (2011) to show that the derived FUV radiation field strength can be affected by the grid resolution near the outflow wall and dust scattering.