Hi all, there is something regarding the x-velocity field I don't really understand. When I analyze a file with yt, I find: ###################################################################### In [1]: filen = 'DD0001/CommonEnvelope0001' In [2]: pf = load(filen) In [3]: region = pf.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0], [1.0,1.0, 1.0]) yt INFO 2010-08-08 14:29:57,375 Getting the binary hierarchy yt INFO 2010-08-08 14:29:57,389 Finished with binary hierarchy reading In [4]: vx = region["x-velocity"] yt INFO 2010-08-08 14:30:06,647 Getting field x-velocity from 1 In [5]: max(vx) Out[5]: 232812.74801568719 ###################################################################### with the latest results being in cm.s-1 However, when I make a slice and plot the x-velocity field, I end up with the figure enclosed, saying that the max velocity is around 1.6e-4 cm.s-1 ! I don't quite understand. Where am I wrong ? Thanks for yoru help, JC
Try making the slice through the position of the peak velocity. I suspect that by some symmetry in your problem you've taken the plane where vx is small. Or, set the zlim on the plot, it looks like you've really clipped the colorbar. d. On Sun, Aug 8, 2010 at 3:09 PM, Jean-Claude Passy <jcpassy@gmail.com> wrote:
Hi all,
there is something regarding the x-velocity field I don't really understand.
When I analyze a file with yt, I find:
###################################################################### In [1]: filen = 'DD0001/CommonEnvelope0001' In [2]: pf = load(filen) In [3]: region = pf.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0], [1.0,1.0, 1.0]) yt INFO 2010-08-08 14:29:57,375 Getting the binary hierarchy yt INFO 2010-08-08 14:29:57,389 Finished with binary hierarchy reading In [4]: vx = region["x-velocity"] yt INFO 2010-08-08 14:30:06,647 Getting field x-velocity from 1 In [5]: max(vx) Out[5]: 232812.74801568719 ######################################################################
with the latest results being in cm.s-1 However, when I make a slice and plot the x-velocity field, I end up with the figure enclosed, saying that the max velocity is around 1.6e-4 cm.s-1 !
I don't quite understand. Where am I wrong ? Thanks for yoru help,
JC
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-- Sent from my Stone Tablet and carried by my Pterodactyl.
Hi David, you were right. There was something wrong with how I set up the zlim. Thanks, JC On 08/08/10 17:40, David Collins wrote:
Try making the slice through the position of the peak velocity. I suspect that by some symmetry in your problem you've taken the plane where vx is small. Or, set the zlim on the plot, it looks like you've really clipped the colorbar.
d.
On Sun, Aug 8, 2010 at 3:09 PM, Jean-Claude Passy<jcpassy@gmail.com> wrote:
Hi all,
there is something regarding the x-velocity field I don't really understand.
When I analyze a file with yt, I find:
###################################################################### In [1]: filen = 'DD0001/CommonEnvelope0001' In [2]: pf = load(filen) In [3]: region = pf.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0], [1.0,1.0, 1.0]) yt INFO 2010-08-08 14:29:57,375 Getting the binary hierarchy yt INFO 2010-08-08 14:29:57,389 Finished with binary hierarchy reading In [4]: vx = region["x-velocity"] yt INFO 2010-08-08 14:30:06,647 Getting field x-velocity from 1 In [5]: max(vx) Out[5]: 232812.74801568719 ######################################################################
with the latest results being in cm.s-1 However, when I make a slice and plot the x-velocity field, I end up with the figure enclosed, saying that the max velocity is around 1.6e-4 cm.s-1 !
I don't quite understand. Where am I wrong ? Thanks for yoru help,
JC
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
David Collins -
Jean-Claude Passy