Skip to content

Commit

Permalink
Partly port old Lookup calculations (#854)
Browse files Browse the repository at this point in the history
  • Loading branch information
riclarsson authored Nov 6, 2024
2 parents ebce23a + 5dd0bf2 commit 54de28f
Show file tree
Hide file tree
Showing 51 changed files with 2,884 additions and 609 deletions.
5 changes: 3 additions & 2 deletions doc/arts/guide.theory.absorption.rst
Original file line number Diff line number Diff line change
Expand Up @@ -151,7 +151,7 @@ and the added emissions are given by
.. math::
K_{NLTE} = \rho \frac{c^2\nu}{8\pi} \left\{r_u\left[
1 - \exp\left(-\frac{h\nu_0}{kT}\right)\right] - \left(r_l \frac{g_u}{g_l} - r_u\right)
1 - \exp\left(\frac{h\nu_0}{kT}\right)\right] - \left(r_l \frac{g_u}{g_l} - r_u\right)
\right\} \frac{ A_{lu}}{\nu_0^3},
where :math:`r_l` and :math:`r_u` are the ratios of the populations of the lower and upper states, respectively.
Expand Down Expand Up @@ -196,12 +196,13 @@ This is seen by putting the above RHS and the expression for :math:`r_u` into th
.. math::
K_{NLTE} = \rho \frac{c^2\nu}{8\pi} \left\{\frac{g_u\exp\left(-\frac{E_u}{kT}\right)}{Q(T)}\left[
1 - \exp\left(-\frac{h\nu_0}{kT}\right)\right] - \left[1 - \exp\left(-\frac{h\nu_0}{kT}\right)\right]\frac{g_u\exp\left(-\frac{E_l}{kT}\right)}{Q(T)}
1 - \exp\left(\frac{h\nu_0}{kT}\right)\right] - \left[1 - \exp\left(-\frac{h\nu_0}{kT}\right)\right]\frac{g_u\exp\left(-\frac{E_l}{kT}\right)}{Q(T)}
\right\} \frac{ A_{lu}}{\nu_0^3} = 0.
The ratio between LTE and non-LTE line strength remaining is:

.. math::
\frac{S_{NLTE}}{S_{LTE}} = \frac{1 - \exp\left(-\frac{h\nu_0}{kT}\right)}{1 - \exp\left(-\frac{h\nu}{kT}\right)}.
It is clear that the non-LTE expression is the one that is incorrect here.
Expand Down
7 changes: 4 additions & 3 deletions examples/recipes/AtmosphericFlux/atmospheric_flux.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -37,7 +37,8 @@
"metadata": {},
"outputs": [],
"source": [
"fop = pyarts.recipe.AtmosphericFlux()"
"fop = pyarts.recipe.AtmosphericFlux(species=[\"H2O-161\", \"O2-66\", \"N2-44\", \"CO2-626\", \"O3-XFIT\"],\n",
" remove_lines_percentile={\"H2O\": 70})"
]
},
{
Expand Down Expand Up @@ -120,7 +121,7 @@
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
Expand All @@ -134,7 +135,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.6"
"version": "3.12.7"
}
},
"nbformat": 4,
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,8 @@
"metadata": {},
"outputs": [],
"source": [
"fop = pyarts.recipe.SpectralAtmosphericFlux()"
"fop = pyarts.recipe.SpectralAtmosphericFlux(species=[\"H2O-161\", \"O2-66\", \"N2-44\", \"CO2-626\", \"O3-XFIT\"],\n",
" remove_lines_percentile={\"H2O\": 70})"
]
},
{
Expand Down Expand Up @@ -130,7 +131,7 @@
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
Expand All @@ -144,7 +145,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.6"
"version": "3.12.7"
}
},
"nbformat": 4,
Expand Down
Loading

0 comments on commit 54de28f

Please sign in to comment.