<div dir="ltr">Hey Team YT,<div><br></div><div style>I have a couple of (I think super-easy) questions:</div><div style><br></div><div style>a. How do you paste a script to the yt pastebin?  I've done it successfully once before, but can't for the life of me figure out how I did it.</div>

<div style><br></div><div style>b. Relatedly, there's a (like 10 line) script appended below where I projection plot gas densities from a Gadget simulation.    I'd like to be able to do the following:</div><div style>

<br></div><div style>    b.i - I'd like to be able to project other fields.  For example, if I look at pf.h.field_list, I get that for particle type 0, the fields:  ('PartType0', 'Masses'), ('PartType0', 'Metallicity') exist.  However, if I try to project:</div>

<div style><br></div><div style>>p = ProjectionPlot(pf,'z',(("deposit", "PartType0_metallicity"))) </div><div style><br></div><div style>(instead of the successful density line below), I get the error :</div>

<div style><br></div><div style><div>>YTFieldNotFound: Could not find field '('deposit', 'PartType0_metallicity')' in snapshot_006.</div><div><br></div></div><div style><br></div><div style>    b.ii.  Assuming the parts from b.i are possible, is it then possible to multiply two fields together in the projection (say, I wanted metal mass and wanted to multiply the particle masses by the particle metallicities).  Similarly, is it possible to multiply a field by a constant for the particle deposition into a projection plot?</div>

<div style><br></div><div style><br></div><div style>Thanks,</div><div style>Desika</div><div style><br></div><div style><br></div><div style>-------------------------</div><div style><div>from yt.mods import *</div><div>

<br></div><div><br></div><div>snap = "/Users/desika/powderday/grid_data/snapshot_006.hdf5"</div><div>bbox = [[-220.0, 220.0],</div><div>        [-220.0, 220.0],</div><div>        [-220.0, 220.0]]</div><div><br>
</div>
<div><br></div><div>unit_base = {'UnitLength_in_cm'         : 3.08568e+21,</div><div>             'UnitMass_in_g'            :   1.989e+43,</div><div>             'UnitVelocity_in_cm_per_s' :      100000}</div>

<div><br></div><div>pf = load(snap,unit_base = unit_base,bounding_box=bbox)</div><div>pf.h</div><div><br></div><div>p = ProjectionPlot(pf,'z',(("deposit", "PartType0_density")))</div><div>p.set_width(200,'kpc')</div>

<div>p.save()</div><div><br></div></div></div>