outline for kymograph: step 1 - kymograph flag added in post_process

Author: mcarcana

src/post_process/m_global_parameters.fpp

@@ -154,6 +154,7 @@ module m_global_parameters
     logical :: parallel_io    !< Format of the data files
     logical :: sim_data
     logical :: file_per_process !< output format
+    logical :: kymograph
 
     integer, allocatable, dimension(:) :: proc_coords !<
     !! Processor coordinates in MPI_CART_COMM
@@ -392,6 +393,7 @@ contains
         qm_wrt = .false.
         schlieren_wrt = .false.
         sim_data = .false.
+        kymograph = .false.
         cf_wrt = .false.
         ib = .false.
 

src/post_process/m_mpi_proxy.fpp

@@ -197,7 +197,7 @@ contains
 
         #:for VAR in [ 'cyl_coord', 'mpp_lim', 'mixture_err',                  &
             & 'alt_soundspeed', 'hypoelasticity', 'parallel_io', 'rho_wrt',    &
-            & 'E_wrt', 'pres_wrt', 'gamma_wrt', 'sim_data',                    &
+            & 'E_wrt', 'pres_wrt', 'gamma_wrt', 'sim_data', 'kymograph',        &
             & 'heat_ratio_wrt', 'pi_inf_wrt', 'pres_inf_wrt', 'cons_vars_wrt', &
             & 'prim_vars_wrt', 'c_wrt', 'qm_wrt','schlieren_wrt', 'bubbles', 'qbmm',   &
             & 'polytropic', 'polydisperse', 'file_per_process', 'relax', 'cf_wrt',     &

src/post_process/m_start_up.f90

@@ -83,7 +83,7 @@ subroutine s_read_input_file
             relax_model, cf_wrt, sigma, adv_n, ib, num_ibs, &
             cfl_adap_dt, cfl_const_dt, t_save, t_stop, n_start, &
             cfl_target, surface_tension &
-            sim_data, hyperelasticity
+            sim_data, hyperelasticity, kymograph
 
         ! Inquiring the status of the post_process.inp file
         file_loc = 'post_process.inp'
@@ -205,6 +205,14 @@ subroutine s_save_data(t_step, varname, pres, c, H)
             call s_write_energy_data_file(q_prim_vf, q_cons_vf)
         end if
 
+        if (kymograph .and. proc_rank == 0) then
+            call s_open_kymo_data_file()
+        end if
+       
+        if (kymograph) then
+            call s_write_kymo_data_file(q_prim_vf)
+        end if
+
         ! Adding the grid to the formatted database file
         call s_write_grid_to_formatted_database_file(t_step)
 
@@ -709,6 +717,10 @@ subroutine s_save_data(t_step, varname, pres, c, H)
             call s_close_energy_data_file()
         end if
 
+        if (kymograph .and. proc_rank == 0) then
+            call s_close_kymo_data_file()
+        end if
+
         ! Closing the formatted database file
         call s_close_formatted_database_file()