Hi guys,

I want to use the nparticles rather than the particles callback so that I only plot PARTICLE_TYPE_MUST_REFINE (4). I am trying (full script at end of email):

size=4.0, col='g', ptype='4')

I tried with and without ' around the 4 but both times I get:

Traceback (most recent call last):
  File "Density_movie.py", line 29, in <module>
    p.modify["nparticles"](1.0,p_size=4.0, col='g', ptype='4')
  File "/scratch/YT-x86_64/lib/python2.5/site-packages/yt-1.6dev-py2.5-linux-x86_64.egg/yt/raven/PlotTypes.py", line 52, in __getitem__
    raise KeyError(item)
KeyError: 'nparticles'

Any idea what I am doing wrong?

Many thanks,



from yt.mods import *
import pylab
import sys
import cmath
import matplotlib.ticker

# Code to make the frames for a movie showing density projection

rho_min = 1.e-22
rho_max = 1.e-32

frame_template = "aaWoC/%s/densityframe_%01i_%04i"

for n in range(12,200): # make density projection with particles on
    #pf = lagos.EnzoStaticOutput("WOCbin/%s/DD%04i/data%04i" % (sys.argv[-1],n, n))
    pf = lagos.EnzoStaticOutput("WOCbin2/%s/DD%04i/data%04i" % (sys.argv[-1],n, n))
    #print pf.h.field_list
    pc = raven.PlotCollection(pf, center=(0.5,0.5,0.5))
    time = pf["InitialTime"]*pf["years"]
    for ax in range(1):
        p = pc.add_projection("Density", ax, weight_field="Density")
        #pc.set_zlim(rho_min, rho_max)
        sp = pf.h.all_data()
        colours = sp["ParticleMassMsun"]
        if sp["ParticleMassMsun"].size > 0:
            p.modify["nparticles"](1.0,p_size=4.0, col='g', ptype='4')
        p.modify["text"]((0.1,0.9),"time = %i years" % (pf["InitialTime"]*pf["years"]), text_args = {'color':'w'})
        #p.set_width(0.1, 'pc')
        pc.save(frame_template % (sys.argv[-1],ax,n), override=True)

Elizabeth Harper-Clark MA MSci
PhD Candidate, Astrophysics, UofT

AIM: edphc1
MSN: edphc1@hotmail.com
Skype: eharperclark
Office phone: 416-978-5759