Re: region_flux in ergs?

On Thu, Mar 20, 2003 at 15:09 +0100, Jorge Sanz wrote:
> 	Is there any function that may help to get the results of
> region_flux in ergs instead of photon units? I would like to calculate
> this for a wide wavelength range and the conversion photon-ergs should be
> done for each point...

Hi Jorge, 

In the current release, there is no equivalent to region_flux
which gives results in ergs.  For the next release, I'll
consider providing a mechanism to do this.

I've appended a simple function which will compute the energy
flux in a wavelength interval, but it doesn't provide the other
statistics computed by region_flux.

-John

-- 
John C. Houck                       MIT Center for Space Research
NE80-6005: 617-253-3849                   77 Massachusetts Avenue
42:21:55.105N, 71:05:28.122W             Cambridge, MA 02139-4307




define flux_ergs (y_lo, y_hi, x)
{
   % input
   %  y_lo, y_hi = wavelength interval
   %  x = struct {bin_lo, bin_hi, value}
   
   % bin-center wavelength, energy
   variable y = 0.5 * (x.bin_lo + x.bin_hi);
   variable e_erg = (Const_hc / y) * 1.e3 * Const_eV;
   
   % convert photons ==> ergs
   x.value *= e_erg;
   
   % integrate over [y_lo, y_hi]
   return rebin (y_lo, y_hi, x.bin_lo, x.bin_hi, x.value);
}
----
You received this message because you are
subscribed to the isis-users list.
To unsubscribe, send a message to
isis-users-request_at_email.domain.hiddenwith the first line of the message as:
unsubscribe