Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aaej_v1_0.gz(35 Kbytes)
Manuscript Title: DUSTCD: a radiative transport code for spherically symmetric dust clouds.
Authors: G.F. Spagna Jr., C.M. Leung
Program title: DUSTCD
Catalogue identifier: AAEJ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 28(1983)337
Programming language: Fortran.
Computer: IBM 3033.
Operating system: MICHIGAN TERMINAL SYSTEM.
RAM: 50K words
Word size: 64
Peripherals: disc.
Keywords: Astrophysics, Radiative transfer, Interstellar matter.
Classification: 1.3, 21.2.

Revision history:
Type Tit le Reference
correction 000A CORRECTION 31/7/84 See below

Nature of problem:
Dust grains are an important component of the interstellar medium. Their presence in the vicinity of stars, in nebulae and in interstellar molecular clouds is confirmed by IR observations. These, coupled with accurate theoretical models, provide a key to understanding the characteristics of the thermally radiating material as well as the physics of the underlying sources which are obscured by the dust. Realistic models require the solution of the problem of photons scattered and absorbed by dust grains.

Solution method:
The radiative transfer problem is cast as a two-point boundary value problem. The resulting system of equations is solved by a Newton- Raphson method.

Restrictions:
The code will handle up to five dust components. A minimum of 60 frequency points, 55 radial points and 60 impact parameters is recommended.

Running time:
For a radial grid of N points, and NF frequency points, the asymptotic estimate of running time is approximated by N**2 * NF+N**3). For a typical model as above, with two dust components, the computation took about 7.6 s per iteration on our IBM 3033, using FORTRAN and double precision. Hence, a typical model should require only one to two minutes of computer time. On a CDC 7600, the program should run about twice as fast.

CORRECTION SUMMARY
Manuscript Title: DUSTCD: a radiative transport code for spherically symmetric dust clouds. (C.P.C. 28(1983)337).
Authors: G.F. Spagna Jr, C.M. Leung
Program title: 000A CORRECTION 31/7/84
Catalogue identifier: AAEJ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 34(1984)223
Classification: 1.3, 21.2.