try out swapping back to old VMEC version

jonathanschilling · May 10, 2024 · e04fc4c · e04fc4c
1 parent 49a47d3
commit e04fc4c
Showing 1 changed file with 16 additions and 8 deletions.
diff --git a/src/NESTOR/analysum2.f90 b/src/NESTOR/analysum2.f90
@@ -61,12 +61,16 @@ 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) + 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 ! flipped
+     bvec (m,-n,  1) = bvec (m,-n,  1) + tlp(i)*bexni(i)*sinm ! flipped
 
      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) + 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 ! flipped
+     grpmn(m,-n,i,1) = grpmn(m,-n,i,1) + slp(i)         *sinm ! flipped
      END IF
 
      IF (lasym) THEN
@@ -79,12 +83,16 @@ 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) + 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 ! flipped
+        bvec(m,-n,2) = bvec(m,-n,2) + tlp(i)*bexni(i)*cosm ! flipped
 
         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) + 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 ! flipped
+           grpmn(m,-n,i,2) = grpmn(m,-n,i,2) + slp(i)*cosm ! flipped
         END IF
      END IF
   END DO