ExoMolHR compatibility

[HajimeKawahara]

With the implementation of ExoMolHR, users can start calculations in the required wavelength range without the need to download the large raw molecular datasets. This improvement should enhance the overall usability of ExoJAX.

https://arxiv.org/abs/2504.08731

[HajimeKawahara]

working on exomolhr branch

[HajimeKawahara]

exojax/tests/endtoend/identify_lines

from exojax.spec.exomolhr import XdbExomolHR
from exojax.spec.exomolhr import list_exomolhr_molecules
from exojax.spec.exomolhr import list_isotopologues
from exojax.utils.grids import wavenumber_grid
import matplotlib.pyplot as plt

nus, _, _ = wavenumber_grid(22800.0, 23600.0, 10, xsmode="premodit", unit="AA")
temperature = 1300.0
molecules = list_exomolhr_molecules()
iso_dict = list_isotopologues(molecules)


lslist = ["-", "--", "-.", ":", "-", "--", "-.", ":"]
lwlist = [1.0, 1.0, 1.0, 1.0, 2.0, 2.0, 2.0, 2.0]
markers_list = [".", "o", "s", "D", "^", "v", "<", ">"]
k=0
for molecule in iso_dict:
    isos = iso_dict[molecule]
    for j, iso in enumerate(isos):
        try:
            mdb = XdbExomolHR(iso, nus, temperature)
            plt.plot(
                1.e8/mdb.nu_lines, mdb.line_strength, markers_list[j], label=iso, ls=lslist[j], lw=lwlist[j]
            )
        except:
            k=k+1
print(k)
plt.yscale("log")
plt.xlabel("wavelength (AA)")
plt.ylabel("Line Strength (cm/molecule)")
plt.legend()
plt.savefig("exomolhr_lines.png")
plt.show()

[HajimeKawahara]

close with #600