Commit c00689d3 authored by Andrea Giannetti's avatar Andrea Giannetti
Browse files

Renamed RADMC files folder.

parent fabc232d
Loading
Loading
Loading
Loading
+1 −1
Original line number Diff line number Diff line
@@ -10,7 +10,7 @@ def main():
        frequency=float(config['radmc']['central_frequency']) * u.Unit(config['radmc']['frequency_units']),
        output_units=u.Unit("micron"))
    with open(os.path.join('mdl', 'radmc3d_postprocessing.sh'), 'w') as outfile:
        outfile.write('cd radmc_input_files\n')
        outfile.write('cd radmc_files\n')
        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"]}')
+4 −4
Original line number Diff line number Diff line
@@ -131,7 +131,7 @@ def write_grid_input_files(profiles, grid_metadata):
    for filename in quantity_mapping:
        write_radmc_input(filename=filename,
                          quantity=quantity_mapping[filename]['quantity'],
                          path=os.path.join('mdl', 'radmc_input_files'),
                          path=os.path.join('mdl', 'radmc_files'),
                          grid_metadata=grid_metadata)


@@ -141,7 +141,7 @@ def write_molecular_number_density_profiles(profiles: dict,
    for species in line_config['species_to_include']:
        write_radmc_input(filename=f'numberdens_{species}.inp',
                          quantity=profiles['gas_number_density'] * float(line_config['molecular_abundances'][species]),
                          path=os.path.join('mdl', 'radmc_input_files'),
                          path=os.path.join('mdl', 'radmc_files'),
                          grid_metadata=grid_metadata)


@@ -164,12 +164,12 @@ def main():
                           profiles=profiles)
    get_moldata(species_names=config['lines']['species_to_include'] + config['lines']['collision_partners'])
    write_radmc_lines_input(line_config=config['lines'],
                            path=os.path.join('mdl', 'radmc_input_files'))
                            path=os.path.join('mdl', 'radmc_files'))
    write_molecular_number_density_profiles(profiles=profiles,
                                            grid_metadata=grid_metadata,
                                            line_config=config['lines'])
    write_radmc_main_input_file(config=config,
                                path=os.path.join('mdl', 'radmc_input_files'))
                                path=os.path.join('mdl', 'radmc_files'))


if __name__ == '__main__':