Skip to content

Commit

Permalink
continuing to make sphinx changes
Browse files Browse the repository at this point in the history
  • Loading branch information
jhmatthews committed Oct 25, 2024
1 parent f84b1d3 commit 6dc9928
Show file tree
Hide file tree
Showing 56 changed files with 69 additions and 783 deletions.
2 changes: 1 addition & 1 deletion docs/sphinx/requirements.txt
Original file line number Diff line number Diff line change
Expand Up @@ -12,4 +12,4 @@ sqlalchemy
jinja2==3.1.4
sphinx_gallery
pyhdf
sphinx-autoapi==3.1.1
sphinx-autoapi==3.1.1
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/bound_bound.rst
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@ We have also used a line list from Verner in the past


Translation to SIROCCO format
============================
======================================
There are several steps to creating the data used in SIROCCO from that in gfall.dat, that are carried out by py_read_kurucz and py_link. The first routine reads the gfall.dat file and creates two output files, a file containing the lines and the associated such as the effective oscillatory strength and a file which contains information about the ion levels. py_read_kurucz chooses only a portion of the Kurucz lines, namely those associated with ions with ionization potentials in a certain range and lines with gf factors exceeding a certain value. The second program py_link attempts to create a model ion with links between the levels and the ions. Both of these routines are driven by .pf files, similar to what are used in SIROCCO. Examples of the .pf files are in the directory py_kurucz

In practice we have not used these data for any SIROCCO publications. At some point early in the AGN project, NSH increased the number of lines, and generated lines\_linked\_ver\_2.dat and levels\_ver\_2.dat. I think this was because there was a small bug which meant the oscillator strength cut that was stated was not that which was applied.
Expand Down Expand Up @@ -69,7 +69,7 @@ one needs to index everything self-consistentl


SIROCCO structure
================
================================

Line data is stored in the data structure **lines**

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,7 @@ These values of :math:`\Upsilon` simply replace :math:`\Omega`.
In the asbsence of data in this format, the Van Regemorter approximation is utilized.

Translation to SIROCCO format
============================
======================================

It is necessary to link each line in our line list with the relevant electron collision strength. This is achieved using the python script "coll_stren_lookup.py" which first reads in the "lines_linked_ver_2.py" line list, then attempts to work out which lines are which by comparing the energy and the oscillator strength of the line. If these match to within a factor of 10% then the code logs this as a possible match. If better matches come along, then the code adopts those instead.

Expand Down Expand Up @@ -67,7 +67,7 @@ and
So, to get :math:`\Upsilon` for a given T, one converts T to x via the correct equation, then linearly interpolate between values of :math:`y(x)`, then convert back to :math:`\Upsilon`.

SIROCCO structure
================
==========================

The data is stored in SIROCCO in the Coll\_stren structure which has memebers

Expand Down
6 changes: 3 additions & 3 deletions docs/sphinx/source/atomic/bound_free_topbase.rst
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ Obtained from The `Opacity project <http://cdsweb.u-strasbg.fr/topbase/topbase.h


Translation to SIROCCO format
============================
======================================

ksl - It's not clear that we are now making use of the topbase data in this way but my original attempt to incorporate topbase photoinization data into SIROCCO is contained in the directory topbase. Processing of these files was done by py_top_phot. My feeling is that we can replace these remarks with those that are more up to date, once Nick and James discuss this section, and delete any mention of my original attempt on this in the data-gen archive.

Expand Down Expand Up @@ -64,12 +64,12 @@ For the macro atom case, the indices relate to the energy levels, that is these
The TopBase energies are inaccurate and so generally adjustments are made using Chianti or some other source to fix up the energy levels.

SIROCCO structure
================
==========================

The data are stored in the Topbase_phot stucture which can be found in atomic.h

Criteria for usage in SIROCCO run
================================
==========================================

Data has to be read into SIROCCO in a logical order. For a set of phototionization x-sections to be accepted, the energy levels (or configuratios) must already have been defiend. See :doc:`levels`

Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/bound_free_verner.rst
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@ There are three sources for this data


Translation to SIROCCO format
============================
======================================

**Tabulation Process**

Expand Down Expand Up @@ -62,7 +62,7 @@ This data is linked to the relevant ion via z and state, islp and level are not
This data is linked to the relevant ion via z and state. the n_shell and l_subshell numbers are used to cross reference to the electron yield records. As above, the last record shows how many points are in the fit, and the data pairs making up the fit follow with the keyword InnerVY.

SIROCCO structure
================
==========================

Where the data is stored internally in SIROCCO

Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/direct_ionization.rst
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ The data comes directly from `Dere 2006, A&A, 466, 771 <https://www.aanda.org/ar


