Code related to "Irregular grids for 3D NLTE radiative transfer in stellar atmospheres" (https://arxiv.org/abs/2306.01041), which is based on my master thesis. I simulate radiation in the solar atmosphere on an irregular grid using a Voronoi tesselation.
The voronoi tesselation is calculated with the open source voro++ library. C files for computing and outputting grid statistics for the Voronoi tesselation are located in the rt_preprocessing/ folder.
The radiative transfer simulation is written in julia, using its native thread parallelism to speed up computations. Code for the RT simulations are in the src/ directory. Results are written to file in an HDF5 format. Analysing results and making plots is done with python, located in the python/ folder.
To run the code, an atmospheric model is needed. The model has to be stored in an HDF5 format. For my calculations, I used a Bifrost model astmosphere, provided by the Institute of Theoretical Astrophysics at the University of Oslo.
Rasiative transer on irregular grids is performed and tested against regular grids with:
- searchlight beam test with the
compare_searchlight.jlscript - Continuum calculations in
compare_continuum.jl - NLTE spectral synthesis are performed in
compare_line.jl
Data are written to file with plot_utils.jl, and plotted/analysed with python. Figures generated from the results are located in the image folder img/.
Any questions regarding the repository can be directed to: [email protected]