first I want to thank for every time I got an answer when I write
here, I really appreciate it.
Now I am writing because I am using PlotCollection (As you can see my
code) and I have some questions
1) the plot is coming in cm^3 and I want to change it to kpc^3 is there
a way to change it?
2) I am attaching the plot can I change to dot or any different symbol
3) why some plot come "empty"? and also I get this error message
/Data has no positive values, and therefore can not be log-scaled./
sphere = pf.h.sphere(c, (250., 'kpc')) # 250.
levels = 8 #8
mi, ma = sphere.quantities["Extrema"]("Density")
contour_values, connected_sets =
sphere.extract_connected_sets("Density", levels, mi, ma)
for i in connected_sets:
eb = connected_sets[i]
print "level " + str(i) + " clmsp# " + str(len(eb))
for j in eb:
obj = eb[j]
com = obj.quantities["CenterOfMass"]()
pc = PlotCollection(pf, com)
"Temperature"],weight="CellMassMsun") # average T(rho)
pc.add_profile_object(obj, ["Density", "CellVolume"],
pc.add_profile_object(obj, ["Density", "CellMassMsun"],
except Exception, e:
Thank you in advance!!