It is currently Tue Aug 06, 2024 1:50 am

 All times are UTC - 8 hours [ DST ]

 Page 1 of 1 [ 1 post ]
 Print view Previous topic | Next topic
Author Message
 Post subject: CR from DM annihilation with a cuspy DM profile.Posted: Fri Sep 25, 2009 4:38 am

Joined: Tue Aug 11, 2009 2:40 pm
Posts: 2
Location: Seoul National University
I used galprop v.50p for my recent paper about CR from decaying DM.

While I used it, I found some bug in the code averaging DM density.
( "DM_profile_av()" in "gen_DM_source.cc")
The calculation with 3-d coordinate or non-cuspy profile is OK. But
the calculation with cuspy profile like NFW in 2-d coordinate make
unacceptably large CR fluxes. It maybe the problem in averaging DM
near the origin (Jacobi factor become significant in this region.).

By simply changing "int nuse" into "double ruse" and "nuse++" into "ruse+=rr",
code produce appropriate result.

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
double Galprop::DM_profile_av(double r,double z,double dr,double dz,double dzz)
{
double DM_profile_av_=0.0;
int nuse=0;
double ruse=0;

for (double zz=z-dz/2.; zz<=z+dz/2.; zz+=dzz)
for (double rr=r-dr/2.; rr<=r+dr/2.; rr+=dr/10.)
{
if (rr<0.) continue;
DM_profile_av_+=rr*DM_profile(rr,0,zz);
ruse+=rr;
}
return (DM_profile_av_/ruse);
}
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

_________________
Ji-Haeng Huh (jhhuh@phya.snu.ac.kr)
High Energy Phenomenology Group
Dept. of Physics/Astronomy
Seoul National University, Seoul, Korea

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 1 post ]

 All times are UTC - 8 hours [ DST ]

#### Who is online

Users browsing this forum: No registered users and 1 guest

 You cannot post new topics in this forumYou cannot reply to topics in this forumYou cannot edit your posts in this forumYou cannot delete your posts in this forumYou cannot post attachments in this forum

Search for:
 Jump to:  Select a forum ------------------ general astrophysics    cosmic rays    gamma rays    interstellar medium    heliospheric modulation    cosmic-ray data galprop general discussion    galprop applications    galprop compilation    galprop debugging    galprop tests    galprop developer    galprop website    galprop features    galprop help    galprop WebRun    2011 GALPROP Workshop: abstracts and discussion other codes    analytical models    semi-analytical models    numerical models    comparisons with galprop    heliospheric modulation