Hi all,

I've used the

yt.analysis_modules.absorption_spectrum.absorption_spectrum_fit

script to fit an Lya absorption spectrum generated by yt.

But there are some abnormal peaks when I did statistics of line width b and column density N. Two figures are attached. I've tried with different redshifts and all of them had peaks in b-distribution at init_b (which is 40km/s here) and maxb(which is 300km/s here). For the N-distribution I got peaks around 10^12 cm^-2 at every redshift. The parameters I used were:

HI_parameters = {'name':'HI',

                 'f': [.4164],

                 'Gamma':[6.265E8],

                 'wavelength':[1215.67],

                 'numLines':1,

                 'maxN': 1E22, 'minN':1E11,

                 'maxb': 300, 'minb':1,

                 'maxz': 5, 'minz':0,

                 'init_b':40,

                 'init_N':1E14}

Could someone please tell me how to solve this? All the analysis were done with the latest development version of yt and I made no changes to the source code. 

Please tell me if you need any datasets.

Thank you!

Pengfei

_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org