Translation to SIROCCO format
============================
======================================


The data table is downloaded in its entirety from the data table associated with the paper. All that happens is that the table is saved to a text file, and the keyword DI_DERE is just prepended to each row.
Expand Down Expand Up @@ -41,7 +41,7 @@ The rate coefficient R(T) is recovered from the scaled rate coefficient in the t
where :math:`E_{1}` is the first exponential integral. In python we use the gsl_sf_expint_E1 routine in gsl.

SIROCCO structure
================
==========================

This data is stored in the dere_di_rate structure with members

Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/electron_yields.rst
Original file line number Diff line number Diff line change
Expand Up @@ -10,7 +10,7 @@ Source
This data comes from `Kaastra and Mewe 1993, A&A, 97, 443 <http://articles.adsabs.harvard.edu/full/1993A%26AS...97..443K>`_ . The data is downloaded from the vizier site linked and put into a file called "electron_yield.data"

Translation to SIROCCO format
============================
======================================

The translation takes place using the python script "kaastra_2_py.py" which takes the saved raw data file "electron_yield.data" and compares it line by line to the inner shell cross section data in "vy_innershell_tab.data"(see above). The n shell and l subshell to which each record applies is coded in the KM data and needs to be decoded. This is what the script does, and all the script then does is output the yield data into a new file "kaastra_electron_yield.data" which contains the n and l cross reference.

