kymograph tentatively done, compiles, need to test next

Author: mcarcana

src/post_process/m_data_output.fpp

@@ -1194,7 +1194,9 @@ contains
     subroutine s_write_kymo_data_file(q_prim_vf)
         type(scalar_field), dimension(sys_size), intent(IN) :: q_prim_vf
         integer :: j, k, l, t !< Generic loop iterators
-        real(kind(0d0)) :: vonMises !< selected planes for kymograph comparison
+        real(kind(0d0)) :: vonMises_d, vonMises_h1  !< selected planes for kymograph comparison
+        real(kind(0d0)) :: vonMises_h2, vonMises_h3 !<kymograph plane comparison continued
+        real(kind(0d0)) :: i !< counting integer
         
         do t = 0, t_stop
            do j = 0, m
@@ -1204,27 +1206,23 @@ contains
               vonMises_d = q_prim_vf(xiend+1)%sf(j, 0, 0)
            end do 
         end do   
-        
+
+        i = m/4d0
         do t = 0, t_stop
            do k = 0, n
-              ! for bubliq, sf(j,0,0) for d and sf(0.5*xcen_bub,k,0) for h
-              ! for bubgel, sf(j,0,0) for d, sf(-xcen_bub,k,0) for h_1,
-              !                              sf(-0.5*xcen_bub,k,0) for h_2
-               if (xcen_bub < 0) then !liq
-                   vonMises_h1 = q_prim_vf(xiend+1)%sf(xcen_bub, k, 0)
-                   vonMises_h2 = q_prim_vf(xiend+1)%sf(0.25*xcen_bub, k, 0)
-               else
-                   vonMises_h1 = q_prim_vf(xiend+1)%sf(-xcen_bub, k, 0)
-                   vonMises_h2 = q_prim_vf(xiend+1)%sf(-0.25*xcen_bub, k, 0)
-               end if
+              !xcen_bub is roughly +/- m/4
+              vonMises_h1 = q_prim_vf(xiend+1)%sf(0, k, 0)
+              vonMises_h2 = q_prim_vf(xiend+1)%sf(i, k, 0)
+              vonMises_h3 = q_prim_vf(xiend+1)%sf(-i, k, 0)
            end do
         end do
 
         if (proc_rank == 0) then 
             write (251, '(10X, 8F24.8)') &
                  vonMises_d, &
                  vonMises_h1, &
-                 vonMises_h2
+                 vonMises_h2, &
+                 vonMises_h3
         end if
 
     end subroutine s_write_kymo_data_file