def iliev_test_5(N=10000,
L=15. | units.kpc,
source_luminosity=5.e48 | units.s**-1,
rhoinit=0.001 | (units.amu / units.cm**3),
Tinit=100. | units.K,
hydro_parameters=dict(),
rad_parameters=dict()):
gamma = 5. / 3.
mu = 1. | units.amu
x, y, z = body_centered_grid_unit_cube(N).make_xyz()
N = len(x)
print N
gas = Particles(N)
# set particles homogeneously in space
# gas.x=L*numpy.random.uniform(-1.,1.,N)
# gas.y=L*numpy.random.uniform(-1.,1.,N)
# gas.z=L*numpy.random.uniform(-1.,1.,N)
gas.x = L * x
gas.y = L * y
gas.z = L * z
gas.x += L * numpy.random.uniform(-1., 1., N) / 1000
gas.y += L * numpy.random.uniform(-1., 1., N) / 1000
gas.z += L * numpy.random.uniform(-1., 1., N) / 1000
# set other properties
gas.h_smooth = 0. | units.parsec
gas.vx = 0. | (units.km / units.s)
gas.vy = 0. | (units.km / units.s)
gas.vz = 0. | (units.km / units.s)
gas.u = 1 / (gamma - 1) * constants.kB * Tinit / mu
gas.rho = rhoinit
gas.mass = rhoinit * (2 * L)**3 / N
gas.xion = 0.1
gas.metallicity = 0.02
# set source in the center of the box
sources = Particles(1)
sources.x = 0. * L
sources.y = 0. * L
sources.z = 0. * L
sources.luminosity = source_luminosity
return gas, sources
