bbody temp in=? xstart=? xstep=? xend=? xlabel=? xunit=? out=?
This routine calculates for a given (vacuum) wavelength or frequency axis the intensity of a black body at given temperature. The intensity is the energy per unit time, per unit area, per unit solid angle, and per unit frequency (and for all polarisations):
2 h nu^3 1 B_nu = ----- --------- c^2 exp(h nu/kT) - 1
where c is the velocity of light, and h and k are the Planck and Boltzmann constants.
bbody 5500 in=in out=out This calculates the black-body spectrum for 5500 K. The spectrum is written to file OUT. The routine tries to find all necessary information for the 1st (and only) axis in OUT from the spectroscopic axis of the file IN. Since LOGAR is left at its default value of FALSE, the data are intensity in Jy/sr. bbody 2.7 logar=true in=! xstart=0 xstep=0.05 xend=6 xlabel=wavelength xunit=log(micron) out=out This calculates the black-body spectrum for 2.7 K. The spectrum is written to OUT. No input file is specified. The axis contains the logarithms of wavelengths in micron, which run from 0 (1 micron) to 6 (1 metre). Since LOGAR=TRUE, the data are the logarithms of intensity in Jy/sr. bbody 1e6 logar=true in=! xstart=-1 xstep=0.05 xend=2 xlabel=frequency xunit=log10(10**15*Hz) out=out This calculates the black-body spectrum for 1 million K. This time the axis is logarithms of frequency, the units used are 10^15 Hz. The frequency range covered is from 10^14 Hz to 10^17 Hz.
FIGARO A general data reduction system