Hi Guido,

Those velocities *should* correspond to the *domain-axis-aligned* velocities.  To get the disk velocities, I think you want fields like "velocity_cylindrical_theta" and "velocity_cylindrical_z" and "velocity_cylindrical_radius".

On Tue, Sep 22, 2020 at 11:07 PM guido <guidogranda@gmail.com> wrote:
Hi Matthew,

Thanks for your reply. Then if there is no translation for the spatial coordinates what is happening to the velocities of the cells belonging to the disk?
I'm asking that because what I see on my code is that the velocities and the disk that is on the negative X region show a distribution that corresponds to the middle of the computational domain and the same happens to the distribution of the velocities on the disk that is on the positive X region. However, since the coordinate system of both the cylinder and the computational domain are at rest there shouldn't be a any difference. So, I wonder what are the velocities that I get when I do something like this:

disk['velocity_x']

The fact that the velocities that I'm getting correspond to another region makes me believe that there is something wrong on the disk structure code and it is always taking cells locate on the center that correspond to a center c=[0.5,0.5,0.5]. Is there any way to verify where the cylinder is located in the computational domain?

Regards,



El sáb., 19 sept. 2020 a las 11:02, <yt-users-request@python.org> escribió:
Send yt-users mailing list submissions to
        yt-users@python.org

To subscribe or unsubscribe via the World Wide Web, visit
        https://mail.python.org/mailman3/lists/yt-users.python.org/
or, via email, send a message with subject or body 'help' to
        yt-users-request@python.org

You can reach the person managing the list at
        yt-users-owner@python.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of yt-users digest..."

Today's Topics:

   1. yt.disk (Guido granda muñoz)
   2. Re: yt.disk (Matthew Turk)


----------------------------------------------------------------------

Message: 1
Date: Fri, 18 Sep 2020 16:05:21 -0500
From: Guido granda muñoz <guidogranda@gmail.com>
Subject: [yt-users] yt.disk
To: yt-users@python.org
Message-ID:
        <CALOHZ_0=sr94xUKaxvBuRki99xZtKQtezhx3nEXXqDV9Ly4S1g@mail.gmail.com>
Content-Type: multipart/alternative;
        boundary="000000000000855e3c05af9cdca7"

--000000000000855e3c05af9cdca7
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

Hello yt-users,

I have a question about the disk data object. The first input it needs is
center, which according to your website is the coordinates of the
geometrical center of the disk on the simulation length units. So, the
center of the computational domain is c=3D[0.5,0.5,0.5].

I set the center of two different cylinders in a way that they don=E2=80=99=
t
overlap. The axis of these cylinders are parallel to the x axis and one is
on the negative X coordinate region while the other one on the positive
region. But, I found that they actually have exactly the same cell centers.
Did I do something wrong?

On the following code the box side is 256 pc, so the center of the first
disk c_1=3D[0.25,0.5,0.5] with a height of 32 pc does not overlap with the
second disk of center c_2=3D[0.75,0.5,0.5] also with height 32 pc.

http://paste.yt-project.org/show/275/

I=E2=80=99m using :

---

Version =3D 3.6.dev0

And the snapshot used is here:

http://use.yt/upload/c58a4c0f


--=20
Guido

--000000000000855e3c05af9cdca7
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr"><br clear=3D"all"><div>
=09
=09


<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">Hello
yt-users,</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">I have
a question about the disk data object. The first input it needs is
center, which according to your website is the coordinates of the
geometrical center of the disk on the simulation length units. So,
the center of the computational domain is c=3D[0.5,0.5,0.5].<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">I set
the center of two different cylinders in a way that they don=E2=80=99t
overlap. The axis of these cylinders are parallel to the x axis and
one is on the negative X coordinate region while the other one on the
positive region. But, I found that they actually have exactly the
same cell centers. Did I do something wrong?<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">On the
following code the box side is 256 pc, so the center of the first
disk c_1=3D[0.25,0.5,0.5] with a height of 32 pc does not overlap with
the second disk of center c_2=3D[0.75,0.5,0.5] also with height 32 pc.=20
</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%"><a href=3D"h=
ttp://paste.yt-project.org/show/275/">http://paste.yt-project.org/show/275/=
</a><br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">I=E2=80=99m
using :<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">---</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;border-top:none;border-bottom:=
1px solid rgb(0,0,0);border-left:none;border-right:none;padding:0in 0in 0.0=
3in;line-height:100%">
Version =3D 3.6.dev0</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">And the
snapshot used is here:<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%"><a href=3D"h=
ttp://use.yt/upload/c58a4c0f">http://use.yt/upload/c58a4c0f</a></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%"><br>

</p></div>-- <br><div dir=3D"ltr" class=3D"gmail_signature" data-smartmail=
=3D"gmail_signature">Guido</div></div>

--000000000000855e3c05af9cdca7--

------------------------------

Message: 2
Date: Fri, 18 Sep 2020 17:13:40 -0500
From: Matthew Turk <matthewturk@gmail.com>
Subject: [yt-users] Re: yt.disk
To: Discussion of the yt analysis package <yt-users@python.org>
Message-ID:
        <CALO3=5GozUpChB+gLs-F1toed22XjDM6F4zEG4j4k=9Nu1HAJw@mail.gmail.com>
Content-Type: multipart/alternative;
        boundary="000000000000d95c4e05af9dd04e"

--000000000000d95c4e05af9dd04e
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

Hi Guido,

