From 13380d479105b7022bee8461ce582e815a7279c7 Mon Sep 17 00:00:00 2001
From: CharlesCNorton <135471798+CharlesCNorton@users.noreply.github.com>
Date: Wed, 9 Sep 2026 10:54:33 -0400
Subject: [PATCH] Pair each analytic vacuum coefficient with the mode the
 caller supplies it for

---
 src/NESTOR/analysum2.f90 | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/src/NESTOR/analysum2.f90 b/src/NESTOR/analysum2.f90
index 22be15eb2..194c879ab 100644
--- a/src/NESTOR/analysum2.f90
+++ b/src/NESTOR/analysum2.f90
@@ -61,12 +61,12 @@ SUBROUTINE analysum2(grpmn, bvec, m, n, l, ivacskip, lasym, m_map, n_map, &
      ! SIN(mu - |n|v) * cmns (l,m,|n|)
      sinp = sinp + temp
 
-     bvec (m, n,  1) = bvec (m, n,  1) + tlp(i)*bexni(i)*sinp
-     bvec (m,-n,  1) = bvec (m,-n,  1) + tlm(i)*bexni(i)*sinm
+     bvec (m, n,  1) = bvec (m, n,  1) + tlm(i)*bexni(i)*sinp
+     bvec (m,-n,  1) = bvec (m,-n,  1) + tlp(i)*bexni(i)*sinm
 
      IF (ivacskip .EQ. 0) THEN
-     grpmn(m, n,i,1) = grpmn(m, n,i,1) + slp(i)         *sinp
-     grpmn(m,-n,i,1) = grpmn(m,-n,i,1) + slm(i)         *sinm
+     grpmn(m, n,i,1) = grpmn(m, n,i,1) + slm(i)         *sinp
+     grpmn(m,-n,i,1) = grpmn(m,-n,i,1) + slp(i)         *sinm
      END IF
 
      IF (lasym) THEN
@@ -79,12 +79,12 @@ SUBROUTINE analysum2(grpmn, bvec, m, n, l, ivacskip, lasym, m_map, n_map, &
         ! COS(mu - |n|v) * cmns (l,m,|n|)
         cosp = cosp + temp
 
-        bvec(m, n,2) = bvec(m, n,2) + tlp(i)*bexni(i)*cosp
-        bvec(m,-n,2) = bvec(m,-n,2) + tlm(i)*bexni(i)*cosm
+        bvec(m, n,2) = bvec(m, n,2) + tlm(i)*bexni(i)*cosp
+        bvec(m,-n,2) = bvec(m,-n,2) + tlp(i)*bexni(i)*cosm
 
         IF (ivacskip .EQ. 0) THEN
-           grpmn(m, n,i,2) = grpmn(m, n,i,2) + slp(i)*cosp
-           grpmn(m,-n,i,2) = grpmn(m,-n,i,2) + slm(i)*cosm
+           grpmn(m, n,i,2) = grpmn(m, n,i,2) + slm(i)*cosp
+           grpmn(m,-n,i,2) = grpmn(m,-n,i,2) + slp(i)*cosm
         END IF
      END IF
   END DO