syllabus example code: stellar_minimal.pyΒΆ

[source code]

"""
   Minimalistic routine for running a stellar evolution code
"""
from amuse.lab import *

def main(M, z, model_time):
#    stellar = MESA()
    stellar = SSE()
    stellar.parameters.metallicity = z
    stellar.particles.add_particle(Particle(mass=M))
    stellar.commit_particles()

    initial_luminosity = stellar.particles.luminosity
    dt = 0.1 | units.Myr
    time = 0 | units.Myr
    while stellar.particles[0].age<model_time:
        time+=dt
        stellar.evolve_model(time)

    final_luminosity = stellar.particles.luminosity
    print "L(t=0)=", initial_luminosity, \
        ", L (t=", stellar.particles.age, ")=", \
        final_luminosity, stellar.particles.radius, stellar.particles.mass.in_(units.MSun)

    stellar.stop()

def new_option_parser():
    from amuse.units.optparse import OptionParser
    result = OptionParser()
    result.add_option("-M", unit= units.MSun,
                      dest="M", type="float",default = 1.0 | units.MSun,
                      help="stellar mass [%default]")
    result.add_option("-t", unit = units.Myr,
                      dest="model_time", type="float",
                      default = 4700.0|units.Myr,
                      help="end time of the simulation [%default]")
    result.add_option("-z", dest="z", type="float",
                      default = 0.02, help="metalicity [%default]")
    return result

if __name__ in ('__main__', '__plot__'):
    o, arguments  = new_option_parser().parse_args()
    main(**o.__dict__)

Keywords: python, amuse, astrophysics, matplotlib, pylab, example, codex (see how-to-search-examples)