From aa71e84c4214859f193b33b977721eddb186d12b Mon Sep 17 00:00:00 2001 From: Spencer Bryngelson Date: Fri, 11 Sep 2026 10:39:43 -0400 Subject: [PATCH] Write vel(2) to 2D probe files and drop the per-step print For a plain 2D case (no bubbles, no hypoelasticity) `s_write_probe_files` wrote only (time, rho, vel(1), pres), silently dropping the second velocity component, and also `print`ed a line to stdout on every time step for every probe. The 3D branch writes all three velocity components, and the 2D bubble and hypoelastic branches write vel(2), so this was the plain-2D branch being incomplete rather than a deliberate choice. The stray `print` floods the job log of any long run. Note this changes the 2D probe file from four columns to five, so 2D probe golden files need regenerating. Fixes #1841 Claude-Session: https://claude.ai/code/session_01HMJ7cycfo7kTFSFq5yhHLG --- src/simulation/m_data_output.fpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/simulation/m_data_output.fpp b/src/simulation/m_data_output.fpp index 57f7eaf9f..394fdb3da 100644 --- a/src/simulation/m_data_output.fpp +++ b/src/simulation/m_data_output.fpp @@ -1590,8 +1590,9 @@ contains & vel(1), vel(2), pres, tau_e(1), tau_e(2), tau_e(3) #:endif else - write (i + 30, '(6X,F12.6,F24.8,F24.8,F24.8)') nondim_time, rho, vel(1), pres - print *, 'time =', nondim_time, 'rho =', rho, 'pres =', pres + #:if not MFC_CASE_OPTIMIZATION or num_dims > 1 + write (i + 30, '(6X,F12.6,F24.8,F24.8,F24.8,F24.8)') nondim_time, rho, vel(1), vel(2), pres + #:endif end if else #:if not MFC_CASE_OPTIMIZATION or num_dims > 2