diff --git a/environment-dev-linux-clang.yml b/environment-dev-linux-clang.yml
index 67c59e177e..6a16528b5f 100644
--- a/environment-dev-linux-clang.yml
+++ b/environment-dev-linux-clang.yml
@@ -23,7 +23,7 @@ dependencies:
         - llvm-openmp
         - nanobind=2.9
         - nbsphinx
-        - matplotlib
+        - matplotlib<3.11
         - netcdf4
         - ninja
         - numpy>=2
diff --git a/environment-dev-linux.yml b/environment-dev-linux.yml
index 640895060d..aadccbec8c 100644
--- a/environment-dev-linux.yml
+++ b/environment-dev-linux.yml
@@ -15,7 +15,7 @@ dependencies:
         - lark-parser
         - libboost-devel
         - make
-        - matplotlib
+        - matplotlib<3.11
         - nanobind=2.9
         - nbsphinx
         - netcdf4
diff --git a/environment-dev-mac.yml b/environment-dev-mac.yml
index a4a7f114b7..49a6c25804 100644
--- a/environment-dev-mac.yml
+++ b/environment-dev-mac.yml
@@ -19,7 +19,7 @@ dependencies:
         - libcxx-devel
         - libmicrohttpd
         - llvm-openmp
-        - matplotlib
+        - matplotlib<3.11
         - nanobind=2.9
         - nbsphinx
         - netcdf4
diff --git a/environment-dev-win.yml b/environment-dev-win.yml
index 6cc3108147..e89e56c7fe 100644
--- a/environment-dev-win.yml
+++ b/environment-dev-win.yml
@@ -10,7 +10,7 @@ dependencies:
         - libboost-devel
         - libcxx
         - llvm-openmp < 21
-        - matplotlib
+        - matplotlib<3.11
         - nanobind=2.9
         - nbsphinx
         - netcdf4
diff --git a/src/core/lbl/lbl_lineshape_voigt_ecs_hartmann.cpp b/src/core/lbl/lbl_lineshape_voigt_ecs_hartmann.cpp
index 7daa2040c7..e747f85311 100644
--- a/src/core/lbl/lbl_lineshape_voigt_ecs_hartmann.cpp
+++ b/src/core/lbl/lbl_lineshape_voigt_ecs_hartmann.cpp
@@ -54,8 +54,8 @@ Numeric reduced_dipole(const Rational Jf,
                        const Rational li,
                        const Rational k) {
   if (not iseven(Jf + lf + 1))
-    return -sqrtr(2 * Jf + 1) * wigner3j(Jf, k, Ji, li, lf - li, -lf);
-  return +sqrtr(2 * Jf + 1) * wigner3j(Jf, k, Ji, li, lf - li, -lf);
+    return -sqrtr(2 * Jf + 1) * wigner3j(Jf, k, Ji, lf, li - lf, -li);
+  return +sqrtr(2 * Jf + 1) * wigner3j(Jf, k, Ji, lf, li - lf, -li);
 }
 
 void relaxation_matrix_offdiagonal(MatrixView& W,