functions in multi.i - _
gb= _multi_bins(nfinal, gb1, gb2, ...)
returns NFINAL+1 boundaries of NFINAL bins constructed by combining
the input bin structures GB1, GB2, etc.
Use NFINAL=0 to get at least the resolution in the finest GBi in
every region of the spectrum.
This is done by constructing a total bin density function
(#bins/energy width), as the maximum of the bin density of each
component. This total bin density function is integrated, and
the integral is divided into NFINAL equal parts; the points in
energy at which this division must be made are the returned bin
In the returned bin structure, the density of bins is everywhere
proportional to the densest bins in any of the GBi.
interpreted function, defined at i/multi.i line 1155
atten_emit= _multi_integrate(f, mesh, time, irays, slimits)
is the default drat_integrate routine.
On entry, file F is positioned at TIME, from which MESH has already
been read. IRAYS and SLIMITS are the rays coordinates (in internal
format) and integration limits.
The result should be ngroup-by-2-by-raydims, where the second index
is 1 for the attenuation factor, 2 for the self-emission (specific
intensity due to emission along the ray).
OPTIONS: drat_linear, drat_ocompute, drat_oadjust,
drat_emult, drat_amult, drat_omult, drat_nomilne,
drat_ekap, drat_akap, drat_glist
interpreted function, defined at i/multi.i line 540