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])))
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])))