def run_chain_parameters():
os.system("rm -rf SPECPAR.log")
logging.basicConfig(filename='SPECPAR.log', format='%(asctime)s : %(message)s', datefmt='%m/%d/%Y %I:%M:%S %p', level=logging.INFO)
logging.info('Started')
path = './'
run_path = 'running_dir/'
save_path = 'save_folder/'
# CREATE folders
# Running ARES
linelist_ew = rp.get_install_dir()+rc.read_config_param('linelists', 'iron_parameters').replace("'", "")
rp.ares_make_mine_opt(run_path, linelist_ew)
logging.info('ARES Started')
rp.run_ares(run_path)
logging.info('ARES Finished')
# Creating moog linelist
filename_lines = rp.get_install_dir()+rc.read_config_param('linelists', 'iron_parameters').replace("'", "")
filename_ares = rc.read_config_param('ares', 'fileout').replace("'", "")
filename_out = 'lines.' + filename_ares
isp.ares_to_lines(run_path+filename_ares, filename_lines, run_path+filename_out, 4000, 9999, 5, 150)
logging.info('Starting AMEBSA for %s', filename_out)
# find_iron_parameters(path,run_path,save_path,filename_out)
find_iron_parameters_tmcalc_prior(path, run_path, save_path, filename_out, filename_ares)
logging.info('Finished')
def run_chain_get_element_abund(moogfile, element):
run_path = 'running_dir/'
save_path = 'save_folder/'
teff, logg, feh, vtur = isp.read_parameters_moogfile(moogfile)
linelist_element = rp.get_install_dir()+rc.read_config_param('linelists', element+'_abund').replace("'", "")
rp.ares_make_mine_opt(run_path, linelist_element)
rp.run_ares(run_path)
filename_ares = rc.read_config_param('ares', 'fileout').replace("'", "")
filename_out = 'lines.' + filename_ares
isp.ares_to_lines(run_path+filename_ares, linelist_element, run_path+filename_out, 4000, 9999, 5, 150)
rp.create_abfind_par(run_path, filename_out)
rp.create_model_kurucz(run_path, teff, logg, feh, vtur)
rp.run_MOOG(run_path, 'abfind.par')
(ele1, ele1_sig, nele1, ele2, ele2_sig, nele2) = rmoog.read_moog_ele_sigma(run_path+'abund_plan_tspec.test', element, 2.)
return (ele1, ele1_sig, nele1, ele2, ele2_sig, nele2)
