def test4(self):
print("test4: evolve test (10 part)")
instance = self.new_instance_of_an_optional_code(
Krome, **default_options)
parts = Particles(10)
parts.number_density = 1.e5 | units.cm**-3
parts.temperature = 50 | units.K
parts.ionrate = 2.e-17 | units.s**-1
Ns = len(instance.species)
parts.abundances = numpy.zeros((10, Ns))
instance.particles.add_particles(parts)
instance.evolve_model(1. | units.Myr)
f = 2 * instance.particles[0].abundances[instance.species["H2"]]
self.assertTrue(f > 0.95) # not much of a test..
#~ for x,i in instance.species.items():
#~ print x, instance.particles[0].abundances[i]
instance.cleanup_code()
instance.stop()
def test4(self):
print "test4: evolve test (10 part)"
instance=self.new_instance_of_an_optional_code(Krome,**default_options)
parts=Particles(10)
parts.number_density=1.e5 | units.cm**-3
parts.temperature=50 | units.K
parts.ionrate=2.e-17 | units.s**-1
Ns=len(instance.species)
parts.abundances=numpy.zeros((1,Ns))
instance.particles.add_particles(parts)
instance.evolve_model( 1. | units.Myr )
f=2*instance.particles[0].abundances[instance.species["H2"]]
self.assertTrue(f> 0.95) # not much of a test..
#~ for x,i in instance.species.items():
#~ print x, instance.particles[0].abundances[i]
instance.cleanup_code()
instance.stop()
def makeparts(self,N):
parts=Particles(N)
numpy.random.seed(1234567)
parts.number_density=(numpy.random.random(N)*1.e5+1.e5)| units.cm**-3
parts.temperature=(numpy.random.random(N)*500+100)| units.K
parts.ionrate=(numpy.random.random(N)*1.e-11+1.e-17)| units.s**-1
return parts
def makeparts(self,N):
parts=Particles(N)
numpy.random.seed(1234567)
parts.number_density=(numpy.random.random(N)*1.e5+1.e5)| units.cm**-3
parts.temperature=(numpy.random.random(N)*500+100)| units.K
parts.ionrate=(numpy.random.random(N)*1.e-11+1.e-17)| units.s**-1
return parts
