From de9ff3a2bb8c7d6b9246d4020dba941112238fd3 Mon Sep 17 00:00:00 2001 From: JHendrikx Date: Wed, 5 Feb 2025 16:07:21 +0100 Subject: [PATCH] Reorder print statements and boundary condition output --- vlaplex.f90 | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/vlaplex.f90 b/vlaplex.f90 index 736c5de..dc21882 100644 --- a/vlaplex.f90 +++ b/vlaplex.f90 @@ -240,10 +240,13 @@ program VlaPlEx do t = 1, nt time = t * dt + t0 call boundaryConditions%get(time, n_bc, u_bc, Temp_bc, Zave_bc) - call writeOutputBoundary(t, dt, n_bc, u_bc, Temp_bc, Zave_bc) + print *, 'Time: ', time * t_ref + print *, 'Temp_bc: ', Temp_bc + print *, 'Zave_bc: ', Zave_bc z_inj = minloc(abs(Z_list - T_to_Z(Temp_bc)),1) Zave_bc = Z_list(z_inj) u_bc = sqrt(Zave_bc * Temp_bc) + call writeOutputBoundary(t, dt, n_bc, u_bc, Temp_bc, Zave_bc) ! f0(j0:nv) = v(j0:nv)**2 / sqrt(PI*Temp_bc**3) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc) f0(j0:nv) = 1.0_dp / sqrt(PI*Temp_bc) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc)