[yt-users] x-velocity

Jean-Claude Passy jcpassy at gmail.com
Sun Aug 8 18:37:48 PDT 2010


  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 at 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 at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>
>



More information about the yt-users mailing list