This code solves essentially the same problem as RADEX written by Van der Tak, except that we take a different approach to solve the statistical equilibrium problem. Given an initial distribution, what we do is to evolve the system towards equilibrium using an ODE solver.

Usage

To use this code, you first need to compile it using the makefile (the executable is named my_radex) by running make in the command line, then edit the configuration file (configure.dat) to meet your needs. After this execute the command

./my_radex configure.dat

and you will get the results.

A python wrapper is also included. To make the wrapper, run make wrapper in the command line. For its usage, see this ipython notebook. The wrapper is preliminary; not all the functionalities in the Fortran source code are included in the wrapper (though usually they are not needed).

We use the LAMDA format for the input energy levels and transition rates.

This code has not been thoroughly tested, and it is possible that some input files won't load properly.

Note

• Works with python 3.7