analyt dbgout

jonathanschilling · Oct 30, 2023 · 86a37ed · 86a37ed
1 parent d08092e
commit 86a37ed
Showing 1 changed file with 9 additions and 5 deletions.
diff --git a/src/NESTOR/analyt.f90 b/src/NESTOR/analyt.f90
@@ -115,17 +115,20 @@ SUBROUTINE analyt(grpmn, bvec, ivacskip, lasym, m_map, n_map, grpmn_m_map, grpmn
   LLOOP: DO l = 0, mf + nf
      fl = fl1
 
+     ! here, tlp/m are available for the current value of l
+     ! --> save into matrix for debugging
+     all_tlp(l,:) = tlp
+     all_tlm(l,:) = tlm
+
      ! COMPUTE SL+ and SL- , Eq (A17)
      ! SLP(M): SL+(-)
      IF (ivacskip .eq. 0) THEN
         slp = (r1p*fl + ra1p)*tlp + r0p*fl*tlp1 - (r1p + r0p)/sqrtc + sign1*(r0p - r1p)/sqrta
         slm = (r1m*fl + ra1m)*tlm + r0m*fl*tlm1 - (r1m + r0m)/sqrtc + sign1*(r0m - r1m)/sqrta
         slpm = slp + slm
 
-        ! here, tlp/m and slp/m are available for the current value of l
+        ! here, slp/m are available for the current value of l
         ! --> save into matrix for debugging
-        all_tlp(l,:) = tlp
-        all_tlm(l,:) = tlm
         all_slp(l,:) = slp
         all_slm(l,:) = slm
      ENDIF
@@ -169,11 +172,12 @@ SUBROUTINE analyt(grpmn, bvec, ivacskip, lasym, m_map, n_map, grpmn_m_map, grpmn
 
     call add_real_3d("all_tlp", mf+nf+1, nv, nu3, all_tlp)
     call add_real_3d("all_tlm", mf+nf+1, nv, nu3, all_tlm)
-    call add_real_2d("bvec", mf1, nf1, bvec)
+    call add_real_2d("bvec", mf1, nf1, bvec) ! (mf+1)x(2*nf+1)x(ndim: 1 or 2)
+
     if (ivacskip .eq. 0) then
       call add_real_3d("all_slp", mf+nf+1, nv, nu3, all_slp)
       call add_real_3d("all_slm", mf+nf+1, nv, nu3, all_slm)
-      call add_real_4d("grpmn", mf1, nf1, nv, nu3, grpmn)
+      call add_real_4d("grpmn", mf1, nf1, nv, nu3, grpmn) ! missing dim: (ndim: 1 or 2)
     else
       call add_null("all_slp")
       call add_null("all_slm")