From 9ba62914f1b3f1f86fad632df1c650d8527c63ee Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Cl=C3=A9ment=20Robert?= Date: Tue, 11 Jul 2023 09:50:49 +0200 Subject: [PATCH 1/2] DOC: fix RadMC3DWriter's docstring (compat with modern yt) --- yt_astro_analysis/radmc3d_export/RadMC3DInterface.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py b/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py index 83ed0762..29d66828 100644 --- a/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py +++ b/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py @@ -110,13 +110,19 @@ class RadMC3DWriter: included in the radiative transfer calculation by radmc3d: >>> import yt + >>> import unyt as un >>> from yt.extensions.astro_analysis.radmc3d_export.api import RadMC3DWriter >>> x_co = 1.0e-4 - >>> mu_h = yt.Quantity(2.34e-24, 'g') + >>> mu_h = un.unyt_quantity(2.34e-24, 'g') >>> def _NumberDensityCO(field, data): ... return (x_co/mu_h)*data["Density"] - >>> yt.add_field("NumberDensityCO", function=_NumberDensityCO) + >>> yt.add_field( + ... ("gas", "NumberDensityCO"), + ... function=_NumberDensityCO, + ... sampling_type="cell", + ... units="cm**-3", + ... ) >>> ds = yt.load("galaxy0030/galaxy0030") >>> writer = RadMC3DWriter(ds) From 88a63634c6c1dd44c4a4ef14691dcdbf8ca3477c Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Cl=C3=A9ment=20Robert?= Date: Thu, 20 Jul 2023 12:30:51 +0200 Subject: [PATCH 2/2] Update yt_astro_analysis/radmc3d_export/RadMC3DInterface.py Co-authored-by: Britton Smith --- yt_astro_analysis/radmc3d_export/RadMC3DInterface.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py b/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py index 29d66828..d44c69c9 100644 --- a/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py +++ b/yt_astro_analysis/radmc3d_export/RadMC3DInterface.py @@ -116,7 +116,7 @@ class RadMC3DWriter: >>> x_co = 1.0e-4 >>> mu_h = un.unyt_quantity(2.34e-24, 'g') >>> def _NumberDensityCO(field, data): - ... return (x_co/mu_h)*data["Density"] + ... return (x_co/mu_h)*data["gas", "density"] >>> yt.add_field( ... ("gas", "NumberDensityCO"), ... function=_NumberDensityCO,