alternative to saving d, h_1, h_2 von Mises stress

Author: mcarcana

src/post_process/m_data_output.fpp

@@ -1198,28 +1198,35 @@ contains
         
         do t = 0, t_stop
            do j = 0, m
-             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
-                 vonMises_d = q_prim_vf(xiend+1)%sf(j, 0, 0)
-                 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
-                 if (proc_rank == 0) then 
-                     write (251, '(10X, 8F24.8)') &
-                     vonMises_d, &
-                     vonMises_h1, &
-                     vonMises_h2
-                 end if
-              end do 
-            end do   
+              vonMises_d = q_prim_vf(xiend+1)%sf(j, 0, 0)
+           end do 
+        end do   
+        
+        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
+           end do
         end do
 
+        if (proc_rank == 0) then 
+            write (251, '(10X, 8F24.8)') &
+                 vonMises_d, &
+                 vonMises_h1, &
+                 vonMises_h2
+        end if
+
     end subroutine s_write_kymo_data_file
 
     subroutine s_close_formatted_database_file() ! -------------------------