or perhaps experimental_spectrum would not be needed at all (other than legacy perhaps)
That could be an Issue too, btw!
Did you also see this: https://radis.readthedocs.io/en/latest/spectrum/spectrum.html#initialize-from-text-file-or-arrays ?
apply_slit(..norm_by="max") method. When it does so, the code tries to process a string "I/I0" as a unit. I have been looking all over the pint documentation but I still haven't been able to identify how this worked. Any leads on how this worked would be appreciated!
or perhaps experimental_spectrum would not be needed at all (other than legacy perhaps)
That could be an Issue too, btw!
Did you also see this: https://radis.readthedocs.io/en/latest/spectrum/spectrum.html#initialize-from-text-file-or-arrays ?
Not yet, looks useful!
I/IO, -ln(I/IO) and eps in conjunction with Astropy? I have been at it for quite some time and it appears simply changing all occurrences of these strings to 1 won't resolve the issue. I have been forced to add a lot of if-else blocks in a lot of places to handle issues arising from the way astropy parses strings. For instance, while we could simply multiply two units represented as strings in pint by concatenating them with a * in between, in Astropy we need to first check whether both the strings are pure units or is one of them 1 (like in our case, following which we need to ensure that this string appears before the other string, nm in this case).
test_slit() locally on my laptop. The traceback seems to be hinting at something going wrong with matplotlib but I still not clear about it.
Also, regarding the many print() comments in your code: I usually debug with the "%debug" feature of ipython, going through the different functions with "u" [up] and "d" [down], and evaluating the variables there. Of course it's up to you!
This message is replying to a Matrix event but we were unable to find associated bridged Gitter message to put it in the appropriate threaded conversation.
Oh! cool stuff
PR https://github.com/radis/radis/pull/271|#271 by Anand Kumar was just merged on the develop branch. It introduces an automatic scaling of the wavenumber grid, which adapts to the pressure, temperature and species considered. It is the best trade-off between accuracy and performance for a single spectrum, so I encourage you to try it. wstep='auto'
Will be available in the public 0.9.30 version incoming next week. And may become the default parameter in future versions if feedbacks are good.
4 replies
How do you display the
linesurvey module in a normal python code? As opposed to a jupyter notebookfrom radis import calc_spectrum
s = calc_spectrum(
wavenum_min=2007,
wavenum_max=2009,
mole_fraction=0.1,
path_length=10, # cm
Tgas=700,
molecule="CO",
isotope=[1,2,3],
databank="hitran",
export_lines=True,
)
s.line_survey(overlay="absorbance", barwidth=0.01)
s.plot('absorbance')
1 reply
• direct link to https://radis.app
• update links to https://radis.readthedocs.io/en/latest/auto_examples/index.html|Examples, and https://radis.github.io/radis-lab|Radis-Lab (online with pre-installed databases)
• mention the <https://www.sciencedirect.com/science/article/pii/S0022407320310049|LDM article>
• update features with link to https://user-images.githubusercontent.com/16088743/121789972-f34ab280-cbda-11eb-99b0-d7ba019d80d5.mp4|GPU-video ; mention of HITEMP auto download
Feedback welcome !
1 reply