Yael-II/MSc1-Internship-Stellar-Cluster
1
0
Fork 0
mirror of https://codeberg.org/Yael-II/MSc1-Internship-Stellar-Cluster.git synced 2026-03-15 03:16:26 +01:00
Code Issues Packages Projects Releases Wiki Activity

Delete COSMIC/COSMIC_v3_emission.py

Browse Source
This commit is contained in:
Yaël M
2024-06-21 16:22:05 +02:00
committed by GitHub
parent abf12651ff
commit d07dd2ed88
1 changed files with 0 additions and 53 deletions
Download Patch File Download Diff File Expand all files Collapse all files

53
COSMIC/COSMIC_v3_emission.py
View File

@@ -1,53 +0,0 @@
"""
* COSMIC - EMISSION
* Version 3
@ Yaël Moussouni
@ Unistra, P&E, MSc1-MoFP
@ Observatory of Strasbourg (Intern)
"""
import numpy as np
import amuse.units.units as u
def rotation(stars, params):
N = params["N_stars"] + np.int32(np.round(params["N_stars"] * params["binary_fraction"]))
stars.X_v = np.random.uniform(0, 10, N) | u.km * u.s**(-1)
p = np.random.uniform(0, 1, N)
stars.X_vsini = stars.X_v * np.sqrt(2*p-p**2)
return stars
def X_emission(stars, params):
for i in range(len(stars)):
if stars[i].temperature.value_in(u.K) > 7300: # O B A
stars[i].X_luminosity = params["X_coefficient_OBA"] * stars[i].luminosity.in_(u.erg * u.s**(-1))
elif stars[i].temperature.value_in(u.K) < 6000: # G K M
stars[i].X_luminosity = params["X_coefficient_GKM"] * (stars[i].X_vsini.value_in(u.km * u.s**(-1)))**(1.9) | u.erg * u.s**(-1)
else: # F
stars[i].X_luminosity = 0 | u.erg * u.s**(-1)
return stars
def plasma_temperature(stars, params):
for i in range(len(stars)):
if stars[i].temperature.value_in(u.K) > 7300 : # O B A
if stars[i].age.value_in(u.Myr) < 150: # young
stars[i].X_temperature_0 = 0.6 | u.kilo(u.eV)
stars[i].X_temperature_1 = 0.0 | u.kilo(u.eV)
else: # old
stars[i].X_temperature_0 = 0.5 | u.kilo(u.eV)
stars[i].X_temperature_1 = 0.0 | u.kilo(u.eV)
elif stars[i].temperature.value_in(u.K) < 6000: # G K M
if stars[i].age.value_in(u.Myr) < 150: # young
stars[i].X_temperature_0 = 0.4 | u.kilo(u.eV)
stars[i].X_temperature_1 = 1.0 | u.kilo(u.eV)
else: # old
stars[i].X_temperature_0 = 0.2 | u.kilo(u.eV)
stars[i].X_temperature_1 = 0.8 | u.kilo(u.eV)
else: # F
if stars[i].age.value_in(u.Myr) < 150: # young
stars[i].X_temperature_0 = 0.6 | u.kilo(u.eV)
stars[i].X_temperature_1 = 0.0 | u.kilo(u.eV)
else: # old
stars[i].X_temperature_0 = 0.5 | u.kilo(u.eV)
stars[i].X_temperature_1 = 0.0 | u.kilo(u.eV)
return stars