-
Notifications
You must be signed in to change notification settings - Fork 271
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
ENH: allow boxlib frontend to read 1D cylindrical datasets #4758
Conversation
I'm triaging this as a bug fix for now because it removes a blocking problem to anyone who would need the functionality, and the change is currently very small, but let's re-assess once we see test go green again. |
Actually, failures are clearly unrelated. |
For what it's worth I've issued a fix for those failures at #4759 |
@yut23 can you rebase to current main so we can get CI green ? |
f21f5ab
to
1f4d8fe
Compare
It might make more sense to split out the |
Works for me, but I think this is also fine as is. |
We should probably add a sample data set to test with. I'm not sure what the process for that is, though. |
You're right. And storage basically free for a 1D dataset so let's do it. |
@yt-fido test this please. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looks good to me :)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
there appears to be a problem with the new data sample:
import yt
ds = yt.load_sample("castro_sedov_1d_cyl_plt00150")
ds = yt.load_sample("castro_sedov_1d_cyl_plt00150")
fails with FileNotFoundError
(the exception is raised on the first call if data is already installed, and on the second otherwise). I'm looking into it.
…lindrical datasets
PR Summary
Removes the check that disallows cylindrical coordinates in 1D. It looks like yt learned how to support these sometime since this code was added in 2014. It also fixes an error in
_parse_header_file()
introduced by #4641.I tested against an example dataset output by https://github.com/AMReX-Astro/Castro/tree/main/Exec/hydro_tests/Sedov using
inputs.1d.cyl
, and the cell volumes checked out (proportional to 1/r).PR Checklist