def slowtest11(self): print "Test for importing new stellar models, also check long-term evolution" instance = EVtwin() instance.parameters.verbosity = True instance.particles.add_particles(Particles(mass = [0.8] | units.MSun)) instance.evolve_model() copy = Particles(1) copy.internal_structure = instance.particles[0].get_internal_structure() instance.particles.add_particles(copy) number_of_zones = instance.particles[0].get_number_of_zones() self.assertEqual(len(instance.particles), 2) self.assertEqual(instance.particles[1].get_number_of_zones(), number_of_zones) self.assertIsOfOrder(instance.particles[1].get_radius_profile()[-1], 1.0 | units.RSun) self.assertIsOfOrder(instance.particles[1].get_temperature_profile()[0], 1.0e7 | units.K) self.assertIsOfOrder(instance.particles[1].get_pressure_profile()[0], 1.0e17 | units.barye) t1, l1 = simulate_evolution_tracks(instance.particles[0], end_time=26.5|units.Gyr) t2, l2 = simulate_evolution_tracks(instance.particles[1], end_time=26.5|units.Gyr) instance.stop() i1 = t1.argmax() # Maximum temperature ~ turnoff main sequence i2 = t2.argmax() self.assertAlmostRelativeEqual(t1[i1], t2[i2], 2) self.assertAlmostRelativeEqual(l1[i1], l2[i2], 2)