It is currently Tue Nov 12, 2024 7:21 am

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 5 posts ] 
Author Message
 Post subject: galprop gas files
PostPosted: Thu May 24, 2007 6:27 am 
Offline

Joined: Tue Feb 20, 2007 1:10 am
Posts: 5
Location: Padova
Which quantity is described in the input file rbands_co4.fits? Is it the number density of hydrogen MOLECULES or the number density of PROTONS in form of molecules?


Top
 Profile  
 
 Post subject: rbands co
PostPosted: Fri May 25, 2007 12:14 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
it's the integrated Tco dv, in units of K km s-1
which has to be converted to molecular hydrogen
column density using the famous Xco factor (e.g. as given in the galdef file).

_________________
Andy Strong, MPE


Top
 Profile  
 
 Post subject: galprop gas files
PostPosted: Fri May 25, 2007 12:36 am 
Offline

Joined: Tue Feb 20, 2007 1:10 am
Posts: 5
Location: Padova
I found some misleading information, both in the file header, where you wrote as unit cm^-2, and in the GALPROP v.50 manual (September 6, 2006), where (page 4 bottom) you can read "the H2 gas number density is defined in the form of table ... and the conversion factor is taken as X= ..."


Top
 Profile  
 
 Post subject: H2 and rbands
PostPosted: Fri May 25, 2007 12:58 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
there are two things:
1 the galprop internally modelled H_2 which is indeed as described in the manual on page 4, using a fixed Xco=1.9E20. This is used for cosmic-ray propagation but not for gamma-rays

2 the Wco rings in rbands, which are used to compute the gamma-rays.
You are apparently right about the error in the header, I informed Seth Digel !


Clearly 1 & 2 are not consistent: we know this and are working on a new version where they will be consistent. However the effect is not great.
But consistency is very desirable of course.

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Sat May 26, 2007 10:30 pm 
Offline

Joined: Wed May 17, 2006 11:30 pm
Posts: 2
Location: SLAC
Luigi is correct that the units are incorrectly stated in the header of rbands_co4.fits. They are K km s^-1, as Andy says. It was just a mistake to list the units as '10E20 cm-2'. This came about because I used the same BUNIT string for the H I rings (rbands_hi7.fits).

10E20 cm-2 now looks ambiguous (or wrong) to me for the H I rings. I meant 10^20 cm^-2. I'll fix the headers of both files in their next update.

Thanks, Luigi, for pointing out the error.

_________________
Seth Digel


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 5 posts ] 

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 forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group