From d87a37bdc5f49f93690d2c40d4452b0f7747e39c 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:22:13 +0200 Subject: [PATCH] Delete COSMIC/COSMIC_v3_evolve.py --- COSMIC/COSMIC_v3_evolve.py | 47 -------------------------------------- 1 file changed, 47 deletions(-) delete mode 100644 COSMIC/COSMIC_v3_evolve.py diff --git a/COSMIC/COSMIC_v3_evolve.py b/COSMIC/COSMIC_v3_evolve.py deleted file mode 100644 index dec35fa..0000000 --- a/COSMIC/COSMIC_v3_evolve.py +++ /dev/null @@ -1,47 +0,0 @@ -""" -* COSMIC - EVOLVE -* 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 COSMIC_v3_output as output -import COSMIC_v3_emission as emission -import COSMIC_v3_coordinates as coords - -def stellar_gravity(codes, stars, channels, params, save_stars = False): - """ - Evolve a model with both gravity and stellar evolution. - """ - for t in params["output_times"]: - print("\033[36m"+"t: "+"\033[0m"+"{:02.2f} Myr/{:02.2f} Myr".format(t.value_in(u.Myr),np.max(params["output_times"].value_in(u.Myr)))) - codes["g"].evolve_model(t) - channels.copy() - if save_stars: output.save_stars(stars, params, t.value_in(u.Myr)) - - print("\033[36m"+"Simulation completed"+"\033[0m") - return 0 - -def stellar_gravity_binaries(codes, stars, binaries, channels, params, save_stars = False, save_binaries = False, compute_X_emission = False): - """ - Evolve a model with both gravity and stellar evolution. - """ - for t in params["output_times"]: - print("\033[36m"+"t: "+"\033[0m"+"{:02.2f} Myr/{:02.2f} Myr".format(t.value_in(u.Myr),np.max(params["output_times"].value_in(u.Myr)))) - codes["b"].evolve_model(t) - channels.copy() - stars = coords.xyz2radecdist(stars) - if compute_X_emission: - stars = emission.rotation(stars, params) - stars = emission.X_emission(stars, params) - stars = emission.plasma_temperature(stars, params) - if save_stars: output.save_stars(stars, params, t.value_in(u.Myr)) - if save_binaries: output.save_binaries(binaries, params, t.value_in(u.Myr)) - - print("\033[36m"+"Simulation completed"+"\033[0m") - return codes, stars, channels - -