A dust continuum radiative transfer module

We have developed a module for the dust continuum radiative transfer calculation as part of "Packages of Unified modeling for Radiative transfer, gas Energetics, and Chemistry (PUREC)". PUREC will be applied to interprete observations of protoplanetary disks. When a disk is under the hydrostatic equilibrium condition, the dust temperature and the vertical density structure should be calculated simultaneously. This module calculates the dust temperature by using the method of mean intensity (Lucy et al. 1999). In the very optically thick mid-palne, the Monte-carlo method is not efficient, thus, we apply "modified random walk" and "Partial Diffusion Approximation" to the module. The module has been verified by bechmark tests.