From c37f37bf09524eaa048af2f6d390bae8b37beb28 Mon Sep 17 00:00:00 2001 From: JRChreim Date: Mon, 24 Nov 2025 15:52:44 -0800 Subject: [PATCH 1/4] Update m_start_up.fpp --- src/post_process/m_start_up.fpp | 15 +++++++++++++++ 1 file changed, 15 insertions(+) diff --git a/src/post_process/m_start_up.fpp b/src/post_process/m_start_up.fpp index bfa0b1ae9e..3554633a04 100644 --- a/src/post_process/m_start_up.fpp +++ b/src/post_process/m_start_up.fpp @@ -419,6 +419,21 @@ contains end if + ! Adding the individual energies to the formatted database file -------- + if (model_eqns == 3 .and. cons_vars_wrt) then + do i = 1, num_fluids + + q_sf = q_cons_vf(i + intxb - 1)%sf(-offset_x%beg:m + offset_x%end, & + -offset_y%beg:n + offset_y%end, & + -offset_z%beg:p + offset_z%end) + + write (varname, '(A,I0)') 'alpha_rho_e', i + call s_write_variable_to_formatted_database_file(varname, t_step) + + varname(:) = ' ' + end do + end if + !Adding Energy cascade FFT if (fft_wrt) then From 1e62531593de82517be62e9d2771477c7c4c2505 Mon Sep 17 00:00:00 2001 From: JRChreim Date: Mon, 24 Nov 2025 16:25:59 -0800 Subject: [PATCH 2/4] Update m_start_up.fpp --- src/post_process/m_start_up.fpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/post_process/m_start_up.fpp b/src/post_process/m_start_up.fpp index 3554633a04..3b4a6e291e 100644 --- a/src/post_process/m_start_up.fpp +++ b/src/post_process/m_start_up.fpp @@ -424,8 +424,8 @@ contains do i = 1, num_fluids q_sf = q_cons_vf(i + intxb - 1)%sf(-offset_x%beg:m + offset_x%end, & - -offset_y%beg:n + offset_y%end, & - -offset_z%beg:p + offset_z%end) + -offset_y%beg:n + offset_y%end, & + -offset_z%beg:p + offset_z%end) write (varname, '(A,I0)') 'alpha_rho_e', i call s_write_variable_to_formatted_database_file(varname, t_step) @@ -433,7 +433,7 @@ contains varname(:) = ' ' end do end if - + !Adding Energy cascade FFT if (fft_wrt) then From f7dbb7df76bb789e19e21f1ca37d348cbec662ed Mon Sep 17 00:00:00 2001 From: JRChreim Date: Mon, 24 Nov 2025 16:37:10 -0800 Subject: [PATCH 3/4] Update m_start_up.fpp --- src/post_process/m_start_up.fpp | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/src/post_process/m_start_up.fpp b/src/post_process/m_start_up.fpp index 3b4a6e291e..fc19c2288f 100644 --- a/src/post_process/m_start_up.fpp +++ b/src/post_process/m_start_up.fpp @@ -423,10 +423,7 @@ contains if (model_eqns == 3 .and. cons_vars_wrt) then do i = 1, num_fluids - q_sf = q_cons_vf(i + intxb - 1)%sf(-offset_x%beg:m + offset_x%end, & - -offset_y%beg:n + offset_y%end, & - -offset_z%beg:p + offset_z%end) - + q_sf = q_cons_vf(i + intxb - 1)%sf(x_beg:x_end, y_beg:y_end, z_beg:z_end) write (varname, '(A,I0)') 'alpha_rho_e', i call s_write_variable_to_formatted_database_file(varname, t_step) From e0c576cff048de523c1bc38317aece0a860d3c3c Mon Sep 17 00:00:00 2001 From: Ben Wilfong <48168887+wilfonba@users.noreply.github.com> Date: Tue, 25 Nov 2025 20:19:58 -0600 Subject: [PATCH 4/4] add alpha_rho_e_wrt and fix junk code failure --- src/post_process/m_global_parameters.fpp | 2 ++ src/post_process/m_mpi_proxy.fpp | 1 + src/post_process/m_start_up.fpp | 18 ++++++++++-------- toolchain/mfc/run/case_dicts.py | 6 +++--- 4 files changed, 16 insertions(+), 11 deletions(-) diff --git a/src/post_process/m_global_parameters.fpp b/src/post_process/m_global_parameters.fpp index 2274e1152f..6ca4626269 100644 --- a/src/post_process/m_global_parameters.fpp +++ b/src/post_process/m_global_parameters.fpp @@ -240,6 +240,7 @@ module m_global_parameters integer :: flux_lim logical, dimension(3) :: flux_wrt logical :: E_wrt + logical, dimension(num_fluids_max) :: alpha_rho_e_wrt logical :: fft_wrt logical :: pres_wrt logical, dimension(num_fluids_max) :: alpha_wrt @@ -432,6 +433,7 @@ contains down_sample = .false. alpha_rho_wrt = .false. + alpha_rho_e_wrt = .false. rho_wrt = .false. mom_wrt = .false. vel_wrt = .false. diff --git a/src/post_process/m_mpi_proxy.fpp b/src/post_process/m_mpi_proxy.fpp index 8d72568d3e..9f9c3f0c42 100644 --- a/src/post_process/m_mpi_proxy.fpp +++ b/src/post_process/m_mpi_proxy.fpp @@ -124,6 +124,7 @@ contains call MPI_BCAST(mom_wrt(1), 3, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr) call MPI_BCAST(vel_wrt(1), 3, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr) call MPI_BCAST(alpha_rho_wrt(1), num_fluids_max, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr) + call MPI_BCAST(alpha_rho_e_wrt(1), num_fluids_max, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr) call MPI_BCAST(alpha_wrt(1), num_fluids_max, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr) do i = 1, num_fluids_max diff --git a/src/post_process/m_start_up.fpp b/src/post_process/m_start_up.fpp index fc19c2288f..41516305e2 100644 --- a/src/post_process/m_start_up.fpp +++ b/src/post_process/m_start_up.fpp @@ -118,7 +118,8 @@ contains lag_id_wrt, lag_pos_wrt, lag_pos_prev_wrt, lag_vel_wrt, & lag_rad_wrt, lag_rvel_wrt, lag_r0_wrt, lag_rmax_wrt, & lag_rmin_wrt, lag_dphidt_wrt, lag_pres_wrt, lag_mv_wrt, & - lag_mg_wrt, lag_betaT_wrt, lag_betaC_wrt + lag_mg_wrt, lag_betaT_wrt, lag_betaC_wrt, & + alpha_rho_e_wrt ! Inquiring the status of the post_process.inp file file_loc = 'post_process.inp' @@ -419,15 +420,16 @@ contains end if - ! Adding the individual energies to the formatted database file -------- - if (model_eqns == 3 .and. cons_vars_wrt) then + ! Adding the individual energies to the formatted database file + if (model_eqns == 3) then do i = 1, num_fluids + if (alpha_rho_e_wrt(i) .or. cons_vars_wrt) then + q_sf = q_cons_vf(i + intxb - 1)%sf(x_beg:x_end, y_beg:y_end, z_beg:z_end) + write (varname, '(A,I0)') 'alpha_rho_e', i + call s_write_variable_to_formatted_database_file(varname, t_step) - q_sf = q_cons_vf(i + intxb - 1)%sf(x_beg:x_end, y_beg:y_end, z_beg:z_end) - write (varname, '(A,I0)') 'alpha_rho_e', i - call s_write_variable_to_formatted_database_file(varname, t_step) - - varname(:) = ' ' + varname(:) = ' ' + end if end do end if diff --git a/toolchain/mfc/run/case_dicts.py b/toolchain/mfc/run/case_dicts.py index 556ee91896..b94d638d91 100644 --- a/toolchain/mfc/run/case_dicts.py +++ b/toolchain/mfc/run/case_dicts.py @@ -528,9 +528,9 @@ def analytic(self): POST_PROCESS['chem_wrt_T'] = ParamType.LOG for fl_id in range(1,10+1): - for append, ty in [("schlieren_alpha", ParamType.REAL), - ("alpha_rho_wrt", ParamType.LOG), - ("alpha_wrt", ParamType.LOG), ("kappa_wrt", ParamType.LOG)]: + for append, ty in [("schlieren_alpha", ParamType.REAL), ("alpha_rho_wrt", ParamType.LOG), + ("alpha_wrt", ParamType.LOG), ("kappa_wrt", ParamType.LOG), + ("alpha_rho_e_wrt", ParamType.LOG)]: POST_PROCESS[f'{append}({fl_id})'] = ty for real_attr in ["gamma", "pi_inf", "ss", "pv", "gamma_v", "M_v", "mu_v", "k_v", "cp_v",