Hi Dave,
thanks a lot, setting y_log=False does the trick!
Wolfram
Message: 3 Date: Fri, 30 Sep 2011 16:32:58 -0600 From: David Collins dcollins@physics.ucsd.edu To: Discussion of the yt analysis package yt-users@lists.spacepope.org Subject: Re: [yt-users] log-linear phase plot Message-ID: CAGva0oysAaoEf8cPXxCevTtJqjcQX55n15Bixsx2WPMFD3TzXQ@mail.gmail.com Content-Type: text/plain; charset=ISO-8859-1
I've run into this as well. I believe you need to give add_phase_object the log value by hand, with
x_log=False
Cheers, d.
On Fri, Sep 30, 2011 at 3:52 PM, Geoffrey So gsiisg@gmail.com wrote:
Just wondering, if FieldOne is already logged, do you still want the limits to encompass 1 to 1e10, have you tried something like (1e9,1e10) ? From G.S.
On Fri, Sep 30, 2011 at 2:37 PM, Wolfram Schmidt schmidt@astro.physik.uni-goettingen.de wrote:
Hi everyone,
I defined two fields, say, FieldOne and FieldTwo, to make 2D phase plot. FieldOne is logarithmic by default, FieldTwo is linear (i.e., I set take_log=False in add_field).
I want to produce a phase plot with logarithmic bins for FieldOne and linear bins for FieldTwo, where the range of FieldTwo is [-10,10]
I thought that
pc.add_phase_object(dd, ["FieldOne", "FieldTwo", "CellMass"], weight=None, x_bins=100,y_bins=100, \ ? ? ? ? ? ? ? ? ? ?x_bounds=[1e0,1e10], y_bounds=[-10,10])
might do the job, but ?yt returns the error message:
Warning: invalid value encountered in log10 yt : [ERROR ? ?] 2011-09-30 16:55:05,610 Your min/max values for x, y have given me a nan. yt : [ERROR ? ?] 2011-09-30 16:55:05,610 Usually this means you are asking for log, with a zero bound. Traceback (most recent call last): ?File "stability.py", line 352, in <module> ? ?x_bounds=[1e0,1e10], y_bounds=[-10,10]) ?File "/nics/e/sw/analysis/yt/2.1/sles11.1_gnu4.3.4/src/yt-hg/yt/visualization/plot_collection.py", line 1149, in add_phase_object ? ?lazy_reader) ?File "/nics/e/sw/analysis/yt/2.1/sles11.1_gnu4.3.4/src/yt-hg/yt/data_objects/profiles.py", line 408, in __init__ ? ?raise ValueError ValueError
It appears that add_phase_object treats FieldTwo logarithmically, although it is linear.
So how can I do a log-linear phase plot?
Best, Wolfram _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-
Wolfram Schmidt