JGonzalez/vlaplex
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Reorder print statements and boundary condition output

Browse source
This commit is contained in:
JHendrikx 2025-02-05 16:07:21 +01:00
commit de9ff3a2bb
1 changed files with 4 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

5
vlaplex.f90
View file

@ -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)