Hi, Everybody!
Does anyone out there have a technique for getting the variance out of
a profile object? A profile object is good at getting <X> vs. B, I'd
then like to get < (X - <X>)^2 > vs B. Matt and I had spittballed the
possibility some time ago, but I was wondering if anyone out there had
successfully done it.
Thanks,
d.
--
Sent from my computer.
Hi,
I'm trying to add text to a particular point in my slice. I've got:
position = na.array([20.2, 13.8, 16.0])
pc = PlotCollection(pf, center=position)
p = pc.add_projection("Density", 2)
p.set_width(2, 'kpc')
for c in range(len(t1clouds)):
r = abs(t1clouds[c]['position']-position)
if r[0] < 1.0 and r[1] < 1.0:
p.modify["point"](t1clouds[c]['position'], str(t1clouds[c]['track']))
p.modify["marker"](t1clouds[c]['position'], marker='o')
pc.save("%s_tracking" % (fn))
But python complains with:
TypeError: text() argument after ** must be a mapping, not NoneType
Did I call the modifier incorrectly?
Elizabeth
Dear all
I have recently tried to install Yt, including Enzo, with the installation
script provided in the Yt webpage; setting the Enzo flag=1 and installing
as superuser in /usr/local all the files
I had the following two problems:
1) the enzo clone was not installed... just before finishing the script
said something like Python SSL not found or supported, and the enzo
installation was aborted.
2) After all this, I was not able to (a) open any pdf file either with
evince or xpdf (under linux; either ubuntu 10.04 and linux mint 13 -
running gnome), (b) login .. my "welcome screen" was off, completely gone
or black as in my laptop (a toshiba with linux mint). so I could not login!
In the case of my desktop at work, I had to reinstall the operating system.
Have any of the people of this forum seen this kind of problems? Any ideas
what might have gone wrong or how to fix this?
cheers,
hector
Hi all,
I've been trying to get a script (provided to me by Britton) working to
create a bunch of light rays from an Enzo simulation. Depending on the
initial seed value, it makes some of them correctly but eventually fails
with the following error:
Traceback (most recent call last):
File "make_light_rays.py", line 32, in <module>
get_los_velocity=True, get_nearest_galaxy=False)
File "/home/student/jomo9722/Enzo/
yt-x86_64/src/yt-hg/yt/analysis_modules/light_ray/light_ray.py", line 275,
in make_light_ray
(sub_ray[field])])
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
line 319, in __getitem__
self.get_data(key)
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
line 508, in get_data
if field not in ("dts", "t") and self._generate_field(field):
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
line 349, in _generate_field
self[field] = self.pf.field_info[field](self)
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py",
line 316, in __call__
dd = self._function(self, data)
File
"/data/moloney/cosmology_test/test_new/ion_balance/SD93IonBalance.py", line
178, in _ion_number_density
return solarAbundance[atom] * data[fractionField] * data['Metallicity']
* \
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
line 319, in __getitem__
self.get_data(key)
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
line 508, in get_data
if field not in ("dts", "t") and self._generate_field(field):
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
line 349, in _generate_field
self[field] = self.pf.field_info[field](self)
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py",
line 316, in __call__
dd = self._function(self, data)
File
"/data/moloney/cosmology_test/test_new/ion_balance/SD93IonBalance.py", line
193, in _ion_fraction_field
fraction = na.power(10, interp(data))
File
"/home/student/jomo9722/Enzo/yt-x86_64/src/yt-hg/yt/utilities/linear_interpolators.py",
line 44, in __call__
x_i = (na.digitize(x_vals, self.x_bins) - 1).astype('int32')
ValueError: Both x and bins must have non-zero length
Working through the debugger, I found that for the segment where it failed
the periodic_ray function had produced the following list of segments.
[[array([ 0.16103374, 0.78830872, 0.24846128]), array([ 0.14031834,
0.57367409, 1. ])], [array([ 0.14031834, 0.57367409, 1.
]), array([ 0.14031834, 0.57367409, 1. ])], [array([ 0.14031834,
0.57367409, 0. ]), array([ 0.13268221, 0.49455532, 0.2770327 ])]]
The second segment here has identical start and end points, which is
leading yt to try to interpolate a field with nothing in it (and therefore
throwing the above error). It seems that what's happening is that due to
rounding errors the first segment being produced actually has a value < 1
for the third dimension so that when it checks to see if the segment hit
any of the box walls it thinks that it didn't. It then tries to start a new
segment that's already at the boundary of the box, producing a ray segment
containing 0 grid cells.
Has anyone run into this problem before, and is there a fix for it floating
around out there somewhere? If not I could probably work something up that
checks to see if the distance from the segment endpoint to the wall is <
SmallNumber and set it equal to the value at the wall if it is.
Thanks.
Hi all,
I have been having trouble running yt on kraken. I am using the same script
that I use on trestles (where it sort of works). The issue I have on kraken
is that when I submit a non-interactive job, I get this error:
+ aprun -n 12 ./try1.py --parallel
Traceback (most recent call last):
File "/var/spool/alps/2310447/try1.py", line 3, in <module>
slice(imin=12, imax=47, width=20)
TypeError: slice() does not take keyword arguments
If I run an interactive job (./try1.py --parallel), I get this error:
yuan@krakenpf2:~/scratch/Thermal09> Assertion failed in file
/ptmp/ulib/mpt/nightly/5.2/051011/mpich2/src/mpid/cray/src/adi/mpid_init.c
at line 115: 0
aborting job:
(null)
I can run a non-parallel interactive job (./try1.py) without any problem.
My script looks like this:
#!/lustre/scratch/proj/sw/yt/dev/bin/python
from profile import *
slice(imin=41, imax=47, width=20)
Thanks!
Yuan
Hi all,
Currently I'm trying to analyze Flash4-beta data with YT, and for the data file uniform grid rather than AMR is used. I always get such error when I load the file:
>>> pf=load('relax_hdf5_chk_0000')
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/convenience.py", line 71, in load
return output_type_registry[candidates[0]](*args, **kwargs)
File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/frontends/flash/data_structures.py", line 214, in __init__
StaticOutput.__init__(self, filename, data_style)
File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py", line 92, in __init__
self._parse_parameter_file()
File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/frontends/flash/data_structures.py", line 377, in _parse_parameter_file
self.min_level = self.parameters["lrefine_min"] -1
KeyError: 'lrefine_min'
So, does YT provides a way to deal with the uniform grid?
Thanks for your care and consideration.
Best wishes,
Suoqing
Hi all,
My name is Devan Conroy and I am working at SLAC as a summer intern in the
SULI program with Jeff Oishi. I am trying to figure out how to count
the number of cells in the grid when loading in data. I tried using
(len(DataIn)), but got the error: object of type 'AMRRegion' has no len().
Any suggestions?
Thank you,
Devan Conroy
--
Devan Conroy
Fairfield University 2013
Mathematics Major
Computer Science/French Minor
Hi all,
My name is Devan Conroy and I am working at SLAC as a summer intern in the
SULI program with Jeff Oishi. I am trying to figure out how to count
the number of cells in the grid when loading in data. I tried using
(len(DataIn)), but got the error: object of type 'AMRRegion' has no len().
Any suggestions?
Thank you,
Devan Conroy
--
Devan Conroy
Fairfield University 2013
Mathematics Major
Computer Science/French Minor
Hi YT Users,
I have a computational domain that is not cubical, and I want to plot
a slice through it that is rectangular rather than square. I can make
a slice with PlotCollection, but when I write the image to a file
using save, I find that it has been padded with whitesapce on the left
and right to make it square. I would like the image to have the same
aspect ratio as the data, not be white-padded. Is there a way to do
this with PlotCollection, or should I just give up and use
FixedResolutionBuffer to extract the data and construct the image
myself?
Thanks.
--
Mark Krumholz
krumholz(a)ucolick.org
