star['T'] = ct.cosmology.a2t(star['sft'])
for c, s, cc, ll, v in zip(center, size, C, L, vel):
ct.view(center=c, size=s)
dC, dL = ct.paint('star', 'T', 'mass', 'sml')
dC[dL == 0] = 0
cc += dC * dL
ll += dL
return C/L, L
if dostar:
Cs, Ls = paintstar(center, radius * 0.2, vel)
for i, c, s, cs, ls in zip(range(len(Cs)), center, radius * 0.2, Cs, Ls):
figure = Figure(dpi=200)
ax = ct.newaxes(figure)
ct.view(center=c, size=s)
now = ct.cosmology.z2t(ct.C['redshift'])
past = ct.cosmology.z2t(11.0)
ct.imshow(n_(cs, vmin=past, vmax=now), nl_(ls, vmin='25db'),
cmap=cm.coolwarm_r)
ct.frame(axis=False, scale={'color':'w'})
ct.print_svg('%03d/star-%05d.svg' % (snapid, int(mass[i])))
else:
Cg, Lg = paintgas(center, radius * 10)
for i, c, s, cg, lg in zip(range(len(Cg)), center, radius * 10, Cg, Lg):
figure = Figure(dpi=200)
ax = ct.newaxes(figure)
ct.view(center=c, size=s)
ct.imshow(nl_(cg, vmin=3.5, vmax=7.5), nl_(lg, vmin='40db'))
ct.frame(axis=False, scale={'color':'w'})
ct.print_svg('%03d/gas-%05d.svg' % (snapid, int(mass[i])))