Expand All @@ -36,7 +36,7 @@ The data is linked to the correct inner shell photoionization cross section (and


SIROCCO structure
================
==========================

The data is stored in SIROCCO in the inner_elec_yield structure which contains

Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/elem_ions.rst
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ This data comes from `Verner, Barthel & Tytler, 1994, ApJ 108, 287. <http://arti


Translation to SIROCCO:
======================
================================

The original data and the translation can be found in py\_verner. A simple awkscript converts the downloaded data to SIROCCO format.

Expand Down Expand Up @@ -93,7 +93,7 @@ Note that only evident changed is the label, but in this case the number of nlte


SIROCCO structure:
=================
===========================
This data is held in SIROCCO in various fields in structures **elements** and **ions**.

Comments:
Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/free-free.rst
Original file line number Diff line number Diff line change
Expand Up @@ -15,7 +15,7 @@ The free-free Gaunt factors are taken from `Sutherland 1998, MNRAS, 300, 321. <
The last file is the one we use to calculate free-free emission, since this in integrated gaunt factor over a population of electrons with a Boltzmann distribution of velocities. The other two files could be of use in the future should we wish to have gaunt factor corrections for the heating rates,in which case we should use the gffgu.dat data file. However generally speaking free-free heating is never important and there would be significant overhead in calculating a gaunt factor for each photon.

Translation to SIROCCO
=====================
===============================
The file is simply modified by hand to put a label "FF\_GAUNT" at the start of each data line and a hash at the start of each comment line.

Datafile - gffint.dat:
Expand Down Expand Up @@ -49,7 +49,7 @@ where
:math:`\Delta=\log(x)-\log(\gamma^2)`

SIROCCO structure
================
==========================
This data is held internally in SIROCCO in the structure **gaunt_total** which has members

- log_gsqrd
Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/levels.rst
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,7 @@ Level information can be derived from a variety of sources, by:


Translation to SIROCCO:
======================
================================



Expand Down Expand Up @@ -94,7 +94,7 @@ Since they quote J, one converts to g = 2J+1
The ionization potential is not used, as it is redundant with the excitation energy which is, and the last column giving the configuration is also for information only.

SIROCCO structure:
=================
===========================
This data is held in SIROCCO in various fields in structure **config**.

Comments:
Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/atomic/photon_yields.rst
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ Source
This data comes from `Kaastra and Mewe 1993, A&A, 97, 443 <http://articles.adsabs.harvard.edu/full/1993A%26AS...97..443K>`_ . The data is downloaded from the vizier site linked and put into a file called "fluorescent\_yield.data"

Translation to SIROCCO format
============================
======================================

The translation takes place using the python script "kaastra_2_py.py". All identical to electron yield, but input file is "fluorescent_yield.data" and output is "kaastra_fluorescent_yield.data"

Expand All @@ -35,7 +35,7 @@ The data is linked to the correct inner shell photoionization cross section (and


SIROCCO structure
================
==========================

The data is stored in SIROCCO in the inner_fluor_yield structure which contains

Expand Down
4 changes: 2 additions & 2 deletions docs/sphinx/source/conf.py
Original file line number Diff line number Diff line change
Expand Up @@ -87,9 +87,9 @@
# built documents.
#
# The short X.Y version.
version = '88'
version = '0.1'
# The full version, including alpha/beta/rc tags.
release = '88a'
release = '0.1'

# The language for content autogenerated by Sphinx. Refer to documentation
# for a list of supported languages.
Expand Down
2 changes: 1 addition & 1 deletion docs/sphinx/source/developer/cuda.rst
Original file line number Diff line number Diff line change
Expand Up @@ -411,6 +411,6 @@ shown below.
# So to compile SIROCCO, we have something which looks vaguely like this. Note that
# we use the CUDA_OBJECTS recipe as a requirement for the sirocco recipe. This CUSOLVER_STATUS_SUCCESS
# the CUDA source to be compiled to object code *if* NVCC is defined
sirocco: startup sirocco.o $(sirocco_objects) $(CUDA_OBJECTS)
sirocco: startup sirocco.o $(sirocco_objects) $(CUDA_OBJECTS)
$(CC) $(CFLAGS) sirocco.o $(sirocco_objects) $(CUDA_OBJECTS) $(kpar_objects) $(LDFLAGS) -o sirocco
4 changes: 2 additions & 2 deletions docs/sphinx/source/py_progs.rst
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
SIROCCO Scripts
Python Scripts
################################

There are several SIROCCO scripts written to prepare input for
There are several python scripts written to prepare input for
and analyse the output of SIROCCO.

Some of the more useful scripts/modules are documented below.
Expand Down
19 changes: 1 addition & 18 deletions docs/sphinx/source/py_progs/MakeMacro/MacroCombine.rst
Original file line number Diff line number Diff line change
Expand Up @@ -4,11 +4,6 @@
.. automodule:: MacroCombine







.. rubric:: Functions

.. autosummary::
Expand All @@ -25,16 +20,4 @@
steer
write_phot
xguess














19 changes: 1 addition & 18 deletions docs/sphinx/source/py_progs/MakeMacro/MakeMacro.rst
Original file line number Diff line number Diff line change
Expand Up @@ -4,11 +4,6 @@
.. automodule:: MakeMacro







.. rubric:: Functions

.. autosummary::
Expand All @@ -23,16 +18,4 @@
make_phot
print_elvlc
write_phot














19 changes: 1 addition & 18 deletions docs/sphinx/source/py_progs/MakeMacro/RedoPhot.rst
Original file line number Diff line number Diff line change
Expand Up @@ -4,11 +4,6 @@
.. automodule:: RedoPhot







.. rubric:: Functions

.. autosummary::
Expand All @@ -19,16 +14,4 @@
redo_one
steer
write_phot_tab














19 changes: 1 addition & 18 deletions docs/sphinx/source/py_progs/development/CheckAtomic.rst
Original file line number Diff line number Diff line change
Expand Up @@ -4,11 +4,6 @@
.. automodule:: CheckAtomic







.. rubric:: Functions

.. autosummary::
Expand All @@ -20,16 +15,4 @@
get_level_tab
get_line_tab
xread














19 changes: 1 addition & 18 deletions docs/sphinx/source/py_progs/development/add_param.rst
Original file line number Diff line number Diff line change
Expand Up @@ -4,11 +4,6 @@
.. automodule:: add_param







.. rubric:: Functions

.. autosummary::
Expand All @@ -17,16 +12,4 @@
doit
read_file
read_table














Loading

0 comments on commit 6dc9928

Please sign in to comment.