fix bug in multipole_error

byuflowlab · Feb 1, 2025 · 58c116f · 58c116f
1 parent 2541d3e
commit 58c116f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/error.jl b/src/error.jl
@@ -160,7 +160,7 @@ function multipole_error(t⃗, r_mp, multipole_expansion, P, error_method, lamb_
     rotate_multipole_y!(weights_tmp_2, weights_tmp_1, Ts, Hs_π2, ζs_mag, θ, nmax, lamb_helmholtz)
 
     # multipole error
-    in0 = (P*(P+1))>>1 + 1
+    in0 = (nmax*(nmax+1))>>1 + 1
     ϕn0 = weights_tmp_2[1,1,in0]
     ϕn1_real = weights_tmp_2[1,1,in0+1]
     ϕn1_imag = weights_tmp_2[2,1,in0+1]