Generalize axes for m

abhro · Aug 2, 2024 · 06f52d7 · 06f52d7
1 parent 6f29f62
commit 06f52d7
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/MagneticFieldSHTCExpander.jl b/src/MagneticFieldSHTCExpander.jl
@@ -99,7 +99,8 @@ function magneticfield(
     end
 
     # set up array of allowed ℓ values
-    ℓ_axes = axes(g)[1]
+    ℓ_axes = axes(g, 1)
+    m_axes = axes(g, 2)
 
     # PHYSICAL THEORY (see Physical-theory.md)
     # Each ℓ,m term of the potential can be written as a product of three
@@ -148,8 +149,7 @@ function magneticfield(
                    /
                    surface_to_surface_scale_denom)
 
-        # FIXME Assuming that m starts at 0. how do we rectify?
-        for m in 0:ℓ
+        for m in m_axes
             # XXX can move this out of the loop?
             G = plm[ℓ,m]
             dG_dθ = -sinθ * dplm[ℓ,m]