The cell centers will indeed be the same -- there's no coordinate
translation that is done when you access just the 'x' 'y' and 'z' fields.
I believe that the fields you might be looking for are "cylindrical_radius"
"cylindrical_z" and "cylindrical_theta".

On Fri, Sep 18, 2020 at 4:06 PM Guido granda mu=C3=B1oz <guidogranda@gmail.=
com>
wrote:

>
> Hello yt-users,
>
> I have a question about the disk data object. The first input it needs is
> center, which according to your website is the coordinates of the
> geometrical center of the disk on the simulation length units. So, the
> center of the computational domain is c=3D[0.5,0.5,0.5].
>
> I set the center of two different cylinders in a way that they don=E2=80=
=99t
> overlap. The axis of these cylinders are parallel to the x axis and one i=
s
> on the negative X coordinate region while the other one on the positive
> region. But, I found that they actually have exactly the same cell center=
s.
> Did I do something wrong?
>
> On the following code the box side is 256 pc, so the center of the first
> disk c_1=3D[0.25,0.5,0.5] with a height of 32 pc does not overlap with th=
e
> second disk of center c_2=3D[0.75,0.5,0.5] also with height 32 pc.
>
> http://paste.yt-project.org/show/275/
>
> I=E2=80=99m using :
>
> ---
>
> Version =3D 3.6.dev0
>
> And the snapshot used is here:
>
> http://use.yt/upload/c58a4c0f
>
>
> --
> Guido
> _______________________________________________
> yt-users mailing list -- yt-users@python.org
> To unsubscribe send an email to yt-users-leave@python.org
> https://mail.python.org/mailman3/lists/yt-users.python.org/
> Member address: matthewturk@gmail.com
>

--000000000000d95c4e05af9dd04e
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr">Hi Guido,<div><br></div><div>The cell centers will indeed =
be the same -- there&#39;s no coordinate translation that is done when you =
access just the &#39;x&#39; &#39;y&#39; and &#39;z&#39; fields.=C2=A0 I bel=
ieve that the fields you might be looking for are &quot;cylindrical_radius&=
quot; &quot;cylindrical_z&quot; and &quot;cylindrical_theta&quot;.</div></d=
iv><br><div class=3D"gmail_quote"><div dir=3D"ltr" class=3D"gmail_attr">On =
Fri, Sep 18, 2020 at 4:06 PM Guido granda mu=C3=B1oz &lt;<a href=3D"mailto:=
guidogranda@gmail.com">guidogranda@gmail.com</a>&gt; wrote:<br></div><block=
quote class=3D"gmail_quote" style=3D"margin:0px 0px 0px 0.8ex;border-left:1=
px solid rgb(204,204,204);padding-left:1ex"><div dir=3D"ltr"><br clear=3D"a=
ll"><div>
=09
=09


<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">Hello
yt-users,</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">I have
a question about the disk data object. The first input it needs is
center, which according to your website is the coordinates of the
geometrical center of the disk on the simulation length units. So,
the center of the computational domain is c=3D[0.5,0.5,0.5].<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">I set
the center of two different cylinders in a way that they don=E2=80=99t
overlap. The axis of these cylinders are parallel to the x axis and
one is on the negative X coordinate region while the other one on the
positive region. But, I found that they actually have exactly the
same cell centers. Did I do something wrong?<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">On the
following code the box side is 256 pc, so the center of the first
disk c_1=3D[0.25,0.5,0.5] with a height of 32 pc does not overlap with
the second disk of center c_2=3D[0.75,0.5,0.5] also with height 32 pc.=20
</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%"><a href=3D"h=
ttp://paste.yt-project.org/show/275/" target=3D"_blank">http://paste.yt-pro=
ject.org/show/275/</a><br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">I=E2=80=99m
using :<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">---</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;border-top:none;border-bottom:=
1px solid rgb(0,0,0);border-left:none;border-right:none;padding:0in 0in 0.0=
3in;line-height:100%">
Version =3D 3.6.dev0</p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%">And the
snapshot used is here:<br></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%"><a href=3D"h=
ttp://use.yt/upload/c58a4c0f" target=3D"_blank">http://use.yt/upload/c58a4c=
0f
</a></p>
<p lang=3D"en-US" style=3D"margin-bottom:0in;line-height:100%"><br>

</p></div>-- <br><div dir=3D"ltr">Guido</div></div>
_______________________________________________<br>
yt-users mailing list -- <a href=3D"mailto:yt-users@python.org" target=3D"_=
blank">yt-users@python.org</a><br>
To unsubscribe send an email to <a href=3D"mailto:yt-users-leave@python.org=
" target=3D"_blank">yt-users-leave@python.org</a><br>
<a href=3D"https://mail.python.org/mailman3/lists/yt-users.python.org/" rel=
=3D"noreferrer" target=3D"_blank">https://mail.python.org/mailman3/lists/yt=
-users.python.org/
</a><br>
Member address: <a href=3D"mailto:matthewturk@gmail.com" target=3D"_blank">=
matthewturk@gmail.com</a><br>
</blockquote></div>

--000000000000d95c4e05af9dd04e--

------------------------------

Subject: Digest Footer

_______________________________________________
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org
https://mail.python.org/mailman3/lists/yt-users.python.org/


------------------------------

End of yt-users Digest, Vol 151, Issue 6
****************************************


--
Guido
_______________________________________________
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org
https://mail.python.org/mailman3/lists/yt-users.python.org/
Member address: matthewturk@gmail.com