# We want to find clumps over the entire dataset, so we'll just grab the whole
# thing! This is a convenience parameter that prepares an object that covers
# the whole domain. Note, though, that it will load on demand and not before!
-data_source = pf.h.all_data()
+data_source = pf.h.disk([0.5, 0.5, 0.5], [0., 0., 1.],
+ 8./pf.units['kpc'], 1./pf.units['kpc'])
# Now we set some sane min/max values between which we want to find contours.
# This is how we tell the clump finder what to look for -- it won't look for
c_min = 10**na.floor(na.log10(data_source[field]).min() )
c_max = 10**na.floor(na.log10(data_source[field]).max()+1)
+# keep only clumps with at least 20 cells
+function = 'self.data[field].size > 20'
# Now find get our 'base' clump -- this one just covers the whole domain.
amods.level_sets.Clump(data_source, None, field)
+master_clump = Clump(data_source, None, field)
# This next command accepts our base clump and we say the range between which
# we want to contour. It recursively finds clumps within the master clump, at
# intervals defined by the step size we feed it. The current value is
# *multiplied* by step size, rather than added to it -- so this means if you
# want to look in log10 space intervals, you would supply step = 10.0.
amods.level_sets.find_clumps(master_clump, c_min, c_max, step)
+find_clumps(master_clump, c_min, c_max, step)
# As it goes, it appends the information about all the sub-clumps to the
# master-clump. Among different ways we can examine it, there's a convenience
f = open('%s_clumps.txt' % pf,'w')
+# We can traverse the clump hierarchy to get a list of all of the 'leaf' clumps
+leaf_clumps = get_lowest_clumps(master_clump)
# If you'd like to visualize these clumps, a list of clumps can be supplied to
-# the "clumps" callback on a plot.
+# the "clumps" callback on a plot. First, we create a projection plot:
+prj = ProjectionPlot(pf, 2, field, center='c', width=(20,'kpc'))
+# Next we annotate the plot with contours on the borders of the clumps
+# Lastly, we write the plot to disk.
+# We can also save the clump object to disk to read in later so we don't have
+# to spend a lot of time regenerating the clump objects.
+# Later, we can read in the clump object like so,
+master_clump = pf.h.load_object('My_clumps')