Commit 84a56d9b authored by Andrea Giannetti's avatar Andrea Giannetti
Browse files

Mapped line modes.

parent 02ed6d93

etl/assets/constants/__init__.py

+7 −0
Original line number Diff line number Diff line
@@ -39,3 +39,10 @@ leiden_url_mapping = { 
    'e-ch3oh': 'https://home.strw.leidenuniv.nl/~moldata/datafiles/e-ch3oh.dat',

    'a-ch3oh': 'https://home.strw.leidenuniv.nl/~moldata/datafiles/ch3oh_a.dat'

}



radmc_lines_mode_mapping = {

    'lte': 1,

    'user_defined_populations': 2,

    'lvg': 3,

    'optically_thin_non_lte': 4,

}

etl/stg/config/config.yml

+1 −1
Original line number Diff line number Diff line
@@ -23,7 +23,7 @@ lines: 
                              "e-ch3oh": 1e-8,

                              "p-h2": 0.25,

    }

    lines_mode: 3

    lines_mode: 'lvg'

    collision_partners: ['p-h2']



radmc:

etl/stg/stg_radmc_input_generator.py

+4 −3
Original line number Diff line number Diff line
@@ -8,7 +8,8 @@ from assets.commons import (compute_power_law_radial_profile, 
                            validate_parameter,

                            setup_logger)

from assets.constants import (mean_molecular_mass,

                              radmc_input_headers)

                              radmc_input_headers,

                              radmc_lines_mode_mapping)

from astropy import units as u

from astropy import constants as cst
@@ -46,7 +47,7 @@ def write_radmc_lines_input(line_config: dict, 
                            path: Union[None, str] = None):

    _path = validate_parameter(path, default='.')



    if line_config['lines_mode'] == 1:

    if line_config['lines_mode'] == 'lte':

        assert len(line_config['collision_partners']) == 0

    else:

        assert len(line_config['collision_partners']) != 0
@@ -212,7 +213,7 @@ def write_radmc_main_input_file(config: dict, 
        outfile.write(f'scattering_mode_max = {config["radmc"]["scattering_mode_max"]}\n')

        outfile.write(f'iranfreqmode = {config["radmc"]["iranfreqmode"]}\n')

        outfile.write(f'tgas_eq_tdust = {config["radmc"]["tgas_eq_tdust"]}\n')

        outfile.write(f'lines_mode = {config["lines"]["lines_mode"]}\n')

        outfile.write(f'lines_mode = {radmc_lines_mode_mapping[config["lines"]["lines_mode"]]}\n')





def main(override_config: Union[dict, None] = None):