Mailman 3 July 2011 - yt-users

derived field and TotalQuantity from two sources possible?
by gso＠physics.ucsd.edu July 19, 2011

July 19, 2011

Hi, I was wondering if it is possible to specify two different sources of input to create a derived field. Currently with a single source, I have a derived field defined as: def _HII_HFraction(field, data): return data["HII_Density"]/(data["HI_Density"]+data["HII_Density"]) add_field("HII_HFraction", function=_HII_HFraction, units=r"\frac{\rho_{HII}}{\rho_H}") And I can access the total quantity in parallel pf=load(file) dd = pf.h.all_data() dd.quantities["TotalQuantity"]("… [View More]

5 5

new fields in output
by Elizabeth Tasker July 17, 2011

July 17, 2011

Hi, I've added PotentialField to the data that gets outputted. I'd like to read this into yt, but how does this work when yt won't be expecting PotentialField to be one of the data sets? Elizabeth

3 3

Erroneous values from 'TotalQuantity' derived quantity
by Anthony Harness July 14, 2011

July 14, 2011

I am trying to use the derived quantity, 'TotalQuantity', to count the number of cells included in my cut_region as follows: cut = dd.cut_region(["grid['x'] > 0.5"]) #Example cut numberCells = cut.quantities['TotalQuantity']('Ones')[0] #Totals a 'one' for each cell, which gives the total number of cells being used This works fine when I have plenty of cells to split up between 84 processors (or when there are exactly zero cells in the cut), but when there is a small number of cells … [View More]

2 1

functions within functions
by Elizabeth Tasker July 14, 2011

July 14, 2011

Hi, I'm not sure where this is actually a yt question or a python one... I'm trying to create a slice of a new field that consists of the PotentialField (in the raw Enzo output) minus the background value, created by the profile tool. I have two functions that I created to do this: def _DiskRadius(field, data): ..... return dradius add_field("DiskRadius", function=_DiskRadius) That returns the value of the radius in the disk plane and def _SubtractBackgroundPotential(field, data):… [View More]

2 2

rendering FLASH 2D simulations
by Florian Bürzle July 11, 2011

July 11, 2011

Hi! While I'm at asking stupid questions - probably I missed it in the documentation (if so, I apologize), but is it possible to make rendered plots from FLASH 2D simulations? Thanks! Florian

2 2

ImportError: No module named amr_utils
by Florian Bürzle July 11, 2011

July 11, 2011

Hi! I've tried to get the hg version of yt working on my Mac Book Pro (OS X 10.6.8) using EPD Python 7.1., however unsuccessfully. I chose the development version hg up -C yt and performed the installation by using 'python setup.py install', and everything seemed to work fine (I captured the screen output in install_log.txt, attached). However, trying to import yt.mods gives the following error: Enthought Python Distribution -- www.enthought.com Version: 7.1-1 (32-bit) Python 2.7.2 |EPD … [View More]

3 4

disk Profiles
by Elizabeth Tasker July 11, 2011

July 11, 2011

Hi all, I'm trying to create a profile object of the density in my disk as a function of radius. I don't actually want to plot it; I just need to grab an average to use as a background value elsewhere. I think I've almost got it: source = pf.h. disk((16,16,16), (0,0,1), 12, 15.6e-3/pf['kpc']) profile = BinnedProfile1D(source, 128, "RadiusCode", 0.5, 12.0) profile.add_fields("CellMass", weight="CellVolume") print profile["CellMass"] I just have a couple of questions: 1. Is there a way to … [View More]

2 2

Re: [yt-users] disk Profiles
by Elizabeth Tasker July 11, 2011

July 11, 2011

In a fit of impatience born of being in the wrong time zone, can I wave this at people again? :) On 2011-07-09, at 6:23 PM, Elizabeth Tasker wrote: > Hi all, > > I'm trying to create a profile object of the density in my disk as a function of radius. I don't actually want to plot it; I just need to grab an average to use as a background value elsewhere. > > I think I've almost got it: > > source = pf.h. disk((16,16,16), (0,0,1), 12, 15.6e-3/pf['kpc']) > profile = … [View More]

1 0

changing the name of saved files
by Elizabeth Tasker July 8, 2011

July 8, 2011

Hi, I'd like to create a bunch of slices of the same field, but with different options (width, modifiers etc). Is it possible to change the file name yt saves the plots under so: GT_BTAccel_256AMR4_PeHeat_sf5_SNe_0301_Slice_z_PotentialField.png Doesn't get over-written with each one? Thanks you~ Elizabeth

3 3

Phase plot question
by david collins July 7, 2011

July 7, 2011

Hi -- I have a dumb question about the contents of the _data member of a BinnedProfile2D. I feel like I may have asked this question before, but I didn't see it in my archive, sorry if this is a dumb one. Why is the size of the Z array the same as Xbins \cross Ybins? For instance, I asked for 512 bins in both X and Y between X1,X2 and Y1,Y2. This gives me 513 bins in the Xbins array, as I expect (same for Y). But the Z array is also 513x513. My assumption has been that in Zarrray[I,J] is … [View More]

3 2