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
