From 137e73b07b7991da5fba8e78db8ac365e70d3175 Mon Sep 17 00:00:00 2001 From: anand43 Date: Wed, 15 Nov 2023 15:24:55 -0500 Subject: [PATCH] fixed a typo --- swiftest/sph_harmonics.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/swiftest/sph_harmonics.py b/swiftest/sph_harmonics.py index 0f1f0b8c9..2f192e86d 100644 --- a/swiftest/sph_harmonics.py +++ b/swiftest/sph_harmonics.py @@ -83,7 +83,8 @@ def clm_from_ellipsoid(mass, density, a, b = None, c = None, lmax = 6, lref_radi ref_radius = shape_SH.expand(normalization = '4pi').coeffs[0, 0, 0] return clm, ref_radius elif(lref_radius == True and ref_radius is not None): - clm = clm.change_ref(r0 = ref_radius) + clm_class = clm_class.change_ref(r0 = ref_radius) + clm = clm_class.to_array(normalization = '4pi') return clm, ref_radius else: return clm @@ -140,7 +141,8 @@ def clm_from_relief(mass, density, grid, lmax = 6, lref_radius = False, ref_radi ref_radius = shape_SH.expand(normalization = '4pi').coeffs[0, 0, 0] return clm, ref_radius elif(lref_radius == True and ref_radius is not None): - clm = clm.change_ref(r0 = ref_radius) + clm_class = clm_class.change_ref(r0 = ref_radius) + clm = clm_class.to_array(normalization = '4pi') return clm, ref_radius else: return clm