  In this paper we investigate how the inclusion of scattering of the stellar
radiation into a passive flaring disk model affects its structure and spectral
energy distribution, and whether neglecting it could significantly decrease the
model reliability. In order to address these questions we construct a detailed
1+1D vertical structure model in which the scattering properties of the dust
can be varied. Models are presented with and without dust scattering, and for
different albedos and phase functions. It is found that scattering has the
effect of reducing the disk temperature at all heights, so that the disk
"shrinks", i.e., the the density at all intermediate heights decreases.
However, this effect in most cases is more than compensated by the increase of
the total extinction (absorption + scattering) cross section, so that the
surface scale height increases, and images in scattered light will see a
slightly thicker disk. The integrated infrared emission decreases as the albedo
increases, because an increasing part of the flux captured by the disk is
reflected away instead of absorbed and reprocessed. The reduction of the
infrared thermal emission of the disk is stronger at short wavelengths (near
infrared) and practically negligible at millimeter wavelengths. For relatively
low albedo (alb <~ 0.5), or for strongly forward-peaked scattering (g roughly
>0.8), the infrared flux reduction is relatively small.
