类hitran数据库物种的非平衡辐射解算器
radis的Python项目详细描述
RADIS
非平衡红外发射与吸收逐行编码 后处理库比较实验光谱和计算光谱。
用户指南、安装过程和示例可在RADIS Website:
开始
安装
假设您安装了python Anaconda分发仅使用:
pip install radis
就这样!现在可以运行下面的第一个示例。如果你 遇到任何问题,或稍后要升级包,请参阅 那detailed installation procedure 是的。
快速启动
使用 calc_spectrum 功能。输出是一个Spectrum object: :
from radis import calc_spectrum
s = calc_spectrum(1900, 2300, # cm-1
molecule='CO',
isotope='1,2,3',
pressure=1.01325, # bar
Tgas=700, # K
mole_fraction=0.1,
path_length=1, # cm
)
s.apply_slit(0.5, 'nm') # simulate an experimental slit
s.plot('radiance')
从hitran数据库(on)计算co非平衡谱 你的第一个电话,这将计算并缓存co(x)rovibrational 能量)::
s2 = calc_spectrum(1900, 2300, # cm-1
molecule='CO',
isotope='1,2,3',
pressure=1.01325, # bar
Tvib=700, # K
Trot=300, # K
mole_fraction=0.1,
path_length=1, # cm
)
s2.apply_slit(0.5, 'nm')
s2.plot('radiance', nfig='same') # compare with previous
快速启动示例自动从
HITRAN-2016,
对于低于700 K的温度有效。对于高温
情况下,可能需要使用其他行数据库,如
HITEMP-2010
(通常T<;2000 K)或
CDSD-4000
(5000千吨)。必须在~/.radis
中描述这些数据库
Configuration
file。
更复杂 examples 需要使用 SpectrumFactory 类,这是Radis逐行计算的核心。 calc_spectrum 是一个包装 SpectrumFactory 对于简单的情况。
可以使用 experimental_spectrum 函数并与 plot_diff 功能。例如:
from numpy import loadtxt
from radis import experimental_spectrum, plot_diff
w, I = loadtxt('my_file.txt').T # assuming 2 columns
sexp = experimental_spectrum(w, I, Iunit='mW/cm2/sr/nm')
plot_diff(sexp, s) # comparing with previously spectrum 's' calculated previously
典型输出 plot_diff:
请参阅 Examples 节获取更多示例,并转到Spectrum page 更多的后处理功能。
在浏览器中(不需要安装!)
或者,您也可以使用 RADIS Interactive Examples 项目。例如,在下面的链接上运行快速入门示例:
或者开始一个光秃秃的萝卜在线会话:
开发人员指南
架构
萝卜的内部结构在Developer Guide:
许可证
此代码位于GNU LESSER GENERAL PUBLIC LICENSE (v3)下的此存储库中
支持
如果遇到任何问题,请首先参考已知列表 Issues 在Github上。感谢您的反馈和建议!
如有任何问题,请加入Gitter上的讨论频道:
参考文献
链接
萝卜:
而且:
PYPI存储库:https://pypi.org/project/radis/