There is a bug in latest versions of GALPROP in file gen_DM_source.cc:
Row 136, we have:
Code:
pow(DM_profile_av(galaxy.x[ix], galaxy.y[ix], galaxy.z[iz], galaxy.dx, galaxy.dy, galaxy.dz, dzz)/DMmass,2)
This coordinate:
Code:
galaxy.y[ix]
is wrong.
Evidently, it must be:
Code:
pow(DM_profile_av(galaxy.x[ix], galaxy.y[iy], galaxy.z[iz], galaxy.dx, galaxy.dy, galaxy.dz, dzz)/DMmass,2)
If someone have used this file, please fix a bug and check your results again.
File with bug fixed is attached.
| Attachments: |
gen_DM_source.cc [10.27 KiB]
Downloaded 1551 times |