[yt-users] x-velocity

Jean-Claude Passy jcpassy at gmail.com
Sun Aug 8 15:09:43 PDT 2010


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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20100808/a92edcd3/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 0001_Slice_z_x-velocity.png
Type: image/png
Size: 27425 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20100808/a92edcd3/attachment.png>


More information about the yt-users mailing list