EVtwin throws for stars with M>20 Msun

[eichryso]

Hello all!

It seems like I cannot use EVtwin for stars more massive than 20 Msun. I run the following script:

from amuse.units import units
from amuse.community.evtwin.interface import EVtwin
from amuse.datamodel import Particle
from amuse.ext.star_to_sph import convert_stellar_model_to_SPH

mass = 30 | units.MSun
age = 3.5 | units.Myr
n_particles = 1000

stellar_evolution = EVtwin()
star = Particle(mass=mass)
se_star = stellar_evolution.particles.add_particle(star)
se_star.evolve_for(age)

sph_particles = convert_stellar_model_to_SPH(se_star, n_particles).gas_particles
stellar_evolution.stop()

and i get an amuse exception error: Error when calling 'evolve_for' of a 'EVtwin', errorcode is -1002.
When I tried the same code using a smaller mass of 10Msun it run just fine. Is this a problem of the EVtwin code?

Thanks in advance,
Eirini

[LourensVeen]

Hi Eirini,

The developers of EVTwin inform me that it's simply not designed to deal with such heavy stars, so it's not you. It would be nice if it gave a better error message though, so let's leave this issue open and take it as a request for better documentation.