Export standard deviation to pandas DataFrames and AstroPy QTables #4410
Conversation
| ) | ||
| assert ( | ||
| prof.standard_deviation[("gas", "density")].units | ||
| == YTArray.from_astropy(at3["density_stddev"]).units |
There was a problem hiding this comment.
shouldn't the field type also be included in the key ("gas_density_stddev" instead of "density_stddev") ?
There was a problem hiding this comment.
@neutrinoceros so the field type is not included already even when we do not include the standard deviation (I don't remember why). So you would access the mean of the profile using at3["density"]
There was a problem hiding this comment.
ok, at least it's consistent.
There was a problem hiding this comment.
We could add the field type, but I'm a bit leery of using _ to connect them since _ could appear in field type strings (and does appear in field name strings), making parsing a bit difficult in some cases.
Is there another sensible character to use to indicate the breaking of the field type and name? Maybe even , is ok? I don't know (cc @matthewturk)
There was a problem hiding this comment.
I guess that would be a breaking change, so I'm actually -1 on this. Let's keep it as is for now.
There was a problem hiding this comment.
Ok then, this is ready to go.
Co-authored-by: Chris Havlin <chris.havlin@gmail.com>
Co-authored-by: Chris Havlin <chris.havlin@gmail.com>
|
Thanks @chrishavlin! |
PR Summary
Profiles in yt can already export data to pandas DataFrames and AstroPy QTables. This PR adds the capability to optionally export the standard deviation of the profile to the data frame or table. Docs and tests are added.
PR Checklist