feat: Add stellar temperature

Author: tgross03

src/radiosim/ppdisks/simulation.py

@@ -36,6 +36,7 @@ def get_default_sampling_config():
             "binary_period": [6.04800e5, 3e9],  # Seconds (logarithmic sampling)
             "binary_eccentricity": [0.0, 0.4],  # 0 = Circle, 0 < e < 1 = Ellipse
             "stellar_mass": [0.5, 5],  # Solar Masses
+            "stellar_temperature": [2000.0, 3000.0],  # Kelvin
             "num_planets": [1, 5],
             "planet_mass": [1.0e-6, 5.0e-3],  # Solar Masses
             "planet_orbit_radius": [7.0, 30.0],  # Astronomical Units
@@ -264,6 +265,11 @@ def simulate(
 
             distances = self._planet_config.get_distances()
 
+            # If binary: distance -> binary eccentricity || binary period != distance
+            # for the two stars in the planet file
+            if planet_parameters["binary_system"]:
+                distances = distances[2:]
+
             param_config["mesh_parameters.ymin"] = (
                 (np.min([mesh_parameters["y_min"], distances.min()]) * un.AU)
                 .to(self._unit_system.length)
@@ -593,6 +599,12 @@ def sample_dict(read_dict):
             size=num_stars,
         )
 
+        planet_parameters["stellar_temperature"] = rng.uniform(
+            low=planet_sampling["stellar_temperature"][0],
+            high=planet_sampling["stellar_temperature"][1],
+            size=num_stars,
+        )
+
         num_planets = rng.integers(
             low=planet_sampling["num_planets"][0],
             high=planet_sampling["num_planets"][1],