fix issue in which H2 collider files rejected themselves because RADEX internally converts H2->OH2+PH2

keflavich / pyradex

Python interface to RADEX

BSD 3-Clause "New" or "Revised" License

19 stars 12 forks source link

fix issue in which H2 collider files rejected themselves because RADEX internally converts H2->OH2+PH2 #19

Open keflavich opened 6 years ago

keflavich commented 6 years ago

Currently fails with:

In [3]: pyradex.Radex(species='c2h_h2_e', collider_densities={'H2': 1e3, 'E':1e-3}, abundance=1e-8, temperature=30)()
 *** Warning: Assuming thermal o/p ratio for H2 of    3.0517887472151187E-002
 *** Warning: Assuming thermal o/p ratio for H2 of    3.0517887472151187E-002
 *** Warning: No rates found for any collision partner