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:
--=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
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
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
is center= s. 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:
-- 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's no coordinate translation that is done when you = access just the 'x' 'y' and 'z' fields.=C2=A0 I bel= ieve that the fields you might be looking for are "cylindrical_radius&= quot; "cylindrical_z" and "cylindrical_theta".</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 <<a href=3D" mailto:= <=> guidogranda@gmail.com">guidogranda@gmail.com</a>> 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 <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/ <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