Skip to content

New ions formulation, a new version of rascas with arbitrary number of scatterers/lines

Leo Michel-Dansac requested to merge merge_new_ions into develop

New version of RASCAS with arbitrary number of scatterers/lines, parametrised with files in ./rascas/ions_parameters/

This new feature has been developed by Valentin Mauerhofer (see his branch new_ions). It allows to set up the gas composition mix in the parameter file, using for instance:

nscatterer = 2

scatterer_names = HI-1216, D-1215

There is only one module_gas_composition. It is not necessary to recompile the codes when changing the gas composition mix. WARNING: old domain dumps and mesh files are not usable anymore. They have to be generated again.

The following developments have then been done by Jeremy Blaizot:

  • made core-skipping work again
  • replaced Lya collisional emission from Goerdt by the new computation from Harley in module_ramses.f90:ramses_get_LyaEmiss_HIDopwidth
  • added option spec_type=Table_lowmem in PhotonsFromStars to process large numbers of particles without using too much memory. Also a little bit of openmp.
  • minor optimisation in module_ramses:read_stars_in_domain.

Changelog:

  • bugfix of the sampling method of stars SEDs in PhotonsFromStars, with option spec_type=Table.
  • update PhotonsFromSourceModel according to changes in PhotonsFromStars
  • update values of the physical constants in module_constants
  • add many ions parameters
  • update examples of parameter files
  • update python directory
  • add a new parameter in module_gas_composition to run RT without dust.
  • encode the dust model string into an integer key to improve the performance of the case statement.
  • debug compilation error with OpenMP in subroutine ramses_read_stars_in_domain (module_ramses.f90)
Edited by Leo Michel-Dansac

Merge request reports