From fd843010fa60c5080af8f79244c7e2b9eb1aa5bd Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ya=C3=ABl=20M?= <154762804+Yael-II@users.noreply.github.com> Date: Fri, 21 Jun 2024 16:21:50 +0200 Subject: [PATCH] Delete COSMIC/COSMIC_v3_config.py --- COSMIC/COSMIC_v3_config.py | 102 ------------------------------------- 1 file changed, 102 deletions(-) delete mode 100644 COSMIC/COSMIC_v3_config.py diff --git a/COSMIC/COSMIC_v3_config.py b/COSMIC/COSMIC_v3_config.py deleted file mode 100644 index b0c9e16..0000000 --- a/COSMIC/COSMIC_v3_config.py +++ /dev/null @@ -1,102 +0,0 @@ -""" -* COSMIC - CONFIG -* Version 3 -@ Yaël Moussouni -@ Unistra, P&E, MSc1-MoFP -@ Observatory of Strasbourg (Intern) -""" - -import numpy as np -import amuse.units.units as u -import datetime - -def params_dic(): - params_list = ["N_stars", "timestep", "output_times", "R0_king", "W0_king", "metallicity_stellar", "M_min_salpeter", "M_max_salpeter", "a_salpeter", "binary_fraction", "mean_period", "std_period", "X_coefficient_OBA", "X_coefficient_GKM", "cluster_position_x", "cluster_position_y", "cluster_position_z", "out_directory", "in_directory", "config", "filename", "workers_stellar", "workers_gravity", "format_type"] - params = {i:0 for i in params_list} - return params, params_list - -def verbatim(params): - print("\033[94m"+"Info: Simulation time of "+"\033[0m"+"{:.2f} Myr".format(np.max(params["output_times"].value_in(u.Myr)))) - print("\033[94m"+"Info: Simulation of "+"\033[0m"+"{:d} stars".format(params["N_stars"])) - print("\033[94m"+"Info: Simulation will be saved with the prefix "+"\033[0m"+"{}".format(params["filename"])) - print("\033[94m"+"Info: Gravity: "+"\033[0m"+"{:d} workers".format(params["workers_gravity"])+"\033[94m"+", Stellar: "+"\033[0m"+"{:d} workers".format(params["workers_stellar"])) - return 0 - -def default(): - params, params_list = params_dic() - - date = datetime.datetime.now() - date = "{:04d}-{:02d}-{:02d}_{:02d}h{:02d}".format(date.year, date.month, date.day, date.hour, date.minute) - - params["N_stars"] = np.int32(10) # Number of particle - params["timestep"] = np.float32(1)|u.Myr # Time increment - params["output_times"] = np.float32([1,2,5,10,20,50,100,200,500,1000,2000,5000])|u.Myr # Output intants - params["R0_king"] = np.float32(1)|u.pc # Cluster virial radius - params["W0_king"] = np.float32(9) # Cluster King's parameter - params["metallicity_stellar"] = np.float32(0.02) # Metallicity - params["M_min_salpeter"] = np.float32(0.1)|u.MSun # Minimum mass for Salpeter's law - params["M_max_salpeter"] = np.float32(125)|u.MSun # Maximum mass for Salpeter's law - params["a_salpeter"] = np.float32(-2.35) # Salpeter's law coefficient - params["binary_fraction"] = np.float32(0.13) # Fraction of binary stars - params["mean_period"] = np.float32(4.54) # Mean orbital period of binary systems (log10(T/days)) - params["std_period"] = np.float32(2.4) # Orbital period deviation of binary systems - params["X_coefficient_OBA"] = np.float32(1.4e-7) # X-ray luminosity coefficient for O, B and A type stars - params["X_coefficient_GKM"] = np.float32(1.4e27) # X-ray luminosity coefficient for G, K and M type stars - params["cluster_position_x"] = np.float32(7000)|u.pc # Position of the cluster in the disk along the sun-center axis - params["cluster_position_y"] = np.float32(0)|u.pc # Position of the cluster in the disk in the perpendicular direction - params["cluster_position_z"] = np.float32(0)|u.pc # Position of the cluster above or below the disk - params["out_directory"] = "./output/" - params["in_directory"] = "./input/" - params["config"] = "default.cfg" - params["filename"] = "COSMIC_{}".format(date) - params["workers_stellar"] = np.int32(1) - params["workers_gravity"] = np.int32(1) - params["format_type"] = "csv" - return params, params_list - -def change_params(params, cfg_param, cfg_values): - if cfg_param == "": - return params - elif cfg_param == "N_stars": params["N_stars"] = np.int32(cfg_values[0]) - elif cfg_param == "timestep": params["timestep"] = np.float32(cfg_values[0])|u.Myr - elif cfg_param == "output_times": params["output_times"] = np.float32(cfg_values)|u.Myr - elif cfg_param == "R0_king": params["R0_king"] = np.float32(cfg_values[0])|u.pc - elif cfg_param == "W0_king": params["W0_king"] = np.float32(cfg_values[0]) - elif cfg_param == "metallicity_stellar": params["metallicity_stellar"] = np.float32(cfg_values[0]) - elif cfg_param == "M_min_salpeter": params["M_min_salpeter"] = np.float32(cfg_values[0])|u.MSun - elif cfg_param == "M_max_salpeter": params["M_max_salpeter"] = np.float32(cfg_values[0])|u.MSun - elif cfg_param == "a_salpeter": params["a_salpeter"] = np.float32(cfg_values[0]) - elif cfg_param == "binary_fraction": params["binary_fraction"] = np.float32(cfg_values[0]) - elif cfg_param == "mean_period": params["mean_period"] = np.float32(cfg_values[0]) - elif cfg_param == "std_period": params["std_period"] = np.float32(cfg_values[0]) - elif cfg_param == "X_coefficient_OBA": params["X_coefficient_OBA"] = np.float32(cfg_values[0]) - elif cfg_param == "X_coefficient_GKM": params["X_coefficient_GKM"] = np.float32(cfg_values[0]) - elif cfg_param == "cluster_position_x": params["cluster_position_x"] = np.float32(cfg_values[0])|u.pc - elif cfg_param == "cluster_position_y": params["cluster_position_y"] = np.float32(cfg_values[0])|u.pc - elif cfg_param == "cluster_position_z": params["cluster_position_z"] = np.float32(cfg_values[0])|u.pc - elif cfg_param == "out_directory": params["out_directory"] = cfg_values[0].replace("\n", "") - elif cfg_param == "in_directory": params["in_directory"] = cfg_values[0].replace("\n", "") - elif cfg_param == "config": params["config"] = cfg_values[0].replace("\n", "") - elif cfg_param == "filename": params["filename"] = cfg_values[0].replace("\n", "") - elif cfg_param == "workers_stellar": params["workers_stellar"] = np.int32(cfg_values[0]) - elif cfg_param == "workers_gravity": params["workers_gravity"] = np.int32(cfg_values[0]) - elif cfg_param == "format_type": params["format_type"] = cfg_values[0].replace("\n", "") - else: - print("\033[93m"+"Warning: cannot understand the parameter in config: "+"\033[0m" + cfg_param) - return params - - -def from_cfg_file(params): - cfg_path = params["in_directory"] + params["config"] - with open(cfg_path, "r") as cfg_file: - for cfg_line in cfg_file.readlines(): - cfg_split = cfg_line.split(" ") - if len(cfg_split) == 1: - print("\033[93m"+"Warning: cannot understand the line in config: "+"\033[0m" + cfg_line) - continue - cfg_param = cfg_split[0] - cfg_values = cfg_split[1:] - params = change_params(params, cfg_param, cfg_values) - verbatim(params) - return params -