Commit 1053d4f6 authored by Andrea Giannetti's avatar Andrea Giannetti
Browse files

Added bash command generator.

parent e617c145

etl/assets/commons/__init__.py

+7 −0
Original line number Diff line number Diff line 
import astropy.units

import yaml

import numpy as np

import os
@@ -106,3 +107,9 @@ def get_moldata(species_names: list): 
        data = urllib.request.urlopen(leiden_url_mapping[species]).read().decode()

        with open(os.path.join('mdl', 'radmc_input_files', f'molecule_{species}.inp'), 'w') as outfile:

            outfile.writelines(data)





def convert_frequency_to_wavelength(frequency: astropy.units.Quantity,

                                    output_units: astropy.units.Unit):

    return frequency.to(output_units, equivalencies=u.spectral())

etl/mdl/config/config.yml

0 → 100644
+9 −0
Original line number Diff line number Diff line 
radmc:

    inclination: 0

    position_angle: 0

    iline: 2

    central_frequency: 230.538

    frequency_units: GHz

    width_kms: 30

    nchannels: 100

    npix: 200

etl/mdl/mdl_prepare_radmc_command.py

0 → 100644
+19 −0
Original line number Diff line number Diff line 
import os

from astropy import units as u

from assets.commons import (load_config_file,

                            convert_frequency_to_wavelength)





def main():

    config = load_config_file(os.path.join('mdl', 'config', 'config.yml'))

    central_frequency = convert_frequency_to_wavelength(

        frequency=float(config['radmc']['central_frequency']) * u.Unit(config['radmc']['frequency_units']),

        output_units=u.Unit("micron"))

    with open('radmc3d_postprocessing.sh', 'w') as outfile:

        outfile.write(f'radmc3d image incl {config["radmc"]["inclination"]} phi {config["radmc"]["position_angle"]}'

                      f' iline {config["radmc"]["iline"]} widthkms {config["radmc"]["width_kms"]} '

                      f'linenlam {config["radmc"]["nchannels"]} npix {config["radmc"]["npix"]}')





if __name__ == '__main__':

    main()

etl/mdl/mdl_temperature_simulation.py

deleted100644 → 0
+0 −0

Empty file deleted.

etl/stg/config/config.yml

+0 −6
Original line number Diff line number Diff line
@@ -21,9 +21,3 @@ radmc: 
    scattering_mode_max: 0

    iranfreqmode: 1

    tgas_eq_tdust: 1
 
    inclination: 0

    position_angle: 0

    central_frequency: 230.538

    frequency_units: GHz

    width_kms: 30

    nchannels: 100