karllark / DIRTY_dustrt

DIRTY: 3D dust radiative transfer for dusty astrophysical sources
https://dirty-dustrt.readthedocs.io
5 stars 1 forks source link

Fix for missing rad field unc and adding nphotons in rad field calculation #33

Closed karllark closed 5 months ago

karllark commented 5 months ago

The radiation field uncertainties were missing from single_wavelength runs. This was due to the sed_lum being set to 1.0, a very small value. Hence the radiation field was below 1e-36, which is the check to make sure the radiation field can be stored as a float. The solution is to set the radiation field to 1e45 and this provides radiation field values that are well with the range that can be stored in a float.

In addition, the number of photons used to compute the radiation field (e.g, the number of photons crossing each cell) is now output.