I am trying to reproduce the results in 1202.4039. I tried the webrun and the match between the two are very well. However the code I compiled on my machine generate very different results. The model I am testing is SNR_z4kpc_R20kpc_Ts150K_EBV5mag, I downloaded the galdef file from galprop website. Because I cannot find rbands_co10mm_v2_2001_qdeg.fits.gz, I use rbands_co10mm_v3_2001_hdeg.fits.gz instead, this is the only change I made. You can see from the log files I uploaded there is no other obvious difference between the webrun and my own run. But pi0 skymap is very different. In my own run there seems to be very few structure/contribution from the high latitude. Did I miss something obvious? I cannot do further testing since I do not have the v2 COR file.

The problem is that the two gas maps you are using are not compatible. The co map is split into 9 annuli while the hi is split into 17. There is unfortunately no check for this inconsistency. If you really need to reproduce the numbers of the papers locally you could merge some of the annuli in the HI map to make them consistent. Otherwise running with the HI map provided with the code is advised.

_________________
Gudlaugur Johannesson, GALPROP developer

Thanks a lot for the explanation, I tested with rbands_hi12_v5_hdeg_zmax1_Ts125.fits.gz and now the output of galprop looks a lot better. I do not know the differences between this and the HI files in the addon package, since the galprop-54_data_111612_addon.tar.gz is released much later I assumed they are more "correct"/better. Is there any plan to release a qdeg CO map to be used with them, or how can I get one? Otherwise the HI maps from the addon package are not really useful.

The new maps are better in the sense that they include the dark gas component derived as the residual dust map after linearly fitting the hi and co maps to the SFD E(B-V) dust map. The problem is that the dark gas correction is then assumed to be distributed similarly to the HI map and used to scale the HI annuli up or down, depending on the residuals.

The maps were released because we published the results in a paper. The co maps are identical apart from the resolution degradation (which is minor) and the different annuli (which is major). The release of the co map has not been discussed, but it makes little sense to have incompatible maps in the data package. I will bring this up with the other developers.

_________________
Gudlaugur Johannesson, GALPROP developer

Sorry for the late reply. The CO map used in 1202.4039 is now available in the download section.

_________________
Gudlaugur Johannesson, GALPROP developer

That is great news. Thank you very much.