def get_stellar_radius(star, SEVCode = None): if SEVCode == None: sev_code = SSE() else: sev_code = SEVCode temp_star = Particle() temp_star.mass = star.mass temp_star.age = star.time sev_code.particles.add_particle(temp_star) sev_code.model_time = star.time sev_code.evolve_model(star.time) radius = sev_code.particles[0].radius print(radius, sev_code.particles[0].age) if SEVCode == None: sev_code.stop() else: sev_code.particles.remove_particle(temp_star) return radius