[yt-users] Trouble plotting 1D scatter data from 3D Athena vtk data.

[Matthew Turk]

Hi Tim,

Thanks for writing, and welcome to yt!

On Sun, Aug 18, 2013 at 8:53 PM, Timothy Waters
<waterst3 at unlv.nevada.edu> wrote:
> Hi,
>
> Here is the data and simple script to make a scatter or profile plot of my
> simulation output,
> https://www.dropbox.com/sh/ibpda938ltsxl0j/2AA2FYDUH8#/
>
> This is a spherically symmetric (cartesian [-0.5,0.5]^3, 32cells^3) hydro
> simulation with initial conditions
> density = 1.0 if radius <= 1/3;
> density falls off as 1/r^2 if r > 1/3;
> velocity is everywhere 0
>
> The output with these ICs is WP12.0000.vtk.  When I run my script the first
> time using iyt I get the first image with the expected ICs.  Running it
> again gives the second image and I don't know why.  I'm hoping those other
> points don't actually exist in my data.  Maybe I am just missing something
> simple?

I'm not sure.  I haven't been able to replicate the problem here, but
I have a handful of suspicions.  I was able to generate these two
images, one for the 0000 dataset and the second for 0030.

http://i.imgur.com/cxnDVjO.png
http://i.imgur.com/MPUSQVH.png

Can you tell me a bit more specifically how you're running the script?
 You start up iyt, then do %run , and then %run again?  Does it happen
if you do:

%run scatter3DAthenaData.py
%run scatter3DAthenaData.py

and look at the resultant plot?  (It does not for me, FWIW.)  Does it
happen if you do multiple datasets in a row?

What's the output of:

yt instinfo

specifically, the changeset hash it provides.

>
> The evolved data, WP12.0030.vtk, is provided as well.  Running this data
> 'cleans the slate' so that I can then go back to see the problem with my
> initial conditions again.  Of course, the problem shows up in this data too.

Interesting.  I do not believe that there is a problem with your data,
but we'll get to the bottom of it!

Best,

Matt

>
> Thanks in advance!
> Tim
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>